proc iml;

start main;



start main;

y={ 10,

    20,

    17,

    12,

    11 };

x={ 1 6  28,

    1 12 40,

    1 10 32,

    1  8 36,

    1  9 34 };



b=inv(x`*x)*x`*y;

Yhat=x*b; 

e=y-x*b;

sse=e`*e;

df=1/2;

rms=sse#df;

vb=inv(x`*x)#rms;

Vyhat=x*inv(x`*x)*x`#rms;

print b vb y yhat e sse rms vyhat;

finish main;

run;quit;


data one;

input fatty_acid over_wt age;

cards;

10  6  28

20 12 40

17 10 32

12  8 36

11  9 34

;

proc glm;

model fatty_acid=over_wt age / solution;

run;

quit;