请教R怎么分组做glm回归

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请教一个问题啊，我想按照某个变量来分组，分别在每个组下面跑回归，请问程序应该怎么写呢？

model1 <- glm(m~g,family=binomial(link='logit'),data=data,subset(group))  这样好像不对，请问应该怎么写呢

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：model mode ODE Mod del

1. dat <- data.frame(group = rep(c('A', 'B'), each=100), x=rnorm(200), y=rnorm(200))
   
2. 
   
3. reg_by <- function(dat_in) {
   
4.     lm(y~x+1, data =dat_in)
   
5. }
   
6. 
   
7. by(dat[,2:3], dat[,1], reg_by)

复制代码

谢谢您！希望再次得到您的指点

还有，希望在模型中显示出or值、95%CI、P值，应该怎么修改程序呢

我自己修改了一下
用model来接收glm的结果，然后在function里面写入
exp(coef(model1))
exp(confint(model1))
但是要么得到or值，要么得到95%ci，不能在这个函数里面直接得到or值和95%ci，请问是怎么回事呢

分开写了两个function来解决这个问题了，等待更好的方法

请问可以分享程序吗

pingguzh 发表于 2017-8-22 16:40
分开写了两个function来解决这个问题了，等待更好的方法

请问大神可以分享一下怎么写的吗