方差分析实验报告

《方差分析实验报告》由会员分享，可在线阅读，更多相关《方差分析实验报告（15页珍藏版）》请在金锄头文库上搜索。

1、数据分析方法课程实验报告实验名称 方差分析 一、实验目的(1)通过 SAS程序针对相关问题做单因素方差分析以及两因素方差分析;(2)熟悉 SAS中的 PROC ANOVA过程，通过此过程建立方差分析表，对因素进行显著性分析以及各水平均值以及两两均值之差的比较，并能分析结果；(3)解决两因素方差分析问题，并作出与单因素相同的分析；(4)解决课后的相关问题二、用文字或图表记录实验过程和结果3.5 解:程序为：DATA THREE_5;INPUT A $ Y;CARDS;a1 7.6a1 8.2a1 6.8a1 5.8a1 6.9a1 6.6a1 6.3a1 7.7a1 6.0a2 6.7a2 8.

2、1a2 9.4a2 8.6a2 7.8a2 7.7a2 8.9a2 7.9a2 8.3a2 8.7a2 7.1a2 8.4a3 8.5a3 9.7a3 10.1a3 7.8a3 9.6a3 9.5;RUN;PROC ANOVA DATA=THREE_5;CLASS A;MODEL Y=A;MEANS A/ BON CLM;MEANS A/ BON CLDIFF;RUN;由题可知，此题为单因素方差分析问题，其中 =9； =12； =6；所以有 = + + =27；运1n23n1n23用 proc anova 过程可以建立方差分析表，如下：Sum ofSource DF Squares Mean

3、Square F Value Pr FModel 2 20.12518519 10.06259259 15.72 Simultaneous 95%A N Mean Confidence Limitsa3 6 9.2000 8.3594 10.0406a2 12 8.1333 7.5389 8.7277a1 9 6.8778 6.1914 7.5641通过中数据可得： =6.8778；同理 =8.1333； =9.2000， 并且得到 、 、1.y2.y3. LM的置信度为 95%的 Bonferroni置信区间分别为：（6.1914,7.5641） ， （7.5389,8.7277） ，H（8

4、.3594,10.0406） ；在通过 SAS程序可得：DifferenceA Between 95% ConfidenceComparison Means Limitsa3 - a2 1.0667 0.2410 1.8923 *a3 - a1 2.3222 1.4519 3.1925 *a2 - a3 -1.0667 -1.8923 -0.2410 *a2 - a1 1.2556 0.5274 1.9837 *a1 - a3 -2.3222 -3.1925 -1.4519 * a1 - a2 -1.2556 -1.9837 -0.5274 *DifferenceA Between Simul

5、taneous 95%Comparison Means Confidence Limitsa3 - a2 1.0667 0.0371 2.0962 *a3 - a1 2.3222 1.2370 3.4074 *a2 - a3 -1.0667 -2.0962 -0.0371 *a2 - a1 1.2556 0.3476 2.1635 *a1 - a3 -2.3222 -3.4074 -1.2370 *a1 - a2 -1.2556 -2.1635 -0.3476 *通过表中数据求的 的置信度为 95%的 Bonferroni同时置信区间为（-2.1635,-0.3476） ，LM的置信度为 95

6、%的 Bonferroni同时置信区间为（-3.4074,-1.2370 ） ， 的置信度为 95%的LH MHBonferroni同时置信区间为（-2.0962,-0.0371） ，从两两均值之差的置信度为 95%的同时置信区间可知： 显著小于 ， 显著小于 ， 显著小于 ，就是说，低档投入获得的生产能力得提高得均值显著小MLHMH于中档并且小于高档；同理，中档投入获得的生产能力得提高得均值显著小于高档，所以得到结论：过去三年科研经费投入越高，当年生产能力的改善越显著。3.6解：(1)（3）的程序如下：DATA THREE_6;INPUT A $ B $ Y;CARDS;.;RUN;PROC

7、 ANOVA DATA=THREE_6;CLASS A B;MODEL Y=A B A*B;MEANS A B A*B;RUN;DATA THREE_61;SET THREE_6;Y1=LOG(Y);RUN;PROC ANOVA DATA=THREE_61;CLASS A B;MODEL Y1=A B A*B;MEANS A B A*B;RUN;PROC ANOVA DATA=THREE_61;CLASS A B;MODEL Y1=A B;MEANS A B/T CLM;RUN;（4）的程序如下：DATA THREE_6;INPUT B $ Y;CARDS;b1 0.71 b1 1.66 b1

