octave mann-whitney/u_test p 值混淆

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 是有意义的，因为在该 p 值处，原假设被拒绝，而替代方案，即 P(x>y) != 1/2 将被接受。

但是，第二个 u_test 表明原假设再次被拒绝，因此替代 P(x>y) > 1/2 被接受，这对我来说没有任何意义。

我哪里做错了？

Answer 1

我查看了u_test.m，该功能似乎有几个错误，感谢您发现并报告它。

你是对的，作者混淆了测试的顺序。根据他的定义，应该是相反的（即你得到的结果应该是'<'）。

此外，z 的计算似乎也包含错误，这似乎在使用在线计算器时通过不同的 p 值得到证实。

请随时在 octave bug tracker 上将此错误报告为错误（并将生成的错误页面链接回 cmets 中的此处；否则，如果您愿意，我会报告错误）。

与此同时，我将开发一个补丁。再次感谢。