资源描述：
《第三章 多元正态总体参数假设检验》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、§3.2例3.2.1x=[3.748.59.3;5.765.18;3.847.210.9;3.253.212;3.155.59.7;4.636.17.9;2.424.814;7.233.17.6;6.747.48.5;5.454.111.3;3.936.912.7;4.558.812.3;3.527.89.8;4.540.28.4;1.513.510.1;8.556.47.1;4.571.68.2;6.552.810.9;4.144.111.2;5.540.99.4]u0=[45010]';n=20;T2=n*(n-1)*(mean(x)'-u0)'*inv(19*cov(x)
2、)*(mean(x)'-u0),p=3;F=(n-p)*T2/((n-1)*p),p=1-fcdf(F,3,17)T2=9.7388F=2.9045p=0.0649在显著性水平0.05下，接受原假设。第二类错误的计算，用非中心的F分布计算，非中心的参数为p67页，中间的参数。ncfcdf(3.2,3,17,20*(mean(x)-[45010])*inv(cov(x))*(mean(x)-[45010])')ans=0.36218248472391例3.2.2x=[3.748.59.3;5.765.18;3.847.210.9;3.253.212;3.155.59.7;4.63
3、6.17.9;2.424.814;7.233.17.6;6.747.48.5;5.454.111.3;3.936.912.7;4.558.812.3;3.527.89.8;4.540.28.4;1.513.510.1;8.556.47.1;4.571.68.2;6.552.810.9;4.144.111.2;5.540.99.4]S=cov(x),[v,d]=eig(S),n=20;p=3;c2=(n-1)*p*3.2/(n*(n-p)),d123=diag(sqrt(d))*sqrt(c2)S=2.879410.0100-1.809110.0100199.7884-5.640
4、0-1.8091-5.64003.6277v=-0.81750.57370.05080.0249-0.05300.9983-0.5754-0.8173-0.0291d=1.30140004.5316000200.4625c2=0.5365d123=0.83561.559210.37032、联立置信区间由3.2.4式计算的T2区间为：x=[3.748.59.3;5.765.18;3.847.210.9;3.253.212;3.155.59.7;4.636.17.9;2.424.814;7.233.17.6;6.747.48.5;5.454.111.3;3.936.912.7;4.
5、558.812.3;3.527.89.8;4.540.28.4;1.513.510.1;8.556.47.1;4.571.68.2;6.552.810.9;4.144.111.2;5.540.99.4]S=cov(x),[v,d]=eig(S),n=20;p=3;c2=(n-1)*p*3.2/(n*(n-p)),d123=diag(sqrt(d))*sqrt(c2)[mean(x)'-sqrt((n-1)*p*3.2/(n-p))*sqrt(diag(S)/n)mean(x)'+sqrt((n-1)*p*3.2/(n-p))*sqrt(diag(S)/n)]ans=3.3971
6、5.882935.047255.75288.570011.3600按3.2.5计算的区间为：[mean(x)'-tinv(0.975,n-1)*sqrt(diag(S)/n)mean(x)'+tinv(0.975,n-1)*sqrt(diag(S)/n)]ans=3.84585.434238.784852.01529.073610.8564§3.3例3.3.1一、假定两总体方差相同x=[65352560;75502055;60453565;75404070;70303050;55403565;60453060;65402560;60503070;55553575]y=[5555
7、4065;50604570;45453575;50505070;55503075;60404560;65554575;50603580;40453065;45504570]n=10;m=10;p=4;mx=mean(x)',my=mean(y)',A1=(n-1)*cov(x),A2=(m-1)*cov(y),D2=(n+m-2)*(mx-my)'*inv(A1+A2)*(mx-my),T2=n*m*D2/(m+n),F=(n+m-p-1)*T2/((n+m-2)*p),p=1-fcdf(F
