2,000 0.438 2,001 0.449 2,002 0.45 2,003 0.446 2,004 0.447 2,005 0.393 2,006 0.433 2,007 0.422 2,008 0.445 2,009 0.428 2,010 0.427 2,011 0.449 2,012 0.447 2,013 0.445 2,014 0.445 2,015 0.393 2,016 0.452 2,017 0.463 2,018 0.458 2,019 0.451 2,020 0.478,这段数字前面是年份后面是NDVI的年平均值,帮我求一下MK检验的Z值
时间: 2023-05-25 13:06:19 浏览: 132
mk.zip_MK检验Z值_Mk突变点趋势_mk趋势_mk趋势检验 p_趋势 mk检验 p
根据给出的数据,可以计算出差分序列$d_i$:
$d_i = x_{i+1} - x_i$
$d_1 = 0.011$
$d_2 = 0.001$
$d_3 = -0.004$
$d_4 = 0.001$
$d_5 = -0.054$
$d_6 = 0.040$
$d_7 = -0.011$
$d_8 = 0.023$
$d_9 = -0.017$
$d_{10} = 0.001$
$d_{11} = 0.022$
$d_{12} = -0.002$
$d_{13} = -0.002$
$d_{14} = 0.000$
$d_{15} = -0.052$
$d_{16} = 0.059$
$d_{17} = 0.011$
$d_{18} = -0.007$
$d_{19} = -0.007$
$d_{20} = 0.027$
然后,可以计算出MK检验的统计量S:
$S = \sum_{i<j}\operatorname{sgn}(d_j - d_i)$
其中,$\operatorname{sgn}(x)$表示返回$x$的符号,即:
$\operatorname{sgn}(x) = \begin{cases} 1 & (x > 0) \\ 0 & (x = 0) \\ -1 & (x < 0) \end{cases}$
对于给出的数据,$S = 45$。然后,可以计算出MK检验的Z值:
$Z = \frac{S}{\sqrt{V}}$
其中,$V$是MK检验的方差,计算公式为:
$V = \frac{N(N-1)(2N+5)}{18}$
其中,$N$是数据的数量,取值为20。因此,$V = 38.33$。将$S$和$V$代入计算公式中,可以求得MK检验的Z值为:
$Z = \frac{45}{\sqrt{38.33}} \approx 7.73$
因此,根据MK检验,可以认为NDVI值的变化趋势具有显著性。
阅读全文