| snk {mutoss} | R Documentation |
Student - Newman - Keuls rejective test procedure.
Description
Student - Newman - Keuls rejective test procedure. The procedure controls the FWER in the WEAK sense.
Usage
snk(formula, data, alpha, MSE=NULL, df=NULL, silent=FALSE)
snk.wrapper(model, data, alpha, silent=FALSE)Arguments
formula |
Formula defining the statistical model containing the response and the factor levels. |
model |
Model with formula, containing the response and the factor levels |
data |
dataset containing the response and the grouping factor. |
alpha |
The level at which the error should be controlled. By default it is alpha=0.05. |
MSE |
Optional for a given variance of the data. |
df |
Optional for a given degree of freedom. |
silent |
If true any output on the console will be suppressed. |
Details
This function computes the Student-Newman-Keuls test for given
data including p samples. The Newman-Keuls procedure is based on a stepwise or
layer approach to significance testing. Sample means are
ordered from the smallest to the largest. The largest
difference, which involves means that are r = p steps apart,
is tested first at \alpha level of significance; if significant,
means that are r = p - 1 steps apart are tested at \alpha level
of significance and so on. The Newman-Keuls procedure provides an
r-mean significance level equal to \alpha for each group of
r ordered means, that is, the probability of falsely rejecting
the hypothesis that all means in an ordered group are equal to
\alpha. It follows that the concept of error rate applies
neither on an experimentwise nor on a per comparison basis-the
actual error rate falls somewhere between the two.
The Newman-Keuls procedure, like Tukey's procedure, requires
equal sample n's. However, in this algorithm, the procedure is
adapted to unequal sample sized which can lead to still
conservative test decisions.
It should be noted that the Newman-Keuls and Tukey procedures
require the same critical difference for the first comparison
that is tested. The Tukey procedure uses this critical difference
for all the remaining tests, whereas the Newman-Keuls procedure
reduces the size of the critical difference, depending on the
number of steps separating the ordered means. As a result, the
Newman-Keuls test is more powerful than Tukey's test. Remember,
however, that the Newman-Keuls procedure does not control the
experimentwise error rate at \alpha.
Value
A list containing:
adjPValues |
A numeric vector containing the adjusted pValues |
rejected |
A logical vector indicating which hypotheses are rejected |
statistics |
A numeric vector containing the test-statistics |
confIntervals |
A matrix containing only the estimates |
errorControl |
A Mutoss S4 class of type |
Author(s)
Frank Konietschke
References
Keuls M (1952). "The use of the studentized range in connection with an analysis of variance". Euphytica 1: 112-122
Examples
x = rnorm(50)
grp = c(rep(1:5,10))
dataframe <- data.frame(x,grp)
result <- snk(x~grp, data=dataframe, alpha=0.05,MSE=NULL, df=NULL, silent = TRUE)
result <- snk(x~grp, data=dataframe,alpha=0.05,MSE=NULL, df=NULL, silent = FALSE)
result <- snk(x~grp, data=dataframe,alpha=0.05,MSE=1, df=Inf, silent = FALSE) # known variance
result <- snk(x~grp, data=dataframe,alpha=0.05,MSE=1, df=1000, silent = FALSE) # known variance