作者celestialgod (攸藍)
看板R_Language
標題Re: [問題] 迴歸參數的擷取
時間Wed Aug 7 14:07:08 2013
# 我造的範例
n = 4; p = 13
logit_trans = function(x) exp(x) / (1+exp(x))
data = matrix(rnorm(n*p, 0, 0.1), n)
data = data.frame(y = sapply(logit_trans(cbind(1, data) %*% 1:(p+1)),
function(p) rbinom(1,1,p)), data)
comb = t(combn(4:12,2))
beta = mapply(function(i, j) matrix(coef(glm( data[,1] ~ data[,2] + data[,3] +
data[,i] + data[,j] + data[,i]:data[,j], family=binomial(link="logit"),
data = data))), comb[,1], comb[,2])
beta = cbind(comb, t(beta))
colnames(beta) <- c("i", "j","Intercept", "Gender", "Age", "CAR1", "CAR2",
"CAR1*CAR2")
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.152.213
※ 編輯: celestialgod 來自: 140.116.152.213 (08/07 14:10)
1F：推 diplazium:雖然不知為何會有NA(singularities是啥?),不過如果只是 08/07 17:44
2F：→ diplazium:要抽出參數值的話,就直接coef(mylogit)即可~ 08/07 17:45
3F：→ celestialgod:singularities就是n<p的情況 造成無法估計 08/07 18:48
4F：→ celestialgod:原po想要完整的TABLE 08/07 18:48
5F：→ celestialgod:我查看summary.glm的code 實在無法解決QQ 08/07 18:49
6F：推 andrew43:mylogit$coefficients 08/07 23:32
※ 編輯: celestialgod 來自: 61.227.245.110 (08/08 13:08)
※ 編輯: celestialgod 來自: 61.227.245.110 (08/08 13:09)