[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