《重复测量资料方差分析.ppt》由会员分享，可在线阅读，更多相关《重复测量资料方差分析.ppt（74页珍藏版）》请在金锄头文库上搜索。
1、中国疾病预防控制中心中国疾病预防控制中心重复测量设计方差分析重复测量设计方差分析中国疾病预防控制中心中国疾病预防控制中心方方差差分分析析单因素：完全随机设计单因素：完全随机设计两因素：随机区组设计两因素：随机区组设计多因素：析因、拉丁方、正交设计多因素：析因、拉丁方、正交设计重复测量设计重复测量设计单组单组多组多组中国疾病预防控制中心中国疾病预防控制中心第一节 重复测量资料数据特征第二节 重复测量资料的方差分析第三节 重复测量方差分析案例及SAS实现3中国疾病预防控制中心中国疾病预防控制中心第一节 重复测量资料数据特征中国疾病预防控制中心中国疾病预防控制中心重复测量的定义重复测量的定义 重重复
2、复测测量量(repeated measure)是是指指对对同同一一研研究究对对象象的的某某一一观观察察指指标标在在不不同同场场合合（occasion，如时间点）进行的多次测量。如时间点）进行的多次测量。中国疾病预防控制中心中国疾病预防控制中心（1）随）随时间变化化趋势，如生，如生长发育模型育模型（2）时间上的稳定性）时间上的稳定性 分析，如环境监测中对同一地区分析，如环境监测中对同一地区 不同时间点的污染情况进行动态监测不同时间点的污染情况进行动态监测（3）结果变量的可重复性，如实验室中某仪器或实验的）结果变量的可重复性，如实验室中某仪器或实验的 准确性准确性 重复测量资料的研究目的重复测量资
3、料的研究目的中国疾病预防控制中心中国疾病预防控制中心重复测量重复测量设计设计重复测量重复测量设计设计（repeated measurement design）受试者内设计（受试者内设计（within-subject design）是指同一观察对象的某观察指标在不同时间点上进行多次观是指同一观察对象的某观察指标在不同时间点上进行多次观察，医学研究中十分常见。察，医学研究中十分常见。如：儿童生长发育中多次测量如：儿童生长发育中多次测量 临床试验中患者接受治疗后在不同时间测量结果临床试验中患者接受治疗后在不同时间测量结果 流行病学研究中接受干预后多次随访观测结果流行病学研究中接受干预后多次随访观测结
4、果中国疾病预防控制中心中国疾病预防控制中心重复测量资料重复测量资料重复测量资料，是重复测量设计中所获取的多次测量资料。重复测量资料，是重复测量设计中所获取的多次测量资料。不同称谓：不同称谓：纵向数据纵向数据(Longitudinal data),如在生长发育研究中如在生长发育研究中重复测量资料重复测量资料(Repeated measures data),如在临床试验研究中如在临床试验研究中面板数据面板数据(Panel data),如在社会学研究中如在社会学研究中中国疾病预防控制中心中国疾病预防控制中心（1）平衡）平衡资料：如果每一个体重复料：如果每一个体重复测量的量的时间和次数相同，和次数相同
5、，称称为平衡平衡资料。料。（2）非平衡）非平衡资料：不同个体重复料：不同个体重复测量的量的时间点不同，或重点不同，或重复次数不同，称复次数不同，称为非平衡非平衡资料。料。重复测量资料重复测量资料中国疾病预防控制中心中国疾病预防控制中心 同一受试对象的同一观察指标在不同时间点上进行多次测量所得的资料，常用来分析该观察指标在不同时间点上的变化。有时是从同一个体的不同部位（或组织）上重复测量获得的指标的观测值。目的：就是比较不同时间点动态变化趋势的特征目的：就是比较不同时间点动态变化趋势的特征重复重复测量测量资料资料中国疾病预防控制中心中国疾病预防控制中心(1)同一同一观察察对象的重复象的重复测量量
6、值之之间是非独立的，后一次是非独立的，后一次测量的数据可能受前一次量的数据可能受前一次测量量结果的影响。果的影响。这不符合不符合传统的的统计分析方法中关于独立性的假分析方法中关于独立性的假设。(2)观察指察指标在所在所测量的量的时间范范围内可能成内可能成趋势性性变化。化。(3)观察察值的的变异来源异来源较多，有来自个体内的多，有来自个体内的变异，个体异，个体间的的变异，或更高水平上的异，或更高水平上的变异；异；变异可能与异可能与时间有关，有关，或其他或其他协变量有关。量有关。重复测量资料的特点重复测量资料的特点中国疾病预防控制中心中国疾病预防控制中心重复测量资料的数据特征重复测量资料的数据特征
7、 当对同一受试对象在不同时间重复测量次数p3时，称为重复测量设计或重复测量数据。测测 量量 时时 间间 点点受试者受试者 1 2 p1 y11 y12 y1p2 y21 y22 y2p:n yn1 yn2 y n p中国疾病预防控制中心中国疾病预防控制中心中国疾病预防控制中心中国疾病预防控制中心实例举例1每一根线代表每一根线代表1只兔子只兔子中国疾病预防控制中心中国疾病预防控制中心实例举例2每一根线代表每一根线代表1位病人位病人中国疾病预防控制中心中国疾病预防控制中心典型的重复测量资料典型的重复测量资料研究对象研究对象时间点时间点1时间点时间点2时间点时间点m1y11y12y1m2y21y22
8、y2m3y31y32y3mnyn1yn2ynm中国疾病预防控制中心中国疾病预防控制中心典型的重复测量资料典型的重复测量资料分组因素分组因素研究对象研究对象时间点时间点1时间点时间点2时间点时间点m11y111y112y11m12y121y122y12m1n121y211y212y21m22y221y222y22m2n2p1yp11yp12yp1mp2yp21yp22yp2mpn3中国疾病预防控制中心中国疾病预防控制中心实际中：重复测量资料比独立资料更多见实际中：重复测量资料比独立资料更多见临床研究中，需要观察病人在不同时间的某些生理、生化或病理指标的变化趋势，研究不同时间或疗程的治疗效果。流行
9、病学研究中，观察队列人群在不同时间上的发病情况。研究不同职业、性别人群实施某种控制后，不同时间的多次效果考察。卫生学研究中，纵向观察儿童生长发育规律等，不同地区和环境营养状况。中国疾病预防控制中心中国疾病预防控制中心p简单处理简单处理: 对平衡的重复测量资料，分别在各时间点上进对平衡的重复测量资料，分别在各时间点上进行分析行分析 p集合分析集合分析: 将各个体的几次不同观察值相加，得到该个体将各个体的几次不同观察值相加，得到该个体的一个综合值，再进行比较分析的一个综合值，再进行比较分析 p两阶段分析两阶段分析: 为每个研究对象独立拟合时间函数，再对拟为每个研究对象独立拟合时间函数，再对拟合结果
10、求均数。合结果求均数。 重复测量资料的常见误用重复测量资料的常见误用中国疾病预防控制中心中国疾病预防控制中心减少样本含量减少样本含量控制个体变异控制个体变异非实验因素（干扰因素）非实验因素（干扰因素）滞留效应（滞留效应（carry-over effectcarry-over effect）潜隐效应（潜隐效应（latent effectlatent effect）学习效应（学习效应（learning effectlearning effect）重复测量设计优缺点重复测量设计优缺点中国疾病预防控制中心中国疾病预防控制中心21重复测量设计与配对设计的区别？重复测量设计与配对设计的区别？中国疾病预防控
11、制中心中国疾病预防控制中心1. 随机区组设计要求每个区组内实验单位彼此独立区组号区组号 A营养素营养素 B营养素营养素 C营养素营养素 1 50.10 58.20 64.50 2 47.80 48.50 62.40 3 53.10 53.80 58.60 4 63.50 64.20 72.50 5 71.20 68.40 79.30 6 41.40 45.70 38.40 7 61.90 53.00 51.20 8 42.20 39.80 46.20 表表 A、B、C 3种营养素喂养种营养素喂养 小白鼠所增体重（克）小白鼠所增体重（克）处理因素只能在区组内随机分配每个实验单位接受处理是不同的见
12、右表：重复测量设计与随机区组设计的区别？重复测量设计与随机区组设计的区别？中国疾病预防控制中心中国疾病预防控制中心2. 重复测量设计区组内实验单位彼此不独立中国疾病预防控制中心中国疾病预防控制中心同一受试的血样重复测量结果是高度相关提示：分析存在一定的复杂性。提示：分析存在一定的复杂性。中国疾病预防控制中心中国疾病预防控制中心小结 重复测量资料p数据：对同一受试对象的某个观察指标进行连续观测所得数据：对同一受试对象的某个观察指标进行连续观测所得p至少有两个因素：处理因素、时间因素；至少有两个因素：处理因素、时间因素； p重复测量的试验结果重复测量的试验结果按时间顺序按时间顺序固定排列，不能像随
13、机区固定排列，不能像随机区组设计的处理那样随机排列。组设计的处理那样随机排列。观察指察指标在所在所测量的量的时间范范围内可能成内可能成趋势性性变化；化；p同一个体不同时间测量值之间同一个体不同时间测量值之间高度相关；高度相关；p观察察值的的变异来源异来源较多，有来自个体内的多，有来自个体内的变异，个体异，个体间的的变异，或更高水平上的异，或更高水平上的变异；异；变异可能与异可能与时间有关，或其他有关，或其他协变量有关。量有关。中国疾病预防控制中心中国疾病预防控制中心第二节 重复测量资料方差分析中国疾病预防控制中心中国疾病预防控制中心想一想？想一想？同一观察单位具有多个观察值，而这些观察值来自同
14、一受试对象的不同时点（部位等）,这类数据间往往有相关性存在，违背了方差分析要求数据满足独立性基本条件。中国疾病预防控制中心中国疾病预防控制中心将总变异总变异分解为： 个体间（个体间（between subjects）变异）变异 个体内个体内（within subject）变异变异，其中个体内变异是与重复因素有关的变量。重复测量资料的方差分析总思想重复测量资料的方差分析总思想中国疾病预防控制中心中国疾病预防控制中心平均值之间的多重比较 先采用配对t检验方法，计算需比较的两两均数的t统计量，然后将这些样本统计量t值与Bonferroni临界t值进行比较。确定P值是否大于中国疾病预防控制中心中国疾病
15、预防控制中心重复测量资料重复测量资料多变量重复测量方差分析单变量重复测量方差分析重复测量资料分类repeated measurement data中国疾病预防控制中心中国疾病预防控制中心多变量重复测量方差分析单变量重复测量方差分析重复测量资料分类重复测量资料分类repeated measurement data两组单组重复测量资料重复测量资料中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析模型重复测量方差分析模型 式中，为总均值，g表示第g个处理水平的效应，j表示第j个测量时间点的效应，表示第g组第i个受试者的效应，该效应为随机效应，表示处理因素第g个水平在第j个测量时间点上的效应，也
16、就是处理因素与时间点的交互效应，为随机误差。中国疾病预防控制中心中国疾病预防控制中心单变量重复测量方差分析单变量重复测量方差分析1. 单组重复测量 指同一组内（或接受同一种处理）的多个受试者，在多个时间点上的反应变量所作的测量，又称为单变量重复测量。测测 量量 时时 间间 点点受试者受试者 1 2 p1 y11 y12 y1p2 y21 y22 y2p:n yn1 yn2 y n p中国疾病预防控制中心中国疾病预防控制中心2.多组重复测量（多组并不等于多因素） 指将受试者按处理的不同水平分为几个组，对这些组内的每一受试者，都在不同时间点对他们的反应变量进行测量。中国疾病预防控制中心中国疾病预防
17、控制中心p重复测量资料方差分析的条件： 1 1. . 正态性正态性 处理因素的各处理水平的样本个体之间是相互处理因素的各处理水平的样本个体之间是相互独立独立的的随机随机样样本，其总体均数服从本，其总体均数服从正态正态分布；分布； 2. 2. 方差齐性方差齐性 相互比较的各处理水平的总体方差相等，即具有方差齐同相互比较的各处理水平的总体方差相等，即具有方差齐同 3. 3. 各各时时间间点点组组成成的的协协方方差差阵阵(covariance (covariance matrix)matrix)具有球形性(sphericity)特征。Box（1954）指出，若球形性质得不到满足，则方差分析的）指出，
18、若球形性质得不到满足，则方差分析的F值是有偏的，这值是有偏的，这会造成过多的拒绝本来是真的无效假设会造成过多的拒绝本来是真的无效假设(增加增加型错误型错误)。 重复测量资料方差分析的前提条件重复测量资料方差分析的前提条件中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析的球性假定重复测量方差分析的球性假定p方差是指在某一时点上测定值方差是指在某一时点上测定值变异性的大小，而异性的大小，而协方差是指在两个不同方差是指在两个不同时点上点上测定定值相互相互变异性异性的大小。的大小。p如果在某个如果在某个时点上的取点上的取值不影响其他不影响其他时点上的取点上的取值，则协方差方差为0，反之，反之，
19、则不不为0。p由方差由方差协方差构成的矩方差构成的矩阵称称协方差方差阵。中国疾病预防控制中心中国疾病预防控制中心p设k、l为两个两个测定定时点，点，代表代表协方差方差阵中的元素。当中的元素。当kl时为方差，方差，kl时为协方差。共有方差。共有a个个测定定时点，将点，将这a个方差和（个方差和（a-1）/2个个协方差方差排成排成协方差方差阵V为：协方差阵的球形假定是指该矩阵主对角线元素（方差）协方差阵的球形假定是指该矩阵主对角线元素（方差）相等、非主对角线元素（协方差）为零。相等、非主对角线元素（协方差）为零。重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制
20、中心所有两两时间点变量间差值对应的方差相等所有两两时间点变量间差值对应的方差相等对于yi与yj两时间点变量间差值对应的方差可采用协方差矩阵计算为：重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心协方差阵协方差阵 A1A2A3A4A11051015A25201520A310153025A415202540s1-22 = 10 + 20 - 2(5) = 20s1-32 = 10 + 30 - 2(10) = 20s1-42 = 10 + 40 - 2(15) = 20s2-32 = 20 + 30 - 2(15) = 20s2-42 = 20 + 40
21、 - 2(20) = 20s3-42 = 30 + 40 - 2(25) = 20所有两两差值对应的方差相等，说明满足球形假定。重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心p球性假定常用球性假定常用Mauchly法法进行行检验pMauchly法法检验的的P值若大于研究者所若大于研究者所选择的的显著性水准著性水准时，说明明协方差方差阵的球形性的球形性质得到得到满足。否足。否则，必，必须对与与时间有关的有关的F统计量的分子、分母自由度量的分子、分母自由度进行行调整，以整，以便减少犯便减少犯I类错误的概率。的概率。p注意：注意：该该法法在小样本时很难检
22、出差异在小样本时很难检出差异。重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心对角线元素（方差）相等非主对角线元素（协方差）为零00中国疾病预防控制中心中国疾病预防控制中心球形对称的检验球形对称的检验用Mauchly法检验协方差阵是否为球形 H0：资料符合球形要求， H1：资料不满足球形要求检验的P值若大于研究者所选择的显著性水准时，说明协方差阵的球形性质得到满足。中国疾病预防控制中心中国疾病预防控制中心球形条件不满足怎么办球形条件不满足怎么办常有两种方法可供选择：1.对重复测量ANOVA检验结果中与时间有关的F值的自由度进行调整（调小）2.采用MA
23、NOVA（多变量方差分析方法）(1) Geenhouse-Geisser调整系数(G-G)(2) Huynh-Feldt调整系数(H-F)分子自由度分子自由度分母自由度分母自由度中国疾病预防控制中心中国疾病预防控制中心自由度调整方法1中国疾病预防控制中心中国疾病预防控制中心自由度调整方法2中国疾病预防控制中心中国疾病预防控制中心调整规则 调整后的F临界值较原先大，提高了拒绝H0的门槛。减少了犯I类错误的概率。中国疾病预防控制中心中国疾病预防控制中心多变量方差分析统计量目前多数统计软件会给出四个统计量目前多数统计软件会给出四个统计量：pWilks lambdapPillais trace、pHo
24、telling-Lawley tracepRoys largest root。中国疾病预防控制中心中国疾病预防控制中心多变量方差分析统计量pWilkslambda指指标的含的含义是最明确的，其是最明确的，其计算公式算公式为|H|/(|H|+|E|)。反映。反映了了组内方差占内方差占总方差（方差（组间方差方差+组内内方差）的比例方差）的比例。该值越小，越小，组间方差越大，意味着方差越大，意味着组间差异越差异越大大。p其它三个指其它三个指标的含的含义并不是很并不是很直直观：Hotelling-Lawleytrace是是HE-1矩矩阵的的对角角线方差之和。方差之和。Pillaistrace是是H(H
25、+E)-1矩矩阵的的对角角线方差之和方差之和。Royslargestroot是是HE-1矩矩阵的最大的最大特征特征值。中国疾病预防控制中心中国疾病预防控制中心第三节 重复测量方差分析案例及SAS实现中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析的SAS语句proc glm ;class 分组变量;model 分析变量=分组变量;lsmeans 分组变量;means 分组变量;repeated 重复测量变量名 /;中国疾病预防控制中心中国疾病预防控制中心 repeated 重复测量变量名重复测量变量名 /;【repeated】后必】后必须指定重复指定重复测量量变量名，量名，测量水平可量
26、水平可选，转换项也是可也是可选的。如：的。如：repeatedtime4;表示指定重复表示指定重复测量量变量名量名为time，测量量水平（即重复次数）水平（即重复次数）为4次，中次，中间均以空格隔开。其中均以空格隔开。其中time由由读者自行命名，也可者自行命名，也可为其它名称，如其它名称，如repeatedt4，repeatedmeasure4等均可。等均可。重复测量方差分析的SAS语句中国疾病预防控制中心中国疾病预防控制中心repeated语句后的句后的主要用于各次主要用于各次测量之量之间的比的比较，常用的有两，常用的有两种种转换方式：方式：contrast和和polynomial。con
27、trast表示指定一个参照水平，其它均与参照水平比表示指定一个参照水平，其它均与参照水平比较。如如repeatedtime4contrast(1);表示共有表示共有4次重复次重复测量，以第量，以第1次次测量量为参参照，第照，第2、3、4次次测量分量分别与第与第1次次测量水平量水平进行比行比较。polynomial表示正交多表示正交多项式式对比，即正交比，即正交转换后判断曲后判断曲线的的n阶水平水平是否有是否有统计学意学意义。如：。如：repeatedtime4polynomial;表示表示4次重复次重复测量，量，分分别对一次、二次和三次曲一次、二次和三次曲线分析是否有分析是否有统计学意学意义，
28、也就是判断，也就是判断测量量值随随时间变化大致呈何种化大致呈何种趋势。重复测量方差分析的SAS语句中国疾病预防控制中心中国疾病预防控制中心printe该该选选项项可可输输出出一一些些相相关关矩矩阵阵（不不是是我我们们所所关关心心的），更重要的是给出球性检验的），更重要的是给出球性检验结果。结果。summary 输输出出重重复复测测量量各各水水平平间间的的比比较较结结果果。如如果果没没有有该该 选选 项项 ， 即即 使使 前前 面面 指指 定定 了了 contrast或或polynomial，也不会显示比较结果，也不会显示比较结果重复测量方差分析的SAS语句中国疾病预防控制中心中国疾病预防控制中
29、心例某研究例某研究评价一医价一医疗仪器在治器在治疗颈椎病中的椎病中的疗效，效，以以VASVAS评分作分作为评价指价指标，共治，共治疗2 2个个疗程。从程。从该研研究中随机究中随机选择2020例研究例研究对象，表象，表1 1列出了列出了2020例研究例研究对象的治象的治疗前、治前、治疗1 1个个疗程、治程、治疗2 2个个疗程的程的VASVAS评分分结果，果，试分析分析VASVAS评分的分的变化化趋势。重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心ID疗前前1个个疗程后程后2个个疗程后程后ID疗前前1个个疗程后程后2个个疗程后程后17861177424431276635541376
30、54665146545441157666665165457666176538877186659444197661087320875重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心dataexample1;inputvas1-vas3;cards;procglm;modelvas1-vas3=/nouni;repeatedtime3contrast(1)/printesummary;run;重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量
31、方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心p从本例分析从本例分析结果可以看出，三次果可以看出，三次VAS评分之分之间存在一定的存在一定的相关性，不相关性，不应采用采用单变量方差分析。重复量方差分析。重复测量方差分析量方差分析结果提示，三次果提示，三次评分之分之间存在一定的存在一定的趋势。p如果我如果我们想想进一步了解存在什么一步了解存在什么样的的趋势，可以通，可以通过在在repeated语句中指定句中指定polynomial来来实现。procglm;model
32、vas1-vas3=/nouni;repeatedtime3polynomial/printesummary;run;重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心图1观测时间VAS评分分疗前第1疗程后第2疗程后01234567重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心例例为探索磁共振探索磁共振扩散加散加权成像在胃成像在胃肠癌肝癌肝转移化移化疗疗效中的价效中的价值，某研究者，某研究者观察了察了20例化例化疗反反应良良好和化好和化疗反反应
33、不良的胃不良的胃肠癌患者分癌患者分别在化在化疗前、化前、化疗后后3天、天、7天、天、42天的成像，并天的成像，并计算表算表观扩散系数散系数ADC值。表中。表中group=0代表反代表反应良好良好组，group=1代代表反表反应不良不良组。试比比较两两组人群的人群的ADC值差异是否差异是否有有统计学意学意义。重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心dataexample2;inputgadc1-adc4;cards;01.86 1.821.831.8601.94 1.771.881.7311.13 1.211.2
34、41.2510.91 1.051.191.16;procglm;classg;modeladc1-adc4=g/nouni;repeatedtime4(03742)polynomial/printesummary;run;注意，当注意，当repeated语句指定语句指定polynomial变换时，如果变换时，如果观测时间间隔不等，必须在观测时间间隔不等，必须在括号中加入具体的时间点。括号中加入具体的时间点。因为正交多项式系数在等间因为正交多项式系数在等间隔和不等间隔时是截然不同隔和不等间隔时是截然不同的。如果不指定，则默认为的。如果不指定，则默认为等间隔，计算的正交多项式等间隔，计算的正交多项式结果会有偏差。结果会有偏差。重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心图2观测时间点点ADC值d0d3d7d420.00.20.40.60.81.01.21.41.6反应良好反应不良重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心
