[問題] 迴歸參數的擷取
[問題類型]:
程式諮詢(我想用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)
推
08/06 23:14, , 1F
08/06 23:14, 1F
→
08/06 23:15, , 2F
08/06 23:15, 2F
→
08/07 10:22, , 3F
08/07 10:22, 3F
推
08/07 20:38, , 4F
08/07 20:38, 4F
推
08/07 22:24, , 5F
08/07 22:24, 5F
推
08/08 01:33, , 6F
08/08 01:33, 6F
→
08/08 01:47, , 7F
08/08 01:47, 7F
推
08/08 02:18, , 8F
08/08 02:18, 8F
→
08/08 02:18, , 9F
08/08 02:18, 9F
→
08/08 02:20, , 10F
08/08 02:20, 10F
→
08/08 10:37, , 11F
08/08 10:37, 11F
附上修改後的程式，謝謝各位的協助
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)
→
08/08 13:09, , 12F
08/08 13:09, 12F
→
08/09 14:07, , 13F
08/09 14:07, 13F
討論串 (同標題文章)