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

[Dieter Menne]

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"))


Dieter