重复测量资料的统计分析方法Word文档格式.docx
- 文档编号:19825336
- 上传时间:2023-01-10
- 格式:DOCX
- 页数:12
- 大小:100.31KB
重复测量资料的统计分析方法Word文档格式.docx
《重复测量资料的统计分析方法Word文档格式.docx》由会员分享,可在线阅读,更多相关《重复测量资料的统计分析方法Word文档格式.docx(12页珍藏版)》请在冰豆网上搜索。
因而本节将借助统计软件Stata,介绍应用混合模型(MixedModel)对重复测量资料进行统计分析。
设观察对象体重的总体均数为μ0,服药后体重总体均数为μ1,即服药前后的体重改变量的总体均数为β=μ1-μ0。
若β=0说明服药前后的体重平均变化为0,即无疗效;
若β<
0,说明服药后的人群平均体重低于服药前的平均体重,即该药物减肥是有效的;
若β>
0,说明服药后的平均体重高于服药前的平均体重,即该药对减肥有不利的作用。
针对本例服药前后的体重总体均数的变化关系,引入自变量t,建立下列服药前后的体重总体均数表达式(即混合模型的确定性部分表达式)。
(12-1)
t=0时,μ为服药前的体重总体均数μ0;
t=1时,μ为服药后的体重总体均数μ1。
应用混合模型可以对本例资料进行统计分析,其中β和μ0的参数估计一般采用限制的最大似然法,然而计算相当复杂,故我们将借助Stata软件对上述资料用混合模型进行统计分析,相应的Stata软件的数据格式如下。
t
y
no
其中y为体重测量值,t为服药时间的自变量,no为观察对象的编号,相应的Stata操作命令如下:
Random-effectsGLSregressionNumberofobs=10
Groupvariable(i):
noNumberofgroups=5
R-sq:
within=0.6533Obspergroup:
min=2
between=.avg=2.0
overall=0.0612max=2
Randomeffectsu_i~GaussianWaldchi2
(1)=7.54
corr(u_i,X)=0(assumed)Prob>
chi2=0.0060
------------------------------------------------------------------------------
y|Coef.Std.Err.zP>
|z|[95%Conf.Interval]
-------------+----------------------------------------------------------------
t|-1.4.509902-2.750.006-2.399389-.4006105
_cons|50.41.37113136.760.00047.7126353.08737
sigma_u|2.9580399
sigma_e|.80622577
rho|.93085106(fractionofvarianceduetou_i)
β估计值为-1.4,μ0估计值为50.4,而μ1的估计值=50.4-1.4=49。
H0:
β=0即无减肥疗效
H1:
β≠0即服药前后的人群平均体重不同
α=0.05
相应的P值=0.006,因此服药前后平均体重的差异有统计学意义,故可以认为该药物有减肥疗效。
例2为了考察某药物在疗程为6个月中的持续减肥作用,现考察5个服用该药的女性肥胖者并且身高为162cm的,这5名女性肥胖者在服用该药前、服药3个月和服药6个月的体重测量值(kg)如下:
服药前
3个月
6个月
42
53
47
这是一组观察对象的多个测量时间点的重复观察测量资料,因此同一对象的不同观察时间点的观察资料是相关的。
(也可以视为配伍区组设计的观察资料,用随机区组设计的方差分析或Friedman秩检验的统计方法检验该药物的减肥作用),因此可用混合模型进行统计分析。
设观察对象在服药前的体重总体均数为μ0、服药3个月时的体重总体均数μ0+β1,服药6个月时的体重总体均数为μ0+β2,即:
β1为服药3个月时的体重平均改变量,β2为服药6个月时的体重平均改变量。
针对本例服药前后的体重总体均数的变化关系,引入自变量t1和t2,建立下列服药前后的体重总体均数表达式
(12-2)
若t1=t2=0时,μ为服药前的体重总体均数μ0;
t1=1,t2=0时,μ为服药3个月时的体重总体均数μ0+β1。
若β1<
0,说明服药3个月时的服药人群平均体重低于服药前的平均体重,即该减肥药有效,反之无疗效;
t1=0,t2=1时,μ为服药6个月时的体重总体均数μ0+β2,而β2<
0和β2>
0同样反映该减肥药有效或无效。
若β2<
1,说明服药6个月时的服药人群平均体重低于3个月时的服药人群平均体重到期间。
我们同样借助Stata软件对上述资料用混合模型进行统计分析,相应的Stata软件的数据格式如下。
t1
t2
其中t1和t2为服药时间的自变量,其他与例12-1相同,相应的Stata操作命令如下:
Stata命令为xtregyt1t2,i(no)
相应输出结果如下:
Random-effectsGLSregressionNumberofobs=15
within=0.9551Obspergroup:
min=3
between=.avg=3.0
overall=0.4602max=3
Randomeffectsu_i~GaussianWaldchi2
(2)=170.18
chi2=0.0000
t1|-1.2.3829708-3.130.002-1.950609-.4493909
t2|-4.8.3829708-12.530.000-5.550609-4.049391
_cons|51.61.10453646.720.00049.4351553.76485
sigma_u|2.394438
sigma_e|.60553007
rho|.93989071(fractionofvarianceduetou_i)
3个月时的体重与6个月时的体重比较的Stata命令和输出结果如下:
testt1=t2(H0:
β2=β1)
(1)t1-t2=0.0
chi2
(1)=88.36
Prob>
β1估计值为-1.2(kg),β2估计值为-4.8(kg),服药前体重总体均数μ0的估计值为51.6(kg);
服药3个月时的体重总体均数μ0+β1的估计值为51.6-1.2=50.4(kg);
服药6个月时的体重总体均数μ0+β2的估计值为51.6-4.8=46.8(kg)。
β1=0即服药3个月时减肥无效
β≠0即服药3个月时与服药前的人群平均体重不同
相应的P值=0.002,因此差异有统计学意义,故可以认为该药物在服药3个月时有减肥疗效。
β2=0即服药6个月时减肥无效
β2≠0即服药6个月时与服药前的人群平均体重不同
相应的P值<
0.001,因此差异有统计学意义,故可以认为该药物在服药6个月时有减肥疗效。
β2=β1即从服药3个月至6个月时,没有继续减肥
β2≠β1即服药6个月时与服药3个月的人群平均体重不同
0.001,因此差异有统计学意义,故可以认为服药3个月至6个月期间,继续有减肥疗效。
多个样本多个时间点重复观察资料
例12-3为了比较A药和B药在疗程为6个月中的持续减肥的疗效,现有10个身高为160cm的女性肥胖者志愿参加这项研究。
随机分成2组,每组各5人。
分别考察这2组肥胖者在服药前、3个月和服药6个月的体重变化。
这2组肥胖者在服用该药前、服药3个月和的体重测量值(kg)如下:
组别和肥胖者编号
A药组1号
A药组2号
A药组3号
41
A药组4号
44
A药组5号
40
B药组1号
54
B药组2号
B药组3号
B药组4号
B药组5号
这是两组观察对象的多个测量时间点的重复观察测量资料,同样对于同一对象的不同观察时间点的观察资料是相关的,但由于需要比较两个药的减肥疗效,所以两因素方差分析,随机区组设计的方差分析或Friedman秩检验的统计方法都不适用于本例的数据统计分析,但仍可用混合模型对本例资料进行统计分析。
由于这是两组观察对象的多个测量时间点的重复观察测量资料,所以仍可以借用上例的总体均数表达式(12-2)分别描述每一组体重变化规律,因此可以得到下列总体均数表达式:
A组
(12-3)
B组
(12-4)
由于应用混合模型进行统计分析需要建立两组统一的总体均数表达式,因此引入统一参数β3=μ10-μ0,β4=β11-β1,β5=β12-β2,代入(12-4)式,得到B组总体均数表达式
(12-5)
引入分组变量g=0表示A组,g=1表示B组,因此两组的总体均数表达式均可表示为
(12-6)
用g=0,1;
t1=0,1和t2=0,1代入(12-6)式,得到相应两组各个时间点的总体均数:
组别
服药前(t1=0,t2=0)
服药3个月(t1=1,t2=0)
服药6个月(t1=0,t2=1)
A组(g=0)
μ0
μ0+β1
μ0+β2
B组(g=1)
μ0+β3
μ0+β1+β3+β4
μ0+β2+β3+β5
两组差异为
β3
β3+β4
β3+β5
若β4和β5不全为0,则称两种药物与服药时间对疗效有交互作用。
两组在3个时间点的总体均数差异分别为β3,β3+β4和β3+β5,因此只需检验H0:
β3=0、H0:
β3+β4=0和H0:
β3+β5=0就可以推断两组总体均数差异。
反之若β4和β5全为0,则称两种药物与服药时间对疗效无交互作用,并且两组各个时间点的总体均数差异均为β3,因此只需检验H0:
β3=0就可以推断两组的总体均数差异。
g
续左侧底部数据
6
7
8
9
10
Stata操作命令如下:
gengt1=g*t1产生交互作用项变量g⨯t1
gengt2=g*t2产生交互作用项变量g⨯t2
xtregyt1t2ggt1gt2,i(no)
Random-effectsGLSregressionNumberofobs=30
noNumberofgroups=10
within=0.8288Obspergroup:
between=0.0973avg=3.0
overall=0.5927max=3
Randomeffectsu_i~GaussianWaldchi2(5)=78.32
t1|-1.81.053565-1.710.088-3.86495.2649502
t2|-81.053565-7.590.000-10.06495-5.93505
g|-.41.612452-0.250.804-3.5603472.760347
gt1|.81.4899660.540.591-2.1202813.720281
gt2|4.21.4899662.820.0051.2797197.120281
_cons|50.61.14017544.380.00048.365352.8347
sigma_u|1.9300259
sigma_e|1.6658331
rho|.57307692(fractionofvarianceduetou_i)
由此得到(12-6)式中的μ0估计值为50.6,β1的估计值为-1.8,β2的估计值为-8,β3的估计值为-0.4,β4的估计值为0.8和4.2。
两组各个时间的总体均数估计如下
A组(g=0)
B组(g=1)
B组-A组
两组差异检验
总体
均数
总体均数
估计值
P值
(t1=0,t2=0)
50.6
50.2
-0.4
0.804
服药3个月时
(t1=1,t2=0)
48.8
49.2
0.4
服药6个月时
(t1=0,t2=1)
42.6
46.4
3.8
0.018
注:
表中均数估计值是参数估计值和总体均数参数表达式计算所得。
如:
服药3个月时A组的总体均数估计值=50.6-1.8=48.8。
3个时间点的两组平均体重比较的Stata统计检验命令和输出结果如下
设α=0.05
服药前两组平均体重比较就是检验H0:
β3=0,相应的P值=0.804>
α,差别无统计学意义,故没有充足证据推断两组在服药前的体重总体均数不等。
服药3个月时的两组平均体重比较的Stata命令和输出结果如下:
testg+gt1=0(H0:
β3+β4=0即服药3个月时的两组体重总体均数相等)
(1)g+gt1=0.0
chi2
(1)=0.06
chi2=0.8041
相应的P值=0.8041>
α,差异无统计学意义,故无证据显示两组总体均数不等。
服药6个月时的两组平均体重比较的Stata命令和输出结果如下:
testg+gt2=0(H0:
β3+β5=0)
(1)g+gt2=0.0
chi2
(1)=5.55
chi2=0.0184
相应的P值=0.0184<
α,差异有统计学意义,故可以认为服A药6个月时的人群平均体重低于服B药6个月的人群平均体重。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 重复 测量 资料 统计分析 方法