《重复测量设计和交叉设计资料.ppt》由会员分享,可在线阅读,更多相关《重复测量设计和交叉设计资料.ppt(58页珍藏版)》请在第壹文秘上搜索。
1、17.1 重复测量资料的数据特征重复测量资料的数据特征一、前后测量设计 最为常见,是重复测量设计的特例,亦称单组前后测量设计,即g=1,m=2,如表12-1。编 号 治疗前 治疗后 差值 1 2 3 4 5 6 7 8 9 10 130 124 136 128 122 118 116 138 126 124 114 110 126 116 102 100 98 122 108 106 16 14 10 12 20 18 18 16 18 18 X 126.2 110.2 16.0 S 7.08 9.31 3.13 表12-1 高血压患者治疗前后的舒张压(mmHg)编 号 哥特里罗紫法 脂肪酸水
2、解法 差值 d 1 0.840 0.580 0.260 2 0.591 0.509 0.082 3 0.674 0.500 0.174 4 0.632 0.316 0.316 5 0.687 0.337 0.350 6 0.978 0.517 0.461 7 0.750 0.454 0.296 8 0.730 0.512 0.218 9 1.200 0.997 0.203 10 0.870 0.506 0.364 表表3-3 两种方法对乳酸饮料中脂肪含量的测定结果两种方法对乳酸饮料中脂肪含量的测定结果(%)比较 1.配对设计中同一对子的两个实验单位可以随机分配处理,两个实验单位同期观察试验结果
3、,可以比较处理组间差别。前后测量设计不能同期观察试验结果不能同期观察试验结果,虽然可以在前后测量之间安排处理,但本质上比较的是前后差别,推论处理是否有效是有条件的,即假定测量时间对观察结果没有影响假定测量时间对观察结果没有影响。与配对设计设计的区别 2.配对 t 检验要求同一对子的两个实验单位的观察结果分别与差值相互独立,差值服从正态分布。前后测量设计前后两次观察结果通常与差值不独立,大多数情况第一次观察结果与差值存在负相关的关系,如表12-1中,治疗前舒张压与差值的相关系数为-0.602。如由表 12-1 计算,治疗前后舒张压的相关系数为 0.963,P0.01,用治疗前舒张压()X推论治疗
4、后舒张压()Y的回归方程为:49.5341.266YX,截距检验 P=0.014,回归系数检验0.01P。3.配对设计用平均差值推论处理的作用,而前后测量设计除了分析平均差值外,还可进行相关回归分析。区别点 配对设计 单组前后测量设计 两实验单位 可随机分配 N 观测时间 同期 两时间点 试验数据与差值关系 独立 N 分析指标 平均差值 平均差值、相关回归 推断 组间差别 前后差别 单组前后测量设计与配对设计的区别区别 二、设立对照的前后测量设计表 12-1 中高血压患者治疗后的舒张压平均下 降了 16 mmHg,虽 然经配 对 t 检 验:16.18,0.01tP,也未必能说明治疗有效,因为
5、住院休息、环境和情绪的改变同样可以使血压恢复平稳。因此,确定疗效的前后测量设计必须增确定疗效的前后测量设计必须增加平行对照加平行对照,如将 20 位轻度高血压患者随机分配到处理组和对照组,试验结果见表 12-2。表12-2 高血压患者治疗前后的舒张压(mmHg)处 理 组 对 照 组 顺序号 治疗前 治疗后 差值()d 顺序号 治疗前 治疗后 差值 ()d 1 2 3 4 5 6 7 8 9 10 130 124 136 128 122 118 116 138 126 124 114 110 126 116 102 100 98 122 108 106 11 12 13 14 15 16 17
6、 18 19 20 118 132 134 114 118 128 118 132 120 134 124 122 132 96 124 118 116 122 124 128 合 计 1262 1102 合 计 1248 1206 均 数 126.2 110.2 16.0 均 数 124.8 120.6 4.2 标准差 7.08 9.31 3.13 标准差 7.90 9.75 8.02 经检验处理组与对照组的差值 d 方差不齐(2212/6.58FSS,0.01P),不符合两均数比较 t 检验的前提条件。三、重复测量设计 当前后测量设计的重复测量次数m3时,称重复测量设计或重复测量数据。放
7、置 时 间(分)编 号 0 45 90 135 1 5.32 5.32 4.98 4.65 2 5.32 5.26 4.93 4.70 3 5.94 5.88 5.43 5.04 4 5.49 5.43 5.32 5.04 5 5.71 5.49 5.43 4.93 6 6.27 6.27 5.66 5.26 7 5.88 5.77 5.43 4.93 8 5.32 5.15 5.04 4.48 表12-3受试者血糖浓度(mmol/L)(g=1)球对称检验:215.44,5,0.01P 重复测量资料的反应变量可以为连续型、离散型或二分类,其中离散型与二分类反应变量资料的分析需要较复杂的统计方法
8、,如广义估计方程(generalized estimating equations,GEEs)等,有兴趣者可参考其它文献。实践中以连续型较常见,本章针对连续型重复测量资料给出若干常用的统计分析方法。在实际研究中,重复测量的时点选择是根据不同在实际研究中,重复测量的时点选择是根据不同的干预措施的背景、受试对象的特点和研究的问的干预措施的背景、受试对象的特点和研究的问题等多个因素确定的。题等多个因素确定的。一般而言,在同样的样本量情况下,采用重复测一般而言,在同样的样本量情况下,采用重复测量的研究设计可以获得更详细的信息,并且可以量的研究设计可以获得更详细的信息,并且可以分析受试对象在观察期间的动
9、态过程,因此在评分析受试对象在观察期间的动态过程,因此在评价指标是无创伤的情况下,应尽可能地采用重复价指标是无创伤的情况下,应尽可能地采用重复测量设计收集资料。测量设计收集资料。与随机区组设计的区别:1重复测量设计中“处理”是在区组(受试者)间随机分配,区组内的各时间点是固定的,不能随机分配,如表12-5,A、B两种处理随机分配给各个患者后,每个患者测量的时间相同。随机区组设计则要求每个区组内实验单位彼此独立,处理只能在区组内随机分配,每个实验单位接受的处理是不相同的,如表4-9。变异来源 自由度 SS MS F P 总变异 31 5.751 区组(受试者)7 2.828 0.361 7.77
10、 0.01 放置时间 3 2.959 0.986 5.85 0.01 误差 21 0.264 0.013 表12-4 表 12-3数据随机区组方差分析表 重复测量数据与随机区组设计数据很相似,如表12-3,而且同样可以计算出随机区组设计的方差分析表(表12-4)。2 重复测量设计区组内实验单位彼此不独立,如表 12-3,即同一受试者的血样重复测量结果是高度相关的,其相关系数见表 12-6。重复测量数据若用随机区组方差分析比较处理组间差异,前提条件是满足“球对称”假设(通过2检验)。放置时间 放置时间(分)(分)0 45 90 135 0 1 0.978*0.936*0.860*45 1 0.8
11、79*0.876*90 1 0.896*135 1 表12-6 表12-3各放置时间点血糖浓度的相关系数*P0.05,处理间时间趋势的差异无统计学意义;处理主效应的P0.56630.05,处理间差异也无统计学意义。故没有足够的证据可以推断两种治疗方案的疗效有差异。时间主效应的P0.0001,可以认为各个时间点的总体均数不全相同。例 12-3 将手术要求基本相同的 15 名患者随机分 3组,在手术过程中分别采用 A,B,C 三种麻醉诱导方法,在0T(诱导前)、1T、2T、3T,4T 五个时相测量患者的收缩压,数据记录见表 12-16。试进行方差分析。麻 醉 诱 导 时 相 诱导 方法 患者 序号
12、 0T 1T 2T 3T 4T A 1 120 108 112 120 117 A 2 118 109 115 126 123 A 3 119 112 119 124 118 A 4 121 112 119 126 120 A 5 127 121 127 133 126 B 6 121 120 118 131 137 B 7 122 121 119 129 133 B 8 128 129 126 135 142 B 9 117 115 111 123 131 B 10 118 114 116 123 133 C 11 131 119 118 135 129 C 12 129 128 121 1
13、48 132 C 13 123 123 120 143 136 C 14 123 121 116 145 126 C 15 125 124 118 142 130 表12-16 不同麻醉诱导时相患者的收缩压(mmHg)麻 醉 诱 导 时 相 诱导方法 0T 1T 2T 3T 4T 合 计(iA)A 605 562 592 629 604 2992 B 606 599 590 641 676 3112 C 631 615 593 713 653 3205 合计(iB)1842 1776 1775 1983 1933 9309 表12-17 不同麻醉诱导、不同时相患者收缩压合计值()(5)ijTn
14、 分解ASS,BSS,ABSS 分组计算不同麻醉诱导、不同时相患者的收缩压的合计值(ijT),见表 12-17。表12-20 不同麻醉诱导、不同时相患者的收缩压(mmHg)麻 醉 诱 导 时 相 诱导方法 T0 T1 T2 T3 T4 A X 121.00 112.40 118.40 125.80 120.80 S 3.54 5.13 5.64 4.71 3.70 B X 121.20 119.80 118.00 128.20 135.20 S 4.32 5.97 5.43 5.22 4.38 C X 126.20 123.00 118.60 142.60 130.60 S 3.63 3.39
15、 1.95 4.83 3.71 重复测量数据统计分析常见的误用情况 1重复进行各时间点的 t 检验。如对表12-20的统计结果,每个时相做3次t检验比较A、B、C三种诱导方法的差别,5个时相要做15次t检验,必然增大假阳性错误。2.忽略个体曲线变化特征。重复测量数据的个体差异是每个观察对象的m次测量结果(即横向差异),不能用纵向均数比较差别。例例22.2 某药物有新、旧两种剂型。为了比较这两种剂型的代谢情况,对16例某病患者服药后0、4、8、12小时的血药浓度()作了测量,具体结果见表22.2,问该药新旧两种剂型的血药浓度-时间曲线的差别是否具有统计学意义。连接各时间点测量值均数的方法将掩盖类
16、连接各时间点测量值均数的方法将掩盖类似图似图22.2所表现的个体曲线位置和形状的所表现的个体曲线位置和形状的特点;特点;各均数所对应的点连起来形成的曲线形状各均数所对应的点连起来形成的曲线形状可能与诸个体的曲线形状毫不相干;可能与诸个体的曲线形状毫不相干;无法体现不同时间点的数据来自同一个体,无法体现不同时间点的数据来自同一个体,而每一个体的多次重复测量值间具有相关而每一个体的多次重复测量值间具有相关性性 信度是指在相同条件下,对同一客观事物重复测量若干次,测量结果的相互符合程度,说明数据的可靠性。3.差值比较缺乏效度 因为前后测量转换为差值后,信度降低,且差值一般不符合正态性和方差齐性的条件。效度效度是指测量指标或观测结果在多大程度上反映了是指测量指标或观测结果在多大程度上反映了事物的客观真实性,说明数据的准确性。事物的客观真实性,说明数据的准确性。17.2 交叉设计模型 药物试验或其它的交叉设计通常要考虑三种效应:药物试验或其它的交叉设计通常要考虑三种效应:阶段效应,处理效应和顺序效应。阶段效应,处理效应和顺序效应。例例17-2 为了比较两种体育运动(运动为了比较两种体育运动(运动