Hi all, Say you have two continuous variables, size and number of eggs laid. These are both measured once for each individual in an experiment. Each individual belongs to a family. In addition, individuals in a family are reared on one of two food treatments (each individual is given only one of these treatments for the duration of the experiment). Family is treated as a random effect and food as a fixed effect. I would like to know how food treatment changes the covariance parameter estimates for the two variables, with family as the subject of the covariance parameter estimates. To produce the covariance parameter estimates for the two traits whilst ignoring food treatment my code is as follows. A dummy variable, ‘variable’, details if a measurement is for size or number of eggs laid, with column y as the value of that variable. ‘number’ is the specific ID given to each individual. Proc mixed data=mylib.mydata; Class family variable number; y = variable; random variable /type=un sub=family g gcorr; repeated /type=un sub=family*number; run; Any suggestions would be greatly appreciated.
... View more