SAS 方差分析 作业5.doc

实验七 方差分析
一、 作业题目：
用二氧化硅（SiO2，原著中译为二氧化矽）50mg对大鼠染尘后，不同时期全肺湿重的变化如表8-1。试比较染尘后1月、3月、6月，三个时期的全肺湿重有无差别（单因素三个水平，完全随机平衡设计的方差分析）？（杨树勤，卫生统计学，第2版，人民卫生出版社，1989，41）
表10-1 不同时期全肺湿重
1月
3月
6月
3.3
4.4
3.6
3.6
4.4
4.4
4.3
3.4
5.1
4.1
4.2
5.0
4.2
4.7
5.5
3.3
4.2
4.7

二、 程序代码：
Data shizhong;
Do rep = 1 to 6;
Do month = 1 to 3;
Input x @;
output;
End;
End;
Cards;
3.3 4.4 3.6
3.6 4.4 4.4
4.3 3.4 5.1
4.1 4.2 5.0
4.2 4.7 5.5
3.3 4.2 4.7
;
proc anova;
class month rep;
model x = month;
MEANS month/DUNCAN;
run;
三、 输出结果：
SAS 系统 09:29 Tuesday, April 16, 2002 1

The ANOVA Procedure

Class Level Information

Class Levels Values

month 3 1 2 3

rep 6 1 2 3 4 5 6

Number of observations 18

SAS 系统 09:29 Tuesday, April 16, 2002 2

The ANOVA Procedure

Dependent Variable: x

Sum of
Source DF Squares Mean Square F Value Pr > F

Model 2 2.52777778 1.26388889 4.52 0.0291

Error 15 4.19666667 0.27977778

Corrected Total 17 6.72444444

R-Square Coeff Var Root MSE x Mean

0.375909 12.46194 0.528940 4.244444

Source DF Anova SS Mean Square F Value Pr > F

month 2 2.52777778 1.26388889 4.52 0.0291

SAS 系统 09:29 Tuesday, April 16, 2002 3

The ANOVA Procedure

Duncan's Multiple Range Test for x

NOTE: This test controls the Type I comparisonwise error rate, not the experimentwise error rate.

Alpha 0.05
Error Degrees of Freedom 15
Error Mean Square 0.279778

Number of Means 2 3
Critical Range .6509 .6823

Means with the same letter are not significantly different.

Duncan Grouping Mean N month

A 4.7167 6 3
A
B A 4.2167 6 2
B
B 3.8000 6 1

四、结果分析
由输出结果知，方差分析模型是显著的，显著水平达0.0291，说明用二氧化硅50mg对大鼠染尘后，不同时期全肺湿重的变化是显著不同的。多重均值测验结果说明，第一种和第二种，第二种和第三种不同时期全肺湿重之间的差异不显著外，其余各种方法之间的差异均显著，显著水平为0.05。