從aov（）提取殘差

Question

我使用以下代碼運行了方差分析：

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

如果使用view(aov2) ，則可以在aov2的結構中看到殘差，但是我希望以不涉及剪切和粘貼的方式提取殘差。 有人可以幫我解決語法問題嗎？

我一直使用的residuals(aov2)各種版本residuals(aov2)僅產生NULL

Answer 1

無法從此模型中提取殘差的原因是，由於鳥鹽比率（???），您指定了隨機效果。 在這里，鳥和鹽的每種唯一組合都被當作具有唯一截距值但具有與鹽和進食量的單位差相關的共同加性效應的隨機簇來對待。

我想不出為什么要在此模型中將這個值指定為隨機效果。 但是為了合理地分析殘差，您可能需要根據擬合模型和最佳截距來計算每個層的擬合差。 我認為這是一件乏味的工作，但信息量不大。

Answer 2

我剛剛得知您可以使用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']