lr.mod <- glm( y ~ x + w, family=binomial ) exp( summary( lr.mod )$coef[,1:2] %*% rbind( c(1,1,1), 1.96*c(0,-1,-1) ) )
should do the job. Pack it in a function if you like, see e.g. the (so far) undocumented function ci.lin in: http://www.biostat.ku.dk/~bxc/R/ci.lin.R (depends on) http://www.biostat.ku.dk/~bxc/R/ci.mat.R Bendix ---------------------- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED] www.biostat.ku.dk/~bxc ---------------------- > -----Original Message----- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of > Cristian Pattaro > Sent: Wednesday, March 17, 2004 11:50 AM > To: [EMAIL PROTECTED] > Subject: [R] Q: Odds Ratio > > > Dear all, > > is there a automatic method to obtain Odds Ratio estimates and their > confidence intervals from a GLM Logistic model? > > Thanks > Cristian > > ============================================= > Cristian Pattaro > ============================================= > Unit of Epidemiology & Medical Statistics > Department of Medicines and Public Health > University of Verona > [EMAIL PROTECTED] > ============================================= > > > [[alternative HTML version deleted]] > > ______________________________________________ > [EMAIL PROTECTED] mailing list > https://www.stat.math.ethz.ch/mailman/listinfo> /r-help > PLEASE > do read the posting guide! > http://www.R-project.org/posting-guide.html > ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
