DoseFinding (version 1.1-1)

critVal: Calculate critical value for multiple contrast test

Description

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.

Usage

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.

Author

Bjoern Bornkamp

See Also

powMCT, optContr, MCTtest

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)