我用GAMS软件解决StoNED模型,运行中出现错误报告(第91行),修改多次仍然没有解决。求高手帮忙!正确解决者奖励100论坛币。
GAMS Rev 240 WIN-VS8 24.0.2 x86/MS Windows 06/06/14 13:36:47 Page 1
G e n e r a l A l g e b r a i c M o d e l i n g S y s t e m
C o m p i l a t i o n
1
2
3
4 SETS
5 i index of industries/i1*i14/
6 t time periods/t1*t7/;
7 alias(i,h)
8 alias(t,s);
9
10 PARAMETERS
11 table Y(i,t) output
.....
28 table K(i,t) capital
.....
45 table L(i,t) labor
.....
62 Parameter Tr(t) /t1 1,t2 2,t3 3,t4 4,t5 5,t6 6,t7 7/;
63
64 Parameter Tr2(t)/t1 1,t2 4,t3 9,t4 16,t5 25,t6 36,t7 49/;
65
66 VARIABLES
67 SS sum of squares of errors
68 V(i,t) residual
69 ThetaK capital trend coefficient
70 PsiK capital trend coefficient 2
71 Thetal labor trend coefficient
72 PsiL labor trend coefficient 2
73 A(i) inefficiency constant
74 B(i) inefficiency slope
75 C(i) inefficiency quadratic;
76
77 POSITIVE VARIABLES
78 BK(i,t) marginal product of capital
79 BL(i,t) marginal product of labor;
80
81 EQUATIONS
82 QSSE objective=sum of squares of errors
83 QREGRESSION(i,t) regression equation
84 QCONCAVITY(i,h,t,s) concavity constraint
85 QKCON(t) nonnegative capital embodied technical progress
86 QLCON(t) nonnegative labor embodied technical progress
87 QEFF(i,t) nonnegative inefficiency;
88
89 QSSE.. SS=e=sum(t,sum(i,V(i,t)*V(i,t)));
90 QREGRESSION(i,t).. Y(i,t) =e= BK(i,t)*K(i,t)+BL(i,t)*L(i,t)+(ThetaK*Tr(t
)+PsiK*Tr2(t))*K(i,t)+(Thetal*Tr(t)+PsiL*Tr2(t))*L(i,t)-(A(i)+B(i)*Tr(t)+C
(i)*Tr2(t))*Y(i,t)+Y(i,t)*V(i,t);
91 QCONCAVITY(i,h,t,s).. BK(i,t)*K(i,t)+BL(i,t)*L(i,t) =1= BK(h,s)*K(i,t)+BL(
**** $37,409
h,s)*L(i,t);
92 QKCON(t).. ThetaK*Tr(t)+PsiK*Tr2(t) =g= 0;
93 QLCON(t).. ThetaL*Tr(t)+PsiL*Tr2(t) =g= 0;
94
95 MODEL StoNED/all/
96
97 SOLVE StoNED using NLP Minimizing SS;
**** $257
GAMS Rev 240 WIN-VS8 24.0.2 x86/MS Windows 06/06/14 13:36:47 Page 2
G e n e r a l A l g e b r a i c M o d e l i n g S y s t e m
Error Messages
37 '=l=' or '=e=' or '=g=' operator expected
257 Solve statement not checked because of previous errors
409 Unrecognizable item - skip to find a new statement
looking for a ';' or a key word to get started again
**** 3 ERROR(S) 0 WARNING(S)
COMPILATION TIME = 0.000 SECONDS 3 Mb WIN240-240 Feb 14, 2013
USER: GAMS Development Corporation, Washington, DC G871201/0000CA-ANY
Free Demo, 202-342-0180,
sales@gams.com,
www.gams.com DC0000
**** FILE SUMMARY
Input C:\Users\wmt\Desktop\Untitled_1234.gms
Output C:\Users\wmt\Documents\gamsdir\projdir\Untitled_1234.lst
**** USER ERROR(S) ENCOUNTERED