方差分析（二）：重复测量资料的方差分析.pptVIP

方差分析（二）：重复测量资料的方差分析 浙江大学医学院公共卫生系 流行病与卫生统计教研室 沈 毅 2005.3 重复测量（repeated measure）是指对同一观察对象的同一观察指标在不同时间点上进行的多次测量，用于分析该观察指标在不同时间上的变化特点。这类测量资料在临床和流行病学研究中比较常见。例如，为研究某种药物对高血压病人的治疗效果，需要定时多次测定受试者的血压，以分析其血压的变动情况。一些传统的统计方法，如t检验、方差分析、线性回归模型等，都要求各次观察是相互独立的。而重复测量资料由于是对同一受试者的某项观察指标进行的多次测量，在同一受试者的多次测量之间可能存在某种相关性，用通常的统计方法就不能充分揭示出其内在的特点，有时甚至会得出错误的结论。因此，有关重复测量资料的分析方法是近代统计学研究的热点之一。 第一节 重复测量资料方差分析对协方差阵的要求 在对重复测量资料进行方差分析时，除要求样本是随机的、在处理的同一个水平上的观察是独立的以及每一水平上的测定值都来自正态总体外，特别强调协方差阵（covariance matrix）的球形性（sphericity）或为符合对称性（compound symmetry）。Box（1954）指出，若球形性质得不到满足，则方差分析的F值是有偏的，这会造成过多的拒绝本来是真的无效假设(增加Ⅰ型错误)。 一、协方差阵的球形性检验 方差是指在某一时点上测定值变异性的大小，而协方差是指在两个不同时点上测定值相互变异性的大小。如果在某个时点上的取值不影响其他时点上的取值，则协方差为0，反之，则不为0。由方差协方差构成的矩阵称协方差阵。 设 k、l为两个测定时点， 代表协方差阵中的元素。当k＝l时为方差，k≠l时为协方差。共有a个测定时点，将这a个方差和（a-1）/2个协方差排成协方差阵V为： 协方差阵的球形性质是指该矩阵主对角线元素（方差）相等、非主对角线元素（协方差）为零。用Mauchly氏法检验协方差阵的球形性质。Mauchly氏检验的 P值若大于研究者所选择的显著性水准α时，说明协方差阵的球形性质得到满足。否则，必须对与时间有关的F统计量的分子、分母自由度进行调整，以便减少犯Ι类错误的概率。调整系数为ε（epsilon）。 二、自由度调整方法 1. 调整系数ε的计算 有两种调整系数。 （1）Greenhouse-Geisser调整系数 为： 式（10-2）中的 是矩阵（10-1）中第k行第l列元素， 是所有元素的总平均值， 是主对角线元素的平均值， 是第k行的平均值。 的取值在1.0与1/（a-1）之间。 （2）Huynh－Feldt调整系数 研究表明，当ε真值在0.7以上时，用 进行自由度凋整后的统计学结论偏于保守，故Huynh和Feldt提出用平均调整值 值进行调整。 值的计算公式为 式（10-3）中的g是对受试对象的某种特征（如性别或年龄）进行分组的组数，n是每组的观察例数。 2．调整规则 只对具有重复测定性质的时间效应的F值的自由度，和处理时间交互作用的F值的自由度进行调整。由于F值有两个自由度v1和v2，调整的分子自由度 ，分母自由度 具体计算时可用 或 代替ε。用调整所得的 及 的值查 F界值表，得 。由于ε≤1，所以调整后的F临界值要大于调整前的F临界值。 第二节 单因素重复测量资料的方差分析 一、单因素重复测量资料的例子 例10-1 一项关于不同药物治疗心律失常效果的对比研究。对9例经常出现心室早搏的病人于用药前测定其心率后进行随机化给药。一部分病人按A药→安慰剂（C药）→B药的顺序给药，另一部分病人按B药→安慰剂（C药）→A药的顺序给药。安慰剂（C药）持续一周作为药物后效的清除期。比较用药前与各种药物及A药与B药之间的心律差别。表10-1列出9名受试病人在用药前、安慰剂（C药）期及用药（A与B）期的心率。 上表的下部及右侧所列的数据均从表体部分的原始测量值派生而来，是为方差分析而准备的。本例药物水平数a＝4，每组观察例数n＝9，观察值总个数N＝a×n＝36。 二、方差分析的步骤 1．提出检验假设