MCTpval {DoseFinding}

Calculate multiplicity adjusted p-values for multiple contrast test

Description

Calculate multiplicity adjusted p-values for a maximum contrast test corresponding to a set of contrasts and given a set of observed test statistics. This function is exported as it may be a useful building block and used in more complex testing situations that are not covered by MCTtest. Most users probably don't need to use this function.

Usage

MCTpval(contMat, corMat, df, tStat,
        alternative = c("one.sided", "two.sided"),
        control = mvtnorm.control())

Arguments

contMat

Contrast matrix to use. The individual contrasts should be saved in the columns of the matrix

corMat

correlation matrix of the contrasts

df

Degrees of freedom to assume in case ‘⁠S⁠’ (a general covariance matrix) is specified. When ‘⁠n⁠’ and ‘⁠sigma⁠’ are specified the ones from the corresponding ANOVA model are calculated.

tStat

Vector of contrast test statistics

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.

Value

Numeric containing the calculated p-values.

Author(s)

Bjoern Bornkamp

References

Pinheiro, J. C., Bornkamp, B., and Bretz, F. (2006). Design and analysis of dose finding studies combining multiple comparisons and modeling procedures, Journal of Biopharmaceutical Statistics, 16, 639–656

See Also

MCTtest, optContr

Examples

## need to add example