/*
** Lesson 2.4: Data Analysis
*/
output file = gpe\output2.4 reset;
load x[17,7]=gpe\longley.txt;
y=x[2:17,7];
x=x[2:17,1:6]~ones(16,1);
call dstat(0,y~x);

b=y/x; @ b=invpd(x'x)*x'y=solvpd(x'x,x'y) @
print "Least Suares Estimator b = " b;

xx=x'*x;
r=eigrs(xx);
cn=sqrt(maxc(r)./minc(r));
print "Condition Number of X = " cn cond(x);
end;