[R] best way to calculate per-parameter differences in across-subject means

Dieter Menne dieter.menne at menne-biomed.de
Wed Aug 2 18:04:52 CEST 2006

René J.V. Bertin <rjvbertin <at> gmail.com> writes:

> I have some data in a data.frame where for each of a number of
> subjects, I have scores for all of a number of symptoms.
> Subjects are subdivided in a number of groups, which have unequal sizes.
> I'd like to plot between-group differences in the scores on the
> various symptoms. Ideally, that would be in a form as would be
> produced by


Maybe it's a bit more than you want, but possibly you are happy with it: see the
example under TukeyHSD.

summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
TukeyHSD(fm1, "tension", ordered = TRUE)
plot(TukeyHSD(fm1, "tension"))


More information about the R-help mailing list