proc iml;
start main;

y={ 2.28460,
    2.06908,
    1.70690,
    2.02030,
    1.40335,
    1.20441,
    1.22212,
    0.66843,
    1.03762,
    1.06027};
V={0.72 0 0 0 0 0 0 0 0 0,
   0 6.00 0 0 0 0 0 0 0 0,
   0 0 12.39 0 0 0 0 0 0 0,
   0 0 0 11.30 0 0 0 0 0 0,
   0 0 0 0 16.97 0 0 0 0 0,
   0 0 0 0 0 19.45 0 0 0 0,
   0 0 0 0 0 0 23.45 0 0 0,
   0 0 0 0 0 0 0 33.52 0 0,
   0 0 0 0 0 0 0 0 40.45 0,
   0 0 0 0 0 0 0 0 0 34.46};

X={1 0.000,
   1 0.157,
   1 0.281,
   1 0.333,
   1 0.396,
   1 0.459,
   1 0.491,
   1 0.584,
   1 0.624,
   1 0.657};

b=inv(x`*x)*x`*y;
Yhat=x*b; 
e=y-x*b;
sse=e`*e;
df=1/8;
rms=sse#df;
vb=inv(x`*x)#rms;
print b vb;

b=inv(x`*inv(v)*x)*x`*inv(v)*y;
Yhat=x*b; 
e=y-x*b;
sse=e`*e;
df=1/8;
rms=sse#df;
vb=inv(x`*inv(v)*x)#rms;
print b vb;

finish main;
run;
quit;