在双因素方差分析中,受试者被分配到两因子的交叉类别组中。以基础安装中的
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