单因素方差分析中,你感兴趣的是比较分类因子定义的两个或多个组别中的因变量均值。以
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