I’d like to fit a mixed model with a bilinear random effect as discussed in (Hoff, P.D. Bilinear mixed-effects models for dyadic data. Journal of the American Statistical Association, 100: 286-295, 2005). In a “normal” mixed modeling situation if we have two random effects, say A and B, then the code random A B; assumes that A~Normal(0,Sigma_A) and B~Normal(0,Sigma_B). However, in the bilinear case the random effect is C = A*B or A T B (if A and B are vectors). Is there a way to incorporate the bilinear effect “C” into any of the mixed procedures in SAS? Thanks.
... View more