Stata在meta分析中的应用.docx
精品word 可编辑资料 - - - - - - - - - - - - -1 定量资料两组比较的meta 分析2 定性资料两组比较的meta 分析实例：分类资料的 meta 分析为了探讨用 Aspirin预防心肌梗塞 (myocardialinfarction，MI) 后死亡的发生；美国在1976 1988 年问进行了 7 个关于 Aspirin预防 MI 后死亡的讨论，具体结果见表1，其中 6 项讨论的结果说明Aspirin组与劝慰剂组的MI 后死亡率的差别无统计学意义；只有1 项结果说明 Aspirin预防 MI 后死亡有效并且差别有统计学意义；现 依据表 1 供应的结果进行meta 分析表 1Aspirin预防心肌梗塞死亡的临床试验结果Aspirin 组劝慰剂组讨论发表年份总例数死亡例数总例数死亡例数MRC-119746154962467CDP19767584477164MRC-21979832102850126GASP19793173230938PARIS19808108540652AMIS198022372462257219ISIS-2操作步骤198885871570860017201 把数据输入 stata软件1 / 10第 1 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -2 变量的说明Study纳入的讨论Year年份Death1Aspirin组的死亡人数Live1Aspirin组的存活人数Death2劝慰剂组的死亡人数Live2劝慰剂组的死亡人数3 进行 meta 分析metandeath1live1dead2live2,orlabel(namevar=study, yearvar=year)结果：2 / 10第 2 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -以上结果分成两部分(1) meta分析的合并统计量合并 OR值 OR=0.897,95% 的可信区间 (0.841,0.957)(2) 给出异质性检验的结果只要异质性检验的P值不小于 0.10( 或者 I-squared小于 50%)就可以认为不存在异质性，可以用应固定效应模型(stata默认的情形 ) ；假如质性检验的P 值小于等于 0.10( 或者 I-squared大于 50%)，就不同的组间存在异质性，应当应用随机效应模型随机效应模型的命令如下：metandeath1 live1 dead2 live2, or label(namevar=study, yearvar=year) random在运行 meta 分析命令的同时stata输出森林图，如下：3 / 10第 3 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -由输出的合并结果和漏斗图可以得出，合并的OR值为 0.90,95%可信区间为 (0.84 ，0.96)4 发表偏倚的检验，命令如下：(1) gen logor=log(_ES)(2) gen selogor=_selogES(3) metabias logor selogor,graph(begg)输出结果如下：4 / 10第 4 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -发表偏倚主要看begg 检验的结果，由上图可以看到发表偏倚假设检验的 z 值为 1.20 ， p 值为 0.230>0.05 ，可以认为没有发表偏倚；进行发表偏倚检验时stata同时给出的漏斗图如下5 结论5 / 10第 5 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -合并 OR小于 1，95%的可信区间不包括1，所以示 Aspirin有预防心肌梗塞后死亡的作用数值变量资料的meta 分析Gotzsche 收集了有关短程小剂量强的松VS劝慰剂或非甾体抗炎药治疗类风湿性关节炎的7 个临床随机对比试验 (RCTs)，观看类风湿性关节炎患者的关节压痛指数(rechieS index)，数据如表 2； 表 27个讨论类风湿性关节炎患者关节压痛指数试验名称发表年份试验组例数均数标准差对比组例数均数标准差Jasni1968916.28.7938.112.8Dick19702417.68.02440.713.0Lee19732130.516.52141.419.8Berry19741213.011.01223.711.1Lee19741814.612.41826.415.1Stenber1992216.31.72111.12.1Geital19952010.84.72016.37.71 在 stata数据编辑器中建立数据文件2 变量名称study纳入的讨论year讨论的年份n1试验组的样本含量6 / 10第 6 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -mean1试验组的均数sd1试验组的标准差n2对比组的样本含量mean2对比组的均数sd2对比组的标准差3 进行 meta 分析metann1mean1sd1n2mean2sd2,label(namevar=study, yearvar=year)输出的结果判定异质性的方法同上面定性资料的meta，有上图可以看到，不同 的组间存在异质性，所以应当应用随机效应模型应用随机效用模型以后输出的结果如下7 / 10第 7 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -stata同时输出森林图如下4 进行发表偏倚的检验Stata命令： metabias_ES _seES,graph(begg)检验结果 ( 结果的说明同分类资料的meta 分析)8 / 10第 8 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -漏斗图如下：由以上分析可以得到没有发表偏倚的存在5 结论9 / 10第 9 页，共 10 页 - - - - - - - - - -精品word 可编辑资料 - - - - - - - - - - - - -合并的效用值为 -1.378 ，95%CI为(-1.951,-0.805)，不包括 0，可以得到短程小剂量强的松有削减类风湿性关节炎患者的关节压痛指数的作用10 / 10第 10 页，共 10 页 - - - - - - - - - -
