完全随机设计多样本资料秩和检验

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