重复测量资料方差分析.docx

重复测量资料方差分析 重复测量(repeated measure)是指对同一观察对象的同一观察指标在不同时 间点上进行的多次测量，用于分析该观察指标在不同时间上的变化特点。这类测 量资料在临床和流行病学研究中比较常见，例如，为研究某种药物对高血压病人 的治疗效果，需要定时多次测量受试者的血压，以分析其血压的变动情况。 1、重复测量资料方差分析中自由度调整方法 1?调整系数的计算 有两个调整系数，第一个是 Greenhouse-Geisser调整系数 G ?)，计算 公式为 (s：s (s： s2) 2 2(s )] 2 2 (s )] (a 1)[ 6) (2a)( (sk)) k l k 2 2 2 2 式中中的Ski是协方差矩阵中的第 k行第I列元素，s = ( Ski)/a是所 k I 2 2 2 2 2 有元素的总平均值，Skk ( si )/a是主对角线元素的平均值，Sk ( Sk」/a I I 是第k行的平均值。？的取值在1.0与1/(a-1)之间。 第2个系数是Huynh-Feldt调整系数 (H F _)。研究表明，当 真值在0.7 以上时，用？进行自由度调整后的统计学结论偏于保守，故 Huynh和Feldt 提出用平均调整值 一值进行调整。一值的计算公式为 - ng(a 1)? 2 (a 1)[( n 1)g (a 1)?] 式中中的g是对受试对象的某种特征(如年龄或性别)进行分组的组数， n 是每组的观察例数。当 —1.0时，取—=1.0。 2.调整规则 只对具有重复测定性质的时间效应的 F值的自由度，和处理时 间交互作用的F值的自由度进行调整。 由于F值的有两个自由度 V1和V2,调整的 分子自由度1 1 分母自由度 2 2 。具体计算时可用或 一代替。用 调整所得的1及2的F值查临界值表，得F()。由于 1.0,所以调整后的 (1,2) F临界值要大于调整前的 F临界值。 2、单因素重复测量资料的方差分析 单因素重复测量资料的例子 一项关于不同药物治疗心律失常效果的对比研 究。对9例经常出现心室早搏的病人于用药前测定其心率后进行随机化给药。一 部分病人按A药 安慰剂（C药） B药的顺序给药，另一部分病人按 B药 安慰剂（C药） A药的顺序给药。安慰剂（C药）持续一周，作为药物后效的 清除期。比较用药前与各种药物及 A药与B药之间的心律差别。图 4-12列出9 名受试病人在用药前、安慰剂（ C药）期及药（A与B）期的心率。 病人号 1 2 3 4 5 6 7 8 9 永药前 94 57 81 82 67 78 87 82 90 A药 67 52 74 59 65 72 75 68 74 B药 90 69 69 71 74 80 106 76 82 C药 67 55 73 72 72 72 74 59 80 图4-12心室早搏病人在用药前后的心率 方差分析的步骤 1.提出检验假设 检验假设为: H0 : pl = JJ2= B= JJ4 ； H1： pi / ph，至少有 个不等式成立。 2?计算离均差平方和、 自由度及均方 有总离均差平方和、 处理因素离均差 平方和、受试对象间离均差平方和及受试对象内离均差平方和等。计算公式为: （1）总离均差平方和SS总及总自由度 总的计算 a n —2 2 SS总 （Yj 丫） S T /N，总 N 1 j 1 i 1 （2）处理因素的离均差平方和a _ （2）处理因素的离均差平方和 a _ —2 SS处理 n （Yj Y） j 1 SS处理及自由度 处理的计算 丄（\2） T2 n j 1 处理 （3）受试对象间离均差平方和SS （3）受试对象间离均差平方和 SS对象间及自由度对象间的计算 ss对象间a (Yi Y)i 1Ti对象间 n 1受试对象内离均差平方和ss对象内 ss对象间 a (Yi Y) i 1 Ti 对象间 n 1 受试对象内离均差平方和 ss对象内及自由度对象内的计算 SS对象内 n a (Yj i 1 2 Y) (Si i 1 Ti2 a)， 对象内n(a 1) 误差的离均差平方和 ss误差与自由度误差的计算 SS误差 SS总 SS处理 SS对象间， 误差 (n 1)(a 1) 根据以上4种离均差平方和与自由度计算所得的均方见表 10-2. 3?计算F值 由于是处理因素的统计学检验，故只计算处理因素的 F值。 F处理 MS处理/ MS误差，F处理服从1 处理与2 误差的F分布 本例，在DPS 本例，在DPS数据处理系统中，按图 4-12方式编辑、定义数据块，然后执 行“试验统计” t 行“试验统计” t “重复测量方差分析” t “单因素分析”功能，得到计算结果 如下。 计算结果 当前日期 02-8-16 8:42:12 处理 样本数