R ANOVA 模型拟合 单因素方差分析

单因素方差分析中，你感兴趣的是比较分类因子定义的两个或多个组别中的因变量均值。以
multcomp包中的cholesterol数据集为例（取自Westfall、Tobia、Rom、Hochberg，1999），50
个患者均接受降低胆固醇药物治疗（trt）五种疗法中的一种疗法。其中三种治疗条件使用药物
相同，分别是20mg一天一次（1time）、10mg一天两次（2times）和5mg一天四次（4times）。剩下
的两种方式（drugD和drugE）代表候选药物。哪种药物疗法降低胆固醇（响应变量）最多呢？分
析过程见代码清单。

> library(multcomp) 
> attach(cholesterol) 
> table(trt)
trt 
 1time 2times 4times drugD drugE 
 10 10 10 10 10 
> aggregate(response, by=list(trt), FUN=mean)
 Group.1 x 
1 1time 5.78 
2 2times 9.22 
3 4times 12.37 
4 drugD 15.36 
5 drugE 20.95 
> aggregate(response, by=list(trt), FUN=sd)
 Group.1 x 
1 1time 2.88 
2 2times 3.48 
3 4times 2.92 
4 drugD 3.45 
5 drugE 3.35 
> fit <- aov(response ~ trt) 
> summary(fit) 
 Df Sum Sq Mean Sq F value Pr(>F)
trt 4 1351 338 32.4 9.8e-13 *** 
Residuals 45 469 10