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]]