重复测量资料的多水平模型(南京).ppt

重复测量资料的多水平模型南京医科大学 陈峰 2007.05.16~20Multilevel Models for Repeated Measurement Data1内容概要o重复测量数据的特点与结构o3个例子n平行组设计n交叉设计n非线性回归o实例分析o重复测量资料的相关结构2重复测量设计o重复测量设计(repeated measurement design)o受试者内设计(within-subject design)o是指同一观察对象的某观察指标在相继的不同 时间点上进行的多次观察 3重复测量资料的特点n在相继的不同时间点上进行的多次观察不是 随机确定的；n重复测量值之间是非独立的4Other terms for repeated measures dataØLongitudinal data, in growth studyØRepeated measures data, in experimental studyØPanel data, in social research5两种设计o3个处理组：A、B、Co4个不同观察时间：1，2，3，4Ø12个试验单元T1T2T3T4AXXXXBXXXXCXXXX6两种设计o3个处理组：A、B、Co4个不同观察时间：1，2，3，4Ø完全随机区组设计n1~12个动物随机分配到12个试验单元 T1T2T3T4 A81311 B2479 C126105Ø重复测量n1~3个动物随机分配到3个组 T1T2T3T4 A1111 B2222 C33337重复测量的数据结构o个体：2水平o重复测量值：1水平 8多中心临床试验中重复测量o中心(医院)：3水平o个体：2水平o重复测量值：1水平9重复测量的数据结构(固定重复时间)t1t2t3t4t5t6 ID1 ××××× ID2 ××××× ID3 ××××× ID4 ××××10重复测量的数据结构(不固定重复时间)ID1 t11 t12 t13 t14 t15 t16 ID2 t21 t22 t23 t24 ID3 t31 t32 ID4 t41 t42 t43 t44 t45 t46 11重复测量概念的推广o重复测量的概念不仅仅是时间上(temporal)的 ，也可以推广到空间(spatial)。

o例如：n同一母鼠所生的仔鼠；n同一家庭的不同成员；n同一患者的两个不同膝关节；n同一肿瘤患者的不同肿块；n同一条河流的不同采样点12重复测量实例1（平行组设计）o观察某溶栓药治疗20名急性脑梗死患者的疗效 ，采用双盲、随机、安慰剂平行对照，每组各 10例观察指标为神经系统体症评分(MDNS) 分别于疗前、疗后每周进行随访，共8周o目的：试验组和对照组的MDNS是否不同？13重复测量实例1：原始资料形式表1 20名急性脑梗死患者治疗不同时期神经系统体症总分id trea tage观察时间(time)0123456781027107 106 106 108 108 112 112 112 1122021107 106 106 106 106 112 112 114 1163021100 100 100 106 109 108 114 116 1164036107 106 106 107 106 111 112 117 1095017110 111 112 112 113 113 113 116 116 ………………………………14重复测量实例1：两组MDNS在时间上的趋势Time(Weeks)0123456788090100110120130140150Time(Weeks)012345678809010011012013014015015重复测量实例2（交叉设计）o研究国产罗格列酮(rosiglitazone maleate tablets)与进口 罗格列酮的生物等效性(bioequivalence)，选择24名志 愿者，进行4×4的交叉设计。

试验安排用下列4×4的拉 丁方： ADBC BACD CBDA DCABnA：进口罗格列酮2mg片剂；nB：国产罗格列酮2mg片剂；nC：进口罗格列酮4mg片剂；nD ：D国产罗格列酮4mg片剂16重复测量实例2：原始资料形式ID顺序第一阶段第二阶段第三阶段第四阶段1DCAB 884.27 905.092330.771936.98 2CBDA 919.502201.98 855.891939.36 3ADBC1738.12 901.701889.72 870.93 4BACD2000.292350.58 952.86 955.46 5CBDA 823.391864.97 710.061372.77 ………………国产罗格列酮与进口罗格列酮4×4交叉试验结果 17重复测量实例3（单组设计，趋势分析）o某种属100只雌性大鼠从出生到第10周的体重 (g)，试建立生长模型18重复测量实例3：原始资料形式idTime（day）04714212835497016.19.6132650799917620526.2111529527810318017436.37.71426457611816819346.3121428567512320520056.3111632548711620120266.3101330517810919221976.4111732618811621222786.41115345484121209198100只雌性大鼠从出生到第10周的体重(g)19Weight(g)Time(day)0 4 7 14 213549701000.0100.0200.0300.0100只雌性大鼠从出生到第10周的体重(g)变化趋势20重复测量实例4（单组设计，趋势分析）o杨珉(1994)对成都市三个妇幼保健院出生的 202 例新生儿的生长发育进行了为期24个月的 随访(最长779天)，其中男婴112例，女婴90例 。

