卫生统计学 方差分析.ppt

先来看一个具体的例题 例8-1 某大学营养与食品卫生研究所将800只条件一致的雌性果蝇随机分配到4种不同浓度的某受试物培养基组，各组200只。经2至3月的培养试验，得各组寿命最高的10只果蝇的生存天数如下： 从本例引伸的问题 属于什么样的研究？设计类型是什么？ 该试验的三要素具体是什么？ 处理因素有几个水平？ 所获得的是什么类型的资料？ 研究目的是什么？ 应该用什么统计方法证实假设？ 第一节 完全随机设计多组 均数比较的方差分析 方差分析的基本思想和应用条件 方差分析（analysis of variance，ANOVA) 是由英国统计学家R.A. Fisher在２０世纪初提出 是用于多个样本均数检验的一种假设检验方法 H0：多个样本总体均数相等 H1：多个样本总体均数不相等或不全等 α=0.05 如何计算统计量 F？ 成组设计方差分析数据满足的条件 独立性 正态性 方差齐性 one-way anova 基本原理 如采用t检验进行两两比较，将增加I类错误 SAS程序 Li8_1 相关SAS程序 Li8_22 SAS 程序： 什么是随机区组设计？ 例8-3 为研究氯化镉CdCl2对V79细胞的毒性作用,以5种不同染毒剂量的氯化镉作为5个区组，以3种不同的染毒时间作为3个处理组，观察V79细胞的克隆率，资料如表8-1。①比较不同剂量氯化镉组的V79细胞克隆率；②比较不同染毒时间组的V79细胞克隆率。 方差分析补充 拉丁方设计方差分析 析因设计方差分析 交叉设计方差分析 重复测量设计的方差分析 拉丁方设计方差分析 SAS程序 析因设计方差分析 概念:是一种多因素多水平交叉分组进行全面实验的设计方法 主要是分析因素间的交互作用 例：为了研究不同氧浓度（因素B）和抗癌药（因素A），对用放射性3H-胸腺嘧啶（3H-TdR）渗入后的人红白血病细胞K562的抑制效果，因素B分为B1（-）、B2（含氧20%），因素A为A1（-），A2（表阿霉素），进行2×2析因设计，数据见下表，试分析A、B两因素对K562细胞的抑制效果。 两因素析因设计的方差分析思想 SS总 =SSA+SSB+SSAB+SS误差 SAS程序 析因分析例题2 SAS程序 交叉设计方差分析 SS总 =SS药物+SS阶段+SS个体+SS误差 SAS程序 阶段不同 二、多组有序变量资料的秩和检验 1.实例及计算分析步骤 例8－3　某大学社会医学与全科医学研究所采用匿名自填式问卷调查了2908名进城农民工对性自慰/手淫的认识，认识程度分为同意、无所谓、反对三个等级，资料如表8－3 表8－3不同文化程度组民工对 性自慰/手淫的认识程度等级比较 174 866 1436 432 合计 2322.0 1736－2908 61 340 587 185 反对 1257.0 779－1735 44 258 522 133 无所谓 389.5 1－778 69 268 327 114 同意 平均秩次 秩次范围 大专及以上 高中或中专 初　中 小　学 例 数 认识等级 试比较小学、初中、高中或中专、大专及以上不同文化程度组的平均等级 ⑴建立检验假设 H0 :4总体分布位置相同 H1 :4总体分布位置全不相同或不全相同　 　α=0.05 ⑵ 编秩求秩和 先计算各等级的合计，再确定各等级的秩次范围及平均秩次，如等级“同意”的合计为778，则其秩次范围为1－778，其平均秩次为（1＋778）/2＝389.5； 计算分析步骤如下： 又如等级“无所谓”的合计为957，则其秩次范围为779－（778＋957），即779－1735，其平均秩次 为1257。以各等级的平均秩次近似代替该等级各调查对象的秩次，计算各组的秩和，如小学组秩和＝389.5×114＋1257×133＋2322×185＝641154。 ⑶ 计算检验统计量 按公式（8－6）计算检验统计量H值，当各样本相同秩较多时，用H除以C的商对H作校正。这里，其中为第j个相同秩次的个数，N为各组例数之和。 ⑷ 确定P值 下结论的方法同多组连续变量资料 的秩和检验。 ⑸ 两两比较 当结论为“各组总体分布位置全不相同或不全相同”时，需要进一步作两两比较。可采用两样本等级资料秩和检验的方法，借助SAS或SPSS软件得到相应的P值（参见第七章），按公式α’=α/k计算各次两两比较的检验水平α’，这里α为总检验水平，k为两两比较的次数。 因变量为等级资料，不呈正态分，不能采用完全