[基础医学]第十一章 多个样本均数比较的方差分析.ppt

第十一章多个样本均数比较的方差分析 南方医科大学生物统计学系 例1： 为了解烫伤后不同时期切痂对肝脏三磷酸腺苷(ATP)的影响，将30只雄性大鼠随机分为3组，每组10只：A组为烫伤组，B组为烫伤后24h（休克期）切痂组，C组为烫伤后96h（非休克期）切痂组。全部动物统一在烫伤后168h处死并测量其肝脏的ATP含量，结果见下表。试问三组的ATP总体均数是否有差别？ 表1 大鼠烫伤后肝脏ATP的测量结果（mg） 均数的比较 主要内容 方差分析的基本思想 （basic thought of ANOVA） 完全随机设计资料的方差分析 (one-way ANOVA) 随机完全区组/配伍组设计资料的方差分析 （tow-way ANOVA） 多个样本均数间的多重比较（基于方差分析的后续性检验） （post hoc test based on ANOVA） 方差分析的基本思想(basic thought of ANOVA) 变异分解 总变异的分解 三种变异之间的关系 基本思想 将总变异分解为组内变异和组间变异； 将平均组间变异与平均组内变异（误差变异）比较，若前者远大于后者，说明处理间的效应不同；若前者与后者接近，甚至小于后者，说明处理间的效应相同，或称处理因素的影响不大。 方差分析应用前提 完全随机设计资料的方差分析(one-way ANOVA) 完全随机设计 变异分解及计算 变异分解及计算 完全随机设计分析步骤 1.建立假设，确定检验水准a： H0：?1=?2=?3（三组的ATP含量相同） H1：三组的ATP含量不同或不全相同 ?=0.05 完全随机设计分析步骤 2.求检验统计量F值及自由度（列方差分析表） 完全随机设计分析步骤 3.求 P 值，下结论。 若F≤1，不必查表，P＞ ?。本例，P＜ ?，拒绝H0 ，接受H1 ，即不同处理的总体均数不同或不全相同（有待多重比较进一步分析）。 例1的方差分析表 随机完全区组设计的方差分析（two-way ANOVA） 随机完全区组设计 随机完全区组设计（randomized complete block design）也称随机区组设计、配伍组设计或单位组设计，是配对设计的扩展。 设计特点： 每个区组的试验单位数与处理组数相等，可减少试验误差，提高统计假设检验的效率。 随机完全区组设计 变异分解 变易分解 变易分解与计算 变易分解与计算 随机完全区组设计的分析步骤 1、建立假设，确定检验水准a H0：?1=?2=?3 =?4 （4个处理组小鼠体内 的尾蚴存活率的总体均数相等）H1：4个处理组小鼠体内的尾蚴存活率的 总体均数?1、?2、?3 、?4不全相等 ?=0.05 随机完全区组设计的分析步骤 2、计算F值 随机完全区组设计的分析步骤 3、结论 查F界值表，确定P值，下结论。按a=0.05水准，拒绝原假设，4个处理组小鼠体内尾蚴存活率的差别有统计学意义，认为各处理组的小鼠体内尾蚴存活率的总体均数不全相等。由各处理组均数可知，对照组小鼠体内的尾蚴存活率最高。 本例中，区组间F=2.33,P0.05,各区组间的总体均数有差别，说明小鼠体重（或各区组的试验条件）对小鼠体内未有存活率有影响。 随机完全区组设计 当随机完全区组设计的区组项（配伍组项）F＜1时，表明所选配伍因素无意义，此时，可将配伍项与误差项合并（变异与自由度同时合并），成为新的误差项，即随机区组设计退化为完全随机设计，以提高实验效率。 Fblock≥1时无须合并。 多个均数间的两两比较（基于方差分析的后续检验）multiple/pairwise comparison/post hoc test 多个均数间的两两比较 当方差分析的结果为Pa时，只说明k个处理组总体均数不相同或不全相同，不能说明各组总体均数间都有差别。 要分析哪两组间均数有差别，需要进行多个均数间的两两比较（multiple comparison） 多个均数间的两两比较 常用方法 SNK-q检验法 (Student-Newman-Keuls’ test, SNK-test) LSD-t检验法 (Least Significant Difference test) Dunnett-t检验法 (Dunnett’s test) SNK- q检验 检验统计量q 例3 （续例1）试比较三个处理组两两之间的差别。 SNK-q检验 建立假设，确定检验水准a H0: ?A＝?B（任两比较组总体均数相同） H1: ?A≠?B ?＝0.05 SNK检验（q检验） 计算q值 SNK检验（q检验）