返回首页  
学科简介 优势特色 教学力量 教学条件 课程设置 科研动态 研究方向 科研资源 学科成果 学习工具
科研资源
怎么做到早期发现肝癌呢
胃炎等慢性炎症是会变成淋巴瘤吗?
肿瘤筛查有哪些指导原则
厨房中怎样防癌?

Content on this page requires a newer version of Adobe Flash Player.

Get Adobe Flash player

科研资源
META分析举例

  为了研究某药物对于治疗急性心肌梗塞的疗效,通过查阅文献,收集到2篇关于该药物的临床随机对照试验的研究并符合文献评价的标准。对照组为常规治疗方案,试验组为该药物+常规治疗的方案。主要结果如下表(括号中的字母表示计算公式中的记号)

试验组

对照组

合计

研究1

发生病危或死亡

182

264

446

(i=1)

未发生病危或死亡

1068

986

2054

合计

1250

1250

2500

研究2

发生病危或死亡

348

204

552

(i=2)

未发生病危或死亡

1273

610

1883

合计

1621

1250

2434

例:为了研究瘦素(Leptin,单位为ng/mol)与肥胖的关系,人的肥胖通常用体重指数BMI(体重除以身高的平方,单位为公斤/2)进行评估和刻划。由于Leptin测试费用较高和麻烦,而LeptinBMI的线性相关程度较高,所以许多研究者希望用BMI近似估计Leptin。通常采用以Leptin为因变量和BMI为自变量作线性回归。为了提高估计精度,特收集了4篇有关LeptinBMI线性回归的文献,希望能通过作Meta分析,提供回归系数b的估计精度。现以以Leptin为因变量和BMI为自变量的回归方程中的回归系数b为观察指标,为了避免年龄的混杂作用,本例限定研究对象均为50岁至65岁的男性正常人。现收集的4个文献的结果如下:

所收集研究结果的编号(no)

回归系数

(b)

回归系数的标准误

SE(b)

1

0.693

0.0848

2

0.556

0.103

3

0.723

0.116

4

0.755

0.085

标化效应指标的Meta分析

有些研究的研究目的是相同的。但由于它们的观察指标往往采用不同的检测方法而使各个研究的结果无法进行直接的比较。因此这些效应指标通过标化后还是有可能进行比较以及进行Meta分析,因此本节介绍标化效应(Standard Effect Size)Meta统计分析方法。

例:为了研究健康教育对2型糖尿病患者饮食等自我控制的作用。现收集了5个关于健康教育对2型糖尿病饮食等自我控制作用的研究。效应指标为糖化血红蛋白HbA1c。资料具体如下:

所收集研究

健康教育组

对照组

P

结果的编号

(no)

样本量

(n1)

均数

(mean1)

标准差

(s1)

样本量

(n2)

均数

(mean2)

标准差

(S2)

1

13

7.699

0.703

10

8.012

0.795

0.3470

2

12

8.29

2.002

19

8.896

2.011

0.4181

3

22

6.411

0.9111

26

7.391

1.112

0.0019

4

20

4.314

3.699

22

6.001

4.097

0.1706

5

27

6.51

0.521

22

6.577

0.801

0.7255

5个研究的结果不一致,而且样本量均较小,所以检验效能1-b都相对较小。因此还是有必要把这5个研究的结果进行Meta分析。

由于各个医院的糖化血红蛋白的检测方法不全相同,因此各个医院的糖化血红蛋白的参考值有一些的差别。因此直接进行Meta分析有可能与背景不符。所以应该进行标化后再进行Meta分析。(括号中的字母表示计算公式中的记号)

合并的标准差公式

" v:shapes="_x0000_i1025">

一般的标化效应公式

" v:shapes="_x0000_i1026">
" v:shapes="_x0000_i1027">
,其中S合并为合并后的标准差,meand为改变量的均数,Sd为改变量的标准差。

标化效应rd的标准误公式

" v:shapes="_x0000_i1028">

由于合并标准差公式在卫生统计学教课书中的t检验章节中均作介绍而且在表中已给出,故不再详细叙述。标化效应指标和标化效应指标的标准误计算都很简单,故不再一一详细进行手算了。表中的合并标准差s、标化的效应指标rd和标化效应指标rd的标准误se均由下列STATA软件操作进行计算的。

STATA软件实现的步骤如下:

1、数据输入结构如下:

健康教育组

对照组

样本量

均数

标准差

样本量

均数

标准差

n1

mean1

s1

n2

mean2

s2

1个研究

13

7.699

0.703

10

8.012

0.795

2个研究

12

8.29

2.002

19

8.896

2.011

3个研究

22

6.411

0.9111

26

7.391

1.112

4个研究

20

4.314

3.699

22

6.001

4.097

5个研究

27

6.51

0.521

22

6.577

0.801

2STATA运算操作命令如下

gen  s=sqrt(((n1-1)*s1^2+(n2-1)*s2^2)/(n1+n2-2))

计算合并标准差S

gen rd=(mean1-mean2)/s

计算标化后的效应指标rd

gen se=sqrt((8+rd^2)/(2*(n1+n2)))

计算标化后的效应指标rd的标准误se

meta  rd  se

进行Meta分析计算

得到下列结果:

Fixed and random effects pooled estimates,

lower and upper 95% confidence limits, and

asymptotic z-test for null hypothesis that true effect=0

Fixed effects estimation

   Est   Lower   Upper  z_va~e  p_va~e

-0.443  -0.731  -0.156  -3.025   0.002 

Test for heterogeneity: Q=  4.418 on 4 degrees of freedom (p= 0.352)

标化要注意的问题:

1.       对于效应指标为线性回归系数,则标化只需对线性回归系数b乘与

" v:shapes="_x0000_i1029">
,其中SX是自变量x的样本标准差,SY是因变量Y的样本标准差。

2.       对于效应指标为Logistic回归模型和Cox模型的回归系数,而相应的协变量是连续性变量,则标化只需对回归系数b乘与SX

3.       对于经典的算法所计算的ORRR等相对指标无需标化,因为这些效应指标与量纲无关。

 
学科简介 | 优势特色 | 教学力量 | 教学条件 | 课程设置 | 科研动态 | 研究方向 | 科研资源 | 学科成果 | 学习工具
陆军军医大学新桥医院肿瘤科 联系电话: 774050
COPYRIGHT @ All Right Reserved XINQIAO hospital.