Kruskal Wallis 多重比较

Kruskal Wallis H 统计量是一个整体检验统计量，支持对所有总体中位数相等情况下的一般假设进行检验。通常情况下，调查者对于这个一般假设不太感兴趣，而会对比较各个组感兴趣。此宏可在非参数设置中执行多重比较。

下载宏

确保 Minitab 知道可在何处找到您下载的宏。选择文件 > 选项 > 常规。在宏位置下，浏览到您保存宏文件的位置。

重要信息

如果您使用较低版本的 Web 浏览器，则在您单击下载按钮时，此文件可能会在 Quicktime 中打开，因为 Quicktime 与 Minitab 宏使用相同的 .mac 文件扩展名。要保存此宏，请右键单击下载按钮并选择目标另存为。

下载 KrusMC.mac

必需输入

一列数值响应数据
一列对应因子水平

注意

您可以通过指定 UNSTACKED 子命令来使用非堆叠数据。

可选输入

UNSTACKED: 指定数据是否为非堆叠数据。
FALPHA: 用于指定所需的整体 alpha 水平（默认值是 .20）。
CONTROL C: 用于指定要控制的列 (C)。

注意

如果数据为堆叠数据，则控制组的响应不得与其他因子水平的响应在同一列中。控制组的响应必须在单独的列中。

运行宏

假设响应数据在 C1 中，因子水平在 C2 中。要运行此宏，请选择查看 > 命令行/历史记录，然后键入以下内容：

%KRUSMC C1 C2

单击运行。

输出

输出内容的第一部分将显示所生成的比较数 (k)

、整体 alpha (α)、Bonferroni 个体 alpha (β)，以及双侧临界 z 值。

下一部分显示了与这些差关联的标准化组均值秩差 (θ) 和 p 值。请注意，这些表格是对称的，所以在其顶部的三角形部分中有星号。同样，在沿表格对角线的位置处还有一些零，这是因为对角线表示组和其自身的比较（这些比较没有意义，在分析中并不考虑在内）。有关如何阅读此表格的示例如下：组 2 和组 4 之间的差值是多少？第三部分显示了中位数的符号置信区间。这些区间的置信水平按照整体 alpha 进行控制。由于这些区间按照全族的整体 alpha 进行控制，因此我们可以采用配对的方式比较这些区间（附录 2）。务必记住，希望得到的置信度可能无法在所有或部分区间上实现。此外，这种通过全族 alpha 实现的“整体”覆盖在样本数量不同的情况下并不准确，但通常还是能达到合理的近似程度。

最后一部分显示了如何才能算作“显著”的差（如果存在的话）。此部分中展示了 z 值、临界 z 值以及和 z 值关联的 p 值。

此图显示了标准化非绝对组均值秩差。在此图中，我们不仅可以看到组差的大小，还可以看到组差的方向，因此非常有用。它还显示了正临界 z 值和负临界 z 值，因此您可以查看差是否“显著”。

数据中的结

如果数据中存在结，则会计算出无偏常量（或校正因子）2。H 统计量随后将调整为 3。标准差 (ξ) 也将根据校正因子 4 进行调整。输出将同时显示调整后和未调整的表格。但是，只会给出调整后的表格的 p 值，因为这些才是应该使用的表格。显示未调整表格的主要原因在于显示结对于 z 值产生的效果。如果结非常广泛，则应当质疑数据的有效性，因为这些检验假定分布是连续的。通常情况下，结只会对您的结论造成很小的影响，或者完全没有影响。

鸣谢

非常感谢 Tom Hettmansperger 博士（宾夕法尼亚州立大学），他对此宏进行了审阅，并针对这项工作进行了多次讨论，付出了宝贵时间和热心帮助。另外，也感谢 Nicholas Bolgiano 先生和 Mike Delozier 先生 (Minitab, Inc.)