matlab滑动T检验源代码分享及出图

sswyj

M=load('wyj.txt');
x=M(:,1);y=M(:,2);t=zeros(30,1);y1=zeros(30,1);
%MMT方法检验%
%先设定子序列的长度（取n1=5;n2=5)%
n1=5;n2=5;
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:25
    averx1=mean(y(i-4:i));
    averx2=mean(y(i+1:i+5));
    s1=var(y(i-4:i));  
    s2=var(y(i+1:i+5));
    s=sqrt((n1*s1+n2*s2)/(n1+n2-2));
    t(i-4)=(averx2-averx1)/(s*m);  
end
plot(x(1:30),t,'r-');
y1(:)=2;                   %n=20-2时，显著性检验标准,99.5%
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1971,2000,5,4]);
xlabel('年');ylabel('t');title('MMT(n=30)检验');

请高人帮我看看源代码有错的地方吗，请帮我指出，谢谢。还有我的这个图1992年以后为什么线条这么直，哪里出现错误了。非常感谢！

小英

最近也在研究，你试试把这两处改一下呢  
第二列 t=zeros(21,1)%21表示步长5到25的序列数%
倒数第六列 plot(x(5:25),t,"r-")

小英

还有你的显著性检验标准应该是3.36， 因为自由度v=n1+n2-2=8,v为8时t0.01=3.36

liang888

谢谢楼主分享

My93131

你好，我现在也在做滑动T检验，运行到倒数第四步时，出现Error using plot   Vectors  must be  the  same lengths 是哪里做错了呢