在线等！

data  ex;
input x1-x3;
cards;
         3.7 48.5 9.3
         5.7 65.1 8.0
         3.8 47.2 10.9
         3.1 55.5 9.7
         4.6 36.1 7.9
         2.4 24.8 14
         7.2 33.1 7.6
         6.7 47.4 8.5
         5.4 54.1 11.3
         3.9 36.9 12.7
         4.5 58.8 12.3
         3.5 27.8 9.8
         4.5 40.2 8.4
         1.5 13.5 10.1
         8.5 56.4 7.1
         4.5 71.6 8.2
         6.5 52.8 10.9
         4.1 44.1 11.2
         5.5 40.9 9.4
;
run;

/*构造统计量，若协方差阵Σ已知n(X平均-u0)'*inv(Σ)(X平均-u0)，服从卡方分布
  若Σ未知，以样本协方差阵代替Σ，上述统计量为霍特林统计量，可以转为F分布*/

proc iml;
use ex;
read all into X;
close ex;
start cov(x);
n=nrow(x);
cov=(x-x[:,])`*(x-x[:,])/(n-1);
return(cov);
finish;
n=nrow(X);
p=ncol(X);
u0={4 50 10};
T2=n*(X[:,]-u0)*inv(cov(X))*(X[:,]-u0)`;
F=(n-p)/(p*(n-1))*T2;
F_alpha=finv(0.95,p,n-p);
prob_F=1-cdf('F',F,p,n-p);
print F F_alpha prob_F;
/*如果F>F_alpha，在0.05置信水平下，拒绝原假设，反之没有足够理由拒绝原假设。
也可以看prob_F值*/

ziyenano好久不来，一看是学妹，居然耐心写这么长代码

这么直白的戳穿我，让我情何以堪~

那就不看