I am comparing an offender group with a control group on a certain task in a generalized linear model with accuracy (0-1) as the DV. The reviewers criticized me by mentioning that one group was imprisoned and the other not. I now want to control for the length they are imprisoned but of course, that is also correlated with age. So, I am doubting about the best way to control for this. Length of imprisonment in the control group is zero (not missing). Let’s say I am interested in the interaction between
Group*A*B. Than the fixed predictors in the model are:
A, B, Group, A*B, A*Group, B*Group, A*B*Group.
Should the new model look like this:
Age, Imprisonment, A, B, Group, A*B, A*Group, B*Group, A*B*Group.
Or like this:
Age, Imprisonment, A, B, Group, A*B, A*Group, B*Group, A*B*Group A*Imprisonment, B*Imprisonment, Group*Imprisonment, A*B*Imprisonment, A*Group*Imprisonment, B*Group*Imprisonment, A*B*Group*Imprisonment.
I also included two random factors:
/RANDOM USE_INTERCEPT=TRUE SUBJECTS=id COVARIANCE_TYPE=VARIANCE_COMPONENTS /RANDOM EFFECTS=trial USE_INTERCEPT=FALSE COVARIANCE_TYPE=VARIANCE_COMPONENTS