From: Addison, Prue <paddison_at_museum.vic.gov.au>

Date: Mon 18 Jul 2005 - 13:59:58 EST

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 Received on Mon Jul 18 14:07:44 2005

Date: Mon 18 Jul 2005 - 13:59:58 EST

Hi,

I am having trouble using the lme function to perform a nested ANOVA with a random nested factor.

My design is as follows:

Site nested within each Location (n=12) (2 Sites nested within each Location) (Random)

Dependent variable: sp (species abundance)

> summary(aov(sp~Location/Site, data=mavric))

Df Sum Sq Mean Sq F value Pr(>F)

Location 4 112366 28092 1.2742 0.2962

Location:Transect 5 121690 24338 1.1039 0.3736

Residuals 40 881875 22047

I have tried the following lme function to specify that Site is random:

> lme1 <- lme(sp~Location, random=~1|Site, data=mavric)

> lme2 <- lme(sp~Location, random=~1|Location/Site, data=mavric)

> anova(lme1)

numDF denDF F-value p-value

(Intercept) 1 40 3.418077 0.0719

> anova(lme2)

Warning in pf(q, df1, df2, lower.tail, log.p) :

Warning: NAs introduced by coercion

(Intercept) 1 40 0.459966 0.5015

? I don't know what this output means

Model df AIC BIC logLik Test L.Ratio p-value

lme1 1 7 603.7534 616.4000 -294.8767

lme2 2 8 605.7534 620.2067 -294.8767 1 vs 2 1.815674e-05 0.9966

? I also don't know what this output means.

Thanks,

Prue

Museum Victoria

+61 3 8341 7777

11 Nicholson St

Carlton

Victoria

www.museum.vic.gov.au

[[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 Received on Mon Jul 18 14:07:44 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:33:46 EST
*