重复测量数据的分析方法.pptVIP

重复测量数据的分析方法 赵耐青 复旦大学卫生统计教研室 重复测量资料简介 为了考察某药物降低血药的作用，现考察5个原发性高血压呀患者，疗程为2个周，这5名原发性高血压患者在服用该药前后的收缩压值(mmHg)如下 单样本重复测量资料的统计分析方法 由于对于每个人不同时刻的血压都在一定范围里随机变化的，所以一般而言，虽然平时临床上测定一个时点的血压，但实际上该时点的血压是作为被测对象在这一时期的平均血压的估计值，故对于一次血压测量而言，可以忽略血压的个体变异，但对于重复测量而言，个体变异是不能被忽略的。 单样本重复测量资料的统计分析方法 设治疗前研究对象的人群收缩压平均值为?0，第i个对象在治疗前的收缩压的平均值为 ，其中 为个体间的随机变异，假定 ，因此 设第i个对象在治疗前实际收缩压测量值为 ，由于个体内收缩压存在随机变异，所以治疗前收缩压测量值可以表示为 单样本重复测量资料的统计分析方法 设治疗后研究对象的人群收缩压平均值为?1，第i个对象在治疗后的收缩压的平均值为 ，其中 为个体间的随机变异，因此 设第i个对象在治疗后实际收缩压测量值为 ，由于个体内收缩压存在随机变异 ，所以治疗前收缩压测量值可以表示为 单样本重复测量资料的统计分析方法 治疗前的收缩压表达式为 治疗后的收缩压表达式为 可以验证：治疗前后的资料是不独立的 所以治疗前后的总体均数改变量为 因此 用t=0和t=1表示治疗前后，则 单样本重复测量资料的统计分析方法 称 为mixed 模型 即混合模型的确定性部分表达式 t=0时，?为服药前收缩压总体均数?0； t=1时，?为服药后收缩压总体均数?1 。 单样本重复测量资料的统计分析方法 由于服药前后的收缩压改变量的总体均数为?＝?1－?0 若?＝0说明服药前后收缩压平均变化为0，即无疗效； 若?0，说明服药后的人群平均收缩压低于服药前，即该药物降低收缩压是有效的； 若?0，说明服药后的人群平均收缩压高于服药前的平均收缩压，即该药对收缩压有不利的作用。 单样本重复测量资料的统计分析方法 不难验证：治疗前后的资料是不独立的。 可以用配对t检验或配对符合秩检验： 其中 为治疗前后的个体内随机变异，故研究问题为?是否为0，如果d服从正态分布，也可以用配对t检验，反之用秩检验 但对于多组重复测量资料要用Mixed模型分析 用Mixed模型进行分析 Mixed模型采用分离两种随机变异后用限制的最大似然估计(REML)拟合模型，解决了资料不独立的情况下对参数进行估计和检验。 Stata数据结构 x=0表示治疗前 x=1表示治疗后 Stata命令 xtreg y x ,i(id) 用Mixed模型进行分析 用Mixed模型进行分析 ?的估计值为-30,P0.001，因此治疗后的收缩压低于治疗前，差异有统计学意义。 ?的95%可信区间为 (-33.87062, -26.12938) 即：有95%可信度可以推断治疗后的人群平均收缩压至少下降了26mmHg 单样本多个时点的重复测量资料 例：为了考察某药物治疗在疗程为6个月中的持续减肥作用，现考察5个服用该药的女性肥胖者并且身高为162cm的，这5名女性肥胖者在服用该药前、服药3个月和服药6个月的体重测量值(kg)如下： 用Mixed模型进行分析 Mixed模型确定部分： 3个月平均体重与基线平均体重之差为?1， 6个月平均体重与基线平均体重之差为?2， 6个月平均体重与3个月平均体重之差为?2-?1 用Mixed模型进行分析 Stata数据格式 Stata命令 xtreg y t1 t2,i(id) 用Mixed模型进行分析 用Mixed模型进行分析 3个月与6个月比较 命令：test t2=t3 ( 1) t2 - t3 = 0 chi2( 1) = 88.36 Prob chi2 = 0.0000 P0.0001，可以推断治疗6个月时平均体重比治疗3个月时平均体重下降，差异有统计学意义。 多组重复测量资料 例 为了比较A药和B药在疗程为6个月中的持续减肥的疗效，现有10个身高为160cm的女性肥胖者志愿参加这项研究。随机分成2组，每组各5人。分别考察这2组肥胖者在服药前、3个月和服药6个月的体重变化。 多组重复测量资料 用Mixed模型进行分析 Mixed模型确定部分，g=0为A药组，g=1为B药组 各个时点的两组总体均数 用Mixed模型进行分析 数据格式