[R] re peated measures
pompon
julien.pompon at agr.gc.ca
Wed Sep 23 14:27:40 CEST 2009
Hi,
I am performing a repeated measures 2-way ANOVA to assess the influence of
plant and leaf on aphid fecundity. Fecundity is measured for each aphid on a
single leaf.
Here is what I typed.
wingless <- reshape(Wingless,
varying =
list(c("d0","d1","d2","d3","d4","d5","d6","d7","d8","d9","d10","d11","d12","d13","d14","d15","d16")),
v.names = c("fecundity"), timevar = "time",
direction = "long")
wingless.aov <- aov(fecundity ~ factor(time) * clip.cage * plant +
Error(factor(id)), data = wingless)
summary(wingless.aov)
and I obtained
Error: factor(id)
Df Sum Sq Mean Sq F value Pr(>F)
factor(time) 4 56.789 14.197 3.0613 0.05925 .
clip.cage 1 14.149 14.149 3.0509 0.10621
plant 1 3.251 3.251 0.7010 0.41880
factor(time):clip.cage 1 0.304 0.304 0.0655 0.80240
clip.cage:plant 1 17.114 17.114 3.6903 0.07880 .
Residuals 12 55.652 4.638
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
factor(time) 16 340.83 21.30 11.5222 < 2e-16 ***
factor(time):clip.cage 16 27.34 1.71 0.9242 0.54195
factor(time):plant 16 46.36 2.90 1.5673 0.07783 .
factor(time):clip.cage:plant 16 24.50 1.53 0.8281 0.65304
Residuals 255 471.44 1.85
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
I don't understand why I have the factor(time) inmy between subject results,
whereas with a similar set of data I don't.
Thank you very much,
Julien Pompon.
--
View this message in context: http://www.nabble.com/repeated-measures-tp25531110p25531110.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list