critVal function

Calculate critical value for multiple contrast test

Calculation of the critical value for a maximum contrast test. This is based on the equicoordinate quantile function of the multivariate normal or t distribution as implemented in the qmvt function from the mvtnorm package.

critVal(
  corMat,
  alpha = 0.025,
  df = NULL,
  alternative = c("one.sided", "two.sided"),
  control = mvtnorm.control()
)

Arguments

• corMat: Correlation matrix of contrasts
• alpha: Significance level for the multiple contrast test
• df: Specify the degrees of freedom to use, if this argument is missing df = Inf is used (which corresponds to the multivariate normal distribution).
• alternative: Character determining the alternative for the multiple contrast trend test.
• control: A list specifying additional control parameters for the qmvt and pmvt calls in the code, see also mvtnorm.control for details.

Examples

R <- matrix(c(1,0.5,0.5,1), nrow=2)
critVal(R, alpha = 0.05, df = 1)
critVal(R, alpha = 0.05, df = 20)
critVal(R, alpha = 0.05, df = Inf)

See Also

powMCT, optContr, MCTtest

Author(s)

Bjoern Bornkamp