重复测量资料的统计分析及SAS程序实例.pdfVIP

，146· 主蕉圣垒瑟j主燃!釜垒墨筮堑鲞蔓2麓 析因设计重复测量资料的统计分析及SAS程序实例 擒伟伟1贺 佳1矗 踩舜杰2 【提要l 目的探讨析因设计重复测量资料的统计分析方法。方法通过构造统计模型，并结合实例给出SAS程 序涞实现析因设计重复测缀资料的统计分析。结粜给出了一般线性模魁和混合线性模型对处理因索和时间因素进行 比较的具体做法。结论畿数据之间相关性较为复杂时。混合线性模型照加适合予析因设计的熬复测量资料的统计分 糖。 【关键词l 辑嚣竣谤霪复瀵量资辩一般线臻摸登混金线毪旗滋 模型可构遣为【i1： 析因设计(factorialexperimentdesign)又称宪全交 叉分组设计，属于多因素、多水平、单效应变量设计方 法，它不仅可以检验每一因素各水平之间的效应差异． +(J9y)舛十(q』9y)曲f+e。甜f 丽且可以检验各因素之阅的交互俸用。重复测璧资料 measurement ≤repeated data)是揍圈一观察对象在不圆分囊代裘各个因素酶不同水平。 时闻点上多次测量同一糖标所得赘的资料，对予同一 其中蛳f为响应变量；∥为总平均值；‰、忍、 观察对象而言，多次测景值之间可能存在相关性。本 (印)曲分别代表因素A、B及其交互作用的效应，属固 文就在一个析因设计的试验中，对观察对象进行不同 时间点上多次测量所获资料的统计分析进行简要探 因素(时间)及其与因素A、B交互作用的效应，也属于 讨。 固定效应；Oi{面)为观察对象的效应，属随枧效应；￡蹦 为隧撬误差。对此模型的限麟条俘尧： 统计方法 l t l l 以最简单的两因素两水平析因设计为例，在5个 时间点进行重复测量，数据结构见表1。 聋而3o，互(口7)巧2o，磊(口)，)可20，蔷(∥)彬。 表1析因设计重复测量资料的数据结构 10o，蚕(缈)西。0，善(秽)nb／2=j； ∑b4互V，0 (口弦)嘶。o．着J_^ N(O，t12)，蠡‘面)与￡删褶互独立。 在上述假定下，响应变摄比“的方差为： Var(y,曲q)=口j+口2 根据此模型进行重复测爨方差分析时，先将总变 异(即总离均差平方和)分解为缀阕合计(观察对象)秘 ’；A、B为处理厦幕，各有掰个承乎(以0、1代表)，％～戳为5个 瘸察瓣瓣焘 缓内会计(霪复溅量)【引，遵蠢文献称之力受试蠹阕因 在重复；燹|量资料的统计分析中，常常可以看到的 素和受试者内因素【¨，然后再进行组闻秘组内的变异 做法有采用随机区组的方差分析，或者以治疗前的基 分解。需要注意的是在计算处理因素和重复测量因素 线数据为协变量做协方差分析。但由于重复测爨资料 (时间)的统计量F时，其分坶均方是不同的，前者使 并不满足独立性的假定，上述方法得到的结论可能是 用组间误藏的均方，后者则使用组内误差的均方。 不当甚至是错误的。所以，对这类资料的分析，应用一 在一般线性模型