关于魏凤英书中MK突变检验发现的问题，大家来讨论一下！

hunteria

今天自己手动进行mk检验，发现了一个问题！
数据就是魏凤英书中举例用的上海1900-1990平均温 data.txt (975 Bytes, 下载次数: 108)

2012-2-18 10:45 上传

点击文件名下载附件
但是我按书中的方法 算出来就是跟书上的图不一样！
倒腾了两天，原来如果把上式中xi>xj改成xi>=xj就跟书上的图一模一样了，难道是魏老师算错了？
求大家论证！

18kingfzh

我来给大家权威介绍下：
情况是这样的：
魏老前辈程序当中，在计算每个样本的累计数有所偏差，其采用的计算方法为：
xi>=xj(1≤j<i)

而M—k检验的原理中关于累计数的计算与此有所差异，尤其是在序列当中存在相同数时，此差异就尤为明显，而且此差异导致的突变点的位置也相差较大，正确的计算方法应该是(修正版采用此计算方法)：
xi>xj(1≤j≤i)

例如序列1，2，2
原版算法得出第3个样本的累计数为2，而正确算法当中得出的累计数为1，当这种情况在长序列中反复出现时，结果就会相差较大。

代码修改如下：
原代码：
42至44行：
        DO 20 J=1,I-1
        IF(Y(I).LT.Y(J))GOTO 20
        M(I)=M(I)+1
应该改为：
        DO 20 J=1,I
        IF(Y(I).GT.Y(J))then       
            M(I)=M(I)+1                       
        end if               
至于j取到I和I-1，这个对结果没从我什么影响，因为当样本和自身比较时，肯定不会满足大于条件，只不过M-K的原理中j是取到I的，我们就统一取到I吧。
PS：因为魏老前辈的程序是错的，所以书上的结果图也自然是错的。

蚊蚊1989

我想用excel算得不得啊？那个程序搞不懂，也算不出来！！！求助！！！！~~~~~

hunteria

什么情况，居然获得小穷神卡，本来就没有金币了

blue00ghost

这个问题去年发现过，但没找到其不一致原因。按照楼主的方法，依然做出的曲线不完全一样，怀疑顺序时>，逆序时>=，可仍然不完全一样。注，所用程序为用另一个fortune程序验证过的vb程序。

hunteria

自己先顶一下，找个沙发等大家来讨论

hunteria

斑竹在不，给个响啊

善人/jw

这个我记得以前也有人说过，不过没算过~~

郭衡

我们老师上课说是书上的印刷错误，正确的如你所说：改为xi>=xj,r=1;同时xi<xj,r=0

言深深

你说的这个是很有可能的，魏老师可能也没有全部自己检测过！
刚翻看了第二版中，貌似还没有改正过来，可以考虑给魏老师发封邮件咨询一下

hunteria

郭衡 发表于 2012-2-18 16:32

                               
登录/注册后可看大图

我们老师上课说是书上的印刷错误，正确的如你所说：改为xi>=xj,r=1;同时xi

那正确的算法应该是>=还是>啊，我看其他文献上也是>号，不是>=号

hunteria

言深深 发表于 2012-2-18 17:03

                               
登录/注册后可看大图

你说的这个是很有可能的，魏老师可能也没有全部自己检测过！
刚翻看了第二版中，貌似还没有改正过来，可以 ...

要不版主代表我们给魏老师去封信问问？

言深深

hunteria 发表于 2012-2-18 19:14

                               
登录/注册后可看大图

要不版主代表我们给魏老师去封信问问？

哈哈，不敢，下次看到魏老师要不说说？