24小时热门版块排行榜    

查看: 17087  |  回复: 18
本帖产生 1 个 博学EPI ,点击这里进行查看

8711193629

金虫 (正式写手)


[交流] 关于Mann-kendall检验方法的计算步骤

各位大侠,帮帮忙吧,哪位会计算Mann-kendall趋势检验的步骤啊,虽然我也看了魏凤英老师编著的那本书及各种文献上的关于计算mann-kendall的计算步骤,可是还是计算不出,配套的Matlab程序一直没有找出,笔算了大半天结果却不对,哪位高手指点一下啊?

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

北拳爸爸

新虫 (初入文坛)


★ ★
8711193629: 金币+2, 博学EPI+1 2012-06-11 08:15:06
引用回帖:
6楼: Originally posted by 8711193629 at 2012-05-13 16:33:30
我是不盼望用程序做了,根本就没有正确的,于是我想Excel手工计算,但还是遇到了问题,感觉魏凤英老师那本《现代气候统计诊断与预测技术》关于MK检验的公式出了好急促错误,都不知到底是哪个对,查阅了不少文献,公 ...

那本书上5.4.3的那个E(sk)吗?一般都是用“-”,不过这个不影响你理解啦~
9楼2012-06-08 21:58:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
2楼: Originally posted by neu234 at 2012-05-09 21:46:11:

不会么?
3楼2012-05-09 21:55:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yangguang788

银虫 (正式写手)



8711193629(金币+1): 谢谢参与
这段时间也想做下这个,同盼交流~
4楼2012-05-10 10:00:25
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yangguang788

银虫 (正式写手)


matlab程序倒是记得在网上见到过,在网上找一找应该能发现的~
5楼2012-05-10 10:01:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
4楼: Originally posted by yangguang788 at 2012-05-10 10:00:25:
这段时间也想做下这个,同盼交流~

我是不盼望用程序做了,根本就没有正确的,于是我想Excel手工计算,但还是遇到了问题,感觉魏凤英老师那本《现代气候统计诊断与预测技术》关于MK检验的公式出了好急促错误,都不知到底是哪个对,查阅了不少文献,公式不一致,鱼的写“+“,有的写”-“,到底信哪个么,???????????
6楼2012-05-13 16:33:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
5楼: Originally posted by yangguang788 at 2012-05-10 10:01:02:
matlab程序倒是记得在网上见到过,在网上找一找应该能发现的~

是啊,我见过,只见过一个,已经下载了,可是套用那个程序,得出的图和书上的图差别不少,哎,该怎么办呢,自己又不是Matlab高手
7楼2012-05-13 16:34:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

北拳爸爸

新虫 (初入文坛)



8711193629(金币+1): 谢谢参与
M-K方法我们专业基本上每个研究生都会哦~一般是用Fortran程序来实现的。你哪里没明白呢?建议先摘10个数字的序列来对着步骤试验一下~这样简单一些方便你手算。Excel还是太麻烦了点。。。
8楼2012-06-08 21:56:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
9楼: Originally posted by 北拳爸爸 at 2012-06-08 21:58:05
那本书上5.4.3的那个E(sk)吗?一般都是用“-”,不过这个不影响你理解啦~...

我这倒是有fortran程序,没接触过,根本不会用啊,我这就30个数据想要做出来怎么那么难
10楼2012-06-11 08:16:25
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
8楼: Originally posted by 北拳爸爸 at 2012-06-08 21:56:00
M-K方法我们专业基本上每个研究生都会哦~一般是用Fortran程序来实现的。你哪里没明白呢?建议先摘10个数字的序列来对着步骤试验一下~这样简单一些方便你手算。Excel还是太麻烦了点。。。

能否加QQ帮忙啊,我的号544215401.。。
11楼2012-06-11 08:20:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ljlyhl0926

铁虫 (小有名气)


最近也想做这方面的,在网上看到两个比较完整的源代码,一个是http://www.qddayu.cn/information.asp?tid=922,另一个是http://blog.sina.com.cn/s/blog_4b4254430100c5sa.html。都可以用,但是不知道是否有错误。方差计算的那一步没弄明白,一个是VarS=n*(n-1)*(2*n+5)/18;另一个是Var=i*(i-1)*(2*i+5)/72; 感觉都不是很符合方差的定义,希望高手指点一二。

[ Last edited by ljlyhl0926 on 2012-7-11 at 20:19 ]
12楼2012-07-11 20:18:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
12楼: Originally posted by ljlyhl0926 at 2012-07-11 20:18:35
最近也想做这方面的,在网上看到两个比较完整的源代码,一个是http://www.qddayu.cn/information.asp?tid=922,另一个是http://blog.sina.com.cn/s/blog_4b4254430100c5sa.html。都可以用,但是不知道是否有错误。方 ...

我看了你的那两个连接,其中方差除以18的那个应该是错误的,我看了好多文献还有魏凤英老师编的《现代气候统计诊断与预测技术》书,书中的公式都是除以72的。我把魏凤英老师编的《现代气候统计诊断与预测技术》第二版买了,书中就有关于这个方差计算的说法,这个方差是对累计数sk的方差,其中在x1,x2,。。。。xn相互独立,且 有相同连续分布时,方差可由Var=k*(k-1)*(2*k+5)/72计算得出,其中k=2,3,4,5,。。。n。
另外网上的这些程序不敢用,因为要发文章,还是确保万无一失就没没用,万一是错误的呢。查阅不少文献很多人的计算方法不一样,没有统一的公式,这使得更没法相信这个mk检验步骤了,想查阅国外的文献,没有关于这种公式的文献,这种公式是国内发明的,悲剧啊。。
13楼2012-07-12 08:30:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ljlyhl0926

