现在滑动t检验的程序发出来，求大家帮我操作一下，不胜感激

龙行天下

matlab程序如下：
function  t=mtt(a,b)

%MTT Moving-t test
%   load ltzhshshl
%   x=z;
%   for i=1:55
%   a=x(1:i);
%   b=x(i:55);
%   c(i)=mtt(a,b);
%   end

n1=length(a);
n2=length(b);
c=(n1+n2)/(n1*n2);
s1=1/n1*sum((a-mean(a)).^2);
s2=1/n2*sum((b-mean(b)).^2);
delta1=(n1-1)*s1+(n2-1)*s2;
delta=delta1/(n1+n2-2);
t=(mean(a)-mean(b))/sqrt(delta*c);
请随便给我举个例子，另外给我讲一下这个程序中的滑动步长是怎么设置的？

kongfeng0824

%MTT Moving-t test
%   load ltzhshshl
%   x=z;
%   for i=1:55
%   a=x(1:i);
%   b=x(i:55);
%   c(i)=mtt(a,b);
%   end
这个就是最好的例子。
看看滑动t的原理，再看看这个例子。

kongfeng0824

这个和滑动F检验都很类似

龙行天下

kongfeng0824 发表于 2013-9-6 17:21

                               
登录/注册后可看大图

%MTT Moving-t test
%   load ltzhshshl
%   x=z;

哦，谢谢啊，已经做出来了。但是还有个小疑问，就是这个例子中的滑动步长是怎么确定的？？望高人赐教，不胜感激！

龙行天下

kongfeng0824 发表于 2013-9-6 17:22

                               
登录/注册后可看大图

这个和滑动F检验都很类似

恩，这两个检验的原理接近，请问你知道这个例子中的滑动步长是怎么确定的？？

kongfeng0824

%   a=x(1:i);
%   b=x(i:55);
这个应该是设置步长的原因所在吧。我觉得。我没看过滑动t的原理，但是我感觉是

龙行天下

kongfeng0824 发表于 2013-9-6 17:43

                               
登录/注册后可看大图

%   a=x(1:i);
%   b=x(i:55);
这个应该是设置步长的原因所在吧。我觉得。我没看过滑动t的原理，但是我感 ...

我也感觉这个是设置步长的，但是不知道这个步长到底是怎么获取的

kongfeng0824

因为没有设置循环，所以你试试吧i取一个特定的值试试。你可以试试excel中华东的结果和这个一样不一样即可以检验

z1985az001

求指教！！！！

风儿飘飘

不错，估计先得懂原理才好应用自如