Question

我想知道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)的结果为我提供了我手动计算的残差（值-个人平均值-小组平均值+总体平均值）。

我的问题是，上面还有其他残差，为什么每个人（感觉上）都使用它们进行正态性检验？

我希望提供一些有用的信息。我已经研究过该功能，但找不到解释。

谢谢。

Answer 1

residual sum of square = total sum of square - Factor sum of squares

在您的情况下，因素是年龄。

残差应为正态分布，这是ANOVA的假设之一。

aov（）的残差如何计算？

1 个答案: