统计专业实验实验.doc

《统计专业实验实验.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、重庆工商大学数学与统计学院《统计专业实验》课程实验报告实验课程：____指导教师：_专业班级：____学生姓名：___学生学号：____实验报告实验项目实验三多元正态总体检验实验日期实验地点实验目的1.掌握单一多元正态总体均值的检验；2.掌握两个多元正态总体均值向量的检验（区分协差阵是否相等）。3.掌握多元方差分析的思想和操作。实验内容1.检验2008年西部9个省区城镇居民大类消费与全国平均水平有无显著差异。2.分析我国上市公司电力、煤气及水生产供应行业和房地产行业在经营绩效（净资产收益率、总资产报酬率、资产负债率和总资产

2、周转率）方面是否存在明显差异，抽样数据见上市公司效绩指标.xls。3.一套生产线同时产出三种产品，分析温度和时间对总体产出率的影响，以及温度和时间对不同产品产出率的影响，数据见三种产品产出率.sav实验思考题解答：1．对协差阵未知时，在相等和不等(n=m)两种情况下，两个正态总体均值向量检验有何差异之处？答：先假设协差相同但未知，n=m和n不等于m的两种情况下，F值会有变化，但是不影响检验结果，最终还是F大于F临界值，拒绝原假设。再假设协差不等，在n=m和n不等于m的两种情况下，F值无变化，F值大于F临界值，所以仍是拒绝原

3、假设。2．对协差阵不等(n=m)的两个正态总体均值向量检验其基本思想如何，与单总体均值向量检验有何联系？答：建立。由于协差不等，且n=m，故检验统计量为，其中，Z=X-Y，S为Z的离差阵。给出检验水平，得到临界值，将计算得到的F值与F临界值比较，如果F大于F临界值，则拒绝原假设。这与单总体均值向量检验的步骤大致相同，只是计算公式，所得到的统计量不同，其基本思想是一致的。实验运行程序、基本步骤及运行结果：实验1、假设检验问题是：H0:西部9个省区城镇居民大类消费与全国平均水平有显著差异H1：西部9个省区城镇居民大类消费与全国

4、平均水平无显著差异首先在所给的表格中选取食品、衣着等8个总体指标，将其他不是西部的城市删掉，计算出西部的地区各个指标的平均值和西部个指标与全国水平的差值，就得到一张关于西部地区在这8个指标上的具体数据，然后用公式s=MMULT(TRANSPOSE(X-Z),X-Z)，就可以求得离差阵s，用公式s^-1==MINVERSE(S)便得到离差阵的逆矩阵，接着计算出T平方统计量，公式：T^2==n*(n-1)*MMULT(MMULT(Z_Z0,S_1),TRANSPOSE(Z_Z0))用公式将T平方统计量转化为F统计量F==(n-

5、p)/((n-1)*p)*G15最后通过查表得到F的临界值，F0.05（8,3）=8.8452>F=5.3618所以拒接原假设，即西部地区的消费与全国相比没有显著性差异。实验2n=13,m=15,n不等于m，设两组样本来自正态总体X和Y，且两组样本相互独立，协差阵相等但未知。假设检验问题：由于有共同未知协差阵，故检验统计量为其中：1分析资料对应的n=11，m=15，p=42利用函数AVERAGE计算X的各指标的平均值，即得到X的平均值向量转置矩阵=Mx,=；同理计算Y的平均值，既得Y的平均值向量转置矩阵=My；计算差值阵，

6、得：3计算X的离差阵，输入公式=MMULT(TRANSPOSE(X-Mx),X-Mx)，计算Y的离差阵，输入公式=MMULT(TRANSPOSE(Y-My),Y-My)计算4计算S的逆矩阵，输入公式=MINVERSE(S)5计算统计量T2=(n+m-2)*(n*m)/(n+m)*MMULT(MMULT(d,S_1),TRANSPOSE(d)，计算统计量F=((n+m-2)-p+1)/((n+m-2)*p)*I6；6计算临界值（取显著水平a=0.05）F0.05=FINV(0.05,p,n+m-p+1)。数据如下：7分析：由

7、上可知，F统计量大于F0.05临界值，所以拒绝原假设，即我国上市公司电力、煤气及水生产供应行业和房地产行业在经营绩效（净资产收益率、总资产报酬率、资产负债率和总资产周转率）方面存在明显差异。实验3：思路：提出零假设―选择检验统计量―计算检验统计量的观测值及概率p值―给出显著性水平做出假设，假设检验问题：H0：不同温度和时间对产品产出率没有显著影响；H1：不同度和时间对产品产出率有显著影响；（1）温度和时间对总体产出率的影响。操作：选择菜单Analyze－GeneralLinearModel－Univariate；再分别选择

8、“总产出率”到随机变量框中，选择“温度”“时间”到控制变量框中，得到结果如下：TestsofBetween-SubjectsEffectsDependentVariable:三种产品和SourceTypeIIISumofSquaresdfMeanSquareFSig.CorrectedModel1891