[R] Dunnett Test in 'multicomp' package

[James F. Henson]

Greetings

Below is my code.

library("multcomp")
viaModel1 <- aov(libido ~ dose, data=viagraData)
dunnettModel <- glht(viaModel1 , linfct = mcp(dose = "Dunnett"), base = "placebo")

The code base="placebo" is ignored.  All treatments are compared to the first treatment in the order, which is "high dose".  It is possible to rearrange the order so that "placebo' is first, but this is inconvenient.

Thanks,
James F. Henson

	[[alternative HTML version deleted]]