作者microcat (微喵)
看板R_Language
标题[问题] 回归参数的撷取
时间Tue Aug 6 12:16:11 2013
[问题类型]:
程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
请把以下不需要的部份删除
使用者
[问题叙述]:
请简略描述你所要做的事情，或是这个程式的目的
执行回归分析，产生下列表单，想撷取第一栏位（Estimate）
因为有出现 NA 的情况，无法完整的将第一栏位撷取出来！
（有很多组需要执行，每组所估算出的 Beta 值长度不一）
希望能保留完整六个 Beta 值
（Intercept,data[,2],data[,3],data[,i],data[,j],data[,i]:data[,j]）
想请教版上先进，可以如何达到这的目标，谢谢各位
Coefficients: (2 not defined because of singularities)
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.779222 1.309672 -0.595 0.552
data[, 2] 1.346329 0.204052 6.598 4.17e-11 ***
data[, 3] 0.035580 0.005365 6.631 3.32e-11 ***
data[, i] -0.771758 0.645828 -1.195 0.232
data[, j] NA NA NA NA
data[, i]:data[, j] NA NA NA NA
[程式范例]:
mylogit <- glm( data[,1] ~ data[,2] + data[,3] + data[,i] +
data[,j] + data[,i]:data[,j],
family=binomial(link="logit"), data = data)
coef <- summary(mylogit)$coefficients
coef[,1]
[关键字]:
感谢各位协助
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.129.5
※ 编辑: microcat 来自: 140.112.129.5 (08/06 14:04)
1F：推 Wush978:能不能提供范例data档案呢？否则无法重现错误 08/06 23:14
2F：→ Wush978:可以参考dump指令 08/06 23:15
3F：→ diplazium:dump...又从版主大学到一招了,谢谢! 08/07 10:22
4F：推 lin15:mylogit$coef 好像可以XD? 08/07 20:38
5F：推 celestialgod:可是 那两行NA出不来XDDD 08/07 22:24
6F：推 lin15:不加summary直接取coef的话就会有NA XD 08/08 01:33
7F：→ celestialgod:原PO意思好像是要完整的表? 08/08 01:47
8F：推 lin15:也可以拉 ifelse(is.na(a$coef),NA,summary(a$coef)) 08/08 02:18
9F：→ lin15:上面好像括号有放错地方... 08/08 02:18
10F：→ lin15:反正就是透过mylogit$coef来判断就可以了! 08/08 02:20
11F：→ microcat:谢谢各位的协助，mylogit$coef 可以达到所需要的目的 08/08 10:37
附上修改後的程式，谢谢各位的协助
Beta_temp <- matrix(0,1,8)
colnames(Beta_temp) <- c("i", "j","Intercept", "Gender", "Age", "CAR1",
"CAR2", "CAR1*CAR2")
for(i in 4:11){
for(j in (i+1):12){
mylogit <- glm( data[,1] ~ data[,2] + data[,3] + data[,i] + data[,j]
+ data[,i]:data[,j], family=binomial(link="logit"),
data = data)
coef <- t(as.matrix(mylogit$coef))
Beta_temp <- rbind(Beta_temp, c(i, j, coef[,1], coef[,2], coef[,3],
coef[,4], coef[,5], coef[,6]))
}
}
Beta_temp[-1,]
※ 编辑: microcat 来自: 140.112.129.5 (08/08 10:40)
12F：→ celestialgod:我按照原PO意思 把程式改好 修改於下面了 08/08 13:09
13F：→ microcat:非常感谢 08/09 14:07