方差分析20071024课件.ppt

* 方差分析 协和医科大学 流行病教研室 王艳红 流行病学资料分析 教学目的 掌握方差分析的基本思想 了解和识别实验设计资料的常见类型 掌握利用SAS软件进行不同设计类型资料的方差分析 正确解读不同设计类型方差分析的结果 一、相关理论的回顾 因变量与自变量 因变量 即：反应变量 受其他因素影响的变量 自变量 即：处理因素 或 效应因子 不受其他因素影响的变量 主效应与交互作用 主效应 效应因子对因变量的独立影响 交互作用 几个效应因子对因变量产生的交叉影响 一、相关理论的回顾 方差分析能解决的问题举例 不同年龄人群治疗前后血压下降水平的比较（主效应：年龄；单向方差分析） 不同药物的降压效果的比较（主效应：药物；单向方差分析） 研究对象按年龄分组后，不同药物降压效果的比较（主效应：年龄、药物；双向方差分析） 将研究对象划分为：ACEI类，钙拮抗剂， ACEI类+钙拮抗剂，对照四组，探讨ACEI类和钙拮抗剂的联合降压效果（主效应：ACEI类、钙拮抗剂、 ACEI类×钙拮抗剂；交互作用） …… 连续变量 （因变量） 例如： 治疗前后收缩压的下降值 1个或多个分类变量 （自变量） 例如： 药物：ACEI类药物、钙拮抗剂、利尿剂 性别：男性与女性 年龄分组：30-39、40-49、50-59、60-69、70-79岁 一、相关理论的回顾 方差分析的用途 两个或多个样本均数间的比较 两个或多个因素的交互作用 回归方程的假设检验 方差齐性检验（单向方差分析） 一、相关理论的回顾 不同处理组样本均数间的差异来源 随机误差 组间处理因素的差异 方差分析的基本思想 全部数据的离均差平方和（ss总）分解为： 处理组间（ss处理组间） 误差（ss误差） 自由度（ν总）也相应分解为： 处理组间（ν处理组间） 误差（ ν 误差） ss总和ν总的分解往往与特定的设计方案有关 二、方差分析中SAS相关语言 二、方差分析中SAS相关语言 anova过程 应用： 方差分析 平衡实验设计资料（也能用于不均衡的完全随机设计资料） 运行速度相对较快 glm过程 应用 方差分析 协方差分析 回归分析 偏相关分析 … 平衡、非平衡实验设计资料 运行时间相对较长 三、常见的实验设计类型 （一）、完全随机设计 （二）、随机区组设计 （三）、拉丁方设计 （四）、析因设计 （一）、完全随机设计 设计方案 实例 某医院妇产科为考察几种卵巢功能异常患者血清中促黄体素的水平，测定其结果如下： 卵巢发育不良：42.50 38.31 35.76 33.60 31.38 丘脑性闭合：6.71 3.32 4.59 1.67 10.51 2.96 11.82 3.86 8.26 2.63 2.20 垂体性闭合：4.50 2.75 11.14 5.98 1.90 5.43 11.05 22.03 问：它们之间的差异是否有显著性意义？ SAS数据集的建立 data aa; if _n_6 then g=1; if _n_=6 and _n_17 then g=2; if _n_=17 then g=3; input y@@; cards; 42.50 38.31 35.76 33.60 31.38 6.71 3.32 4.59 1.67 10.51 2.96 11.82 3.86 8.26 2.63 2.20 4.50 2.75 11.14 5.98 1.90 5.43 11.05 22.03 ; run; proc univariate data=aa normal; var y; class g; run; proc anova / glm data=aa; class g; model y=g; means g; means g / hovtest snk dunnett(1); run; quit; SAS过程步 练习1 为了解烫伤后不同时期切痂对肝脏三磷酸腺苷（ATP）含量的影响，将30只雄性动物随机分成3组，每组10只： A组：烫伤对照组， B组：烫伤后24小时（休克期）切痂组 C组：烫伤后96小时（非休克组）切痂组 全部动物统一在烫伤后168小时处死测量其肝脏ATP含量，结果见右表。 问：三组大鼠肝脏的ATP含量有无差别？ （二）、随机区组设计 设计方案 实例 为研究注射不同剂量雌激素对大白鼠子宫体重的影响，取4窝不同种系的大白鼠，每窝3只，随机分配到每组内接受不同剂量雌激素的注射，然后测定其子宫重量，结果见下表。 问：注射不同剂量的雌激素对大白鼠子宫重量是否有影响？ SAS数据集的建立 data aa; do block=1 to 4; do dose=1 to 3; input weight @@; output; end; end; cards;