R 单因素协方差分析

单因素协方差分析（ANCOVA）扩展了单因素方差分析（ANOVA），包含一个或多个定量的
协变量。下面的例子来自于multcomp包中的litter数据集（见Westfall et al.，1999）。怀孕小鼠
被分为四个小组，每个小组接受不同剂量（0、5、50或500）的药物处理。产下幼崽的体重均值
为因变量，怀孕时间为协变量。分析代码见代码清单

> data(litter, package="multcomp") 
> attach(litter) 
> table(dose) 
dose 
 0 5 50 500 
 20 19 18 17 
> aggregate(weight, by=list(dose), FUN=mean) 
 Group.1 x 
1 0 32.3 
2 5 29.3 
3 50 29.9 
4 500 29.6 
> fit <- aov(weight ~ gesttime + dose) 
> summary(fit) 
 Df Sum Sq Mean Sq F value Pr(>F) 
gesttime 1 134.30 134.30 8.0493 0.005971 ** 
dose 3 137.12 45.71 2.7394 0.049883 * 
Residuals 69 1151.27 16.69 
--- 
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

利用table()函数，可以看到每种剂量下所产的幼崽数并不相同：0剂量时（未用药）产崽
20个，500剂量时产崽17个。再用aggregate()函数获得各组均值&#