r - 按组进行简单回归并显示输出

Question

我想问一些关于 R 中回归可能性的一般性问题。

例如，我有 58 个区域的两个变量之间的数据。我想进行整个回归过程，包括每个区域的假设检查、模型拟合和诊断，但通过一个命令获得整体结果，这意味着没有循环。

我已经知道我可以使用 lmList 函数在一次试验中进行模型拟合。但是，我不知道是否有可能一次性获得所有 58 个回归的 QQ 正态残差图。

有谁知道这是否可行？如果是这样，我可能需要什么样的功能？

Answer 1

取决于“一个命令”的含义，以及为什么要避免循环。怎么样：

library(nlme)
L <- lmList(y~x|region,data=yourData)
lapply(L,plot,which=2)

应该管用; 但是，它会依次吐出 58 个图。如果您尝试将它们全部捕获在一个页面上，您可能会收到关于页边距太小的错误。

根据lmList返回的回归列表，您还有很多其他选择。例如，

library(plyr)
qqDat <- ldply(L,function(x) as.data.frame(qqnorm(residuals(x))))

将为您提供一个数据框，其中包含数据中每个组的 QQ 图信息（预期值和观察值）。