Wilcoxon 秩和 / Mann-Whitney-U 检验的精确统计量答案

【问题标题】：Exact statistic for Wilcoxon rank-sum / Man-Whitney-U testWilcoxon 秩和 / Mann-Whitney-U 检验的精确统计量
【发布时间】：2021-03-17 01:32:31
【问题描述】：

我想计算一个小样本 (n=32) 的秩和检验统计量。 scipy.stats.mannwhitneyu 对检验统计量和相关的 p 值使用正态近似值（请参阅 here）。对于小样本，这可能会产生误导。

在 python 中有没有办法计算 exact 测试统计数据和/或相关的 p 值（如 here 所述，当然不是在 excel 中）？

[R 有这样的选择，请参阅here，但我更愿意留在 python 中。]

【问题讨论】：

Pingouin 统计包或许

标签： python statistics

【解决方案1】：

看起来pingouin 部署了简单的蛮力方法。在小样本上，它确实为您提供了不同的 p 值结果

代码，Python 3.9 Win10 x64

import numpy as np
import pingouin as pg

import scipy.stats as st

x = [20, 22, 19, 20, 22, 18, 24, 20, 19, 24, 26, 13]
y = [38, 37, 33, 29, 14, 12, 20, 22, 17, 25, 26, 16]

print(pg.wilcoxon(x, y, tail='two-sided'))

st.wilcoxon(x, y, zero_method='wilcox', alternative='two-sided')

【讨论】：