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