8、 2.01 b1 2.16 b1 2.42 b1 2.42b1 2.56 b1 2.60 b1 3.31 b1 3.64 b1 3.74 b1 3.74b1 4.39 b1 4.50 b1 5.07 b1 5.26 b1 8.15 b1 8.24b2 2.20 b2 2.93 b2 3.08 b2 3.49 b2 4.11 b2 4.95b2 5.16 b2 5.54 b2 5.68 b2 6.25 b2 7.25 b2 7.90b2 8.85 b2 11.96 b2 15.54 b2 15.89 b2 18.30 b2 18.59b3 2.25 b3 3.93 b3 5.08 b3 5.82

9、 b3 5.84 b3 6.89b3 8.50 b3 8.56 b3 9.44 b3 10.52 b3 13.46 b3 13.57b3 14.76 b3 16.41 b3 16.96 b3 17.56 b3 22.82 b3 29.13b1 2.20 b1 2.69 b1 3.54 b1 3.75 b1 3.83 b1 4.08b1 4.27 b1 4.53 b1 5.32 b1 6.18 b1 6.22 b1 6.33b1 6.97 b1 6.97 b1 7.52 b1 8.36 b1 11.65 b1 12.45b2 4.04 b2 4.16 b2 4.42 b2 4.93 b2 5.4

10、9 b2 5.77b2 5.86 b2 6.28 b2 6.97 b2 7.06 b2 7.78 b2 9.23b2 9.34 b2 9.91 b2 13.46 b2 18.40 b2 23.89 b2 26.39b3 2.71 b3 5.43 b3 6.38 b3 6.38 b3 8.32 b3 9.04b3 9.56 b3 10.01 b3 10.08 b3 10.62 b3 13.80 b3 15.99b3 17.90 b3 18.25 b3 19.32 b3 19.87 b3 21.60 b3 22.25;RUN;DATA THREE_61;SET THREE_6;Y1=LOG(Y);

11、RUN;PROC ANOVA DATA=THREE_61;CLASS B;MODEL Y1=B;MEANS B/BON CLM;MEANS B/BON CLDIFF;RUN;DATA THREE_6;INPUT A $ Y;CARDS;a1 0.71 a1 1.66 a1 2.01 a1 2.16 a1 2.42 a1 2.42a1 2.56 a1 2.60 a1 3.31 a1 3.64 a1 3.74 a1 3.74a1 4.39 a1 4.50 a1 5.07 a1 5.26 a1 8.15 a1 8.24a1 2.20 a1 2.93 a1 3.08 a1 3.49 a1 4.11 a

12、1 4.95a1 5.16 a1 5.54 a1 5.68 a1 6.25 a1 7.25 a1 7.90a1 8.85 a1 11.96 a1 15.54 a1 15.89 a1 18.30 a1 18.59a1 2.25 a1 3.93 a1 5.08 a1 5.82 a1 5.84 a1 6.89a1 8.50 a1 8.56 a1 9.44 a1 10.52 a1 13.46 a1 13.57a1 14.76 a1 16.41 a1 16.96 a1 17.56 a1 22.82 a1 29.13a2 2.20 a2 2.69 a2 3.54 a2 3.75 a2 3.83 a2 4.

13、08a2 4.27 a2 4.53 a2 5.32 a2 6.18 a2 6.22 a2 6.33a2 6.97 a2 6.97 a2 7.52 a2 8.36 a2 11.65 a2 12.45a2 4.04 a2 4.16 a2 4.42 a2 4.93 a2 5.49 a2 5.77a2 5.86 a2 6.28 a2 6.97 a2 7.06 a2 7.78 a2 9.23a2 9.34 a2 9.91 a2 13.46 a2 18.40 a2 23.89 a2 26.39a2 2.71 a2 5.43 a2 6.38 a2 6.38 a2 8.32 a2 9.04a2 9.56 a2 10.01 a2 10.08 a2 10.62 a2 13.80 a2 15.99a2 17.90 a2 18.25 a2 19.32 a2 19.87 a2 21.60 a2 22.25;RUN;DATA THREE_61;SET THREE_6;Y1=LOG(Y);RUN;PROC ANOVA DATA=THREE_61;CLASS A;MODEL Y1=A;MEANS A/BON CLM;MEANS A/BON CLDIFF;RUN;通过 PROC ANOVA过程输出了因素在个组合水平上得观测值得样本均值和标准差，如下表：The ANOVA ProcedureLevel of -