## Paul
## I think you are looking for interaction2wt
y <- rnorm(36)
f1 <- rep(c("after", "before"), 18)
f2 <- rep(1:3, 12)
f3 <- rep(1:2, each=18)
## your definition of ff was faulty. It gave a constant.
f3.f1 <- interaction(f3, f1)
interaction.plot(f3.f1, f2, y)
f2 <- factor(f2)
f3 <- factor(f3)
## lattice would be better
xyplot(y ~ f2 | f3, groups=f2)
## if you don't have HH, you will need to install it and its dependencies
## install.packages("HH")
## You may need to close R and restart R
require(HH)
interaction2wt(y ~ f1 + f2 + f3)
position(f1) <- c(1.5, 2.5)
position(f3) <- c(1.5, 2.5)
interaction2wt(y ~ f1 + f2 + f3, box.ratio=.6)
## Rich
______________________________________________
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.
