我刚开始学习 R，需要一些帮助来找到我的数据的残差平均值和中位数。我计算了 lm 并在总结中得到如下残差：
min 1Q median 3Q Max
-111.86 -34.90 -7.6 33.46 182.58
问题：残差的中位数是 -7.6，但我的平均值是多少？ 或者是否有计算残差的均值和中位数？ 我打算做 mean(resid(trees.lm) 还是应该输入 mean(trees.lm$resid)
请澄清一下，因为我的同学们对同一个数据集都有不同的反应。
【问题讨论】：
这里一个具体问题的答案是：
mean(resid(trees.lm))
您不应该像这样深入研究合适的模型对象并去除任意组件。当你意识到你刚刚通过以下方式提取了工作残差时，在像 GLM 这样更复杂的东西上这样做会咬你的手：
glm.mod.obj$residuals
这对你不太可能有用。
即使对于像 lm() 对象这样的简单事物，使用 resid() 或访问 $residuals 的内容也可能会有所不同，具体取决于模型的拟合方式（例如，na.action 参数的设置是什么？）。
此外，线性模型假设残差是 i.i.d。具有均值 0 和方差 $\hat{\sigma}^2}$ 的高斯（或正态）随机变量，因此均值应该非常接近 0（即非常、非常、非常接近 0，但不完全是因为这是一台计算机和浮点运算正在发挥作用）。
【讨论】：
R 通常报告1e-14 或类似的平均值。这就是为什么我将我的评论设为“接近于零”而不是“完全为零”。