var yF cF xF kF lF zF tlF txF gF pF O; varexo ezF elF exF egF epF; parameters gamma delta beta eta nu teta css lss pss rhoz rhozz rhozl rhozx rhozg rhol rholz rholl rholx rholg rhox rhoxz rhoxl rhoxx rhoxg rhog rhogz rhogl rhogx rhogg rhop rhopz rhopl rhopx rhopg rhopp rhozp rholp rhoxp rhogp; gamma=1.00567267; delta=0.025; nu=1.5; beta=0.99; teta=0.34; eta=2.3094; css=0.57; lss=0.21; pss=1; model; eta*nu*(lF)^(nu - 1)= (1-tlF)*(1-teta)*exp(zF)*kF^(teta)*lF^(-teta); (cF(+1) - eta*(lF(+1))^(nu))*(1+txF)*gamma = beta*((cF - eta*lF^(nu))*teta*exp(zF)*(kF)^(teta - 1)*(lF)^(1 - teta) + (1+txF(+1))*(1-delta)); yF = exp(zF)*kF^(teta)*lF^(1-teta); gamma*kF = xF + (1-delta)*kF(-1); yF = cF + xF + gF; O = ((css)/(css - eta*lss^nu))*((cF - css)/(css)) - ((eta*nu*lss^nu)/(css - eta*lss^nu))*((lF-lss)/(lss)) - (pF - pss)/pss; zF = rhoz + rhozz*zF(-1) + rhozl*tlF(-1) + rhozx*txF(-1) + rhozg*gF(-1) + rhozp*pF(-1) + ezF; tlF = rhol + rholz*zF(-1) + rholl*tlF(-1) + rholx*txF(-1) + rholg*gF(-1) + rholp*pF(-1) + elF; txF = rhox + rhoxz*zF(-1) + rhoxl*tlF(-1) + rhoxx*txF(-1) + rhoxg*gF(-1) + rhoxp*pF(-1) + exF; gF = rhog + rhogz*zF(-1) + rhogl*tlF(-1) + rhogx*txF(-1) + rhogg*gF(-1) + rhogp*pF(-1) + egF; pF = rhop + rhopz*zF(-1) + rhopl*tlF(-1) + rhopx*txF(-1) + rhopg*gF(-1) + rhopp*pF(-1) + epF; end; initval; cF = 0.57; xF = 0.23; kF = 1; gF = 0.2044397; lF = 0.21; zF = 0; pF = 1; tlF = 0; txF = 0; O = 0.001; end; estimated_params_bounds; rhoz, -100, 100; rhozz, -100, 100; rhozl, -100, 100; rhozx, -100, 100; rhozg, -100, 100; rhol, -100, 100; rholz, -100, 100; rholl, -100, 100; rholx, -100, 100; rholg, -100, 100; rhox, -100, 100; rhoxz, -100, 100; rhoxl, -100, 100; rhoxx, -100, 100; rhoxg, -100, 100; rhog, -100, 100; rhogz, -100, 100; rhogl, -100, 100; rhogx, -100, 100; rhogg, -100, 100; rhop, -100, 100; rhopz, -100, 100; rhopl, -100, 100; rhopx, -100, 100; rhopg, -100, 100; rhopp, -100, 100; rhoxp, -100, 100; rhogp, -100, 100; rholp, -100, 100; rhozp, -100, 100; end; varobs yF xF cF lF kF O; estimation(datafile = data2); steady;