Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/82.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181

Warning: file_get_contents(/data/phpspider/zhask/data//catemap/0/xml/12.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 为什么可以';当我提取残差时,我不能找到我的因子名称吗?_R_Regression - Fatal编程技术网

R 为什么可以';当我提取残差时,我不能找到我的因子名称吗?

R 为什么可以';当我提取残差时,我不能找到我的因子名称吗?,r,regression,R,Regression,我正在处理一些选举数据,试图将其按“州”和“选举”分开 我运行了一个对州和年份有固定影响的回归(如下所示),得到了我的汇总数据,并且一直在尝试使用resid()函数来提取残差 m5 <- lm(demVote ~ state*year, data=presidentialElections) plot(resid(m5) ~ fitted(m5)) resid.m5 <- resid(m5) m5正如前面的评论中所说,您必须意识到返回的残差与您在数据集中的观察值的顺序相同 下面是

我正在处理一些选举数据,试图将其按“州”和“选举”分开 我运行了一个对州和年份有固定影响的回归(如下所示),得到了我的汇总数据,并且一直在尝试使用resid()函数来提取残差

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