文档介绍:完全随机设计多样本资料秩和检验
多组比较的秩和检验
Kruskal-Wallis法
例:某泌尿科医生欲比较三种病人的血中前列腺特异性抗原(PSA)水平,随机选择6名良性前列腺肥大(BPH)患者,8名前列腺癌病理活检阳性患者,8名病理活检阴性的非患者,测定其PSA水平。请作统计分析。
多组比较的秩和检验
统计分析步骤
:
H0:三组病人PSA水平总体分布相同
H1:三组病人PSA水平总体分布不全相同
α=
方法同前
各组秩次相加,得到秩和Ri
统计分析步骤
H0(3组总体分布相同)为真时,在大多数情况下, 很小,所以H一般较小,反之H0非真时, 的绝对值增大,所以H一般会比较大。
本例:H=
:
若组数k=3,每组样本例数≤5,可以查附表得到P值。
若组数超过3,或者至少一组的样本例数超过5,则H值近似服从2分布,自由度为组数k-1。
统计分析步骤
统计分析步骤
本例三组的样本例数都超过5,查2界值表,得P<。
因此,拒绝无效假设,接受备择假设,认为三组病人PSA水平总体分布不相同或不全相同。从平均秩可以看出前列腺癌病理活检阳性患者PSA水平最高。
STATA 命令
STATA 命令
kwallis psa,by(group)
组间两两比较可以用Bonferroni方法校正水平。
Thank You !