为什么提取残差时找不到因子名称？答案

【问题标题】：Why can't I find my factor names when I extract residuals?为什么提取残差时找不到因子名称？
【发布时间】：2014-02-01 18:49:59
【问题描述】：

我正在处理一些选举数据，试图将其按“州”和“选举”分开。我对州和年份进行了具有固定效应的回归（如下所示），得到了我的汇总数据，并一直在尝试使用 resid() 函数来提取残差。

m5 <- lm(demVote ~ state*year, data=presidentialElections)

plot(resid(m5) ~ fitted(m5))
resid.m5 <- resid(m5)

我认为上面的一切都很完美。但是，这就是我迷路的地方 - 如果我做 summary(resid.m5) （我把提取的残差放在哪里，或者我想），我似乎再也找不到我的因子名称了。如果我想查看每个州或每年的残差（例如，按州/年的平均值），那么如何使用 resid() 函数访问它？谢谢！

【问题讨论】：

残差的顺序应该与presidentialElections 中的state 和year 的顺序相匹配，因此根据您要对分组残差执行的操作，您可以执行以下操作：@987654325 @，或aggregate(resid.m5, list(presidentialElections$state, presidentialElections$year), mean)。

【解决方案1】：

正如之前在 cmets 中所说，您必须意识到返回的残差与您在数据集中的观察顺序相同。

这是一个使用每个 R 安装附带的 iris 数据集的示例（以及可能非常荒谬的回归）：

data(iris)
m5 <- lm(Sepal.Length ~ Species*Sepal.Width, data=iris)
resid.m5 <- resid(m5)

dta.complete <- data.frame(iris, r.m5=resid.m5)

在这里，残差与原始数据相结合。这可能有点不正统，但为什么不把事情放在一起。现在您可以随心所欲地使用所有经典子集。例如：

with(dta.complete, by(r.m5, Species, mean))

祝你好运！

【讨论】：