我做了一个简单的回归分析,并希望根据拟合值绘制残差,以查看它们的行为。之后我想检查残差是否正常分布。所以我做了:
summary(lm(Gesamt~PopTotal+HDI))
residuen<-summary(lm(Gesamt~PopTotal+HDI))$residuals
fitted<-summary(lm(Gesamt~PopTotal+HDI))$fitted.values
plot(fitted,residuen)
问题是,
summary(lm(Gesamt~PopTotal+HDI))$fitted.values
结果给出NULL
,这是否意味着没有拟合值?我想这是由于一些缺失值,但我不知道。如何计算残差而不是拟合值?
我的数据集可在此处找到:http://www.sendspace.com/file/8e27d0
答案 0 :(得分:4)
data=olympiadaten
来使问题更具可重现性
summary.lm()
,泛型summary()
在其参数为类lm
时调用的方法没有$fitted.values
,但lm()
没有。尝试将summary(lm(Gesamt~PopTotal+HDI))$fitted.values
更改为
lm(Gesamt~PopTotal+HDI, data=olympiadaten)$fitted.values
str()