首先,我们需要确认一下游程检验的公式,公式很难输入大家自己查(心理统计学,胡竹菁,2010,p283)。
这里有两种选择:
1.你直接用matlab里自带的函数runstest
格式:
[h,sig,stats]=runstest(x1,alpha)
h为判断结果
sig为接受概率
stats分4个部分(nruns为游程数;n1为取值1的数目;n0为取值0的数目;Z为检验值)
好了程序如下:
clear all;
x=input('请选择数据文件: ');
alpha=input('请输入显著性水平: ');
[n k]=size(x);
for i=1:k
x1=x(:,i);
[h,sig,stats]=runstest(x1,alpha)
if h(:,1)==1
jl=['结论:在显著性水平',num2str(alpha),'(双侧)下,拒绝虚无假设,即设所检验的数据的序列是非随机的'];
else
jl=['结论:在显著性水平',num2str(alpha),'(双侧)下,接受虚无假设,即设所检验的数据的序列是随机的'];
end
diary('H:\data
file\MATLAB\数据运行结果\游程检验结果.doc') ;
diary on;