Re: [R] Random structure of nested design in lme

From: Doran, Harold <HDoran_at_air.org>
Date: Wed 19 Jul 2006 - 21:53:42 EST


Can you provide an example of what you have done with lme so we might be able to evaluate the issue?

> -----Original Message-----
> From: r-help-bounces@stat.math.ethz.ch
> [mailto:r-help-bounces@stat.math.ethz.ch] On Behalf Of ESCHEN Rene
> Sent: Wednesday, July 19, 2006 7:37 AM
> To: r-help@stat.math.ethz.ch
> Subject: [R] Random structure of nested design in lme
>
> All,
>
> I'm trying to analyze the results of a reciprocal transplant
> experiment using lme(). While I get the error-term right in
> aov(), in lme() it appears impossible to get as expected. I
> would be greatful for any help.
>
> My experiment aimed to identify whether two fixed factors
> (habitat type and soil type) affect the development of
> plants. I took soil from six random sites each of two types
> (arable and grassland) and transplanted them back into the
> sites of origin in such way that in each of the sites there
> were six pots containing arable soil and six pots of
> grassland soil, each containing a seedling.
>
> With aov(), I got the analysis as I expected, with habitat
> type tested against destination site, and soil type tested
> against origin site:
>
> summary(aov(response~soiltype*habitat+Error(destination+origin)))
> #
> #Error: destination
> # Df Sum Sq Mean Sq F value Pr(>F)
> #habitat 1 1.0000 1.0000 0.699 0.4226
> #Residuals 10 14.3056 1.4306
> #
> #Error: origin
> # Df Sum Sq Mean Sq F value Pr(>F)
> #soiltype 1 1.77778 1.77778 11.636 0.006645 **
> #Residuals 10 1.52778 0.15278
> #---
> #Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 #
> #Error: Within
> # Df Sum Sq Mean Sq F value Pr(>F)
> #soiltype:habitat 1 0.2500 0.2500 2.1774 0.1427
> #Residuals 120 13.7778 0.1148
>
> However, when I try to replicate this analysis in lme, I am
> unable to get the structure of the random factors (origin and
> destination) correct. Does anyone have a suggestion how to
> resolve this problem?
>
> Thanks in advance.
>
> René Eschen
>
> CABI Bioscience Centre Switzerland
> Rue des Grillons 1
> 2800 Delémont
> Switzerland
>
> [[alternative HTML version deleted]]
>
>



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu Jul 20 01:20:19 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 20 Jul 2006 - 04:16:46 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.