Мне интересно, как вычисляются невязки в 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) дает мне остатки, которые я вычислил вручную (значение - личное среднее - среднее по группе + общее среднее).
У меня вопрос: каковы остальные остатки выше и почему каждый (так он чувствует) использует их для проверки нормальности?
Я хотел бы получить полезную информацию. Я уже погрузился в эту функцию, но не смог найти объяснения.
Спасибо.