R 双因素方差分析

在双因素方差分析中，受试者被分配到两因子的交叉类别组中。以基础安装中的
ToothGrowth数据集为例，随机分配60只豚鼠，分别采用两种喂食方法（橙汁或维生素C），各
喂食方法中抗坏血酸含量有三种水平（0.5mg、1mg或2mg），每种处理方式组合都被分配10只豚
鼠。牙齿长度为因变量

> attach(ToothGrowth) 
> table(supp, dose) 
 dose 
supp 0.5 1 2 
 OJ 10 10 10 
 VC 10 10 10 
> aggregate(len, by=list(supp, dose), FUN=mean) 
 Group.1 Group.2 x 
1 OJ 0.5 13.23 
2 VC 0.5 7.98 
3 OJ 1.0 22.70 
4 VC 1.0 16.77 
5 OJ 2.0 26.06 
6 VC 2.0 26.14 
> aggregate(len, by=list(supp, dose), FUN=sd) 
 Group.1 Group.2 x 
1 OJ 0.5 4.46 
2 VC 0.5 2.75 
3 OJ 1.0 3.91 
4 VC 1.0 2.52 
5 OJ 2.0 2.66 
6 VC 2.0 4.80 
> dose <- factor(dose) 
> fit <- aov(len ~ supp*dose) 
> summary(fit) 
 Df Sum Sq Mean Sq F value