铁虫 (小有名气)


引用回帖:
13楼: Originally posted by 8711193629 at 2012-07-12 08:30:24
我看了你的那两个连接,其中方差除以18的那个应该是错误的,我看了好多文献还有魏凤英老师编的《现代气候统计诊断与预测技术》书,书中的公式都是除以72的。我把魏凤英老师编的《现代气候统计诊断与预测技术》第二 ...

好像M-K一个是用于趋势检验,另一个是用于突变检验,两者公式不一样。后面那个链接有两个代码,前面的是用于趋势检验,后面的是用于突变检验,魏凤英老师的那本书给出的是用于突变检验的公式,突变检验的公式在计算方差时好像都是除以72,至于趋势检验,就不懂了。哎,数理统计基础不好的伤不起啊。
14楼2012-07-12 12:28:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
14楼: Originally posted by ljlyhl0926 at 2012-07-12 12:28:31
好像M-K一个是用于趋势检验,另一个是用于突变检验,两者公式不一样。后面那个链接有两个代码,前面的是用于趋势检验,后面的是用于突变检验,魏凤英老师的那本书给出的是用于突变检验的公式,突变检验的公式在计算 ...

对,国内发明的好像就是突变检验,但是公式都没有统一的,不知道哪个对,哪个错。。后来我哪个软件不用,我就用excel,然后套用魏凤英书里的公式一点点计算,用的是那个书里的数据,可是做出来的图和书中那个图有偏差,大体的趋势很相像,总是有几个地方不一样,于是我总修改公式,怎么修改也做出来的和书上的图不一样。还有一点让我怀疑的是,魏老师的书第一版中秋期望的公式是E=k(K+1)/4,到了第二版,不一样了,公式改为E=k(k-1)/4..但是第一版和第二版做出来的突变图竟然是一模一样,这也太奇怪了,公式都改了,图怎么不改呢,我怀疑那个图应该是错误的。。。。
15楼2012-07-12 14:19:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

rubbyyang

新虫 (正式写手)


Lz 最后完成了么?
16楼2012-07-19 10:56:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

crystral

至尊木虫 (正式写手)


引用回帖:
15楼: Originally posted by 8711193629 at 2012-07-12 14:19:13
对,国内发明的好像就是突变检验,但是公式都没有统一的,不知道哪个对,哪个错。。后来我哪个软件不用,我就用excel,然后套用魏凤英书里的公式一点点计算,用的是那个书里的数据,可是做出来的图和书中那个图有偏 ...

期望值的公式应该是E=K(K-1)/4,第二版的是对的,LZ做出来了吗?我也在研究这个,还没搞明白,头都大了
17楼2012-08-06 08:06:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
16楼: Originally posted by rubbyyang at 2012-07-19 10:56:20
Lz 最后完成了么?

前天弄完了,用公式一点一点计算的,我们很多网友一块讨论的,应该是对的,借用的Excel做的,利用的魏凤英那本树上的例子
18楼2012-08-17 11:30:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

8711193629

金虫 (正式写手)


引用回帖:
17楼: Originally posted by crystral at 2012-08-06 08:06:41
期望值的公式应该是E=K(K-1)/4,第二版的是对的,LZ做出来了吗?我也在研究这个,还没搞明白,头都大了...

恩,弄出来了,没有用什么软件,跟一个网友交流过,她用的DPS做出来的和书上的一模一样,但是经过讨论,书上的那个图与公式不符合,公式中Xi>Xj是对的,而图上是Xi≥Xj,前者是对的。困惑半年了都,烦死了都要
19楼2012-08-17 11:33:03
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
neu2342楼
2012-05-09 21:46  
8711193629(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 8711193629 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 硕博巨婴,也许才刚刚开始 +30 SNaiL1995 2024-06-12 72/3600 2024-06-16 08:34 by 小龙虾biubiubiu
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +5 babu2015 2024-06-14 5/250 2024-06-15 20:40 by lc231001
[教师之家] 请问事业编制和年薪制冲突吗? +6 ZHONGWU_U 2024-06-14 6/300 2024-06-15 20:16 by Ermito
[基金申请] 博后基金,博管会会提前知道消息吗? +3 yuyiang 2024-06-13 3/150 2024-06-15 19:52 by Lynn212
[教师之家] 我们学院常年位居 各学院 倒数第二。专业撤销的话,在编者有什么补偿? +13 河西夜郎 2024-06-09 14/700 2024-06-15 19:44 by LittleBush
[硕博家园] 博士毕业高校和就业的相关问题 +4 SCITOPPP 2024-06-14 6/300 2024-06-15 18:54 by SCITOPPP
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +8 zju2000 2024-06-12 14/700 2024-06-15 13:59 by chemhua
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[硕博家园] 关于硕博连读的一些疑问? +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[基金申请] 化学会年会改成一年开一次了? +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[论文投稿] 投稿kbs被拒后,系统提供的推荐的其他期刊有用吗? +3 chenrui2015 2024-06-10 3/150 2024-06-14 10:30 by Tr.sjx1997
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[有机交流] 旋蒸请教 +6 办公室摸鱼 2024-06-09 9/450 2024-06-13 21:31 by wuxiaominggo
[基金申请] 连续两年医学口青年项目初审体会 +11 进击的荣耀 2024-06-09 18/900 2024-06-13 17:27 by 进击的荣耀
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
信息提示
请填处理意见