请帮我看看matlab做出的mk图与魏凤英老师书中图不一样 哪里错了

katze

请帮我看看matlab做出的mk图与魏凤英老师书中图不一样 哪里错了 谢谢！

houtian

建议一下  把程序直接提贴出来吧   

二爷名声在外

说实话，看起来差不多，不知道楼主说的不一样在哪里……

katze

交点对应的坐标

katze

houtian 发表于 2014-6-25 00:38
建议一下  把程序直接提贴出来吧

A=xlsread('shiyan.xls')
x=A(:,1);
y=A(:,2);
N=length(y);
n=length(y);
Sk=zeros(size(y));
UFk=zeros(size(y));
s = 0;
for i=2:n
   for j=1:i
         if y(i)>y(j)
           s=s+1;
         else
           s=s+0;
         end;
   end;
   Sk(i)=s;
   E=i*(i-1)/4;
  Var=i*(i-1)*(2*i+5)/72;
  UFk(i)=(Sk(i)-E)/sqrt(Var);
end;
y2=zeros(size(y));
Sk2=zeros(size(y));
UBk=zeros(size(y));
s=0;
for i=1:n
    y2(i)=y(n-i+1);
end;
for i=2:n
   for j=1:i
         if y2(i)>y2(j)
           s=s+1;
         else
           s=s+0;
         end;
   end;
   Sk2(i)=s;
   E=i*(i-1)/4;
  Var=i*(i-1)*(2*i+5)/72;
  UBk(i)=0-(Sk2(i)-E)/sqrt(Var);
end;
UBk2=zeros(size(y));
for i=1:n
   UBk2(i)=UBk(n-i+1);
end;
xlswrite('f:\test2.xls',UFk,'Sheet1','A1');
xlswrite('f:\test2.xls',UBk2,'Sheet1','B1');
figure(3)
plot(x,UFk,'r-','linewidth',1.5);
hold on
plot(x,UBk2,'b-.','linewidth',1.5);
plot(x,1.96*ones(N,1),':','linewidth',1);
axis([min(x),max(x),-4,8]);
legend('UF统计量','UB统计量','0.05显著水平');
xlabel('t (year)','FontName','TimesNewRoman','FontSize',12);
ylabel('统计量','FontName','TimesNewRoman','Fontsize',12);
hold on
plot(x,0*ones(N,1),'-.','linewidth',1);
plot(x,1.96*ones(N,1),':','linewidth',1);
plot(x,-1.96*ones(N,1),':','linewidth',1);

katze

二爷名声在外 发表于 2014-6-25 08:21
说实话，看起来差不多，不知道楼主说的不一样在哪里……

交点对应的时间 还有统计值坐标 都不一样

二爷名声在外

katze 发表于 2014-6-25 08:31
交点对应的时间 还有统计值坐标 都不一样

看出来了……

言深深

看下你的图，UF是对的，UB似乎错了，请检查一这个问题。ps：不懂matlab···

cx6689942

我前不久也碰到这个问题，做详细分析和参考相关文献后证明你在matlab做的是正确的。

映山红

魏凤英老师的那个图可能错了