例题:
使用H检验法:
x=[21,24,27,20;20,18,19,15;22,25,27,22];
x1=reshape(x',numel(x),1);
group=[1,1,1,1,2,2,2,2,3,3,3,3];
[p,table,stats]=kruskalwallis(x1,group);
Source SS df MS Chi-sq Prob>Chi-sq
------------------------------------------------------
Groups 93.875 2 46.9375 7.3 0.026
Error 47.625 9 5.2917
Total 141.5 11
p值为0.026<0.05,故落在拒绝域里,拒绝原假设,认为3组样本来自分布不相同的总体。
EX.两两多重比较
c=multcompare(stats)
c =
1.0000 2.0000 -0.8189 5.1250 11.0689 0.1072
1.0000 3.0000 -7.3189 -1.3750 4.5689 0.8505
2.0000 3.0000 -12.4439 -6.5000 -0.5561 0.0280
可得第2组数据与第3组数据服从的总体分布不相同。