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

看云淡风清

莫莫丶格格子° 发表于 2017-1-8 10:12
研究之后，个人认为，魏老师的程序是正确的，源程序中的代码为
DO 10 I=2,N
        M(I)=0

关键就是xi＞xj还是xi≥xj问题啊，方法介绍和程序中的不一样，到底是怎样的呢？

zfb0817

有没有确定的答案

何佳田力

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

酒儿

下载了，改天试试看

lu'an

20楼正解，@18kingfzh。
以前发过excel版，今日复习，重看了一遍。
代码修改如下：（转录20楼）
原代码：
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  
或者直接改43行.LT.为.LE.也可（有网友已提到）。
书中计算公式无误，但图有误。
https://www.idosi.org/wasj/wasj28(7)13/2.pdf P906 另处佐证。
另附改后for文件及书中数据与改后输出结果。
运行时输入91 1900
即：91个数据点，1900年为第一年。

小英

郭衡 发表于 2012-2-18 16:32
我们老师上课说是书上的印刷错误，正确的如你所说：改为xi>=xj,r=1;同时xi

但做出的图好像和书上也不一样呢

小英

hunteria 发表于 2012-2-18 19:13
那正确的算法应该是>=还是>啊，我看其他文献上也是>号，不是>=号

对啊，我也想问

无剑道

{:5_275:}{:5_275:}{:5_275:}{:5_275:}{:5_275:}

angfou

原来真的是有错误呀

xueqixiang

把1:i;改成1:（i-1）;同时用xi>=xj；然后，好像和教材上差不多了