powered by
mcalinfct(model, focus, mmm.data=model$model, formula.in=terms(model), linfct.Means= multcomp:::meanslinfct(model, focus, mmm.data, formula.in, contrasts.arg=model$contrasts), type="Tukey" )
linfct
glht
MMC
## See the examples in HH/scripts/MMC.cc176.R
Run the code above in your browser using DataLab