[R] confint for glm (general linear model)

casperyc casperyc at hotmail.co.uk
Sun Dec 13 02:19:29 CET 2009


for an example,


counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9); treatment <- gl(3,3)
glm.D93 <- glm(counts ~ outcome + treatment, family=poisson())
confint(glm.D93)
confint.default(glm.D93)  # based on asymptotic normality

to verify the confidence interval (confint.default(glm.D93))  for outcome2

-4.542553e-01 + c(-1,1) * 0.2021708 * qt(0.975,df=4)
-1.0155714  0.1070608

does not give me 
outcome2    -0.8505027 -0.05800787
as in confint.default(glm.D93)

Thanks
-- 
View this message in context: http://n4.nabble.com/confint-for-glm-general-linear-model-tp954071p962790.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list