系统评价中的常用统计分析方法(2013.1.15).pptVIP

结果 Meta分析举例----定性资料 a：干预组戒烟人数；r1:干预组总人数； c：对照组戒烟人数；r2:对照组总人数； State软件操作步骤如下 1、命令方式： ○gen b=r1-a; /*计算出干预组未发生人数*/ ○gen d=r2-c; /*计算出对照组未发生人数*/ metan a b c d, label(namevar=name, yearvar=year) fixed rr xlabel(0.1,0.2,0.5,0.75,1,1.25,2,5,10) Stata软件操作步骤如下 2、菜单方式： 2.1 分别计算干预组和对照组未发生人数。操作如下：Data→Create or Change data →Create new variable,出现如下对话框，如图填入即可。同理可计算d。 Stata软件操作步骤如下（续） 2、菜单方式 2.2 meta分析。操作如下：User菜单→Meta-Analysis →of Binary and Continuous（metan）,弹出如下对话框： Stata软件操作步骤如下（续） 菜单方式 2.3 上图graph opts选项卡设置如下： 设置坐标轴刻度 结果 结果 L’Abbe图的制作过程—异质性评价工具 User菜单→Meta-Analysis→L’abbe，出现如下对话框并设置如下。 L’Abbe图的作用 判断研究结果有无异质性。 L’Abbe图是以研究中的干预组事件发生率相对于对照组事件发生率作图，若研究结果同质，则所有点呈直线分布，偏离线过远，则表明研究结果不同质。 图中圆圈的大小表示效应量的大小。 Galbraith星状图—异质性评价工具 User菜单→Meta-Analysis→galbraith plot for…，出现如下对话框并设置如下。 Galbraith星状图作用及说明 判断研究结果有无异质性。 原理：Galbraith星状图是以标准化估计值（如logOR/lnRR/SE相对于其标准误的倒数做图，若散点图斜率较为接近则说明研究间同质）。 如上图，有一个研究落在了95%CI线以外，其它研究的散点斜率差别不是很大，提示研究间有一定同质性，但可接受。 异质性的统计量—H 异质性评价统计量——H的计算 是通过对统计量Q进行自由度（文献数）的校正，H统计量： 异质性统计量—H的意义 漏斗图的制作过程—发表性偏倚评价工具 User菜单→Meta-Analysis→funnel Graph,ver…，出现如下对话框并设置如下。 metafunnel _ES _selogES, var xtitle(risk ration) ytitle(1/SE(Effect size)) 发表偏倚的评价—Egger、Begg检验及Egger图、Begg图 步骤：User菜单→Meta-Analysis→publication Bias，出现如下对话框并设置如下。 发表偏倚检验结果 本质为旋转的漏斗图 剪补法（metatrim）—发表偏倚识别工具 User菜单→Meta-Analysis→trim and Fill Analysis，出现如下对话框并设置如下。 敏感性分析—metainf 敏感性分析，即通过比较纳入或排除某些研究进行的meta分析结果，探讨被去除的研究对合并效应的影响程度，从而判断结果的稳定性。且通过比较可了解异质性的来源。 操作过程： User菜单→Meta-Analysis→influence analysis，metan-based，出现如下对话框并设置如下。 NNT—metannt 。 Meta分析举例----定量资料 例：比较氟含量对身体发育的影响，以女童掌II皮质厚度为指标，收集11个研究的数据如下，试进行meta分析。 State软件操作步骤如下 操作如下：User菜单→Meta-Analysis →of Binary and Continuous（metan）,弹出如下对话框： 比值比（odds ratio,OR） 干预组的比值比为=121/160=0.76 对照组的比值比为=160/121=1.32 group Common cold 合计 Yes No 治疗组 121 160 281 对照组 160 121 281 合计 281 281 562 比值比=121/160=0.76=0.58 160/121 1.32 OR=1，表明干预组与对照组无差异 效应指标的含义 率，即危险度，用来说明某事发生可能性大小的指标； RR: 相对危险度，队列研究中常用指标。说明暴露组的发生率是非暴露组的多少倍。RR越大，表明暴露与结局关联的强度越大。 相对危险度与关联的强度见下表 OR OR即比值比，也