mann-kendall突变检验奇怪的结果

cashia371081

mann-kendall突变检验法做的突变检验。从原始资料看，数据本身是逐年增加的（增加趋势稳定从0.01一直增加到1.1，由于数据比较保密，不敢放上来请大家见谅）。哪位高手帮忙分析一下这个结果。uk与ub的相交点不在临界线之间。而且统计值uk、ub普遍偏大,我看很多文献都是在-4~6之间。请大神帮忙解答下。程序是用的网上找到的一个程序（针对魏老师书中的错误，改程序作了改动）。
C THIS IS A PROGRAM FOR DETECTING ABRUPT CLIMATIC CHANGE
C BY USING MANN-KENDALL TECHNIQUE
PROGRAM MK
DIMENSION Y(1000),YY(1000),U(1000),UF(1000),UB(1000),
     & M(1000),MD(1000)
WRITE(*,10)
  10 FORMAT(2X,'N=?,NYEAR=?')
READ(*,*)N,NYEAR
C ***************************************************
C * N:     SAMPLE SIZE                              *
C * NYEAR: FIRST YEAR OF THE TIME SERIES            *
C * Y(N):  ORIGINAL TIME SERIES                     *
C * UF(N): ORIGINAL SERIES OF U(LN)                 *
C * UB(N): COUNTER SERIES OF U(LN)                  *
C * A,B:   CRITICAL VALUE 1.96 AND -1.96            *
C ***************************************************   
OPEN(2,FILE='G:\发生面积率.txt')
READ(2,*)(Y(I),I=1,N)
CALL SMK(Y,M,MD,UF,N)
DO 20 I=1,N
  20 YY(I)=Y(N+1-I)
CALL SMK(YY,M,MD,U,N)
DO 30 I=1,N
  30 UB(I)=-U(N+1-I)
OPEN(3,FILE='G:\area_result_recversion.txt',STATUS='NEW')
A=1.96
B=-1.96
DO 40 I=1,N
WRITE(3,50)NYEAR+I-1,UF(I),UB(I),A,B
  50 FORMAT(1X,I4,4F8.2)
  40 CONTINUE
CLOSE(3)
STOP
END
C***********************************************************
SUBROUTINE SMK(Y,M,MD,U,N)
DIMENSION Y(N),M(N),MD(N),U(N)
M(1)=0
DO 10 I=2,N
M(I)=0
MD(I)=0
DO 20 J=1,I     !将原版里的i-1改为i;
IF(Y(I).GT.Y(J))then  !将原版里的LT情况下直接跳到下一循环改
  M(I)=M(I)+1     !为GT情况执行累计数操作，特别强调：
end if       !这两种情况在序列中存在相同数时是不一样的。
  20 CONTINUE
MD(I)=MD(I-1)+M(I)
  10 CONTINUE
U(1)=0.0
DO 30 I=2,N
E=I*(I-1)/4.00
VAR=I*(I-1)*(2*I+5)/72.00
U(I)=(MD(I)-E)/SQRT(VAR)
  30 CONTINUE
RETURN
END


言深深

这个结果是可信的，表面你的原数据序列是持续增加的吧···你不妨给出你的原序列

nuister

临界线是没问题的，好多都是-2~2，包括书上的例子，但是你的交点不在临界线之间，就不能判定它是突变点了。

cashia371081

言深深 发表于 2013-5-13 08:33

                               
登录/注册后可看大图

这个结果是可信的，表面你的原数据序列是持续增加的吧···你不妨给出你的原序列

可是交点不在临界线之间，能说交点是突变点吗

吴艳锋

没有突变点 正序列和反序列的交点不在置信区间线条内

cashia371081

吴艳锋 发表于 2013-5-13 11:33

                               
登录/注册后可看大图

没有突变点 正序列和反序列的交点不在置信区间线条内

魏老师书里说当uk或ub超过临界线时，表明上升或下降趋势显著。那在这张图上可可以说自1959年以后增加趋势明显吗

吴艳锋

cashia371081 发表于 2013-5-13 11:50

                               
登录/注册后可看大图

魏老师书里说当uk或ub超过临界线时，表明上升或下降趋势显著。那在这张图上可可以说自1959年以后增加趋势 ...

有突变点才开始有升降趋势，举个例子，假设你做的楼主的图中纵轴的+2改为-2，纵轴的+10改为+2，蓝色和红色的曲线分布不变，就可以说有突变点，1975年（看不太清，大概是吧）为突变点，1975年后气候开始呈下降趋势，在2003年达到极显著水平（p〈0.01或者Zc小于-1.98），气温下降趋势明显。

言深深

cashia371081 发表于 2013-5-13 10:57

                               
登录/注册后可看大图

可是交点不在临界线之间，能说交点是突变点吗

临界线就是一个显著性判断标准，有交点就可以认为有“突变”，但是可能没那么显著，再者就是这样的突变本来就不是特别的准确。你这个主要就是看看趋势吧···

cashia371081

言深深 发表于 2013-5-14 08:52

                               
登录/注册后可看大图

你好歹用源数据画图噻···光光的数据如何看？

第一次发帖啦，没经验。。。

cashia371081

言深深 发表于 2013-5-13 14:42

                               
登录/注册后可看大图

临界线就是一个显著性判断标准，有交点就可以认为有“突变”，但是可能没那么显著，再者就是这样的突变本 ...

我将数据上传了，那这个突变点可以采用吗