I am a new R user (am using it through the Rcmdr package) and have struggled to find out how to report OR and RR directly when running GLM models (not only reporting coefficients.) Example of the syntax that I have used: GLM.2 <- glm(diarsev ~ treatmentarm +childage +breastfed, family=binomial(logit), data=fieldtrials2) summary(GLM.2) This works well except that I manually have to calculate the OR based on the coefficients. Can I get these directly (with confidence intervals) by just amending the syntax? Will be grateful for advice! -- View this message in context: http://r.789695.n4.nabble.com/Reporting-odds-ratios-or-risk-ratios-from-GLM-tp3357209p3357209.html Sent from the R help mailing list archive at Nabble.com.
OR <- exp(coef(GLM.2)[-1]) OR.ci <- exp(confint(GLM.2)[-1,]) -Aaron On Tue, Mar 15, 2011 at 1:25 PM, lafadnes <ubuntu@fadnes.net> wrote:> I am a new R user (am using it through the Rcmdr package) and have > struggled > to find out how to report OR and RR directly when running GLM models (not > only reporting coefficients.) > > Example of the syntax that I have used: > > GLM.2 <- glm(diarsev ~ treatmentarm +childage +breastfed, > family=binomial(logit), data=fieldtrials2) > summary(GLM.2) > > This works well except that I manually have to calculate the OR based on > the > coefficients. Can I get these directly (with confidence intervals) by just > amending the syntax? > > Will be grateful for advice! > > -- > View this message in context: > http://r.789695.n4.nabble.com/Reporting-odds-ratios-or-risk-ratios-from-GLM-tp3357209p3357209.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]
On Mar 15, 2011, at 1:25 PM, lafadnes wrote:> I am a new R user (am using it through the Rcmdr package) and have > struggled > to find out how to report OR and RR directly when running GLM models > (not > only reporting coefficients.) > > Example of the syntax that I have used: > > GLM.2 <- glm(diarsev ~ treatmentarm +childage +breastfed, > family=binomial(logit), data=fieldtrials2) > summary(GLM.2) > > This works well except that I manually have to calculate the OR > based on the > coefficients. Can I get these directly (with confidence intervals) > by just > amending the syntax?Look at predict.glm {stats}. It gives both facilities and the help page looks very clear. -- David Winsemius, MD West Hartford, CT
Thanks a lot! The following strategy as you suggested worked fine: OR <- exp(coef(GLM.2)[-1]) OR.ci <- exp(confint(GLM.2)[-1,]) OR OR.ci -- View this message in context: http://r.789695.n4.nabble.com/Reporting-odds-ratios-or-risk-ratios-from-GLM-tp3357209p3381020.html Sent from the R help mailing list archive at Nabble.com.