I have a problem with the proc GLIMMIX . I use it for fit a mixed model Y=XB + Zv +e where XB is the fixed effects and Zv is the random effect .
The vectors v and e follow a Gaussian law . the covariance matrix of v is G and the covariance matrix of e is R .
I want to fix the matrix R than i know.
I know i can use " PARMS " with the option "hold=" but my matrix R is 300x300.
I would like set R from a data.
Can you help me to find a solution ?