Odds ratio & 95% CI from generalized estimating equations with robust standard error

[RRuser199]

I am working on logistic regression and was thinking of incorporating robust sandwich standard errors. I was able to get the coefficient but have trouble getting its associated odds ratio and 95% CI.

Here is my R code. All the feedback is very much appreciated.

model <- glm (disease ~ prescore, data = mydata, family=binomial) ## regular logistic regression
library(sandwich)
library(lmtest)
coeftest(model, vcov = sandwich) ## coefficient with sandwich standard errors

[williaml]

There is an example of the 95% CI using sandwich here:

https://data.library.virginia.edu/understanding-robust-standard-errors/

There is a general example of odds ratios for logistic regression here:

https://stats.idre.ucla.edu/r/dae/logit-regression/