[R] Balanced incomplete block analysis

[Murray Jorgensen]

Excuse me everyone, but I don't have to teach this very often!

Has anyone got some R code for doing adjusted treatment means and the
recovery of inter-block information in the analysis of balanced incomplete
block designs?

With data from O. L. Davies, page 207 I've got as far as:

wear <- scan()
238 196 254 238 213 312 279 334 421 308 367 412

treatment <- factor(rep(1:4,rep(3,4)))
tyre <- scan()
1 2 3 1 2 4 1 3 4 2 3 4

tyre <- as.factor(tyre)
tapply(wear,list(tyre,treatment),sum)
wear.lm <- lm(wear ~ tyre + treatment)
anova(wear.lm)
summary(wear.lm)


Which is neat as far as it goes, but adjusted treatment means and their
LSDs would be nice. 


Murray Jorgensen


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._