Hello John,
As you said, I could also take a means model and test linear hypothesis for
the desired effects - would this also be the case for the repeated measure i
did in the first place.
I copied the model from the car model where you first call:
> modx<-lm(cbind(div_h, div_l) ~ site, divrep)
(?Could I test linear hypothesis here, instead of continuing as I did
beneath)
> idat
cover
1 high
2 low
> (av.ok1 <- Anova(modx, idata=idat, idesign=~cover))
Type II Repeated Measures MANOVA Tests: Pillai test statistic
Df test stat approx F num Df den Df Pr(>F)
site 1 0.49908 9.9631 1 10 0.010220 *
cover 1 0.28145 3.9169 1 10 0.075984 .
site:cover 1 0.53963 11.7216 1 10 0.006507 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> divrep
repl. site div_h div_l
1 1 Scrub 4.18 5.23
2 2 Scrub 5.47 7.18
3 3 Scrub 3.74 4.97
4 4 Scrub 2.62 5.17
5 5 Scrub 3.33 6.43
6 6 Scrub 1.62 8.96
7 1 Tall_Forb 4.70 3.88
8 2 Tall_Forb 3.65 1.97
9 3 Tall_Forb 2.50 1.19
10 4 Tall_Forb 1.87 2.37
11 5 Tall_Forb 5.33 3.56
12 6 Tall_Forb 3.06 3.60
Your answers helped a lot -
Thank you very much for the quick reply!
Best wishes,
Kay
--
View this message in context:
http://n4.nabble.com/repeated-measures-anova-car-package-tp1573721p1574747.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
