【问题标题】:Residuals but no fitted values from a regression model in RR中回归模型的残差但没有拟合值
【发布时间】:2012-09-28 07:37:21
【问题描述】:

我做了一个简单的回归分析,并希望根据拟合值绘制残差以查看它们的行为。之后我想检查残差是否是正态分布的。所以我做了:

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,这是否意味着没有拟合值?我想这是由于一些缺失值,但我不知道。 R如何计算残差而不是拟合值?

我的数据集可以在这里找到:http://www.sendspace.com/file/8e27d0

【问题讨论】:

  • 总结没有拟合值,模型有。试试m1 &lt;- lm(Gesamt~PopTotal+HDI) 然后装(m1) 或m1$fitted

标签: r regression


【解决方案1】:
  1. 您可以通过将data=olympiadaten 添加到您的代码中来使您的问题更具重现性
  2. summary.lm(),泛型summary() 在其参数属于lm 类时调用的方法,没有$fitted.values,但lm() 有。尝试将summary(lm(Gesamt~PopTotal+HDI))$fitted.values 更改为 lm(Gesamt~PopTotal+HDI, data=olympiadaten)$fitted.values
  3. 一般情况下,您可以通过str() 在 R 中检查对象的组成部分

【讨论】:

  • 好的,谢谢,所以没关系,如果我用 ks 检验检查残差是否为正态分布,对吗? @miura
  • 不,不建议这样做。您应该看一下残差与拟合图。试试plot(lm(Gesamt~PopTotal+HDI, data=olympiadaten)),这是四个情节中的第一个,有点令人担忧。
  • 第二个图显示了一个正常的残差 QQ 图,它处理了您在评论中提出的残差正态性问题。它看起来也不太好。这里有很多讨论回归诊断的线程,您可以查看。
猜你喜欢
  • 1970-01-01
  • 2018-12-07
  • 2018-10-13
  • 2017-10-03
  • 1970-01-01
  • 2019-05-15
  • 2017-03-17
  • 1970-01-01
  • 2018-04-09
相关资源
最近更新 更多