我使用以下代码运行anova：
aov2 <- aov(amt.eaten ~ salt + Error(bird / salt),data)
如果我使用view(aov2)我可以看到aov2结构中的残差，但我想以不涉及剪切和粘贴的方式提取它们。有人可以用语法帮助我吗？
我一直使用的residuals(aov2)的各种版本只生成NULL
答案 0 :(得分：1)
由于鸟盐比（???），您无法从此模型中提取残差的原因是您指定了随机效应。在这里，鸟和盐的每个独特组合被视为具有独特截距值的随机簇，但是与盐的单位差异和食用量相关联的常见加性效应。
我无法想象为什么我们要将此值指定为此模型中的随机效果。但是为了明智地分析残差，您可能希望根据拟合模型和最优截距计算每个层中的拟合差异。不过，我认为这是繁琐的工作，而且信息量不大。
答案 1 :(得分：1)
我了解到您可以使用proj：
x1 <- gl(8, 4)
block <- gl(2, 16)
y <- as.numeric(x1) + rnorm(length(x1))
d <- data.frame(block, x1, y)
m <- aov(y ~ x1 + Error(block), d)
m.pr <- proj(m)
m.pr[['Within']][,'Residuals']