我想知道如何计算aov()中的残差。我已经看了几个小时,但无法弄清楚。
我使用方差分析进行重复测量。
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)
的结果为我提供了我手动计算的残差(值-个人均值-组均值+总体均值)。我的问题是,上面还有其他残差吗?为什么每个人(感觉上)都使用它们进行正态性检验?
我希望提供一些有用的意见。我已经研究过该功能,但找不到解释。
谢谢。
最佳答案
residual sum of square = total sum of square - Factor sum of squares
在您的情况下,因素是年龄。
残差应为正态分布,这是ANOVA的假设之一。
关于r - aov()的残差如何计算?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/54065304/