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]]


