data test; do time = 1 to 10; do mom = 1 to 50; x3 = rannor(123); do kid = 1 to 3; y = rannor(345); x1 = rannor(3221); x2 = rannor(999); output; end; end; end; run;proc means data=test; run;
proc mixed data=test; class mom kid time; model y = x1 x2 x3 / solution ; repeated time / subject = mom(kid) type=ar(1) ; run;