Question

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

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

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

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

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))))

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