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:>Hello.
>
>I am looking for some help using anova's in RGui.
>
>My experiment ie data, has a fixed burning treatment (Factor A) 2 levels,
>unburnt/burnt.
>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
>aov(dependent
>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
>email?
>
>Thank-you,
>Frith
>
>______________________________________________
>R-help at 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
Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat.
Centre for Resource and Environmental Studies
The Australian National University
Canberra ACT 0200
Australia
T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au
F: +61 2 6125 0757
CRICOS Provider # 00120C