Brunner-Munzel 検定 ッテ,もはや,参照する価値もない記事です。
まったくのところ,ストーカーに成り下がっています
さて,brunnermunzel を 1.3.3 にもどし,さらに permutationtest を t ではなく p に基づくように変更する。つまり,以下のようjに
calc_statistics.f において
stat = my - mx
bm_test.f において
nvx = n1 * vx; nvy = n2 * vy
以上2箇所変更により
brunnermunzel.test は 1.3.3 に戻り,
> brunnermunzel.test(x,y)
Brunner-Munzel Test
data: x and y
Brunner-Munzel Test Statistic = 3.1375, df = 17.683, p-value = 0.005786
95 percent confidence interval:
0.5952169 0.9827052
sample estimates:
P(X<Y)+.5*P(X=Y)
0.788961
brunnermunzel.permutation.test は t ではなく p に基づいて並べ替え検定を行うようになった。
> brunnermunzel.permutation.test(x, y)
permuted Brunner-Munzel Test
data: x and y
p-value = 0.00669
これでよいのかどうか?