应用多元分析第四章---多元正态总体.pptVIP

下载本文档

2
0
约2.95万字
约 152页
2019-07-01 发布于浙江
举报
版权申诉

应用多元分析第四章---多元正态总体.ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 或等价于其中例1 proc iml; n1=20;n2=20;n3=20; n=n1+n2+n3;k=3; p=4; x1={ 260 75 40 18 , 200 72 34 17 , 240 87 45 18 , 170 65 39 17 , 270 110 39 24 , 205 130 34 23 , 190 69 27 15 , 200 46 45 15 , 250 117 21 20 , 200 107 28 20 , 225 130 36 11 , 210 125 26 17 , 170 64 31 14 , 270 76 33 13 , 190 60 34 16 , 280 81 20 18 , 310 119 25 15 , 270 57 31 8 , 250 67 31 14 , 260 135 39 29 } ; x2={ 310 122 30 21 ,310 60 35 18 ,190 40 27 15 ,225 65 34 16 , 170 65 37 16 ,210 82 31 17 ,280 67 37 18 ,210 38 36 17 , 280 65 30 23 ,200 76 40 17 ,200 76 39 20 ,280 94 26 11 , 190 60 33 17 ,295 55 30 16 ,270 125 24 21 ,280 120 32 18 , 240 62 32 20 ,280 69 29 20 ,370 70 30 20 ,280 40 37 17 }; x3={ 320 64 39 17 ,260 59 37 11 ,360 88 28 26 ,295 100 36 12 , 270 65 32 21 ,380 114 36 21 ,240 55 42 10 ,260 55 34 20 , 260 110 29 20 ,295 73 33 21 ,240 114 38 18 ,310 103 32 18 , 330 112 21 11 ,345 127 24 20 ,250 62 22 16 ,260 59 21 19 , 225 100 34 30 ,345 120 36 18 ,360 107 25 23 ,250 117 36 16 }; xx=x1//x2//x3; ln={[20] 1} ; x10=(ln*x1)/n1; print x10; mm1=i(n1)-j(n1,n1,1)/n1; mm=i(n)-j(n,n,1)/n; s1=x1`*mm1*x1/(n1-1); s2=x2`*mm1*x2/(n2-1); s3=x3`*mm1*x3/(n3-1); tt=xx`*mm*xx/(n-1); s=(s1+s2+s3)*(n1-1)/(n-3); ds1=det(s1); ds2=det(s2); ds3=det(s3); ds=det(s); d=(2*p**2+3*p-1)*(k+1)/((6*p+1)*(n-k)); f=p*(p+1)*(k-1)/2; m=(n-k)*log(det(s))-(n1- 1)*(log(det(s1))+ log(det(s2))+log(det(s3))); kc=(1-d)*m; p0=1-probchi(kc,df) ; print s1 s2 s3 s f d m kc; 记其中对例1 中的数据，判断三个组(即三个总体)的均值向量和协方差矩阵是否全相等(α=0.05) proc iml; n1=20;n2=20;n3=20