08-15-2011 10:33 AM
I am using proc glimmix for logistic regression of longitudinal data with dichotonomous outcomes.
The code i am using is this.
proc glimmix data = t1 method = quad;
class Id month feed;
model case = feed month feed*month / s dist = binomial link = logit oddsratio;
random intercept / subject = Id type = un;
The exp(estimates) does not match the odds ratio in the odds ratio table. I wan to know if I am interpreting the tables correctly or am I missing something. Any help in this matter is greatly appreciated. Thx, N
07-24-2012 08:28 PM
I am having exactly the same issue in a similar model. Here's what it says in the users guide under GLIMMIX proc
I did some tests to see if I could figure out what was going on with my model. When I include only the types of variables described above, the exp(estimate) does in fact equal the odds ratio displayed in the OR table. However, as soon as I enter additional variables into my model (for which GLIMMIX can't calculate OR's) , the exp(estimate) and odds ratios no longer match.
Is it possible that the odds ratio table calculates and displays OR's as if only variables it can use (e.g. those listed above) are in the model, rather than displaying correctly adjusted OR's calculated for the complete model (including the interaction terms etc that it can't calculate OR's for.)? If so, the output would be extremely misleading.
Anyone else out there that can shed light on this issue? Which results should we trust - the calculated OR's or the estimated parameters?