Question

我试图了解如何在R中正确运行重复测量或嵌套ANOVA，而不使用混合模型。从咨询教程，单变量重复测量anova的公式是：

aov(Y ~ IV+ Error(SUBJECT/IV) )

其中IV是受试者内部，受试者是受试者的身份。但是，大多数示例显示具有两个层次的输出：错误：主题和错误：主题：WS。同时我得到三个阶层（错误：主题和错误：主题：WS，错误：在内）。当我试图仅指定两个（内部和之间）时，为什么我有三个阶层？

这是一个可重复的例子：

data(beavers)
id = rep(c("beaver1","beaver2"),times=c(nrow(beaver1),nrow(beaver2)))
data = data.frame(id=id,rbind(beaver1,beaver2))
data$activ=factor(data$activ)
aov(temp~activ+Error(id/activ),data=data)

温度是温度的连续测量值，id是海狸激活的身份是活动的二元因子。模型的输出是：

Error: id
      Df Sum Sq Mean Sq
activ  1  28.74   28.74

Error: id:activ
           Df Sum Sq Mean Sq F value Pr(>F)
activ      1 15.313  15.313   18.51  0.145
Residuals  1  0.827   0.827               

Error: Within
           Df Sum Sq Mean Sq F value Pr(>F)
Residuals 210   7.85 0.03738

为什么aov（）中的Error（）给出了三个级别？

0 个答案: