Hi Jil & Mike,
I have created another variable "repA" through "repC" to capture the replications in each treatment cell, e.g., in the first cell I recorded "1, 1, 1, 1, 1" for the five replications in each of the rep variables. Then I had to define the "subject = " statement with repC, the random variable. Is this in order? the data set, code and output are as shown below:
proc mixed data=FFRgroup ; class factorA factorB factorC groupA groupB groupC repA repB repC; model yield= groupA | groupB / ddfm=satterthwaite; random int factorA factorB factorA*factorB/subject=repC v; lsmeans groupA | groupB / pdiff cl alpha=0.05 ; Title "FFR Model"; run;
Looks like the V matrix is still a diagonal matrix by default in PROC mixed, as highlighted by Mike.
... View more