R/mcalinfct.R

In HH: Statistical Analysis and Data Display: Heiberger and Holland

## MCA multiple comparisons analysis (pairwise)
if.R(r=
mcalinfct <-
  function(model, focus,
           mmm.data=model$model,
           formula.in=terms(model),
           linfct.Means=
           ## multcomp:::meanslinfct(model, focus, mmm.data, formula.in,
           multcomp.meanslinfct(model, focus, mmm.data, formula.in,
           contrasts.arg=model$contrasts),
           type="Tukey"
           ) {
    if (inherits(model, "lme")) ## changed rmh
      mmm.data <- model$data    ## changed rmh

    lev <- levels(mmm.data[[focus]])
    names(lev) <- lev
    mca <- contrMat(table(mmm.data[[focus]]), type)
##     result <-
    mca[, lev] %*% linfct.Means[lev,]
##     attr(result, "type") <- "Tukey"
##     attr(result, "class") <- c("contrMat", "matrix")
##     mca <- list(result)
##     names(mca) <- focus
##     attr(mca, "class") <- "mcp"
##     mca
  }
,s={})

## source("~/HH-R.package/HH/R/mcalinfct.R")