Kruskal-Wallis检验

秩和检验测试

描述

p= Kruskal-Wallis检验(X返回p在矩阵的每一列 - 值对的零假设的数据X来自相同的分布,使用秩和检验测试。另一种假设是,并不是所有的样品都来自相同的分布。Kruskal-Wallis检验也返回一个方差分析表和箱形图。

p= Kruskal-Wallis检验(X返回p-VALUE对于零假设的试验是,在每个分类组中的数据,如由分组变量指定的来自相同的分布。另一种假设是,并非所有的组来自相同的分布。

p= Kruskal-Wallis检验(Xdisplayopt返回p测试的 - 值,并可以显示或抑制方差分析表和箱形图。

[pTBL统计] = Kruskal-Wallis检验(___还返回ANOVA表作为单元阵列TBL和结构统计包含有关测试统计信息。

例子

全部收缩

创建两个不同的正态分布的对象。第一个分布亩= 0西格玛= 1和第二分配具有亩= 2西格玛= 1

PD1 = makedist('正常');PD2 = makedist('正常''亩',2,“西格玛”,1);

通过从这两个分布产生随机数创建采样数据的矩阵。

RNG('默认');%,持续重现X = [随机(pd1,20,2),随机(pd2,20,1)];

的前两列X含有从第一分布产生的数据,而第三列包含从第二分布生成的数据。

检验零假设,从每列中的样本数据X来自相同的分布。

P = Kruskal-Wallis检验(x)的

P = 3.6896e-06

的返回值p表明Kruskal-Wallis检验拒绝零假设,这三个数据样本来自在1%的显着性水平相同的分布。方差分析表提供了额外的测试结果,以及箱形图直观地呈现在每列的汇总统计X

创建两个不同的正态分布的对象。第一个分布亩= 0西格玛= 1。第二分布亩= 2西格玛= 1

PD1 = makedist('正常');PD2 = makedist('正常''亩',2,“西格玛”,1);

通过从这两个分布产生随机数创建采样数据的矩阵。

RNG('默认');%,持续重现X = [随机(pd1,20,2),随机(pd2,20,1)];

的前两列X含有从第一分布产生的数据,而第三列包含从第二分布生成的数据。

检验零假设,从每列中的样本数据X来自相同的分布。抑制输出的显示器,并生成结构统计在进一步的测试使用。

[P,TBL,统计] = Kruskal-Wallis检验(X,[],“关”
P = 3.6896e-06
TBL =4×6单元阵列列1至5 { '源'} { 'SS'} { 'DF'} { 'MS'} { '智平方'} { '列'} {[7.6311e + 03]} {[2]} {[3.8155e+03]} {[ 25.0200]} {'Error' } {[1.0364e+04]} {[57]} {[ 181.8228]} {0x0 double} {'Total' } {[ 17995]} {[59]} {0x0 double } {0x0 double} Column 6 {'Prob>Chi-sq'} {[ 3.6896e-06]} {0x0 double } {0x0 double }
统计=同场的结构:gnames:[3X1字符] N:[20 20 20]来源: 'Kruskal-Wallis检验' meanranks:[26.7500 18.9500 45.8000] SUMT:0

的返回值p表明测试拒绝在1%的显着性水平的零假设。您可以使用结构统计以执行附加的后续测试。单元阵列TBL包含相同的数据作为图形ANOVA表,包括列和行的标签。

进行后续试验,以确定哪个数据样本来自不同的分布。

C = multcompare(统计)
注:间隔可用于测试,但不是同时置信区间。

C =3×61.0000 2.0000 -5.1435 7.8000 0.3345 20.7435 1.0000 3.0000 -31.9935 -19.0500 -6.1065 0.0016 2.0000 3.0000 -39.7935 -26.8500 -13.9065 0.0000

结果表明,有组1和组3之间的差异显著,所以检验拒绝零假设在这两个组中的数据来自相同的分布。这同样适用于组2和3真实但是,没有组1和2之间的显著差异,所以测试不拒绝零假设,这两个组来自相同的分布。因此,这些结果表明,在组1和2中的数据来自相同的分布,以及在组3中的数据来自不同的分布。

创建载体,强度含的金属束强度的测量值。创建的第二矢量,合金,指示从中相应光束由金属合金的类型。

强度= [82 86 79 83 84 85 86 87 74 82...78 75 76 77 79 79 77 78 82 79];合金= {'ST''ST''ST''ST''ST''ST''ST''ST'...'AL1''AL1''AL1''AL1''AL1''AL1'...'AL2''AL2''AL2''AL2''AL2''AL2'};

检验零假设,光束强度测量跨所有三种合金相同的分布。

P = Kruskal-Wallis检验(强度,合金,“关”
p值= 0.0018

的返回值p表明测试拒绝在1%的显着性水平的零假设。

输入参数

全部收缩

对于假设检验的样本数据,指定为载体或-通过-ñ矩阵。如果X是一个-通过-ñ矩阵中,每个ñ列表示含有一个独立的样品相互独立意见。

数据类型:|

分组变量,指定为数值或逻辑向量,一个字符或字符串数​​组,或字符向量的单元阵列。

  • 如果X是矢量,然后在每个元件识别组,其在相应的元件X所属,和必须是相同的长度的矢量X。如果某行的包含一个空值,即行和相应的观测X被忽略。为NaN在任一值X要么同样被忽略。

  • 如果X是一个矩阵,则在每一列X代表不同的组,并且可以使用指定这些列的标签。元素的数量和列的数X必须相等。

这些标签包含在还注释箱形图。

例:{ '红', '蓝', '绿', '蓝', '红', '蓝', '绿', '绿色', '红色'}

数据类型:||合乎逻辑|烧焦||细胞

显示选项,指定为'上'要么“关”。如果displayopt'上'Kruskal-Wallis检验显示了以下数字:

  • 含正方形,自由度,以及其他量的总和的ANOVA表来计算基于所述数据的行列X

  • 在数据矩阵中的每一列的数据的箱形图X。箱形图是根据实际的数据值,而不是在行列。

如果displayopt“关”Kruskal-Wallis检验不显示这些数字。

如果指定的值displayopt,你还必须指定一个值。如果你没有一个分组变量,指定[]

例:“关”

输出参数

全部收缩

p-VALUE试验,返回作为在范围[0,1]的标量值的。p是观察检验统计量极端或比更加极端,在零假设下所观察到的值的概率。小值p在零假设的有效性产生怀疑。

测试结果的方差分析表,返回为单元阵列。TBL包括正方形的基于数据的中的行列计算的总和,自由度,以及其他量X,以及列和行标签。

测试数据,返回的结构。您可以通过使用成对样品中位数的后续执行多重比较检验multcompare统计作为输入值。

更多关于

全部收缩

秩和检验测试

克鲁斯卡尔 - 沃利斯检验是经典的单因素方差分析的非参数版本和Wilcoxon秩和检验的扩展到两个以上的组。它比较数据的组的中位数X以确定该样品来自相同的群体(或,等价地,来自不同种群具有相同的分布)。

该秩和检验测试使用的数据的行列,而不是数值,计算检验统计量。它通过从最小的数据排序,以最大的跨所有群体,并采取这种排序的数字索引查找行列。对于捆绑观察的秩等于与它捆绑所有观测的平均评级。该F在古典单向ANOVA使用t-统计通过卡方统计代替,并且p- 值措施的卡方统计的意义。

在Kruskal-Wallis检验假定所有样本来自从由于组效果可能不同的位置具有相同的连续分布,除了人口,并且所有的观察是相互独立的。相比之下,传统的单因素方差分析替换较强的假设总体具有正态分布的第一个假设。

R2006a前推出