Re: Elasticities in Almost Idela demand system

maurizio — Thu, 22 Jun 2017 15:27:28 GMT

Thank yuo.

Here is the log. The SAS version is 9.4. I suppose the IML version varies accordingly

Regards

NOTE: Copyright (c) 2002-2012 by SAS Institute Inc., Cary, NC, USA.
NOTE: SAS (r) Proprietary Software 9.4 (TS1M0)
      Licensed to UNIVERSITA CATTOLICA DEL SACRO CUORE, Site 70123393.
NOTE: This session is executing on the X64_7PRO platform.

NOTE: Updated analytical products:

      SAS/STAT 12.3 (maintenance)
      SAS/ETS 12.3 (maintenance)
      SAS/OR 12.3 (maintenance)
      SAS/IML 12.3 (maintenance)
      SAS/QC 12.3 (maintenance)

NOTE: Additional host information:

X64_7PRO WIN 6.1.7601 Service Pack 1 Workstation

NOTE: SAS initialization used:
      real time           2.39 seconds
      cpu time            1.26 seconds

NOTE: Writing HTML Body file: sashtml.htm
1    Proc model data=bau.consumo3;
2
3    witabread_less= B0+B1*LP_ITABREAD_LESS+ B2*LP_ITACHEESE_LESS + B3*LP_ITADAIRYUH_LESS +
3 ! B4*LP_ITAFORNO_LESS + B5*LP_ITAPASTA_LESS + B6*LP_LATTE_LESS + B7*LP_PFRESCA_LESS +
3 ! B8*LP_SUBITABREAD + B9*LP_SUBITACHEESE + B10*LP_SUBITADAIRYUHT + B11*LP_SUBITAFORNO +
3 ! B12*LP_SUBITALATTE + B13*LP_SUBITAPASTA + B14*LP_SUBITAPFRESCA + B15*LP_OTHERPLUS + B16*LREXP
3 ! ;
4
5
6    WITACHEESE_LESS = C0+C1*LP_ITABREAD_LESS+ C2*LP_ITACHEESE_LESS + C3*LP_ITADAIRYUH_LESS +
6 ! C4*LP_ITAFORNO_LESS + C5*LP_ITAPASTA_LESS + C6*LP_LATTE_LESS + C7*LP_PFRESCA_LESS +
6 ! C8*LP_SUBITABREAD + C9*LP_SUBITACHEESE + C10*LP_SUBITADAIRYUHT + C11*LP_SUBITAFORNO +
6 ! C12*LP_SUBITALATTE + C13*LP_SUBITAPASTA + C14*LP_SUBITAPFRESCA + C15*LP_OTHERPLUS+C16*LREXP ;
7
8
9
10   WITADAIRY_UHTLESS = U0+U1*LP_ITABREAD_LESS+ U2*LP_ITACHEESE_LESS + U3*LP_ITADAIRYUH_LESS +
10 ! U4*LP_ITAFORNO_LESS + U5*LP_ITAPASTA_LESS + U6*LP_LATTE_LESS + U7*LP_PFRESCA_LESS +
10 ! U8*LP_SUBITABREAD + U9*LP_SUBITACHEESE + U10*LP_SUBITADAIRYUHT + U11*LP_SUBITAFORNO +
10 ! U12*LP_SUBITALATTE + U13*LP_SUBITAPASTA + U14*LP_SUBITAPFRESCA + U15*LP_OTHERPLUS+ U16*LREXP
10 ! ;
11
12
13
14   WITAHOVEN_LESS = H0+H1*LP_ITABREAD_LESS+ H2*LP_ITACHEESE_LESS + H3*LP_ITADAIRYUH_LESS +
14 ! H4*LP_ITAFORNO_LESS + H5*LP_ITAPASTA_LESS + H6*LP_LATTE_LESS + H7*LP_PFRESCA_LESS +
14 ! H8*LP_SUBITABREAD + H9*LP_SUBITACHEESE + H10*LP_SUBITADAIRYUHT + H11*LP_SUBITAFORNO +
14 ! H12*LP_SUBITALATTE + H13*LP_SUBITAPASTA + H14*LP_SUBITAPFRESCA + H15*LP_OTHERPLUS+ H16*LREXP
14 ! ;
15
16
17   WITA_PASTALESS = P0+P1*LP_ITABREAD_LESS+ P2*LP_ITACHEESE_LESS + P3*LP_ITADAIRYUH_LESS +
17 ! P4*LP_ITAFORNO_LESS + P5*LP_ITAPASTA_LESS + P6*LP_LATTE_LESS + P7*LP_PFRESCA_LESS +
17 ! P8*LP_SUBITABREAD + P9*LP_SUBITACHEESE + P10*LP_SUBITADAIRYUHT + P11*LP_SUBITAFORNO +
17 ! P12*LP_SUBITALATTE + P13*LP_SUBITAPASTA + P14*LP_SUBITAPFRESCA + P15*LP_OTHERPLUS + P16*LREXP
17 ! ;
18
19
20
21   WITADAIRY_LESS = D0+D1*LP_ITABREAD_LESS+ D2*LP_ITACHEESE_LESS + D3*LP_ITADAIRYUH_LESS +
21 ! D4*LP_ITAFORNO_LESS + D5*LP_ITAPASTA_LESS + D6*LP_LATTE_LESS + D7*LP_PFRESCA_LESS +
21 ! D8*LP_SUBITABREAD + D9*LP_SUBITACHEESE + D10*LP_SUBITADAIRYUHT + D11*LP_SUBITAFORNO +
21 ! D12*LP_SUBITALATTE + D13*LP_SUBITAPASTA + D14*LP_SUBITAPFRESCA + D15*LP_OTHERPLUS+ D16*LREXP
21 ! ;
22
23
24   WITAFRESHP_LESS =F0+ F1*LP_ITABREAD_LESS+ F2*LP_ITACHEESE_LESS + F3*LP_ITADAIRYUH_LESS +
24 ! F4*LP_ITAFORNO_LESS + F5*LP_ITAPASTA_LESS + F6*LP_LATTE_LESS + F7*LP_PFRESCA_LESS +
24 ! F8*LP_SUBITABREAD + F9*LP_SUBITACHEESE + F10*LP_SUBITADAIRYUHT + F11*LP_SUBITAFORNO +
24 ! F12*LP_SUBITALATTE + F13*LP_SUBITAPASTA + F14*LP_SUBITAPFRESCA + F15*LP_OTHERPLUS+F16*LREXP ;
25
26
27
28   WITABREAD_SUB = S10+S11*LP_ITABREAD_LESS+ S12*LP_ITACHEESE_LESS + S13*LP_ITADAIRYUH_LESS +
28 ! S14*LP_ITAFORNO_LESS + S15*LP_ITAPASTA_LESS + S16*LP_LATTE_LESS + S17*LP_PFRESCA_LESS +
28 ! S18*LP_SUBITABREAD + S19*LP_SUBITACHEESE + S110*LP_SUBITADAIRYUHT + S111*LP_SUBITAFORNO +
28 ! S112*LP_SUBITALATTE + S113*LP_SUBITAPASTA + S114*LP_SUBITAPFRESCA +
28 ! S115*LP_OTHERPLUS+S116*LREXP ;
29
30
31
32   WSUB_ITACHEESE= S30+S31*LP_ITABREAD_LESS+ S32*LP_ITACHEESE_LESS + S33*LP_ITADAIRYUH_LESS +
32 ! S34*LP_ITAFORNO_LESS + S35*LP_ITAPASTA_LESS + S36*LP_LATTE_LESS + S37*LP_PFRESCA_LESS +
32 ! S38*LP_SUBITABREAD + S39*LP_SUBITACHEESE + S310*LP_SUBITADAIRYUHT + S311*LP_SUBITAFORNO +
32 ! S312*LP_SUBITALATTE + S313*LP_SUBITAPASTA + S314*LP_SUBITAPFRESCA + S315*LP_OTHERPLUS +
32 ! S316*LREXP;
33
34
35
36
37   WSUB_ITAHOVEN =S50+S51*LP_ITABREAD_LESS+S52*LP_ITACHEESE_LESS + S53*LP_ITADAIRYUH_LESS +
37 ! S54*LP_ITAFORNO_LESS + S55*LP_ITAPASTA_LESS + S56*LP_LATTE_LESS + S57*LP_PFRESCA_LESS +
37 ! S58*LP_SUBITABREAD + S59*LP_SUBITACHEESE + S510*LP_SUBITADAIRYUHT + S511*LP_SUBITAFORNO +
37 ! S512*LP_SUBITALATTE + S513*LP_SUBITAPASTA + S514*LP_SUBITAPFRESCA + S515*LP_OTHERPLUS+
37 ! S516*LREXP ;
38
39
40   WITADAIRY_SUB=S20+ S21*LP_ITABREAD_LESS+ S22*LP_ITACHEESE_LESS + S23*LP_ITADAIRYUH_LESS +
40 ! S24*LP_ITAFORNO_LESS + S25*LP_ITAPASTA_LESS + S26*LP_LATTE_LESS + S27*LP_PFRESCA_LESS +
40 ! S28*LP_SUBITABREAD + S29*LP_SUBITACHEESE + S210*LP_SUBITADAIRYUHT + S211*LP_SUBITAFORNO +
40 ! S212*LP_SUBITALATTE + S213*LP_SUBITAPASTA + S214*LP_SUBITAPFRESCA + S215*LP_OTHERPLUS+
40 ! S216*LREXP ;
41
42   WSUB_ITAPASTA = S60+S61*LP_ITABREAD_LESS+ S62*LP_ITACHEESE_LESS + S63*LP_ITADAIRYUH_LESS +
42 ! S64*LP_ITAFORNO_LESS + S65*LP_ITAPASTA_LESS + S66*LP_LATTE_LESS + S67*LP_PFRESCA_LESS +
42 ! S68*LP_SUBITABREAD + S69*LP_SUBITACHEESE + S610*LP_SUBITADAIRYUHT + S611*LP_SUBITAFORNO +
42 ! S612*LP_SUBITALATTE + S613*LP_SUBITAPASTA + S614*LP_SUBITAPFRESCA + S615*LP_OTHERPLUS+
42 ! S616*LREXP ;
43
44
45
46   WSUB_ITAFRESHP = S40+S41*LP_ITABREAD_LESS+ S42*LP_ITACHEESE_LESS + S43*LP_ITADAIRYUH_LESS +
46 ! S44*LP_ITAFORNO_LESS + S45*LP_ITAPASTA_LESS + S46*LP_LATTE_LESS + S47*LP_PFRESCA_LESS +
46 ! S48*LP_SUBITABREAD + S49*LP_SUBITACHEESE + S410*LP_SUBITADAIRYUHT + S411*LP_SUBITAFORNO +
46 ! S412*LP_SUBITALATTE + S413*LP_SUBITAPASTA + S414*LP_SUBITAPFRESCA + S415*LP_OTHERPLUS +
46 ! S416*LREXP;
47
48   WOTHERPLUS = O0+O1*LP_ITABREAD_LESS+ O2*LP_ITACHEESE_LESS + O3*LP_ITADAIRYUH_LESS +
48 ! O4*LP_ITAFORNO_LESS + O5*LP_ITAPASTA_LESS + O6*LP_LATTE_LESS + O7*LP_PFRESCA_LESS +
48 ! O8*LP_SUBITABREAD + O9*LP_SUBITACHEESE + O10*LP_SUBITADAIRYUHT + O11*LP_SUBITAFORNO +
48 ! O12*LP_SUBITALATTE + O13*LP_SUBITAPASTA + O14*LP_SUBITAPFRESCA + O15*LP_OTHERPLUS+ O16*LREXP;
49
50
51

