[R-sig-eco] manova repeated measures + Error
Marina Vergara Fagundes
marvefa at gmail.com
Fri Feb 13 16:58:55 CET 2015
Dear all,
I conducted a interaction plant experiment, with 20 adult plant species (5
individuals each) and 3 seedling species, and a 100 plots in total. Each
plot consist of two 2m² quadrants, the first contains one adult individual
as a potential “nurse” and 3 seedling species below their canopy, and the
second just with 3 seedling species with no canopy influence (in an open
area).
I measured growth variables of the seedlings (leaves number and height)
over 8 weaks (times).
For each growth variable I calculated a index (rii) with the difference of
growth below and far from canopy.
I'm analyzing this data with MANOVA with repeated measures. Where Y axes
are the index of growth measurements (for leafs and height), and the X axes
are the 20 nurse species and interactions with seedlings species and time.
This is how my data frame is organized:
plot nurse target time height_rii leaves_rii
113 algodao angico 1 0.03703704 0.00000000
113 algodao aroeira 2 0.05660377 0.09090909
113 algodao cating 3 0.03846154 0.00000000
115 algodao angico 1 0.01818182 0.11111111
115 algodao aroeira 2 0.05660377 0.50000000
115 algodao cating 3 0.03703704 0.00000000
This is the code for the MANOVA:
Y <- with(tabela3, cbind(altura_rii, folhas_rii))
m1 <- manova( Y ~ nurse*target*time+Error(tempo/plot),data=tabela3)
summary(m1)
When i ask for anova table i got this specific message:
*Error: timeNo error degrees of freedom Dftempo 7*Error: time:plot
Df Pillai approx F num Df den Df Pr(>F)
nurse 19 0.17070 3.14316 38 1280 8.633e-10 ***
nurse:time 133 0.24253 0.66405 266 1280 1
Residuals 640
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Pillai approx F num Df den Df Pr(>F)
target 2 0.013439 4.3297 4 2560 0.001715 **
nurse:target 38 0.178898 3.3090 76 2560 < 2.2e-16 ***
target:time 14 0.027914 1.2941 28 2560 0.138459
nurse:target:time 266 0.174314 0.4594 532 2560 1.000000
Residuals 1280
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Any idea what is it about?
Thanks a lot for help!
Marina
[[alternative HTML version deleted]]
More information about the R-sig-ecology
mailing list