NPMTest {PMCMRplus} | R Documentation |
Performs Nashimoto and Wright's all-pairs comparison procedure for simply ordered mean ranksums. Their test denoted as NPM test is basically an extension of Nemenyi's procedure for testing increasingly ordered alternatives.
The modified procedure uses the property of a simple order, θ_m' - θ_m ≤ θ_j - θ_i ≤ θ_l' - θ_l \qquad (l ≤ i ≤ m~\mathrm{and}~ m' ≤ j ≤ l'). The null hypothesis H_{ij}: θ_i = θ_j is tested against the alternative A_{ij}: θ_i < θ_j for any 1 ≤ i < j ≤ k.
The p-values are estimated from the studentized range distribution.
If the medians are already increasingly ordered, than the NPM-test simplifies
to the ordinary Nemenyi test (see kwAllPairsNemenyiTest
).
NPMTest(x, ...) ## Default S3 method: NPMTest(x, g, ...) ## S3 method for class 'formula' NPMTest(formula, data, subset, na.action, ...)
x |
a numeric vector of data values, or a list of numeric data vectors. |
... |
further arguments to be passed to or from methods. |
g |
a vector or factor object giving the group for the
corresponding elements of |
formula |
a formula of the form |
data |
an optional matrix or data frame (or similar: see
|
subset |
an optional vector specifying a subset of observations to be used. |
na.action |
a function which indicates what should happen when
the data contain |
A list with class "PMCMR"
containing the following components:
a character string indicating what type of test was performed.
a character string giving the name(s) of the data.
lower-triangle matrix of the estimated quantiles of the pairwise test statistics.
lower-triangle matrix of the p-values for the pairwise tests.
a character string describing the alternative hypothesis.
a character string describing the method for p-value adjustment.
a data frame of the input data.
a string that denotes the test distribution.
## Example from Sachs (1997, p. 402) x <- c(106, 114, 116, 127, 145, 110, 125, 143, 148, 151, 136, 139, 149, 160, 174) g <- gl(3,5) levels(g) <- c("A", "B", "C") NPMTest(x, g)