Code: Select all
system(model=NZMod1)
variables loil FFR lrgdp lgdpdef lm1 day90 lner
lags 1 to nlags
det constant TREND
end(system)
*
estimate(print,resids=varresids)
....
dec frml[rect] afrml
nonlin c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 c12 c13 c14 c15 c16
frml afrml = ||1.0,0.0,0.0,0.0,0.0,0.0,0.0|$
c1,1.0,0.0,0.0,0.0,0.0,0.0,0.0|$
c2,0.0,1.0,0.0,0.0,0.0,0.0|$
c3,0.0,c4,1.0,0.0,0.0,0.0|$
0.0,0.0,c5,c6,1.0,c7,0.0|$
c8,0.0,0.0,0.0,c9,1.0,c10|$
c11,c12,c13,c14,c15,c16,1.0||
compute c1=c2=c3=c4=c5=c6=c7=c8=c9=c10=c11=c12=c13=c14=c15=c16=0.0
cvmodel(a=afrml,title="NZ model 1",method=bfgs,pmethod=simplex,piters=500,factor=afactor) %sigma
Thanks