52   FIT witabread_less WITACHEESE_LESS WITADAIRY_UHTLESS WITAHOVEN_LESS WITA_PASTALESS
52 ! WITADAIRY_LESS WITAFRESHP_LESS WITABREAD_SUB WSUB_ITACHEESE WSUB_ITAHOVEN WITADAIRY_SUB
52 ! WSUB_ITAPASTA WSUB_ITAFRESHP WOTHERPLUS /ITSUR NESTIT OUTS=REST OUTEST=BAU.CONSUMO5
52 ! OUT=RESID2 CONVERGE=.00001 MAXIT=1000;
53   PARMS B0 B1 B2 B3 B4 B5 B6 B7 B8 B9 B10 B11 B12 B13 B14 B15 B16
54
55   C0 C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16
56
57   U0 U1 U2 U3 U4 U5 U6 U7 U8 U9 U10 U11 U12 U13 U14 U15 U16
58   H0 H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11 H12 H13 H14 H15 H16
59
60   P0 P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11 P12 P13 P14 P15 P16
61   D0 D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 D16
62   F0 F1 F2 F3 F4 F5 F6 F7 F8 F9 F10 F11 F12 F13 F14 F15 F16
63   S10 S11 S12 S13 S14 S15 S16 S17 S18 S19 S110 S111 S112 S113 S114 S115 S116
64   S30 S31 S32 S33 S34 S35 S36 S37 S38 S39 S310 S311 S312 S313 S314 S315 S316
65   S50 S51 S52 S53 S54 S55 S56 S57 S58 S59 S510 S511 S512 S513 S514 S515 S516
66   S20 S21 S22 S23 S24 S25 S26 S27 S28 S29 S210 S211 S212 S213 S214 S215 S216
67   S60 S61 S62 S63 S64 S65 S66 S67 S68 S69 S610 S611 S612 S613 S614 S615 S616
68   S40 S41 S42 S43 S44 S45 S46 S47 S48 S49 S410 S411 S412 S413 S414 S415 S416
69   O0 O1 O2 O3 O4 O5 O6 O7 O8 O9 O10 O11 O12 O13 O14 O15 O16;
70   run;

