Demsar's Many-to-One Test for Unreplicated Blocked Data
Performs Demsar's non-parametric many-to-one comparison test for Friedman-type ranked data.
frdManyOneDemsarTest(y, ...)
## Default S3 method:
frdManyOneDemsarTest(
y,
groups,
blocks,
alternative = c("two.sided", "greater", "less"),
p.adjust.method = p.adjust.methods,
...
)y |
a numeric vector of data values, or a list of numeric data vectors. |
groups |
a vector or factor object giving the group for the
corresponding elements of |
blocks |
a vector or factor object giving the block for the
corresponding elements of |
alternative |
the alternative hypothesis. Defaults to |
p.adjust.method |
method for adjusting p values
(see |
... |
further arguments to be passed to or from methods. |
For many-to-one comparisons (pairwise comparisons with one control) in a two factorial unreplicated complete block design with non-normally distributed residuals, Demsar's test can be performed on Friedman-type ranked data.
Let there be k groups including the control, then the number of treatment levels is m = k - 1. A total of m pairwise comparisons can be performed between the i-th treatment level and the control. H_i: θ_0 = θ_i is tested in the two-tailed case against A_i: θ_0 \ne θ_i, ~~ (1 ≤ i ≤ m).
The p-values are computed from the standard normal distribution.
Any of the p-adjustment methods as included in p.adjust
can be used for the adjustment of p-values.
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.
Demsar, J. (2006) Statistical comparisons of classifiers over multiple data sets, Journal of Machine Learning Research 7, 1–30.
## Sachs, 1997, p. 675 ## Six persons (block) received six different diuretics ## (A to F, treatment). ## The responses are the Na-concentration (mval) ## in the urine measured 2 hours after each treatment. ## Assume A is the control. y <- matrix(c( 3.88, 5.64, 5.76, 4.25, 5.91, 4.33, 30.58, 30.14, 16.92, 23.19, 26.74, 10.91, 25.24, 33.52, 25.45, 18.85, 20.45, 26.67, 4.44, 7.94, 4.04, 4.4, 4.23, 4.36, 29.41, 30.72, 32.92, 28.23, 23.35, 12, 38.87, 33.12, 39.15, 28.06, 38.23, 26.65),nrow=6, ncol=6, dimnames=list(1:6, LETTERS[1:6])) ## Global Friedman test friedmanTest(y) ## Demsar's many-one test frdManyOneDemsarTest(y=y, p.adjust = "bonferroni") ## Exact many-one test frdManyOneExactTest(y=y, p.adjust = "bonferroni") ## Nemenyi's many-one test frdManyOneNemenyiTest(y=y)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.