欢迎来到天天文库
浏览记录
ID:47580145
大小:328.00 KB
页数:18页
时间:2019-09-21
《重复测量资料的统计剖析方法》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、重复测量资料的统计分析方法在临床医学研究中,一些干预研究和纵向研究都经常会涉及到同一研究对象的多次观察,而同一个对象的多次观察的记录资料称为重复测量的资料。由于同一对象不同时间点的观察往往存在相关的问题,也就是存在不独立性的问题,而大多数的医学统计方法都要求资料是独立,所以这些资料的统计分析需要比较特殊的统计方法进行分析。本节将先举例介绍常见的重复测量资料,并介绍相应的重复测量资料的统计分析方法。一、单个样本的重复测量资料例1为了考察某药物减肥的作用,现考察5个身高为160cm、服用该药的女性肥胖
2、者,疗程为3个月,这5名女性肥胖者在服用该•药前后的体重测量值(kg)如下:吧胖者编号12345服药前体重5052495546Y01服药后体重4851495245Yn这是一组观察对象的资料,每个观察对象有两个时间点的测量资料,因此这是最简单的重复观察测量资料(也可以认为配对设计的资料)。由于各个观察对象在服药前的体重不全相同,所以其体重含有服药前的体重个体变异成分,而在服药后,各个观察对象的体重下降幅度也不全相同,故存在体重下降幅度的个体变异成分,因此观察对象重复测量资料的统计分析方法在临床医学研
3、究中,一些干预研究和纵向研究都经常会涉及到同一研究对象的多次观察,而同一个对象的多次观察的记录资料称为重复测量的资料。由于同一对象不同时间点的观察往往存在相关的问题,也就是存在不独立性的问题,而大多数的医学统计方法都要求资料是独立,所以这些资料的统计分析需要比较特殊的统计方法进行分析。本节将先举例介绍常见的重复测量资料,并介绍相应的重复测量资料的统计分析方法。一、单个样本的重复测量资料例1为了考察某药物减肥的作用,现考察5个身高为160cm、服用该药的女性肥胖者,疗程为3个月,这5名女性肥胖者在服
4、用该•药前后的体重测量值(kg)如下:吧胖者编号12345服药前体重5052495546Y01服药后体重4851495245Yn这是一组观察对象的资料,每个观察对象有两个时间点的测量资料,因此这是最简单的重复观察测量资料(也可以认为配对设计的资料)。由于各个观察对象在服药前的体重不全相同,所以其体重含有服药前的体重个体变异成分,而在服药后,各个观察对象的体重下降幅度也不全相同,故存在体重下降幅度的个体变异成分,因此观察对象在服药后的体重中不仅含有体重下降幅度的个体变异成分,而且还含有服药前的体重个
5、体变异成分,故服药前后的体重资料不独立。对于这种不独立资料的统计分析一般采用变异成分的分解或消除某一个体变界成分的方法进行统计处理的。如配对t检验和符号秩检验就是采用服药前后资料相减作为统计分析数据,因而消除了服药前体重的个体变异,使进入统计分析的资料仅含有体重下降幅度的个体变异,但这种消除某种不独立的变异成分的统计方法无法对比较复杂的重复测量资料进行统计分析。因而本节将借助统计软件Stata,介绍应用混合模型(MixedModel)对重复测量资料进行统计分析。设观察对象体重的总体均数为0,服药后
6、体重总体均数为1,即服药前后的体重改变量的总体均数为=1—0。若=0说明服药前后的体重平均变化为0,即无疗效;若<0,说明服药后的人群平均体重低于服药前的平均体重,即该药物减肥是有效的;若>0,说明服药后的平均体重高于服药前的平均体重,即该药对减肥有不利的作用。针对本例服药前后的体重总体均数的变化关系,引入自变量t,建立下列服药前后的体重总体均数表达式(即混合模型的确定性部分表达式)。“(12-1)t=0时,为服药前的体重总体均数o;t=l时,为服药后的体重总体均数lo应用混合模型可以对本例资料进
7、行统计分析,其屮和0的参数估计一般采用限制的最大似然法,然而计算相当复杂,故我们将借助Stata软件对上述资料用混合模型进行统计分析,相应的Stata软件的数据格式如下。tyno0501052204930554046514811512149315241455其中y为体重测量值,t为服药时间的口变量,no为观察对象的编号,和应的Stata操作命令如下:Random-effectsGLSregressionNumberofobs二10Groupvariable(i):noNumberofgroups=
8、5R-sq:within二0.6533Obspergroup:min二2between—■avg二2.0overall二0.0612max二2Randomeffectsu_i~GaussianWaldchi2(1)二7.54corr(u_i,X)=0(assumed)Prob>chi2二0.0060yCoef・Std・Err.zP>
9、z[95%Conf.Interval]+t
10、-1.4.509902-2.750.006-2.399389-4006105_cons
11、50.41.371
此文档下载收益归作者所有