前一篇文章我们介绍了秩和检验,这里我们介绍另一种秩和检验方法,Wilcoxon signed rank sum test
文章目录
Wilcoxon signed rank sum test包含两种检验方式,数据成对的情况和数据是独立观察的情形。下面分开讨论:
对于成对的数据,Wilcoxon signed rank sum test检验的流程如下:
- 声明空假设: H 0 H_0 H0,这里就是两个数据之差的中值为0;
- 计算数据对的差 d i = x i − y i d_i = x_i - y_i di=xi−yi
- 对 d i d_i di取绝对值之后进行排序
- 计算正数的秩 W + W^+ W+与 W − W^- W−,
对于单次的观测
- 声明空假设中值为M;
- 计算 d i = x i − M d_i = x_i - M di=xi−M
- 其余步骤类似
空假设下,我们认为
d
i
d_i
di应该是在0两侧对称分布的。
4. 计算
W
=
min
(
W
+
,
W
−
)
W=\min(W^+, W^-)
W=min(W+,W−)作为统计量
5. 查表计算p值,双侧检验的p值是单侧检验p值的两倍;
如果数据足够多,可以利用正态分布进行近似
对于数据集:
A: 2.0, 3.6, 2.6, 2.6, 7.3, 3.4, 14.9, 6.6, 2.3, 2.0, 6.8, 8.5
B: 3.5, 5.7, 2.9, 2.4, 9.9, 3.3, 16.7, 6.0, 3.8, 4.0, 9.1, 20.9
利用scipy.stats
包进行秩和检验的结果如下:
from scipy import stats
_, p = stats.wilcoxon(a, b) #> p=0.012