prints sur estimation results
prtsys(res,out)
* res = a tlist provided by a sur, twosls or threesls function
* out = the file where the results are printed (optional; default: %io(2), i.e the screen)
* nothing, just print the results
load(GROCERDIR+'/data/grun.dat') // grunfeld investment data // see page 650, Green 1997 ige = grun(:,1); // general electric fge = grun(:,2); cge = grun(:,3); iwest = grun(:,4); // westinghouse fwest = grun(:,5); cwest = grun(:,6); igm = grun(:,7); // general motors fgm = grun(:,8); cgm = grun(:,9); ich = grun(:,10); // chrysler fch = grun(:,11); cch = grun(:,12); iuss = grun(:,13); // us steel fuss = grun(:,14); cuss = grun(:,15); // (order follows that in Green, 1997) eq1='igm=a1*fgm+a2*cgm+a3' eq2='ich=a4*fch+a5*cch+a6' eq3='ige=a7*fge+a8*cge+a9' eq4='iwest=a10*fwest+a11*cwest+a12' eq5='iuss=a13*fuss+a14*cuss+a15' r=sur(eq1,eq2,eq3,eq4,eq5,'noprint') prtsys(r,%io(2)) | ![]() | ![]() |