习题4.5.sas

SAS 源码 《应用多元分析》王学民编著的SAS程序

proc iml;
x1={441.16,0.13,-3.36};
x2={505.97,0.09,-4.57};
x3={432.51,0.14,-3.31};
s1={294.74 -0.60 -32.57,-0.60 0.0013 0.073,-32.57 0.073 4.23};
s2={1596.18 -1.19 -91.05,-1.19 0.0009 0.071,-91.05 0.071 5.76};
s3={182.67 -0.42 -22.00,-0.42 0.001 0.056,-22.00 0.056 3.14};
x=(x1+x2+x3)/3;
sse=75*(s1+s2+s3);
sstr=76*((x1-x)*t(x1-x)+(x2-x)*t(x2-x)+(x3-x)*t(x3-x));
sst=sse+sstr;
k=det(sse)/det(sst);
print k;