Mann-Whitney 检验能否用于 Kruskal-Wallis 之后的事后比较?

机器算法验证 假设检验 事后 wilcoxon-mann-whitney 检验 克鲁斯卡尔-沃利斯测试 邓恩测试
2022-03-27 03:36:04

我有一个模拟,动物被放置在一个敌对的环境中,并定时查看它使用某种生存方法可以生存多长时间。它可以使用三种方法来生存。我使用每种生存方法对动物进行了 300 次模拟。所有模拟都在相同的环境中进行,但存在一些随机性,因此每次都不同。我计算动物在每次模拟中存活的秒数。活得越久越好。我的数据如下所示:

Approach 1, Approach 2, Approach 2
45,79,38
48,32,24
85,108,44
... 300 rows of these

我不确定在此之后我所做的一切,所以如果我做的事情愚蠢和错误,请告诉我。我试图找出使用特定方法的寿命是否存在统计差异。

我对每个样本进行了夏皮罗测试,它们返回的 p 值很小,所以我相信数据没有标准化。

行上的数据彼此没有关系。用于每个模拟的随机种子是不同的。因此,我相信数据没有配对。

由于数据未标准化、未配对且样本超过两个,因此我进行了 Kruskal Wallis 检验,其 p 值为 0.048。然后我转而选择了 Mann Whitney。真的不确定是否应该在这里使用 Mann Whitney。

我通过执行曼惠特尼检验将每种生存方法与其他方法进行了比较,即{(方法 1,方法 2),(方法 1,方法 3),(方法 2,方法 3)}。使用双尾检验未发现这对(方法 2,方法 3)之间的统计显着性,但使用单尾检验发现显着性差异。

问题:

  1. 我不知道这样使用 Mann Whitney 是否有意义。
  2. 我不知道我是否应该使用一尾或二尾 Mann Whitney。
4个回答

不,您不应该在这种情况下使用 Mann-Whitney检验。U

原因如下:Dunn 检验是一种适当的事后检验*在拒绝 Kruskal-Wallis 检验之后。如果从拒绝 Kruskal-Wallis 转向执行普通的成对秩和(即 Wilcoxon 或 Mann-Whitney)检验,则得到两个问题:(1)用于成对秩和检验的秩是不是 Kruskal-Wallis 检验使用的等级;(2) 秩和检验不使用 Kruskal-Wallis 原假设隐含的汇总方差。邓恩的测试没有这些问题

拒绝 Kruskal-Wallis 检验(已针对多重比较进行了调整)后的事后检验可能无法拒绝针对给定全族错误率或给定与综合检验的给定相对应的错误发现率的所有成对检验,只是与任何其他多重比较综合/事后测试场景一样。α

除非你有理由相信一个群体的生存时间比另一个群体的生存时间长或短否则你应该使用双边检验。

Dunn 的测试可以在 Stata 中使用dunntest (type net describe dunntest, from(https://www.alexisdinno.com/stata)) 执行,在 R 中使用dunn.test包。

另外,我想知道您是否可以采用生存分析方法来评估动物是否以及何时根据不同的条件死亡?


*一些不太知名的事后成对测试遵循被拒绝的 Kruskal-Wallis,包括 Conover-Iman(类似于 Dunn,但基于t分布,而不是z分布,在conovertest中为 Stata 实现,对于conover.test包中的 R)和 Dwass-Steel-Citchlow-Fligner 测试。

Kruskal-Wallis/Wilcoxon 的一个统一概括是比例优势模型,该模型允许与优势比的逐点或同时置信区间进行一般对比。这是在我的 Rrmsormcontrast.rms函数中实现的。

您也可以使用 Conover 之后的临界差值或 Schaich 和 Hamerle 之后的临界差值。前者更自由,而后者精确但缺乏一点权力。这两种方法都在我的网站brightstat.com 上进行了说明,brightstat 的 webapp 还可以让您计算这些关键差异并立即执行事后测试。 Brightstat.com 上的 Kruskal-Wallis

如果您使用的是 SPSS,请使用 Bonferroni 校正(p 值除以组数)进行事后 Mann-Whitney。