关于M-K突变检验中累计数的疑问？？

虾米的尾巴

      在论坛里看了关于M-K检验的帖子，现在对其累计数还不太清楚：
      在这篇帖子中http://bbs.06climate.com/forum.php?mod=viewthread&tid=8356说到魏凤英老师的书中关于累计数是Xi>Xj,但程序中是Xi>=Xj,因而导致出图结果不同。在查阅文献后，关于累计数的都是Xi>Xj,下面这是符淙斌院士1992年文章的图片：

      其中关于累计数的解释是Xi>Xj,另外上面提到的帖子中还说到魏凤英老师书中E（Sk)=n*(n+1)/4也是应该是错的，应该是k*(k-1)/4,并且用n还会引起歧义，使这项的计算得出的是定值（我看到论坛里的一位童鞋就说犯了这个错误）。在查阅一些文献后，都是k*(k-1)。
      那么是不是应该采用Xi>Xj而不是Xi>=Xj呢？并且用E（Sk)=k*(k-1)呢？不知有没有哪位童鞋能给个明确说法？

虾米的尾巴

有没有人？别沉了

gradfine

应该采用Xi>=Xj,并且用E（Sk)=k*(k-1)/4.我用这个实例做出的结果与魏凤英老师的图一致。

虾米的尾巴

gradfine 发表于 2014-2-13 09:41

                               
登录/注册后可看大图

应该采用Xi>=Xj,并且用E（Sk)=k*(k-1)/4.我用这个实例做出的结果与魏凤英老师的图一致。

用k*(k-1)/4应该是没问题了。但对于累计数魏凤英老师的书中及我找到的文献中都写的是Xi>Xj，没有理由在计算编程的时候又采用Xi>=Xj,而且有童鞋发帖说这两种方法得出的突变点相差10年，这是我疑惑的地方。

Nonego

应该采用Xi>Xj而不是Xi>=Xj，在http://bbs.06climate.com/forum.php?mod=viewthread&tid=15390中，你可以看到，国外学者在碰到有相同变量的时候引用了“结（tie）”的概念，你可以好好看看链接中的内容，虽然讲的是趋势分析，不是突变，但是对于累计量的确定道理是一样的，而E（Sk)=k*(k-1)应该是正确的，我看了好多文献也都是减！希望对你有帮助。

Nonego

应该采用Xi>Xj而不是Xi>=Xj，在http://bbs.06climate.com/forum.php?mod=viewthread&tid=15390中，你可以看到，国外学者在碰到有相同变量的时候引用了“结（tie）”的概念，你可以好好看看链接中的内容，虽然讲的是趋势分析，不是突变，但是对于累计量的确定道理是一样的，而E（Sk)=k*(k-1)应该是正确的，我看了好多文献也都是减！希望对你有帮助。

虾米的尾巴

Nonego 发表于 2014-2-13 16:37

                               
登录/注册后可看大图

应该采用Xi>Xj而不是Xi>=Xj，在http://bbs.06climate.com/forum.php?mod=viewthread&tid=15390中，你可以看 ...

谢谢！综合看来采用Xi>Xj应该是对的。

小小人

是不是啊，正在思考中！