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.
MCTpval(contMat, corMat, df, tStat,
alternative = c("one.sided", "two.sided"),
control = mvtnorm.control())Numeric containing the calculated p-values.
Contrast matrix to use. The individual contrasts should be saved in the columns of the matrix
correlation matrix of the contrasts
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.
Vector of contrast test statistics
Character determining the alternative for the multiple contrast trend test.
A list specifying additional control parameters for the qmvt and pmvt calls in the code, see also mvtnorm.control for details.
Bjoern Bornkamp
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
MCTtest, optContr