[R] Anova's in R

Simon Blomberg simonb at cres10.anu.edu.au
Fri Jul 29 01:39:56 CEST 2005

If I understand you correctly, this looks like a split-plot design.

The anova is:

aov(response~burning*temperature+Error(site), data=dataset)

Alternatively in lme:

lme(response~burning*temperature, random=~1|site, data=dataset)

At 03:09 PM 29/07/2005, you wrote:
>I am looking for some help using anova's in RGui.
>My experiment ie data, has a fixed burning treatment (Factor A) 2 levels,
>Nested within each level of Factor A are 2 random sites (Factor B).
>All sites are crossed with a fixed temperature treatment (Factor C) 2
>levels, 0 degreesC/2 degreesC, with many replicates of these temperature
>treatments randomly located at each site.
>I am trying the following
>variable~burning*temperature*site+Error(replicate),data=dataset) and
>variations on that, however can't get it quite right.... the F ratios are
>not correct. I imagine this is a fairly common experimental design in
>ecology and would ask that anyone who has some advice please reply to this
>R-help at stat.math.ethz.ch mailing list
>PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat.
Centre for Resource and Environmental Studies
The Australian National University
Canberra ACT 0200
T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au
F: +61 2 6125 0757
CRICOS Provider # 00120C

More information about the R-help mailing list