NOTE: At ITSUR Iteration 2 CONVERGE=0.00001 Criteria Met.
NOTE: The data set WORK.RESID2 has 60 observations and 33 variables.
NOTE: The data set BAU.CONSUMO5 has 1 observations and 242 variables.
NOTE: The data set WORK.REST has 14 observations and 17 variables.

NOTE: PROCEDURE MODEL used (Total process time):
      real time           1:17.70
      cpu time            7.36 seconds

71   Proc means data=bau.consumo3 noprint;
72   Var;
73
74   OUTPUT OUT=MEANW MEAN= witabread_less WITACHEESE_LESS WITADAIRY_UHTLESS WITAHOVEN_LESS
74 ! WITA_PASTALESS WITADAIRY_LESS WITAFRESHP_LESS WITABREAD_SUB WSUB_ITACHEESE WSUB_ITAHOVEN
74 ! WITADAIRY_SUB WSUB_ITAPASTA WSUB_ITAFRESHP WOTHERPLUS LREXP;
75   run;

NOTE: There were 60 observations read from the data set BAU.CONSUMO3.
NOTE: The data set WORK.MEANW has 1 observations and 17 variables.
NOTE: PROCEDURE MEANS used (Total process time):
      real time           0.11 seconds
      cpu time            0.09 seconds

