Here is your example. The table you displayed in gigawiz ignored the two-way factor structure and interpreted the data as a single factor with 6 levels. I created the interaction of a and b to get that behavior.
## your example, with data stored in a data.frame tmp <- data.frame(x=c(76, 84, 78, 80, 82, 70, 62, 72, 71, 69, 72, 74, 66, 74, 68, 66, 69, 72, 72, 78, 74, 71, 73, 67, 86, 67, 72, 85, 87, 74, 83, 86, 66, 68, 70, 76, 78, 76, 69, 74, 72, 72, 76, 69, 69, 82, 79, 81), a=factor(rep(c("A1", "A2"), each = 24)), b=factor(rep(c("B1", "B2", "B3"), each=8, times=2))) x.aov <- aov(x ~ a*b, data=tmp) summary(x.aov) ## your request require(multcomp) tmp$ab <- with(tmp, interaction(a, b)) xi.aov <- aov(x ~ ab, data=tmp) summary(xi.aov) xi.glht <- glht(xi.aov, linfct=mcp(ab="Tukey")) confint(xi.glht) ## graphs ## boxplots require(lattice) bwplot(x ~ ab, data=tmp) ## interaction plot ## install.packages("HH") ## if you don't have HH yet require(HH) interaction2wt(x ~ a*b, data=tmp) On Wed, Apr 4, 2012 at 1:10 PM, Jinsong Zhao <jsz...@yeah.net> wrote: > On 2012-04-03 20:03, Rmh wrote: > >> yes. See ?glht in the multcomp package, and the examples using glht in >> ?MMC in the HH package. >> >> Sent from my iPhone >> >> > Thank you very much for the clues. However, I can't figure out how to > construct the linfct in glht. > > I also tried to inverse the computation based on: > http://www.gigawiz.com/**images12/twowayrmposthoc.jpg<http://www.gigawiz.com/images12/twowayrmposthoc.jpg> > However, I can't catch the MSE used in the above figure. > > Regards, > Jinsong > > [[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.