爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: zhouyuanshen

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

  [复制链接]

新浪微博达人勋

发表于 2015-5-23 22:16:59 | 显示全部楼层
zhouyuanshen 发表于 2015-5-23 18:58
论坛的好友都在回复里面了,您可以浏览一下。大部分在第一页

可是看的不是很明白,我想知道你最后是不是用MK做的突变分析,还是用其他的方法做的,求回复
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-7-24 13:46:19 | 显示全部楼层
代码有问题啊。
第23行
应该是这样才对
U[i] <- (s[i]- (i*(i+1)/4))/sqrt(i*(i-1)*(2*i+5)/72)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-10-12 23:35:17 | 显示全部楼层
愚遥 发表于 2015-7-24 13:46
代码有问题啊。
第23行
应该是这样才对

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

魏凤英老师数据

魏凤英老师数据


自己的数据

自己的数据






再次谢谢您
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-10-12 23:36:05 | 显示全部楼层
高地不保 发表于 2015-5-23 22:16
可是看的不是很明白,我想知道你最后是不是用MK做的突变分析,还是用其他的方法做的,求回复

用了,但不是主要的结论。仅仅是辅助的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-18 08:58:04 | 显示全部楼层
超过一个突变点的时候,MK方法就不适用了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-18 21:43:55 | 显示全部楼层
正好需要,先研究下,不明白的请教大神
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-10-19 13:36:11 | 显示全部楼层
wk1984 发表于 2015-10-18 08:58
超过一个突变点的时候,MK方法就不适用了

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

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

但是如果是因为突变点的个数造成的不准确,那是不是可以分段做呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-10-19 13:36:28 | 显示全部楼层
polydandy 发表于 2015-10-18 21:43
正好需要,先研究下,不明白的请教大神

恩,相互学习。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-21 11:33:03 | 显示全部楼层
好复杂啊 {:5_217:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-10-28 10:27:02 | 显示全部楼层
过来学习一下
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表