求助：如何在R中获得回归系数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如题，求问如何得到coefficients里某一解释变量的回归系数？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：回归系数 coefficients coefficient EFFICIENT 解释变量 矩阵 系数 如何

coefficients
是什么

clzu 发表于 2012-10-30 20:23
coefficients
是什么

就是summary里的              
             Estimate  Std. Error  t value     Pr(>|t|)
(Intercept) 13.7799267 0.283330374 48.63554 1.988756e-70
X            0.6757055 0.008076221 83.66605 5.923357e-93
这个，我想用语句去取得X的回归系数就是0.6757055……

merry0312 发表于 2012-10-30 20:29
就是summary里的              
             Estimate  Std. Error  t value     Pr(>|t|)
(Intercept ...

1. ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
   
2. trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
   
3. group <- gl(2,10,20, labels=c("Ctl","Trt"))
   
4. weight <- c(ctl, trt)
   
5. lm.D9 <- lm(weight ~ group)
   
6. lm.D9$coeff
   
7. (Intercept) groupTrt
   
8. 5.032 -0.371
   
9. lm.D9$coeff[2]
   
10. groupTrt
    
11. -0.371

复制代码

data(Journals,package="AER")
reg<-lm(log(subs)~log(price/citations),data=Journals)
coef(reg)[2]

1. data(Journals,package="AER")
   
2. reg<-lm(log(subs)~log(price/citations),data=Journals)
   
3. coef(reg)[2]

复制代码

其实 lm() 进行的回归运算，用符号 “<-” 将其赋给一个对象，楼上的例子中就是reg。
这个reg属于一个lm类（class），其实类就是数据在R中的表现形式，而对应的模式（mode）就是数据在内存中记录形式。可是试试下面两个命令看看结果：

1. class(reg)
   
2. mode(reg)

复制代码

关键是mode命令返回“list”这一结果，那么就是对reg进行list操作了，即提取楼主想要的各种回归结果数据。使用下面的命令看看reg里面有什么：

1. ls(reg)

复制代码

这个时候就可以用list提取数据的操作符“$”来提取数据了，例如提取系数就可以用：

1. reg$coefficients

复制代码

这时返回一个数字向量，包括了所有的系数，在用下标操作符“[ ]”就得到想要的结果了：

1. reg$coefficients[2]

复制代码

建议楼主研究一下R的数据类型和结构，以及操作方法，可以为使用R打个结实点儿的基础。
推荐《R in a Nutshell》

xwlove 发表于 2012-10-30 22:10

好使！感谢！