硕士研究生医学统计学之 秩和检验.pptVIP

秩转换的非参数检验 第一节 配对样本比较的Wilcoxon符号秩检验 第二节 两个独立样本比较的Wilcoxon秩和检验 第三节 完全随机设计多个样本比较的Kruskal-Wallis H检验 第四节 随机区组设计多个样本比较的Friedman M检验 3. 确定P值，作出统计推断 ν=g-1=3-1=2 接种9D与11C比较，按α=0.05的水准，拒绝H0，接受H1，可以认为两组小白鼠存活日数不同；接种9D与DSC1比较，按α=0.05的水准，拒绝H0，接受H1，可以认为两组小白鼠存活日数不同；接种11C与DSC1比较，按α=0.05的水准，不拒绝H0，还不能认为两组小白鼠存活日数不同。 (二)、bonfferoni调整法 改变检验水准后，使用Wilcoxon秩和检验，对任意两样本进行两两比较。 要使多次比较后犯Ⅰ型错误的累计概率α’’保持不变或至少不超过原水准α，则每次两两比较的检验水准α’可利用bonfferoni不等式确定。 1。多个处理组间的两两比较 2。实验组与同一个对照组的比较 其中：k——样本个数(处理组数) α——原检验水准 Bofferoni提出，若每次检验水准为α’，共进行m次比较，当H0为真时，犯Ⅰ型错误的累计概率α’’不超过mα’，即有bonfferoni不等式α’’mα’成立。 4个样本均数比较，拒绝H0，接受H1， 两两比较需6次，m=6， 此时犯Ⅰ型错误的累计概率：α’’=1-(1-0.05)6=0.265 mα=0.05×6=0.30 要使多次比较后犯Ⅰ型错误的累计概率α’’保持不变或至少不超过原水准α，则每次两两比较的检验水准α’可利用bonfferoni不等式确定：即α’=α/m。 一、多个相关样本比较的Friedman M检验 Friedman M检验用于推断随机区组设计的多个相关样本所来自的总体分布是否有差别。 在理论上讲，检验假设H0应为多个总体分布相同，即多个样本来自同一总体。但由于M检验对总体分布的形状差别不敏感，而对总体分布的位置敏感，所以实际应用中H0可以写作多个总体分布的位置相同 例8-9 8名受试对象在相同实验条件下分别接受4种不同频率声音的刺激，他们的反应率（%）资料见表8-12。问4种频率声音刺激的反应率是否有差别？ 29.5 ─ 　 23.5 ─ 　 16 ─ 　 11 ─ 　 4 10.8 3 8.5 2 8.2 1 7.8 8 4 11.4 3 10.6 2 9.0 1 8.9 7 4 10.6 2 9.6 3 9.8 1 8.6 6 3.5 8.6 3.5 8.6 1 8.0 2 8.3 5 4 12.0 3 9.9 1 8.7 2 9.8 4 3 9.8 4 10.4 1 9.1 2 9.4 3 3 12.0 2 11.8 4 12.7 1 11.6 2 4 11.7 　 3 9.8 　 2 9.6 　 1 8.4 　 1 秩 反应率 秩 反应率 秩 反应率 秩 反应率 频率D 频率C 频率B 频率A 受试号 表8-12 8名受试对象对4种不同频率声音刺激的反应率（%）比较 分析：本例实验设计为随机区组设计，资料为百分率资料，WB=0.810，PB=0.037，选择Friedman M检验 1。 建立假设并确定检验水准 H0：四种频率声音刺激的反应率的总体分布位置相同 H1：四种频率声音刺激的反应率的总体分布位置不全相同 α=0.05 2。计算统计量M ① 将每个区组的数据小到大编秩，遇到数据相同的取平均秩次 ② 分别求各样本秩次之和，用Ri表示。 ③ 计算样本的平均秩和 ④ 计算统计量M 其中：n表示区组个数；g表示处理组个数 3. 确定P值，作出统计推断 ① 查表法：当n≤15以及g≤15时，可以查M界值表确定概率值。 本例：n=8，g=4，查界值表M0.05/2=105 P0.05 按α=0.05的水准，拒绝H0，接受H1，可以认为四种频率声音刺激的反应率的总体分布位置不全相同。 ② χ2近似法：当n15或g15时，可以用χ2近似法确定概率值。 其中：tj为按区组而言的第j个相同秩次的个数 查χ2界值表，P0.005，结论同前 在实际工作中g4或者g=4且n5或者g=3且n9时，就可以使用χ2近似法。 二、多个相关样本两两比较 当M检验拒绝H0，接受H1，认为多个总体分布位置不全相同，需要进一步推断那两个总体分布不同，即多重比较。 (一)、q检验 其中： n——区组个数 g——处理组个数 Ri(Rj)——第i(j)个处理组的秩和 tj