Hello list,
I' d very much appreciate some help with a two sample repeated measures
ANOVA.
I did the analysis yielding sign. main effects (between subj.=site, within
subj.=cover) and a sign. interaction:
Univariate Type II Repeated-Measures ANOVA Assuming Sphericity
SS num Df Error SS den Df F Pr(>F)
site 18.7620 1 18.831 10 9.9631 0.010220 *
cover 6.4481 1 16.462 10 3.9169 0.075984 .
site:cover 19.2963 1 16.462 10 11.7216 0.006507 **
---
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
But then I was interested in the significance of the within factor in each
of the levels of the between factor.
Also, I'm a bit confused if the type of SS does play a role here - I choose
the default, the design is balanced.
In despair of a better solution I did two seperate paired t-tests on each
level, but I guess that's not the best solution.
Thanks for any advise,
Kay Cichini
--
View this message in context:
http://n4.nabble.com/repeated-measures-anova-car-package-tp1573721p1573721.html
Sent from the R help mailing list archive at Nabble.com.
Dear Kay,> -----Original Message----- > From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On > Behalf Of Kay Cichini > Sent: March-01-10 10:45 AM > To: r-help at r-project.org > Subject: [R] repeated measures anova, car package > > > Hello list, > > I' d very much appreciate some help with a two sample repeated measures > ANOVA. > I did the analysis yielding sign. main effects (between subj.=site, within > subj.=cover) and a sign. interaction: > > Univariate Type II Repeated-Measures ANOVA Assuming Sphericity > > SS num Df Error SS den Df F Pr(>F) > site 18.7620 1 18.831 10 9.9631 0.010220 * > cover 6.4481 1 16.462 10 3.9169 0.075984 . > site:cover 19.2963 1 16.462 10 11.7216 0.006507 ** > --- > Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 > > But then I was interested in the significance of the within factor in each > of the levels of the between factor. > Also, I'm a bit confused if the type of SS does play a role here - I choose > the default, the design is balanced. > In despair of a better solution I did two seperate paired t-tests on each > level, but I guess that's not the best solution. > > Thanks for any advise, > Kay CichiniIn this case, the repeated-measures ANOVA is equivalent to analyzing the sum and difference of the two measures on cover. If you really want to test for the difference between the two levels of cover in each level of site, then the paired t-tests are perfectly reasonable. You could also these test hypotheses this using the car functions, e.g., by fitting the model mod <- lm(cbind(cover1, cover2) ~ site - 1, data=Data), which would fit a means model for site rather than use an intercept and dummy-coded contrast. Then you can use the linear.hypothesis() function to test each of the two hypotheses -- e.g., linear.hypothesis(mod, "siteA", P=matrix(c(1, -1)). I'm making up the level name for site (A) and the names of the repeated-measures variables (cover1, cover2) because you didn't say what they are. Because there's just one between-subjects factor, there's no distinction between type-II and -III tests, as you would have discovered had you tried both. I hope this helps, John> -- > View this message in context: http://n4.nabble.com/repeated-measures-anova- > car-package-tp1573721p1573721.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help at 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.
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)> idatcover 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> divreprepl. 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.
sorry, it should say "I copied the model from the car manual.." more precise, I did the same as there in doing the lm(cbind()) call first an then the av.ok, etc. -- View this message in context: http://n4.nabble.com/repeated-measures-anova-car-package-tp1573721p1577432.html Sent from the R help mailing list archive at Nabble.com.
gReetings
Running R 2.10.1 with Notepad++ and NpptoR.exe, when I submit a script to
run :-
1. On my work XP machine , it runs fine in an already-open Rgui window
2. On my home Windows 7 machine, it ignores the existing Rgui window and
fires up
a fresh one, which doesn't use the same .RData, so it fails.
Anyone else hit this and, ideally, solved it ... ?
cheers Bob
Robert Kinley
bobk@lilly.com
[[alternative HTML version deleted]]
Maybe Matching Threads
- Lost setting up IAXmodem after drive crash
- Re: HA: Accessing iSCSI disc images from the RHEV Manager using libguestfs
- Kernel build failure
- Re: HA: Accessing iSCSI disc images from the RHEV Manager using libguestfs
- [xen-unstable] Commit 2ca9fbd739b8a72b16dd790d0fff7b75f5488fb8 AMD IOMMU: allocate IRTE entries instead of using a static mapping, makes dom0 boot process stall several times.