我想知道aov()中的残差是如何计算的。我已经找了好几个小时了,但还是弄不明白。
我使用ANOVA进行重复测量。
Data <- data.frame(subject = factor(rep(1:10, 3)),
age = factor(c(rep(4, 10),
rep(10, 10),
rep(35, 10))),
weight = c(20, 9, 16, 14, 30, 26, 26, 27, 13, 15,
27, 18, 30, 26, 43, 48, 38, 38, 22, 47,
50, 44, 52, 46, 64, 70, 73, 57, 54, 63))
ANOVA_MW <- aov(weight ~ age +
Error(subject / age),
data = Data)
summary(ANOVA_MW)
我知道下面的命令给了我一些东西。round(ANOVA_MW$subject:age$residuals, 2)
但是,我只得到20个值,而不是30个值。它从11开始。这可能与主体的残差有关。我不知道。
proj(ANOVA_MW)
的结果给出了我手动计算的残差(值-个人平均值-组平均值+总体平均值)。
我的问题是,上面的其他残差是什么,为什么每个人(所以它感觉)都在使用它们进行正态性测试?
我想要一些有用的意见。我已经深入研究了这个函数,但是找不到一个解释。
谢谢。
发布于 2019-09-29 18:14:46
residual sum of square = total sum of square - Factor sum of squares
在你的例子中,因素是年龄。
残差应该是正态分布的,这是ANOVA的假设之一。
https://stackoverflow.com/questions/54065304
复制