Hello dear R members.
I have been learning the Anova syntax in order to perform an SS type III
Anova with repeated measures designs (thank you Prof. John Fox!)
And another question came up: where/what are the (between/within) residuals
for my model?
############ Play code:
phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),
levels=c("pretest", "posttest", "followup"))
hour <- ordered(rep(1:5, 3))
idata <- data.frame(phase, hour)
idata
mod.ok <- lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5,
post.1, post.2, post.3, post.4, post.5,
fup.1, fup.2, fup.3, fup.4, fup.5) ~
treatment*gender,
data=OBrienKaiser)
av.ok <- Anova(mod.ok, idata=idata, idesign=~phase*hour)
summary(av.ok, multivariate=FALSE)
## Univariate Type II Repeated-Measures ANOVA Assuming Sphericity
##
## SS num Df Error SS den Df F
Pr(>F)
## treatment 211.286 2 228.056 10 4.6323
0.037687
## gender 58.286 1 228.056 10 2.5558
0.140974
## treatment:gender 130.241 2 228.056 10 2.8555
0.104469
## phase 167.500 2 80.278 20 20.8651
1.274e-05
## treatment:phase 78.668 4 80.278 20 4.8997
0.006426
## gender:phase 1.668 2 80.278 20 0.2078
0.814130
## treatment:gender:phase 10.221 4 80.278 20 0.6366
0.642369
## hour 106.292 4 62.500 40 17.0067
3.191e-08
## treatment:hour 1.161 8 62.500 40 0.0929
0.999257
## gender:hour 2.559 4 62.500 40 0.4094
0.800772
## treatment:gender:hour 7.755 8 62.500 40 0.6204
0.755484
## phase:hour 11.083 8 96.167 80 1.1525
0.338317
## treatment:phase:hour 6.262 16 96.167 80 0.3256
0.992814
## gender:phase:hour 6.636 8 96.167 80 0.6900
0.699124
## treatment:gender:phase:hour 14.155 16 96.167 80 0.7359
0.749562
--
----------------------------------------------
My contact information:
Tal Galili
Phone number: 972-50-3373767
FaceBook: Tal Galili
My Blogs:
www.talgalili.com
www.biostatistics.co.il
[[alternative HTML version deleted]]
______________________________________________
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.
