Re: Strip-split-plot analysis with lme

[Prof Brian Ripley <ripley@stats.ox.ac.uk>]

On Wed, 19 Jan 2005, quataepa wrote:

> Dear S-Plusers,
>
> How to perform a strip-split-plot analysis in S+ (version 6.2) with lme?
>
> Specifically I want to perform following analysis
>
> Afit <- aov ( Gas ~ Treatment * Wplot * Splot +
>   Error( Plant / (Wplot + Splot) ) , data=Wset)
>
> where Wplot and Splot are the whole plot and strip plot variables.
I would expect Splot to be nested in Wplot, not in Plant, given 
your description, so how have you encoded the variables?  And what is 
Plant?
> With aov no problem, however I do no succeed to translate this to lme.
>
> I tried:
>
> Afit <- lme ( Gas ~ Treatment * Wplot * Splot,
>   random = ~1 | Plant / (Wplot + Splot) , data=Wset)
>
> but this is not allowed by the syntaxis.
You need to specify the nesting, probably by Plant/Wplot/Splot.

--
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595