MK趋势分析
Mann-Kendall(MK)趋势分析在气象、水文学中已得到广泛应用,该方法不要求变量遵从正太分布,故可以检验水文、气象序列。
计算原理
对于一个时间序列Xt = (x1,x2,…,xn),MK趋势检验的统计量S计算公式如下:
其中,sgn为检验函数:
然后计算方差Var以及Z统计量:
利用双边检验对数据进行检验,当-Z1-a<=Zs<=Z1-a时,即表明序列具有显著变化趋势,反之趋势不显著。当Z1-a给定的显著水平a=0.05时,临界值为1.96。
工作环境
Python3.6、numpy
准备工作
首先导入numpy包,并输入数据。
计算S值及Var
注意索引值和我们平时排列计数时的区别。
s_end = [