scipy.stats: Brunner-Munzel 検定 `brunnermunzel`

Brunner-Munzel 検定を行う。

brunnermunzel(x, y, alternative='two-sided', distribution='t', nan_policy='propagate')

distribution　$p$ 値の計算に $t$ 分布を使う（デフォルト）とき 't'，標準正規分布を使うとき 'normal' を指定する。

from scipy.stats import brunnermunzel
import numpy as np

x = np.array([1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 4, 1, 1])
y = np.array([3, 3, 4, 3, 1, 2, 3, 1, 1, 5, 4])
brunnermunzel(x, y)

BrunnerMunzelResult(statistic=3.1374674823029505, pvalue=0.005786208666151538)

brunnermunzel(x, y, distribution='normal')

BrunnerMunzelResult(statistic=3.1374674823029505, pvalue=0.0017041417600383024)