R 为什么可以';当我提取残差时,我不能找到我的因子名称吗?
我正在处理一些选举数据,试图将其按“州”和“选举”分开 我运行了一个对州和年份有固定影响的回归(如下所示),得到了我的汇总数据,并且一直在尝试使用resid()函数来提取残差R 为什么可以';当我提取残差时,我不能找到我的因子名称吗?,r,regression,R,Regression,我正在处理一些选举数据,试图将其按“州”和“选举”分开 我运行了一个对州和年份有固定影响的回归(如下所示),得到了我的汇总数据,并且一直在尝试使用resid()函数来提取残差 m5 <- lm(demVote ~ state*year, data=presidentialElections) plot(resid(m5) ~ fitted(m5)) resid.m5 <- resid(m5) m5正如前面的评论中所说,您必须意识到返回的残差与您在数据集中的观察值的顺序相同 下面是
m5 <- lm(demVote ~ state*year, data=presidentialElections)
plot(resid(m5) ~ fitted(m5))
resid.m5 <- resid(m5)
m5正如前面的评论中所说,您必须意识到返回的残差与您在数据集中的观察值的顺序相同
下面是一个使用iris数据集的示例,该数据集随每个R安装一起提供(可能是一个非常荒谬的回归):
祝你好运 残差的顺序应该与总统选举中的州
和年
的顺序相匹配,因此根据您想要对分组残差执行的操作,您可以执行以下操作:绘图(resid.m5~presidentialElections$state)
,或聚合(resid.m5,list(presidentialElections$state,presidentialElections$year),mean)
。如果缺少数据,则需要小心一点;由于剩余向量的长度不等于原始数据的行数,顺序已更改-可能使用m5$模型而不是原始df 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))