OLE Error: 80040064. Struttura FORMATETC non valida.
76   PROC IML;
NOTE: IML Ready
77   use bau.consumo5;
78   read all var { B1 B2 B3 B4 B5 B6 B7 B8 B9 B11 B12 B13 B14 B15 B16
79   C1 C2 C3 C4 C5 C6 C7 C8 C9 C11 C12 C13 C14 C15 C16
80
81   U1 U2 U3 U4 U5 U6 U7 U8 U9 U11 U12 U13 U14 U15 U16
82   H1 H2 H3 H4 H5 H6 H7 H8 H9 H11 H12 H13 H14 H15 H16
83
84   P1 P2 P3 P4 P5 P6 P7 P8 P9 P11 P12 P13 P14 P15 P16
85   D1 D2 D3 D4 D5 D6 D7 D8 D9 D11 D12 D13 D14 D15 D16
86   F1 F2 F3 F4 F5 F6 F7 F8 F9 F11 F12 F13 F14 F15 F16
87   S11 S12 S13 S14 S15 S16 S17 S18 S19 S111 S112 S113 S114 S115 S116
88   S31 S32 S33 S34 S35 S36 S37 S38 S39 S311 S312 S313 S314 S315 S316
89   S51 S52 S53 S54 S55 S56 S57 S58 S59 S511 S512 S513 S514 S515 S516
90   S21 S22 S23 S24 S25 S26 S27 S28 S29 S211 S212 S213 S214 S215 S216
91   S61 S62 S63 S64 S65 S66 S67 S68 S69 S611 S612 S613 S614 S615 S616
92   S41 S42 S43 S44 S45 S46 S47 S48 S49 S411 S412 S413 S414 S415 S416
93   O1 O2 O3 O4 O5 O6 O7 O8 O9 O11 O12 O13 O14 O15 O16};
94   close bau.consumo5;
95
96   /*READ MEANS*/
97
98   use meanw;
99   READ ALL VAR{ witabread_less WITACHEESE_LESS WITADAIRY_UHTLESS WITAHOVEN_LESS WITA_PASTALESS
99 ! WITADAIRY_LESS WITAFRESHP_LESS WITABREAD_SUB WSUB_ITACHEESE WSUB_ITAHOVEN WITADAIRY_SUB
99 ! WSUB_ITAPASTA WSUB_ITAFRESHP WOTHERPLUS LREXP };
100
101
102 close meanw;
103 /*BUDGET SHARE VECTOR*/
104 W=witabread_less // WITACHEESE_LESS // WITADAIRY_UHTLESS //WITAHOVEN_LESS//WITA_PASTALESS
104! //WITADAIRY_LESS //WITAFRESHP_LESS //WITABREAD_SUB //WSUB_ITACHEESE //WSUB_ITAHOVEN//
104! //WITADAIRY_SUB//WSUB_ITAPASTA //WSUB_ITAFRESHP// WOTHERPLUS;
     --
     22
     76
ERROR 22-322: Syntax error, expecting one of the following: a name, a quoted string,
              a numeric constant, a datetime constant, a missing value, +.

ERROR 76-322: Syntax error, statement will be ignored.

105 /*GAMMA IJ MATRIX */
106 /* TOLTO (16) PERCHE’ E’ SU REDDITO tolta costante e numero 10*/
107 gij=(B1||B2||B3||B4||B5||B6||B7||B8||B9||B11||B12||B13||B14||B15)//
108
109 (C1||C2||C3||C4||C5||C6||C7||C8||C9||C11||C12||C13||C14||C15)//
110
111 (U1||U2||U3||U4||U5||U6||U7||U8||U9||U11||U12| |U13||U14| U15)//
                                                    -
                                                    22
                                                    76
ERROR 22-322: Syntax error, expecting one of the following: a name, a quoted string,
              a numeric constant, a datetime constant, a missing value, +.

ERROR 76-322: Syntax error, statement will be ignored.

112 (H1||H2||H3||H4||H5||H6||H7||H8| |H9||H11| |H12||H13||H14||H15|)//
113
114 (P1||P2||P3||P4||P5||P6||P7||P8||P9||P11||P12||P13||P14||P15)//
115 (D1||D2||D3||D4||D5||D6||D7||D8||D9||D11||D12||D13||D14||D15)//
116 (F1||F2||F3||F4||F5||F6||F7||F8||F9||F11||F12||F13||F14||F15)//
117 (S11||S12||S13||S14||S15||S16||S17||S18||S19||S111||S112||S113| |S114||S115)//
118 (S30||S31||S32||S33||S34||S35||S36||S37||S38||S39||S310||S311||S312||S313| |S314||S315)//
119 (S51||S52||S53||S54||S55||S56||S57||S58||S59||S511||S512||S513| |S514||S515)//
120 (S21||S22||S23||S24||S25||S26||S27||S28||S29||S210||S211||S212||S213| |S214||S215)//
121 (S61||S62||S63||S64||S65||S66||S67||S68||S69||S611||S612||S613| |S614||S615)//
122 (S41||S42|S43||S44||S45||S46||S47||S48||S49||S410||S411||S412||S413||S414| |S415)//
123 (O1||O2||O3||O4||O5||O6||O7||O8||O9||O11| |O12||O13||O14||O15);
124 /* SUBITADAIRYUHT PARAMETERS*/
125 /* Bisognerebbe mettere una equazione per W SUBITADAIRYUHT
126 =0- B16-C16 U16-H16-P16-D16-F16-S116-S316-S516-S216-S616-   S416-O16;*/
127
128 /*CALCULATE EACH ELEMENT OF THE ELASTICITY MATRIX*/
129 nk=ncol(gi);
130 mi=-1#I(nk);
ERROR: (execution) Invalid operand to operation.

