线性回归函数lm中残差的自由度问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，有如下代码（原出处，薛毅老师的书）

1. > x<-c(0.10,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.20,0.21,0.23)
   
2. > y<-c(42.0,43.5,45.0,45.5,45.0,47.5,49.0,53.0,50.0,55.0,55.0,60.0)
   
3. > lm.sol<-lm(y ~ 1+x)
   
4. > summary(lm.sol)
   
5. 
   
6. Call:
   
7. lm(formula = y ~ 1 + x)
   
8. 
   
9. Residuals:
   
10. Min 1Q Median 3Q Max
    
11. -2.0431 -0.7056 0.1694 0.6633 2.2653
    
12. 
    
13. Coefficients:
    
14. Estimate Std. Error t value Pr(>|t|)
    
15. (Intercept) 28.493 1.580 18.04 5.88e-09 ***
    
16. x 130.835 9.683 13.51 9.50e-08 ***
    
17. ---
    
18. Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
    
19. 
    
20. Residual standard error: 1.319 on 10 degrees of freedom
    
21. Multiple R-squared: 0.9481, Adjusted R-squared: 0.9429
    
22. F-statistic: 182.6 on 1 and 10 DF, p-value: 9.505e-08

复制代码

那么我们如何由线性回归的结果lm.sol直接给出来残差的自由度呢？当然后来他的书里面也提到了用

1. lm.sol$da.residual

复制代码

现在我有点不明白的是如何确定用

1. lm.sol$da.residual

复制代码

即可，从哪里可以看出来？
谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：线性回归 自由度 coefficients coefficient Intercept 自由度

说实话，你这里的中文实在不怎么通顺，我看了之后不太明白你到底要问什么。能不能重新表述一下？

names(lm.sol)

怪了怪 发表于 2014-12-17 21:57
说实话，你这里的中文实在不怎么通顺，我看了之后不太明白你到底要问什么。能不能重新表述一下？

谢谢你的回复，简单来说就一句话：
如何根据lm函数的结果直接给出残差平方和的自由度。

kk22boy 发表于 2014-12-17 22:47
names(lm.sol)

谢谢你，看到了！

kaifengedu 发表于 2014-12-18 10:01
谢谢你，看到了！

残差的自由度就是n-k啊。你的例子中n=12,k=2（截距项一个，x一个）。
我觉得这不需要算啊，看看就知道的。

怪了怪 发表于 2014-12-18 11:06
残差的自由度就是n-k啊。你的例子中n=12,k=2（截距项一个，x一个）。
我觉得这不需要算啊，看看就知道的 ...

您说得对，不需要计算，我是想看下是否可以从lm的计算结果看出来。

也可用str(lm.sol)查看该对象的结构。

kaifengedu 发表于 2014-12-18 11:55
您说得对，不需要计算，我是想看下是否可以从lm的计算结果看出来。

如果是这样，如果你用Rstudio， 可以从右上角的Environment窗口里看到这里提到的全部信息，然后就可以调用。点开回归方程的名称就能看到。

怪了怪 发表于 2014-12-18 21:57
如果是这样，如果你用Rstudio， 可以从右上角的Environment窗口里看到这里提到的全部信息，然后就可以调用 ...

多谢多谢！