Hi,

I have the following dataset (shown only part of it):

WFPS    sampling day    slurry  depth   replicates      N2O (mg/m2/h)
1       1       1       3       1       0.007
1       1       1       3       2       0.022
1       1       1       3       3       0.015
1       1       2       3       1       0.033
1       1       2       3       2       0.010
1       1       2       3       3       0.022
1       1       2       6       1       0.038
1       1       2       6       2       0.048
1       1       2       6       3       0.043
1       1       3       0       1       0.024
1       2       1       3       1       0.017
1       2       1       3       2       0.016
1       2       1       3       3       0.015
1       2       2       3       1       0.047
1       2       2       3       2       0.083
1       2       2       3       3       0.042
1       2       2       6       1       0.036
1       2       2       6       2       0.048
1       2       2       6       3       0.030
1       2       3       0       1       0.020
1       3       1       3       1       0.006
1       3       1       3       2       0.007
1       3       1       3       3       0.019
1       3       2       3       1       0.055
1       3       2       3       2       0.004
1       3       2       3       3       0.025
1       3       2       6       1       0.030
1       3       2       6       2       0.024
1       3       2       6       3       0.042
1       3       3       0       1       0.017


What it is: daily N2O emission rate after the injection of 2 different types of slurry into the soil moistened at 2 different levels.
So I want to test all the factors on the N2O emissions:
1) injection depth
2) slurry type
3) water content into the soil (=WFPS)
Bu I want to take into account the replicates. How can I write the R syntax in this case? Something like:
 model2<-aov(N2O.mg.m2.h. ~ WFPS*sampling_day*slurry*depth + Error(replicates/(WFPS + sampling_day + slurry + depth)))



Thanks in advance,
Joao



	[[alternative HTML version deleted]]


