r - 从 aov() 中提取残差

Question

我使用以下代码运行了方差分析：

aov2 <- aov(amt.eaten ~ salt + Error(bird / salt),data)

如果我使用view(aov2)，我可以看到 aov2 结构中的残差，但我想以不涉及剪切和粘贴的方式提取它们。有人可以帮我解决语法吗？

residuals(aov2)我一直在使用的各种版本只生产NULL

Answer 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']         

Answer 2

您无法从该模型中提取残差的原因是您指定了由于鸟盐比 (???) 而产生的随机效应。在这里，鸟和盐的每个独特组合都被视为具有唯一截距值的随机簇，但具有与盐的单位差异和食用量相关的共同加性效应。

我无法理解为什么我们要在这个模型中指定这个值作为随机效应。但为了明智地分析残差，您可能需要根据拟合模型和最佳截距计算每个层中的拟合差异。然而，我认为这是一项乏味的工作，而且信息量不大。