s-news
[Top] [All Lists]

Strip-split-plot analysis with lme

To: Snews <s-news@lists.biostat.wustl.edu>
Subject: Strip-split-plot analysis with lme
From: quataepa <paul.quataert@lin.vlaanderen.be>
Date: Wed, 19 Jan 2005 13:41:45 +0100
User-agent: Mozilla/5.0 (Windows; U; WinNT4.0; en-US; rv:1.0.2) Gecko/20021120 Netscape/7.01
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="">
 
where Wplot and Splot are the whole plot and strip plot variables.
 
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="">
 
but this is not allowed by the syntaxis.
 
Thanks,
Paul Quataert
Institute for Forestry and Game Management
Belgium
<Prev in Thread] Current Thread [Next in Thread>