o重复测量的次数由8 到17 次不等，共2219个测 量数据o测量的指标有体重(kg)，身长(cm)，头围(cm) ，坐高(cm)，胸围(cm)等指标 21202婴儿从出生到24月的体重(kg)变化趋势图223 资料的分析o传统方法及其弊端n对平衡的重复测量资料，分别在各时间点上进行分析孤 立地看待各时点数据，增加I型误差 n将各个体的几次不同观察值相加，得到该个体的一个综合 值，再进行比较分析(aggregated analysis)人为地减少误差 ，它忽略了对不同来源的变异的分析；未考虑观察值在时 间上的变化规律，也未考虑其他协变量与时间的交互作用 对结果的影响损失了很多信息n将n个患者的几次不同观察均作为因变量，时间以及其他变 量作为自变量，样本含量为 ，拟合线性(或广义线 性)模型将非独立数据当做独立数据看待，增加了power， 增加假阳性23实例1分析：分析资料形式24数据结构o水平2：个体o水平1：重复测量25该研究要回答以下几个问题o平均来说，试验组和对照组的MDNS是否不同 ？个体间是否有差异？o试验组和对照组的MDNS随时间的变化趋势是 否不同？这种趋势在个体间是否有差异？o年龄对结论是否有影响？26传统线性模型与2水平模型27传统模型与2水平模型的拟合结果28传统模型与2水平模型的结果比较效应参数传统模型2水平模型固定效应0(常数项)108.856(1.082)108.856(1.699)2 (treat) 13.744(1.531) 13.744(2.402)随机效应2水平(ID)－19.280(9.187)1水平(时点)105.437(11.114)86.157(9.633)-2LL1349.2791334.99429考虑随时间的变化（两组斜率相同）30考虑随时间的变化（不同组斜率不同）31考虑年龄的影响32考虑方差是否随时间而改变33效应参数系数(估计误差)模型A模型B模型C模型D模型E固定效应0(常数项)108.856(1.69 9)103.722(1.77 3)103.722(1.56 8)98.615(3.878) 98.029(3.444)2 (treat)13.744(2.40 2)1.558(2.507)1.558(2.217) -0.325(2.717) -0.541(2.351)1 (time)－1.283(0.127)1.283(0.207)1.283(0.127)1.283(0.207)4 (time*treat)－3.047(0.180)3.047(0.292)3.047(0.180)3.047(0.292)3 (age)－－－0.216(0.148)0.241(0.132)随机效应2水平(ID)19.280(9.187) 27.779(9.125) 21.792(7.781) 24.986(8.242) 18.217(6.652)－－0.171(0.731)－0.274(0.674)－－0.304(0.136)－0.304(0.136)1水平(时点)86.157(9.633) 9.673(1.081)7.390(0.883)9.673(1.081)7.390(0.883)-2LL1334.994985.093963.517983.058960.46520名急性脑梗塞患者治疗不同时期神经系统体症总分的2水平模型 34不同自变量取值组合时的预测值(未考虑2水平残差)35不同自变量取值组合时的预测值(考虑2水平残差)36o观察值随时间的改变而增加；o观察值随时间的改变幅度(斜率)在两个组间是 不同的，处理组的斜率大；o2水平上 有统计学意义，说明个体的截距是 有差异的；o2水平上 有统计学意义，说明个体的观察值 随时间改变的幅度(斜率)也是有差异的。

例1资料分析结论37重复测量实例2：分析资料形式38数据结构o水平2：个体o水平1：重复测量39主要问题o罗格列酮的生物等效性(bioequivalence)n国产2mg与进口2mg是否相同n国产4mg与进口4mg是否相同n国产4mg与2mg的差别与进口4mg与2mg的 差别是否相同40交叉设计的2水平模型41例2交叉设计的2水平模型估计[log(AUC)]42例2交叉设计的2水平模型估计[log(tmax)]43交叉试验资料的2水平模型效应参数log(AUC)log(Cmax) 系数(估计误差)系数(估计误差) 固定效应0(常数项) 7.669(0.049) 5.896(0.040) 1 (dose)-0.775(0.032) 0.754(0.030) 2 (dome)-0.054(0.032)-0.067(0.030) 3 (dose*dome) 0.033(0.046) 0.042(0.042) 4 (sequence(2/1))-0.042(0.063) 0.075(0.051)5 (sequence(3/1))-0.075(0.063)-0.050(0.051)6 (sequence(4/1))-0.042(0.063) 0.058(0.051)随机效应 2水平(个体)0.009(0.003)0.005(0.002) 1水平(阶段)0.013(0.002)0.011(0.002) －2LL-122.418-145.56744例2分析结论o国产罗格列酮与进口罗格列酮的生物等效性尚 不成立，国产药物需要进一步改进工艺。

45实例3数据格式46主要问题o描述生长规律n平均值的增长规律n变异的变化规律472水平生长模型 其中，H0表示固定效应模型的阶数，H1表示随机效应的阶数例如，H0=4，H1=2时： 48必要的准备o每个time, time2, time3, time4必须中心化 否则，方差部分可能估计不出来aver c4 b1 b2。