Hello- I have run a few models with glimmix and there is something that has confused me and still does. I hope someone can help me figure out what is going on. Right now, I am working on some data whose outcome variable is binary. There are also repeated measurements. So if I run a 'relatively' simple model, say:
proc glimmix data=subset5;
class wave hhid;
model injury(event='1')= arthrit male wave wave*arthrit wave*male / s d=b;
random _residual_/ sub=hhid type=cs;
Where the outcome injury and the covariates arthrit male are coded as 1 or 0.
So I obtain my Type III tests of fixed effects. The p values for these are:
Effect DF DF F Value Pr > F
However, when I look at the solution, the p-values do not correspond at all. For instance, the p-values for simple effects such as arthrit and male are:
Effect wave Estimate Error DF t Value Pr > |t|
The p-values for male reported in the Parameter Estimates table represents a test of the effect of gender in the reference wave only. But the p-value for male reported in the Type III tests of fixed effects table represents a test of the gender effect across all waves.
To better understand this, try adding the following ESTIMATE statements to your GLIMMIX run:
estimate "Male effect in reference group"
male*wave 0 0 1;
estimate "Average male effect"
male*wave 1 1 1 / divisor=3;
The first ESTIMATE statement should give you the effect estimate and p-value returned in the Parameter Estimates table. The second ESTIMATE should give you the p-value returned in the Type III tests of fixed effects table.
Your interest is probably not to test the effect of gender in wave 3 only. Rather, you probably would like to test the effect of gender in the entire data set. That is the test given by the second ESTIMATE statement representing the "Average male effect" across all waves and also the test presented by the Type III tests of fixed effects table.