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())
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.
Numeric containing the calculated p-values.
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
# NOT RUN {
## need to add example
# }
Run the code above in your browser using DataLab