Page 1 of 1

GAMLj mixed model logistic regression

Posted: Tue Mar 31, 2020 4:09 pm
by iainjgallagher
Hi

I wonder if someone could help me with the interpretation of a logistic regression mixed model.

I am trying to model response (either high or low) within control or intervention groups at two timepoints (before or after intervention or placebo). The omv file is attached. I'm using a mixed model just to account for paired observations (i.e. I only have a random intercept assuming the slope would be the same for everyone and the same across the groups).

My problem is in interpreting the Estimated Marginal Means i.e. the probability of answering low.

In the GAMLj output for example the probability of answering low in the control group is 1 (first Est Marginal Means table).

The actual numbers are ctrl:low = 14; not ctrl:low = 60. So P(low|ctrl) = 14/74 = 0.189

The df and CI's are also rather odd.

I'm rather new to mixed models so if someone could gently talk me through the interpretation (or more likely the mistake) here I can probably work out the rest myself.

Thanks.

Re: GAMLj mixed model logistic regression

Posted: Thu Apr 09, 2020 10:19 am
by mcfanda@gmail.com
Hi,
your estimated probabilities are either almost 1 or almost 0. It is indeed a bit weird, but I've checked your models and I get the same results in R glmer() and emmeans(), so the results seems to be correct. On the side, I would not use "dummy" as the coding methods because with "dummy" you obtain simple effects rather than main effects. Use "simple" so you get the main effects that you probably need