你好,想问下这个滑动T检验程序运行完后怎么画图啊?

你好,想问下这个滑动T检验程序运行完后怎么画图啊?,第1张

=[]%读入你需要的矩阵for i=1:36 index1=i:(i+9)index2=(i+10):(i+19)a=x(index1)b=x(index2)n1=length(a)n2=length(b)c=(n1+n2)/(n1*n2)s1=var(a)s2=var(b)delta1=n1*s1+n2*s2delta=delta1/(n1+n2-2)t(1,i)=(mean(a)-mean(b))/sqrt(delta*c)endЩШЬЫбФУТлмйиёнжёлмджж━┼┴┤┬┤

长时间尺度的气候研究中通常涉及到年际变化和年代际变化。文献中也经常出现所谓年代际突变这类的描述,这次便介绍两个用来检验年代际变化的方法, 其一是滑动t检验(sliding t-test),另一个则是曼肯德尔检验(Man-Kendall test)。

滑动t检验是考察两组样本平均值的差异是否显著来检验突变。

滑动t检验的基本思想是:把一气候序列中两段子序列均值有无显著差异看为来自两个总体均值有无显著差异的问题来检验。如果两段子序列的均值差异超过了一定的显著性水平,可以认为均值发生了质变,有突变发生。

原文也提到了该方法的局限性,需要人为设置滑动步长,具有一定主观性,需反复设置不同步长最终确定合适的突变点。

自由度n1+n2-2,根据置信度检验表查找对应显著性阈值。

曼肯检验是一种非参检验方法,避免了滑动t检验的局限性,在年代际变化研究中受到广泛应用。

以下是计算步骤:

对于同一组数据:

分别使用两种方法进行检验,首先是滑动t检验:

存在多个突变点,这时便需要调整滑动补偿,选取合适的步长。

而利用MK检验:

对于该组数据,相比之下,MK检验的效果要优于滑动t检验。

两种检验的分析方法如下:

滑动t检验:

M-K检验:

T检验法:应用t分布理论对正态总体或近似服从正态分布的总体当方差σ2未知时关于平均数的检验方法.

Q检验法:首先把数据按照从大到小排序,找出最大值与最小值,并计算可疑出其与相邻值的差值,并将其与最大值与最小值之差做商.得出Q与题目给出的要求的Q对比,要是大于,则是舍去,应该正确.


欢迎分享,转载请注明来源:内存溢出

原文地址: https://outofmemory.cn/yw/12070265.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-20
下一篇 2023-05-20

发表评论

登录后才能评论

评论列表(0条)

保存