关于matlab的滑动t检验

m0932 · 发表于 2014-5-12 19:50:31

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

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

x

本帖最后由 m0932 于 2014-5-12 19:52 编辑

最近做突变检验，想用滑动t检验来做，于是照着公式写了下，可是画出来的图明显不对啊。我把程序贴上来，希望大家帮忙看看哪里不对了。另外，我把数据也发上来，以帮助修改。图为改数据的滑动t检验结果，步长为100，虚线为99%的检验线。

function [t,P]=sttest(V,step)
windows=step;
data=V;
clear V;
len=length(data);
%j=1;
for ii=windows:len-windows
      x1=data(ii-windows+1:ii);
      x2=data(ii+1:ii+windows);
      s=sqrt((windows*var(x1,1)+windows*var(x2,1))/(windows+windows-2));
      n=sqrt(1/windows+1/windows);
      t(ii)=(mean(x2)-mean(x1))/(s*n);
      %j=j+1;
end
t(1:windows-1)=nan;
tep=zeros(1,length(len-windows+1:len));
tep(:)=nan;
t=[t,tep];
P=tinv(1-0.01,windows+windows-2);

m0932 · 发表于 2014-5-13 18:45:19

有人能用这个数据画一下吗？

m0932 · 发表于 2014-5-15 10:35:54

求高手指点下

依景止水 · 发表于 2014-6-4 14:46:06

其实关于读MK图的这个问题可以专门开个讨论帖

路在何方 · 发表于 2014-6-28 23:22:52

这个程序我直接运行不了，为什么，程序有问题吗

路在何方 · 发表于 2014-6-28 23:23:15

小白期待有人解答

z1985az001 · 发表于 2015-1-2 23:06:57

学习学习！！

逍遥仙境 · 发表于 2015-3-11 15:39:49

明白了，谢谢！

霸道天儿 · 发表于 2015-5-11 10:24:56

家园有很多这样的贴及资料

走路睡觉的猫WD · 发表于 2015-6-15 14:47:05

我也学习一下，正好遇到同样的问题

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

[讨论] 关于matlab的滑动t检验

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