No code, only some thoughts at first. These might make it possible to address what you are looking for. Possession is some sort of factor, what it looks like to me is a random effect. Unfortunately, that doesn't seem to help much. What would happen if you considered the subject to be subj*possession, and possession as a random effect. I think you would get an analysis that is relatively complete then if you tried:
proc glimmix data= ic=q;
class subj possession tmt prime;
model location=prime consistency consistency*prime tmt tmt*prime extrinsicness extrinsicness*tmt/ddfm=kr(firstorder) s;
random possession/subject=subj ;
random tmt/residual subject=subj*possession type=chol;
run;
Note that I really don't have a good idea for the type of covariance structure for the first random statement. If you have the computing power available, the obvious choice would be unstructured, using type=chol.
I really hope you get more answers on this one, as I feel this is just dipping one toe into the pond.
SteveDenham
... View more
