octave mann-whitney/u_test p-value confusion

Question

Я нахожу результат теста mann-whitney запутанным, и документация gnu мне не помогла. https://www.gnu.org/software/octave/doc/interpreter/Tests.html

Вот простой пример, который я попробовал:

octave:1> x=[1,1,1,1,1] 
x = 

    1 1 1 1 1 

octave:2> y=[2,2,2,2,2,2] 
y = 

    2 2 2 2 2 2 

octave:3> [p,z]=u_test(x, y, "<>") 
p = 0.0061699 
z = -2.7386 
octave:4> [p,z]=u_test(x, y, ">") 
p = 0.0030849 
z = -2.7386 

Первый u_test имеет смысл, так как при этом значение р, нулевая гипотеза отвергается, и альтернатива, которая является то, что Р (х> у) ! = 1/2 будет принято.

Однако второй u_test предполагает, что нулевая гипотеза снова отвергается, поэтому принимается альтернатива P (x> y)> 1/2, что для меня не имеет никакого смысла.

Где я ошибся?

Answer 1

Я посмотрел u_test.m, и функция, похоже, содержит несколько ошибок, благодарю вас за обнаружение и отчетность.

Вы правы, автор перепутал порядок испытаний. Это должно быть наоборот в соответствии с его определением (т. Е. Результат, который вы получаете, должен быть равен '<').

Кроме того, кажется, что для вычисления z также имеется ошибка, это, похоже, подтверждается другим значением p при использовании онлайн-калькуляторов.

Не стесняйтесь сообщать об этом как об ошибке в octave bug tracker (и, пожалуйста, ссылку на полученную страницу с ошибкой здесь, в комментариях, в противном случае я сообщу об ошибке, если вы захотите).

Тем временем я буду работать над патчем. Еще раз спасибо.