proc iml;
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;