基礎的すぎるから知らなかったのか confint(ModelName, level = 1-alpha) という関数があったらしい、プロファイル尤度に基づく信頼区間を計算できる (2012/6/19 記述を修正)。 d <- c( 1,1,1,1,0,0,0,0,0,0, 1,1,1,1,1,1,1,1,0,0, 1,0,0,0,0,0,0,0,0,0, 1,1,1,1,1,0,0,0,0,0 ) trt <- gl(2, 20, 40) foo <- gl(2, 10, 40) GLM.1 <- glm(d ~ trt + foo, family=binomial()) summary(GLM.1) ModelName <- GLM.1 alpha <- 0.05 x <- summary(ModelName) y <- confint(ModelName, level=1-alpha)