请围观：关于Mann-Kendall检验结果的分析(附R分析代码)

高地不保 · 发表于 2015-5-23 22:16:59

zhouyuanshen 发表于 2015-5-23 18:58
论坛的好友都在回复里面了，您可以浏览一下。大部分在第一页

可是看的不是很明白，我想知道你最后是不是用MK做的突变分析，还是用其他的方法做的，求回复

愚遥 · 发表于 2015-7-24 13:46:19

代码有问题啊。
第23行
应该是这样才对
U[i] <- (s[i]- (i*(i＋1)/4))/sqrt(i*(i-1)*(2*i+5)/72)

zhouyuanshen · 发表于 2015-10-12 23:35:17

愚遥发表于 2015-7-24 13:46
代码有问题啊。
第23行
应该是这样才对

谢谢您的提醒，我对着魏老师的书看了一遍，好像是输错了。在这里给大家道歉了，希望对大家的影响不大。
我试了一下，用魏老师的数据和我自己的数据重新运行了，错误代码的结果似乎比较保守。。。。。下面两张图，左边是经过您提醒后修改后的代码，右边是原来的错误代码

魏凤英老师数据

自己的数据

再次谢谢您

zhouyuanshen · 发表于 2015-10-12 23:36:05

高地不保发表于 2015-5-23 22:16
可是看的不是很明白，我想知道你最后是不是用MK做的突变分析，还是用其他的方法做的，求回复

用了，但不是主要的结论。仅仅是辅助的

wk1984 · 发表于 2015-10-18 08:58:04

超过一个突变点的时候，MK方法就不适用了

polydandy · 发表于 2015-10-18 21:43:55

正好需要，先研究下，不明白的请教大神

zhouyuanshen · 发表于 2015-10-19 13:36:11

wk1984 发表于 2015-10-18 08:58
超过一个突变点的时候，MK方法就不适用了

这个方法没有深究过，因为我主要不做突变。

好像大家都不完全依靠这个方法，正如您说的，可能不适用。

但是如果是因为突变点的个数造成的不准确，那是不是可以分段做呢

zhouyuanshen · 发表于 2015-10-19 13:36:28

polydandy 发表于 2015-10-18 21:43
正好需要，先研究下，不明白的请教大神

恩，相互学习。

小预报员 · 发表于 2015-10-21 11:33:03

好复杂啊

{:5_217:}

double001 · 发表于 2015-10-28 10:27:02

过来学习一下

		自动登录	找回密码
密码			立即注册