关于魏凤英《现代气候统计诊断与预测技术》书中Mann-Kendall突变检验法

yutian2999 · 发表于 2012-7-14 02:00:58

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

在魏凤英老师的《现代气候统计诊断与预测技术》书中第70页里Mann-Kendall突变检验法中的公式(5.4.3)求解E[sk]=n*(n+1)/4有误！参考唐启义2006年的《DPS数据处理系统——试验设计、统计分析及模型优化》书中第841页。根据书中第71页给的1900-1990年上海年平均气温序列数据，计算得到结果与图5.2有出入，经过分析与计算，最后找到原因：计算ri时采用xi>=xj(j=1,2,...,i-1)，请对比参考第69页公式(5.4.1)。DPS数据处理系统软件也正是采用以上xi>=xj来判断记数，得出上海气温突变始于1921年，但如果采用xi>xj来记数得到正确突变点是1931年。魏凤英根据图5.2得出了“上海年平均气温的增暖趋势及发生突变的时间均与北半球年平均气温完全一致。”（见书72页）。现给出两种不同记数得到不同的突变点，相差达10年！可见MK方法采用累计记数对数据中相同值个数相当敏感。
本贴探讨分析仅供学习之用，特别是新人少走弯路，无对老师不敬，也提醒各老师撰写出书时能多注意算法和细节内容正确，以免怠误后人！另需说明，DPS数据处理软件计算也存不当之处。

上海1900-1990气温MK检验出入.jpg

上海气温Mann-Kendall突变分析对比图。

足无际 · 发表于 2015-5-6 23:25:22

楼主好，既然已知道问题所在，我想知道xi>=xj变为xi>xj在DPS软件中能不能修改？若不能，可不可以推荐一款可以正确计算突变的步骤和方法以及软件？

topmad · 发表于 2012-7-14 08:39:05

欢迎多多讨论...

言深深 · 发表于 2012-7-14 21:43:39

这样的工作需要多做一做的，很有意义，赞

flyingbird2008 · 发表于 2012-7-20 15:56:12

谢谢楼主的分享

wl650915 · 发表于 2012-7-20 16:13:52

楼主高人

honview · 发表于 2012-8-3 15:00:01

大家多讨论才是王道

朝雾 · 发表于 2012-8-3 15:13:54

很有意义，有时使用时确实没有思考过

museum · 发表于 2012-8-21 17:11:22

那。。那这个最终该相信哪个的结果呢？

vanhojer · 发表于 2012-8-23 10:15:23

前辈也有失误的时候，mk挺好的

scaler · 发表于 2012-8-31 20:18:18

楼主看的是哪一版的我看的第二版跟你说的不一致啊

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

关于魏凤英《现代气候统计诊断与预测技术》书中Mann-Kendall突变检验法

登录后查看更多精彩内容~

评分