Here are my data:
Time Treatment Ephemeroptera Trichoptera Plecoptera Before Control 98 76 32 After Control 91 54 26 Before Impact 24 30 12 After Impact 29 21 45
Exp = read.csv("fast.csv", header = T) attach(Exp) names(Exp) m1 = aov(Ephemeroptera ~ Time) summary(m1) #It gives me output with p-value, df etc. (which is expected) m2 = aov(Ephemeroptera ~ Treatment) summary(m2) #It gives me output with p-value, df etc. (which is expected) m3 = aov(Ephemeroptera ~ Time*Treatment) summary(m3) #It gives me the interaction output (Time:Treatment), but it doesn't provide a p-value.
I am not on my home computer, so I can't include what I got for the output right now. However, I will update it in a little bit (when I get home). Until then, is there any reason why it does not give me any p-values for the output when I try the repeated measures, but it provides it for the regular anovas? Should I try lme instead?