好吧 学妹把一道练习题搬上来了 求指导

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

题：x=（x1 x2 x3）'服从三元正态分布，现有20个观测值
         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

问：检验H0：u=u0=(4 50 10)' ; H1:u不等于u0


题很简单  但是初涉sas 程序不会写~求助各位大侠了

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：练习题 求指导 SOSO 各位大侠 正态分布 练习题 指导

在线等！

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 发表于 2013-5-6 14:03
data  ex;
input x1-x3;
cards;

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

Imasasor 发表于 2013-5-6 15:07
ziyenano好久不来，一看是学妹，居然耐心写这么长代码

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

那就不看

ziyenano 发表于 2013-5-6 14:03
data  ex;
input x1-x3;
cards;

十分感谢~不过程序运行不出来耶
？？？？

silvia1992 发表于 2013-5-6 22:11
十分感谢~不过程序运行不出来耶
？？？？

没有问题啊，
最后输出F,F_alpha,prob_F三个值。

ziyenano 发表于 2013-5-6 14:03
data  ex;
input x1-x3;
cards;

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

Σ未知的时候 应该转化成的是T分布吧~C值等于一个常数与Fa的值，而不是F分布

silvia1992 发表于 2013-5-6 22:11
十分感谢~不过程序运行不出来耶
？？？？

没有问题啊，
最后输出F F_alpha prob_F三个值