When I'm running the program below I get the following error "( invalid
name"
I know that something's wrong with the post command because when I don't
include it the program runs OK.
Can anybody help?
Dorit
*************************************
program define test
version 7.0
if "`1'" == "?" {
global S_1 "sigmaCK sigmaCL sigmaKL"
exit
}
nl ces lnva93
scalar B0 = _b[B0]
scalar R = _b[R]
scalar Dc = _b[Dc]
scalar Dk = _b[Dk]
scalar Bcc = _b[Bcc]
scalar Bkk = _b[Bkk]
scalar Bll = _b[Bll]
scalar Bck = _b[Bck]
scalar Bcl = _b[Bcl]
scalar Bkl = _b[Bkl]
scalar C = 14.86249
scalar lnC = 2.69884
scalar K = 576.6429
scalar lnK = 6.357223
scalar L = 486.6226
scalar lnL = 6.187489
scalar Z = Dc*C^(-R)+Dk*K^(-R)+(1-Dc-Dk)*L^(-R)
scalar V = exp(B0 - (1/R)*ln(Dc*(C^(-R)) + Dk*(K^(-R)) +
(1-Dc-Dk)*(L^(-R))) + Bcc*(lnC^2) + Bkk*(lnK^2) + Bll*(lnL^2) + Bck*lnC*lnK
+ Bcl*lnC*lnL + Bkl*lnK*lnL) scalar Fc =
V*(Dc*C^(-R-1)/Z+Bck/C*lnK+Bcl/C*lnL+2*Bcc/C*lnC)
scalar Fk = V*(Dk*K^(-R-1)/Z+Bck/K*lnC+Bkl/K*lnL+2*Bkk/K*lnK)
scalar Fl = V*((1-Dc-Dk)*L^(-R-1)/Z+Bcl/L*lnC+Bkl/L*lnK+2*Bll/L*lnL)
scalar Fcc = Fc^2/V - Fc/C +
V*((-R*Dc*C^(-R-2))/Z+(Dc^2*C^(-2*R-2)*R)/Z^2+2*Bcc/C^2)
scalar Fkk = Fk^2/V - Fk/K +
V*((-R*Dk*K^(-R-2))/Z+(Dk^2*K^(-2*R-2)*R)/Z^2+2*Bkk/K^2)
scalar Fll = Fl^2/V - Fl/L +
V*((-R*(1-Dc-Dk)*L^(-R-2))/Z+((1-Dc-Dk)^2*L^(-2*R-2)*R)/Z^2+2*Bll/L^2)
scalar Fck = (Fc*Fk)/V +
V*((R*Dc*Dk*C^(-R-1)*K^(-R-1)/Z^2)+Bck/(C*K))
scalar Fcl = (Fc*Fl)/V +
V*((R*Dc*(1-Dc-Dk)*C^(-R-1)*L^(-R-1)/Z^2)+Bcl/(C*L))
scalar Fkl = (Fk*Fl)/V +
V*((R*Dk*(1-Dc-Dk)*K^(-R-1)*L^(-R-1)/Z^2)+Bkl/(K*L))
matrix H1 = (0,Fc,Fk,Fl\Fc, Fcc,
Fck,Fcl\Fk,Fck,Fkk,Fkl\Fl,Fcl,Fkl,Fll)
scalar DH = det(H1)
matrix Hck = (0, Fc,Fl\Fk,Fck,Fkl\Fl, Fcl,Fll)
scalar DHck = - det(Hck)
scalar sigma_ck = ((C*Fc+K*Fk+L*Fl)/(C*K))*(DHck/DH)
matrix Hcl = (0,Fc,Fk\ Fk,Fck,Fkk\ Fl,Fcl,Fkl)
scalar DHcl = det(Hcl)
scalar sigma_cl = ((C*Fc+K*Fk+L*Fl)/(C*L))*(DHcl/DH)
matrix Hkl = (0,Fc,Fk\Fc,Fcc,Fck\Fl,Fcl,Fkl)
scalar DHkl = - det(Hkl)
scalar sigma_kl = ((C*Fc+K*Fk+L*Fl)/(K*L))*(DHkl/DH)
post `1' (sigma_ck) (sigma_cl) (sigma_kl)
end
*******************************************************
Dorit Nevo
PhD Candidate, MIS
Faculty of Commerce, University of British Columbia
http://people.commerce.ubc.ca/phd/dnevo
(604) 822-4772
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
