[R] Difficulty getting standard deviation of ALL odds ratios with glm function, logistic regression, need cov of parameters

From: Garrett Fox <gfox_at_stat.rice.edu>
Date: Wed 27 Jul 2005 - 08:37:46 EST


I am trying to do logistic regression with a categorical predictor variable with the glm() function, family=binomial. Using glm() I would like to be able to calculate the confidence intervals of all three possible odds ratios for a factor (the factor has three categories). Three categories imply two columns of 0's and 1's in the design matrix, and two parameter estimates with their SE's. Two confidence intervals for odds ratios can be easily calculated, the third confidence interval can be calculated if I know the SE of B1-B2, but this requires the covariance matrix of parameter estimates, which is not given to my knowledge.

My initial thought was that Cov(B1,B2)=0, but this cannot be true or the variance estimates would depend on how you set up the design matrix (the factor symbolized by 0,0 would always have a higher SD).

Please help me find the covariance of the parameter estimates.

Thank you in advance,

Garrett Fox



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 27 08:42:52 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:34:01 EST