Re: [問題] 迴歸參數的擷取

[celestialgod (攸藍)]

# 我造的範例
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)
推 diplazium:雖然不知為何會有NA(singularities是啥?),不過如果只是
→ diplazium:要抽出參數值的話,就直接coef(mylogit)即可~
→ celestialgod:singularities就是n<p的情況 造成無法估計
→ celestialgod:原po想要完整的TABLE
→ celestialgod:我查看summary.glm的code 實在無法解決QQ
推 andrew43:mylogit$coefficients
※ 編輯: celestialgod    來自: 61.227.245.110       (08/08 13:08)
※ 編輯: celestialgod    來自: 61.227.245.110       (08/08 13:09)