单因素多水平设计定量资料的统计分析及SAS实现.ppt

Kruskal-Wallis H检验 计算统计量H值 本例组数k=3，n1=n2=n3=5，查H界值表，HH0.01=7.98，P0.01，按?=0.01水准拒绝H0，接受H1，差别有统计学意义。认为三种药物杀灭钉螺的死亡率不同 Kruskal-Wallis H检验 比较小白鼠接种三种不同菌型伤寒杆菌9D、11C和DSC1后存活日数，结果如下。问小白鼠接种三种不同菌型伤寒杆菌的存活日数有无差别？ 9D ： 2 2 2 3 4 4 4 5 7 7 11C ： 5 5 6 6 6 7 8 10 12 DSC1 ：3 5 6 6 6 7 7 9 10 11 11 Kruskal-Wallis H检验 本例为时间资料，不符合正态分布，采用Kruskal-Wallis H检验 H0 ：接种三种不同菌型伤寒杆菌的存活日数总体分布位置相同； H1 ：接种三种不同菌型伤寒杆菌的存活日数总体分布位置不同； a=0.05。 Kruskal-Wallis H检验 把三个样本数据混合由小到大编秩，遇数据相等者取平均秩 计算各组秩和，详见下表 Kruskal-Wallis H检验 计算统计量H值 Kruskal-Wallis H检验 本例各组样本数均大于5，查 界值表，可得P0.01，按?=0.05水准拒绝H0，接受H1，差别有统计学意义。可以认为小白鼠接种三种不同菌型伤寒杆菌的存活日数有差别 多个平均秩的两两比较 使用Kruskal-Wallis检验对多个独立样本进行比较，在检验结果为拒绝H0，接受H1时，需要进一步作多个平均秩的两两比较 平均秩的两两比较可以使用Mann-Whitney U检验，也可以采用Dunn法或Nemenyi法 Nemenyi法 Nemenyi法采用 统计量 当各组样本例数较大时，任意两组进行比较的 值为 多个平均秩的两两比较 比较小白鼠接种三种不同菌型伤寒杆菌9D、11C和DSC1后存活日数。通过Kruskal-Wallis检验得知接种三种不同菌型伤寒杆菌的存活日数有差别，试考察任意两者之间的关系 多个平均秩的两两比较 9D与11C 9D与DSC1 11C与DSC1 多个平均秩的两两比较 根据 查 界值表得P0.05，可认为小白鼠接种11C的存活日数高于接种9D的存活日数 根据 查 界值表得P0.025，可认为小白鼠接种DSC1的存活日数高于接种9D的存活日数 根据 查 界值表得P0.99，尚不能认为小白鼠接种11C的存活日数与接种DSC1的存活日数有差别 多总体方差齐性检验 Bartlett的?2检验：也适用于两总体方差齐性检验，要求资料具有正态性 Levene的F检验：也适用于两总体方差齐性检验，不要求资料具有正态性 多总体方差齐性检验 O’Brien的F检验：它是对Levene的F检验法的修正方法 Brown-Forsythe的F检验：模拟的结果说明此法较好 数据变换 对数变换 平方根变换 平方根反正弦变换 倒数变换 完全随机设计资料方差分析表 变异来源 SS 自由度 MS F值 组间变异 k-1 组内变异 N-k 总变异 N-1 单因素多水平设计定量资料的方差分析 某研究者随机选取34例住院患者，采血测定血糖及血总胆固醇。根据血糖水平将患者分为:2型糖尿病组(10例)、糖代谢异常组(12例)和正常对照组(12例)。三组受试对象的血总胆固醇结果如下，试比较他们之间血总胆固醇含量的差异是否有统计学意义？ 2型糖尿病组： 1.64 1.71 2.13 2.28 1.58 1.78 3.24 0.75 3.25 2.05 糖代谢异常组： 1.12 1.85 1.16 0.67 2.12 2.10 1.05 1.84 1.97 0.84 0.93 1.72 正常对照组： 1.55 0.80 1.36 1.06 1.24 0.86 1.19 1.24 0.80 0.72 2.00 0.80 单因素多水平设计定量资料的方差分析 建立假设，确定检验水准 H0 ：m1= m2 = m3（三组患者血总胆固醇含量的总体均数相等）； H1 : 三组患者血总胆固醇含量的总体均数不全相等； a=0.05。 单因素多水平设计定量资料的方差分析 计算检验统计量 单因素多水平设计定量资料的方差分析 单因素多水平设计定量资料的方差分析 * 方差分析表 变异来源 自由度 SS MS F 组间变异 2 4