st: proc mixed translation

[James Valcour <valcourj@uoguelph.ca>]

I've been trying to convert the following SAS code to Stata.  The 
experiment involved a latin square design as follows: for each sex(M/F) 
5 squares were set up with 3 rats/square and 3 periods/square 
(treatments random allocated to each cell).  The outcome was measured in 
triplicate for each cell of the square.

proc mixed noitprint noclprint covtest cl method=reml;
 class square trt rat period sex;
 model resp = trt sex trt*sex / outp=new;
 random square(sex) rat(square sex) period  period*trt*square(rat);
run;

I've been trying to run this through -xtmixed- but I'm not sure if it 
can handle how the analysis has been specified.  I'm particularly 
stumped over how the syntax for the "period*trt*square(rat)" random 
effect would be specified.

Any help is appreciated.

Thanks!


--

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
James Valcour
Ph.D. Candidate
Department of Population Medicine
Ontario Veterinary College
University of Guelph
Guelph, ON, CA
N1G 2W1
phone: (519)824-4120 x54595
fax: (519)763-8621
email: valcourj@uoguelph.ca
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/