Hello, I was wondering how to obtain OR estimates and confidence intervals for a fixed effects model I am fitting with GLIMMIX (GLM with logit link)- so basically a logistic regression. I am able to obtain predicted probabilities at each level of the interaction with the lsmeans statement: lsmeans var_a*var_b/ilink cl, and realize that I can probably exponentiate the parameter estimates for each interaction term to obtain odds ratios (provided I've set my referent group for both var a and var b), but the confidence intervals are tricky. Can I use a command or option in the lsmeans statement to automate this and make it easier? Thanks!
... View more