[R-sig-ME] sas to R

Kevin Wright kw.stat at gmail.com
Mon Jun 25 17:25:29 CEST 2012


This could be similar to a multi-location RCB design were "trial" is
location.  Since no distribution is specified, the distribution is
assumed to be Gaussian.  Make sure that trial, block, trt are factors,
this should be similar to SAS:

lmer(y ~ trt + (1|trial/block/trt), data=df)

> proc glimmix data=df;
> class trial block trt;
> model y=trt;
> random trial block(trial) trt(block*trial);

Kevin Wright



More information about the R-sig-mixed-models mailing list