Dear R Help -

I am analyzing data from an ecological experiment and am having problems
with the ANOVA functions I've tried thus far. The experiment consists of a
blocked/split-plot design, with plant biomass as the response. The following
is an overview of the treatments applied (nitrogen addition, phosphorus
addition, and seeded/not seeded) and at what level (block, main-plot, and
sub-plot):

- 6 experimental blocks divided into 2 main-plots
- each main-plot divided into 8 sub-plots, for a total of 96 sub-plots (6
blocks * 2 main-plots * 8 sub-plots)
- 16 experimental treatment conditions in a 4 x 2 x 2 factorial design:
   - N at 4 levels for each sub-plot
   - P at 2 levels for each sub-plot
   - Seed at two levels for each block (one level per main-plot)
- response variable = biomass
Block
   Main Plot 0 : No Seed
       Sub-Plot 0.1 : N0 P0
       Sub-Plot 0.2 : N0 P1
       Sub-Plot 0.3 : N1 P0
       Sub-Plot 0.4 : N1 P1
       Sub-Plot 0.5 : N2 P0
       Sub-Plot 0.6 : N2 P1
       Sub-Plot 0.7 : N3 P0
       Sub-Plot 0.8 : N3 P1
   Main Plot 1 : Seed
       Sub-Plot 1.1 : N0 P0
       Sub-Plot 1.2 : N0 P1
       Sub-Plot 1.3 : N1 P0
       Sub-Plot 1.4 : N1 P1
       Sub-Plot 1.5 : N2 P0
       Sub-Plot 1.6 : N2 P1
       Sub-Plot 1.7 : N3 P0
       Sub-Plot 1.8 : N3 P1

I've tried several different approaches to run an ANOVA (lmer, aov, lme) on
this data, trying to use type III SSs and include a random factor, but am
having problems. Any suggestions?

Thank you for your help!

-Mari Pesek
-- 
Dept of Ecology and Evolutionary Biology
University of Kansas
Haworth Hall, 1200 Sunnyside Ave
Lawrence, KS 66045

	[[alternative HTML version deleted]]