operation : I at line 130 column 8
operands : nk

nk      1 row       1 col     (numeric)

         0

statement : ASSIGN at line 130 column 1
131 ff2=j(nk,nk,0);
ERROR: (execution) Invalid operand to operation.

operation : J at line 131 column 6
operands : nk, nk, *LIT1004

nk      1 row       1 col     (numeric)

         0

nk      1 row       1 col     (numeric)

         0

*LIT1004      1 row       1 col     (numeric)

         0

statement : ASSIGN at line 131 column 1
131!                         /*initialize marshallian elasticity */
132 fic2=j(nk,nk,0);
ERROR: (execution) Invalid operand to operation.

operation : J at line 132 column 7
operands : nk, nk, *LIT1005

nk      1 row       1 col     (numeric)

         0

nk      1 row       1 col     (numeric)

         0

*LIT1005      1 row       1 col     (numeric)

         0

statement : ASSIGN at line 132 column 1
132!                    /*Initialize Hicksian elasticity*/
133 fi2=j(nk,1,0) /*Income elasticity vector*/
134 /*plotting elasticities*/
135 /*initialize index vector */
136
137 x = j(nk*nk,1,0);
     -
     22
     76
ERROR 22-322: Syntax error, expecting one of the following: #, ##, &, (|, *, **, +, -, /, //, :,
              <, <=, <>, =, >, ><, >=, @, [, ^=, |, ||.

ERROR 76-322: Syntax error, statement will be ignored.

138    y = j(nk*nk,1,0);
ERROR: (execution) Invalid operand to operation.

operation : J at line 138 column 8
operands : _TEM1001, *LIT1008, *LIT1009

_TEM1001      1 row       1 col     (numeric)

         0

*LIT1008      1 row       1 col     (numeric)

         1

*LIT1009      1 row       1 col     (numeric)

         0

statement : ASSIGN at line 138 column 3
139
140    /* initialize vector to store elasticity matrices */
141    Helast = j(nk*nk,1,0);
ERROR: (execution) Invalid operand to operation.

operation : J at line 141 column 13
operands : _TEM1001, *LIT1010, *LIT1011

_TEM1001      1 row       1 col     (numeric)

         0

*LIT1010      1 row       1 col     (numeric)

         1

*LIT1011      1 row       1 col     (numeric)

         0

statement : ASSIGN at line 141 column 3
142    Melast = j(nk*nk,1,0);
ERROR: (execution) Invalid operand to operation.

operation : J at line 142 column 13
operands : _TEM1001, *LIT1012, *LIT1013

_TEM1001      1 row       1 col     (numeric)

         0

*LIT1012      1 row       1 col     (numeric)

         1

*LIT1013      1 row       1 col     (numeric)

         0

statement : ASSIGN at line 142 column 3
143
144    i=1;
145    do i=1 to nk;
146       fi2[i,1] = 1 + b[i,]/w[i,];
147       j=1;
148       do j=1 to nk;
149          ff2[i,j] = mi[i,j] + (gij[i,j] - b[i,]#(w[j,]-b[j,]#(lx-p)))/w[i,];
150          fic2[i,j] = ff2[i,j] + w[j,]#fi2[i,];
151          x[(i-1)*nk+j,1] = i ;
152    y[(i-1)*nk+j,1] = j ;
153    Melast[(i-1)*nk+j,1] = ff2[i,j] ;
154    Helast[(i-1)*nk+j,1] = fic2[i,j] ;
155       end;
156    end;
157 run;
NOTE: Module MAIN is undefined in IML; cannot be RUN.
NOTE: Exiting IML.
NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE IML used (Total process time):
      real time           2:22.51
      cpu time            3.54 seconds

158 PROC IML;
NOTE: IML Ready
159 use bau.consumo5;
160 read all var { B1 B2 B3 B4 B5 B6 B7 B8 B9 B11 B12 B13 B14 B15 B16
161 C1 C2 C3 C4 C5 C6 C7 C8 C9 C11 C12 C13 C14 C15 C16
162
163 U1 U2 U3 U4 U5 U6 U7 U8 U9 U11 U12 U13 U14 U15 U16
164 H1 H2 H3 H4 H5 H6 H7 H8 H9 H11 H12 H13 H14 H15 H16
165
166 P1 P2 P3 P4 P5 P6 P7 P8 P9 P11 P12 P13 P14 P15 P16
167 D1 D2 D3 D4 D5 D6 D7 D8 D9 D11 D12 D13 D14 D15 D16
168 F1 F2 F3 F4 F5 F6 F7 F8 F9 F11 F12 F13 F14 F15 F16
169 S11 S12 S13 S14 S15 S16 S17 S18 S19 S111 S112 S113 S114 S115 S116
170 S31 S32 S33 S34 S35 S36 S37 S38 S39 S311 S312 S313 S314 S315 S316
171 S51 S52 S53 S54 S55 S56 S57 S58 S59 S511 S512 S513 S514 S515 S516
172 S21 S22 S23 S24 S25 S26 S27 S28 S29 S211 S212 S213 S214 S215 S216
173 S61 S62 S63 S64 S65 S66 S67 S68 S69 S611 S612 S613 S614 S615 S616
174 S41 S42 S43 S44 S45 S46 S47 S48 S49 S411 S412 S413 S414 S415 S416
175 O1 O2 O3 O4 O5 O6 O7 O8 O9 O11 O12 O13 O14 O15 O16};
176 close bau.consumo5;
177
178 /*READ MEANS*/
179
180 use meanw;
181 READ ALL VAR{ witabread_less WITACHEESE_LESS WITADAIRY_UHTLESS WITAHOVEN_LESS WITA_PASTALESS
181! WITADAIRY_LESS WITAFRESHP_LESS WITABREAD_SUB WSUB_ITACHEESE WSUB_ITAHOVEN WITADAIRY_SUB
181! WSUB_ITAPASTA WSUB_ITAFRESHP WOTHERPLUS LREXP };
182 close meanw;
183 run;
NOTE: Module MAIN is undefined in IML; cannot be RUN.
NOTE: Exiting IML.
NOTE: PROCEDURE IML used (Total process time):
      real time           1:39.34
      cpu time            1.95 seconds

184 PROC IML;
NOTE: IML Ready
185 use bau.consumo5;
186 read all var { B1 B2 B3 B4 B5 B6 B7 B8 B9 B11 B12 B13 B14 B15 B16
187 C1 C2 C3 C4 C5 C6 C7 C8 C9 C11 C12 C13 C14 C15 C16
188
189 U1 U2 U3 U4 U5 U6 U7 U8 U9 U11 U12 U13 U14 U15 U16
190 H1 H2 H3 H4 H5 H6 H7 H8 H9 H11 H12 H13 H14 H15 H16
191
192 P1 P2 P3 P4 P5 P6 P7 P8 P9 P11 P12 P13 P14 P15 P16
193 D1 D2 D3 D4 D5 D6 D7 D8 D9 D11 D12 D13 D14 D15 D16
194 F1 F2 F3 F4 F5 F6 F7 F8 F9 F11 F12 F13 F14 F15 F16
195 S11 S12 S13 S14 S15 S16 S17 S18 S19 S111 S112 S113 S114 S115 S116
196 S31 S32 S33 S34 S35 S36 S37 S38 S39 S311 S312 S313 S314 S315 S316
197 S51 S52 S53 S54 S55 S56 S57 S58 S59 S511 S512 S513 S514 S515 S516
198 S21 S22 S23 S24 S25 S26 S27 S28 S29 S211 S212 S213 S214 S215 S216
199 S61 S62 S63 S64 S65 S66 S67 S68 S69 S611 S612 S613 S614 S615 S616
200 S41 S42 S43 S44 S45 S46 S47 S48 S49 S411 S412 S413 S414 S415 S416
201 O1 O2 O3 O4 O5 O6 O7 O8 O9 O11 O12 O13 O14 O15 O16};
202 close bau.consumo5;
203
204 /*READ MEANS*/
205
206 use meanw;
207 READ ALL VAR{ witabread_less WITACHEESE_LESS WITADAIRY_UHTLESS WITAHOVEN_LESS WITA_PASTALESS
207! WITADAIRY_LESS WITAFRESHP_LESS WITABREAD_SUB WSUB_ITACHEESE WSUB_ITAHOVEN WITADAIRY_SUB
207! WSUB_ITAPASTA WSUB_ITAFRESHP WOTHERPLUS LREXP };
208 close meanw;
209 quit;
NOTE: Exiting IML.
NOTE: PROCEDURE IML used (Total process time):
      real time           0.04 seconds
      cpu time            0.03 seconds

210 /*BUDGET SHARE VECTOR*/
211 W=witabread_less // WITACHEESE_LESS // WITADAIRY_UHTLESS //WITAHOVEN_LESS//WITA_PASTALESS
     -
     180
211! //WITADAIRY_LESS //WITAFRESHP_LESS //WITABREAD_SUB //WSUB_ITACHEESE //WSUB_ITAHOVEN//
211! //WITADAIRY_SUB//WSUB_ITAPASTA //WSUB_ITAFRESHP// WOTHERPLUS;

ERROR 180-322: Statement is not valid or it is used out of proper order.

212 /*GAMMA IJ MATRIX */
213 /* TOLTO (16) PERCHE’ E’ SU REDDITO tolta costante e numero 10*/
214 gij=(B1||B2||B3||B4||B5||B6||B7||B8||B9||B11||B12||B13||B14||B15)//
     ---
     180
ERROR 180-322: Statement is not valid or it is used out of proper order.
215
216 (C1||C2||C3||C4||C5||C6||C7||C8||C9||C11||C12||C13||C14||C15)//
217
218 (U1||U2||U3||U4||U5||U6||U7||U8||U9||U11||U12| |U13||U14| U15)//
219 (H1||H2||H3||H4||H5||H6||H7||H8| |H9||H11| |H12||H13||H14||H15|)//
220
221 (P1||P2||P3||P4||P5||P6||P7||P8||P9||P11||P12||P13||P14||P15)//
222 (D1||D2||D3||D4||D5||D6||D7||D8||D9||D11||D12||D13||D14||D15)//
223 (F1||F2||F3||F4||F5||F6||F7||F8||F9||F11||F12||F13||F14||F15)//
224 (S11||S12||S13||S14||S15||S16||S17||S18||S19||S111||S112||S113| |S114||S115)//
225 (S30||S31||S32||S33||S34||S35||S36||S37||S38||S39||S310||S311||S312||S313| |S314||S315)//
226 (S51||S52||S53||S54||S55||S56||S57||S58||S59||S511||S512||S513| |S514||S515)//
227 (S21||S22||S23||S24||S25||S26||S27||S28||S29||S210||S211||S212||S213| |S214||S215)//
228 (S61||S62||S63||S64||S65||S66||S67||S68||S69||S611||S612||S613| |S614||S615)//
229 (S41||S42|S43||S44||S45||S46||S47||S48||S49||S410||S411||S412||S413||S414| |S415)//
230 (O1||O2||O3||O4||O5||O6||O7||O8||O9||O11| |O12||O13||O14||O15);

231 /* SUBITADAIRYUHT PARAMETERS*/
232 /* Bisognerebbe mettere una equazione per W SUBITADAIRYUHT
233 =0- B16-C16 U16-H16-P16-D16-F16-S116-S316-S516-S216-S616-   S416-O16;*/
234
235 /*CALCULATE EACH ELEMENT OF THE ELASTICITY MATRIX*/
236 nk=ncol(gi);
     --
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

237 mi=-1#I(nk);
     --
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

238 ff2=j(nk,nk,0);
     ---
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

239 fic2=j(nk,nk,0);
     ----
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

240 fi2=j(nk,1,0) /*Income elasticity vector*/
     ---
     180
ERROR 180-322: Statement is not valid or it is used out of proper order.
241 /*plotting elasticities*/
242 /*initialize index vector */
243
244 x = j(nk*nk,1,0);

245    y = j(nk*nk,1,0);
       -
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

246
247    /* initialize vector to store elasticity matrices */
248    Helast = j(nk*nk,1,0);
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

249    Melast = j(nk*nk,1,0);
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

250
251    i=1;
       -
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

252    do i=1 to nk;
       --
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

253       fi2[i,1] = 1 + b[i,]/w[i,];
          ---
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

254       j=1;
          -
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

255       do j=1 to nk;
          --
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

256          ff2[i,j] = mi[i,j] + (gij[i,j] - b[i,]#(w[j,]-b[j,]#(lx-p)))/w[i,];
             ---
             180

ERROR 180-322: Statement is not valid or it is used out of proper order.

257          fic2[i,j] = ff2[i,j] + w[j,]#fi2[i,];
             ----
             180

ERROR 180-322: Statement is not valid or it is used out of proper order.

258          x[(i-1)*nk+j,1] = i ;
             -
             180

ERROR 180-322: Statement is not valid or it is used out of proper order.

259    y[(i-1)*nk+j,1] = j ;
       -
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

260    Melast[(i-1)*nk+j,1] = ff2[i,j] ;
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

261    Helast[(i-1)*nk+j,1] = fic2[i,j] ;
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

262       end;
          ---
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

263    end;
       ---
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

264 W=witabread_less // WITACHEESE_LESS // WITADAIRY_UHTLESS //WITAHOVEN_LESS//WITA_PASTALESS
     -
     180
264! //WITADAIRY_LESS //WITAFRESHP_LESS //WITABREAD_SUB //WSUB_ITACHEESE //WSUB_ITAHOVEN//
264! //WITADAIRY_SUB//WSUB_ITAPASTA //WSUB_ITAFRESHP// WOTHERPLUS;

ERROR 180-322: Statement is not valid or it is used out of proper order.

265 /*GAMMA IJ MATRIX */
266 /* TOLTO (16) PERCHE’ E’ SU REDDITO tolta costante e numero 10*/
267 gij=(B1||B2||B3||B4||B5||B6||B7||B8||B9||B11||B12||B13||B14||B15)//
     ---
     180
ERROR 180-322: Statement is not valid or it is used out of proper order.
268
269 (C1||C2||C3||C4||C5||C6||C7||C8||C9||C11||C12||C13||C14||C15)//
270
271 (U1||U2||U3||U4||U5||U6||U7||U8||U9||U11||U12| |U13||U14| U15)//
272 (H1||H2||H3||H4||H5||H6||H7||H8| |H9||H11| |H12||H13||H14||H15|)//
273
274 (P1||P2||P3||P4||P5||P6||P7||P8||P9||P11||P12||P13||P14||P15)//
275 (D1||D2||D3||D4||D5||D6||D7||D8||D9||D11||D12||D13||D14||D15)//
276 (F1||F2||F3||F4||F5||F6||F7||F8||F9||F11||F12||F13||F14||F15)//
277 (S11||S12||S13||S14||S15||S16||S17||S18||S19||S111||S112||S113| |S114||S115)//
278 (S30||S31||S32||S33||S34||S35||S36||S37||S38||S39||S310||S311||S312||S313| |S314||S315)//
279 (S51||S52||S53||S54||S55||S56||S57||S58||S59||S511||S512||S513| |S514||S515)//
280 (S21||S22||S23||S24||S25||S26||S27||S28||S29||S210||S211||S212||S213| |S214||S215)//
281 (S61||S62||S63||S64||S65||S66||S67||S68||S69||S611||S612||S613| |S614||S615)//
282 (S41||S42|S43||S44||S45||S46||S47||S48||S49||S410||S411||S412||S413||S414| |S415)//
283 (O1||O2||O3||O4||O5||O6||O7||O8||O9||O11| |O12||O13||O14||O15);

284 /* SUBITADAIRYUHT PARAMETERS*/
285 /* Bisognerebbe mettere una equazione per W SUBITADAIRYUHT
286 =0- B16-C16 U16-H16-P16-D16-F16-S116-S316-S516-S216-S616-   S416-O16;*/
287
288 /*CALCULATE EACH ELEMENT OF THE ELASTICITY MATRIX*/
289 nk=ncol(gi);
     --
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

290 mi=-1#I(nk);
     --
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

291 ff2=j(nk,nk,0);
     ---
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

292 fic2=j(nk,nk,0);
     ----
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

293 fi2=j(nk,1,0) /*Income elasticity vector*/
     ---
     180
ERROR 180-322: Statement is not valid or it is used out of proper order.
294 /*plotting elasticities*/
295 /*initialize index vector */
296
297 x = j(nk*nk,1,0);

298    y = j(nk*nk,1,0);
       -
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

299
300    /* initialize vector to store elasticity matrices */
301    Helast = j(nk*nk,1,0);
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

302    Melast = j(nk*nk,1,0);
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

303
304    i=1;
       -
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

305    do i=1 to nk;
       --
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

306       fi2[i,1] = 1 + b[i,]/w[i,];
          ---
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

307       j=1;
          -
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

308       do j=1 to nk;
          --
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

309          ff2[i,j] = mi[i,j] + (gij[i,j] - b[i,]#(w[j,]-b[j,]#(lx-p)))/w[i,];
             ---
             180

ERROR 180-322: Statement is not valid or it is used out of proper order.

310          fic2[i,j] = ff2[i,j] + w[j,]#fi2[i,];
             ----
             180

ERROR 180-322: Statement is not valid or it is used out of proper order.

311          x[(i-1)*nk+j,1] = i ;
             -
             180

ERROR 180-322: Statement is not valid or it is used out of proper order.

312    y[(i-1)*nk+j,1] = j ;
       -
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

313    Melast[(i-1)*nk+j,1] = ff2[i,j] ;
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

314    Helast[(i-1)*nk+j,1] = fic2[i,j] ;
       ------
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

315       end;
          ---
          180

ERROR 180-322: Statement is not valid or it is used out of proper order.

316    end;
       ---
       180

ERROR 180-322: Statement is not valid or it is used out of proper order.

317 /*BUDGET SHARE VECTOR*/
318 W = witabread_less // WITACHEESE_LESS // WITADAIRY_UHTLESS //WITAHOVEN_LESS//WITA_PASTALESS
     -
     180
318! //WITADAIRY_LESS //WITAFRESHP_LESS //WITABREAD_SUB //WSUB_ITACHEESE //WSUB_ITAHOVEN//
318! //WITADAIRY_SUB//WSUB_ITAPASTA //WSUB_ITAFRESHP// WOTHERPLUS;

ERROR 180-322: Statement is not valid or it is used out of proper order.

319 W = witabread_less // WITACHEESE_LESS // WITADAIRY_UHTLESS ;
     -
     180

ERROR 180-322: Statement is not valid or it is used out of proper order.

topic Re: Elasticities in Almost Idela demand system in SAS/IML Software and Matrix Computations

Elasticities in Almost Idela demand system

Re: Elasticities in Almost Idela demand system

Re: Elasticities in Almost Idela demand system

Re: Elasticities in Almost Idela demand system

Re: Elasticities in Almost Idela demand system

Re: Elasticities in Almost Idela demand system

Re: Elasticities in Almost Idela demand system