R kruskal.test Kruskal-Wallis 秩和检验

说明

执行 Kruskal-Wallis 秩和测试。

用法

kruskal.test(x, ...)

## Default S3 method:
kruskal.test(x, g, ...)

## S3 method for class 'formula'
kruskal.test(formula, data, subset, na.action, ...)

参数

x	数据值的数值向量，或数值数据向量的列表。列表中的非数字元素将被强制，并带有警告。
g	一个向量或因子对象，为 x 的相应元素提供组。如果 x 是列表，则忽略并显示警告。
formula	response ~ group 形式的公式，其中 response 给出数据值，group 给出相应组的向量或因子。
data	包含公式 formula 中的变量的可选矩阵或 DataFrame (或类似的：请参阅 model.frame )。默认情况下，变量取自environment(formula)。
subset	一个可选向量，指定要使用的观测子集。
na.action	一个函数，指示当数据包含 NA 时应该发生什么。默认为 getOption("na.action") 。
...	要传递给方法或从方法传递的更多参数。

细节

kruskal.test 对每个组(样本)中 x 的分布位置参数相同的零值执行 Kruskal-Wallis 秩和检验。另一种选择是它们至少在一个方面有所不同。

如果x是一个列表，则其元素将作为要比较的样本，因此必须是数值数据向量。在这种情况下，g被忽略，可以简单地使用kruskal.test(x)来执行测试。如果示例尚未包含在列表中，请使用 kruskal.test(list(x, ...)) 。

否则， x 必须是数值数据向量，并且 g 必须是与 x 长度相同的向量或因子对象，为 x 的相应元素提供组。

值

类"htest" 的列表包含以下组件：

statistic	Kruskal-Wallis 排名总和统计。
parameter	检验统计量的近似卡方分布的自由度。
p.value	检验的 p 值。
method	字符串 "Kruskal-Wallis rank sum test" 。
data.name	给出数据名称的字符串。

例子

## Hollander & Wolfe (1973), 116.
## Mucociliary efficiency from the rate of removal of dust in normal
##  subjects, subjects with obstructive airway disease, and subjects
##  with asbestosis.
x <- c(2.9, 3.0, 2.5, 2.6, 3.2) # normal subjects
y <- c(3.8, 2.7, 4.0, 2.4)      # with obstructive airway disease
z <- c(2.8, 3.4, 3.7, 2.2, 2.0) # with asbestosis
kruskal.test(list(x, y, z))
## Equivalently,
x <- c(x, y, z)
g <- factor(rep(1:3, c(5, 4, 5)),
            labels = c("Normal subjects",
                       "Subjects with obstructive airway disease",
                       "Subjects with asbestosis"))
kruskal.test(x, g)

## Formula interface.
require(graphics)
boxplot(Ozone ~ Month, data = airquality)
kruskal.test(Ozone ~ Month, data = airquality)

参考

Myles Hollander and Douglas A. Wolfe (1973), Nonparametric Statistical Methods. New York: John Wiley & Sons. Pages 115-120.

也可以看看

Wilcoxon 秩和检验 (wilcox.test) 作为两个样本的特殊情况； lm 与anova 一起用于在正态性假设下执行one-way 位置分析；以学生 t 检验 (t.test) 作为两个样本的特殊情况。

coin 包中的 wilcox_test 用于精确、渐近和蒙特卡罗条件 p 值，包括存在联系的情况。