資源簡介
M-K檢驗法Matlab函數式M文件代碼,含計算、繪圖部分,能很好符合魏鳳英老師的例子
代碼片段和文件信息
function?pp=MK(x)
l=length(x);
r=zeros(l1);
s=zeros(l1);
E=zeros(l1);
V=zeros(l1);
UF=zeros(l1);
UB=zeros(l1);
%…………………………………………………………………初始化……………………………………………………………………………
%計算統計量s
for?i=1:l
????for?j=1:i
????????if?x(i)>x(j)
???????????r(i)=r(i)+1;
????????end
????end
s(i)=sum(r);
end
%計算UF
for?i=2:l
????E(i)=i*(i-1)/4;
????V(i)=i*(i-1)*(2*i+5)/72;
????UF(i)=(s(i)-E(i))/sqrt(V(i));
end
UF(1)=0;
%寫入文件
????????%DLMWRITE(‘MK1.txt‘UF);
%…………………………………………………………第一部分完成…………………………………………………………………………
%逆序
x=flipud(x);
%計算統計量s
r=zeros(l1);
for?i=1:l
????for?j=1:i
????????if?x(i)>x(j)
???????????r(i)=r(i)+1;
????????end
????end
s(i)=sum(r);
end
%計算UB
for?i=2:l
????
評論
共有 條評論