2022年Stata在meta分析中的應(yīng)用

1、1 定量資料兩組比較的 meta分析2 定性資料兩組比較的 meta 分析實(shí)例：分類資料的 meta 分析為了探討用 Aspirin預(yù)防心肌梗塞 myocardial infarction,MI后死亡的發(fā)生；美國(guó)在19761988 年問(wèn)進(jìn)行了 7 個(gè)關(guān)于 Aspirin預(yù)防 MI 后死亡的討論,具體結(jié)果見(jiàn)表1,其中 6 項(xiàng)討論的結(jié)果說(shuō)明1Aspirin組與勸慰劑組的MI 后死亡率的差別無(wú)統(tǒng)計(jì)學(xué)意義；只有項(xiàng)結(jié)果說(shuō)明 Aspirin預(yù)防 MI 后死亡有效并且差別有統(tǒng)計(jì)學(xué)意義；現(xiàn)依據(jù)表 1 供應(yīng)的結(jié)果進(jìn)行meta 分析表 1 Aspirin預(yù)防心肌梗塞死亡的臨床試驗(yàn)結(jié)果討論發(fā)表年份Aspirin 組

2、勸慰劑組總例數(shù)死亡例數(shù)總例數(shù)死亡例數(shù)MRC-1 1974 615 49 624 67 CDP 1976 758 44 771 64 MRC-2 1979 832 102 850 126 GASP 1979 317 32 309 38 PARIS 1980 810 85 406 52 AMIS 1980 2237 246 2257 219 ISIS-2 1988 8587 1570 8600 1720 操作步驟1 把數(shù)據(jù)輸入 stata 軟件1 / 10 2 變量的說(shuō)明Study 納入的討論Year 年份Death1 Aspirin 組的死亡人數(shù)Live1 Aspirin 組的存活人數(shù)Death

3、2 勸慰劑組的死亡人數(shù)Live2 勸慰劑組的死亡人數(shù)3 進(jìn)行 meta 分析metan death1 live1 dead2 live2, or labelnamevar=study, yearvar=year結(jié)果：2 / 10 以上結(jié)果分成兩部分 1meta 分析的合并統(tǒng)計(jì)量合并 OR值 OR=0.897,95%的可信區(qū)間 0.841,0.957 2 給出異質(zhì)性檢驗(yàn)的結(jié)果只要異質(zhì)性檢驗(yàn)的P值不小于 0.10 或者 I-squared小于 50%就可以認(rèn)為不存在異質(zhì)性,可以用應(yīng)固定效應(yīng)模型 stata 默認(rèn)的情形 ；假如質(zhì)性檢驗(yàn)的 P 值小于等于 0.10 或者 I-squared 大于 50

4、%,就不 同的組間存在異質(zhì)性,應(yīng)當(dāng)應(yīng)用隨機(jī)效應(yīng)模型隨機(jī)效應(yīng)模型的命令如下：metan death1 live1 dead2 live2, or labelnamevar=study, yearvar=year random 在運(yùn)行 meta 分析命令的同時(shí) stata 輸出森林圖,如下：3 / 10 由輸出的合并結(jié)果和漏斗圖可以得出,合并的 區(qū)間為 0.84 ,0.96 4 發(fā)表偏倚的檢驗(yàn),命令如下：1 gen logor=log_ES 2 gen selogor=_selogES 3 metabias logor selogor,graphbegg 輸出結(jié)果如下：4 / 10 OR值為 0.

5、90,95%可信發(fā)表偏倚主要看 begg 檢驗(yàn)的結(jié)果,由上圖可以看到發(fā)表偏倚假設(shè)檢驗(yàn)的 z 值為 1.20 ,p 值為 0.2300.05 ,可以認(rèn)為沒(méi)有發(fā)表偏倚；進(jìn)行發(fā)表偏倚檢驗(yàn)時(shí)stata同時(shí)給出的漏斗圖如下5 結(jié)論5 / 10 合并 OR小于 1,95%的可信區(qū)間不包括1,所以示 Aspirin有預(yù)防心肌梗塞后死亡的作用數(shù)值變量資料的 meta 分析Gotzsche 收集了有關(guān)短程小劑量強(qiáng)的松VS勸慰劑或非甾體抗炎藥治療類風(fēng)濕性關(guān)節(jié)炎的 7 個(gè)臨床隨機(jī)對(duì)比試驗(yàn) RCTs,觀看類風(fēng)濕性關(guān)節(jié)炎患者的關(guān)節(jié)壓痛指數(shù) rechie S index ,數(shù)據(jù)如表 2；表 2 7 個(gè)討論類風(fēng)濕性關(guān)節(jié)炎患

6、者關(guān)節(jié)壓痛指數(shù)試驗(yàn)發(fā)表試驗(yàn)組均數(shù)標(biāo)準(zhǔn)差對(duì)比組均數(shù)標(biāo)準(zhǔn)差名稱年份例數(shù)例數(shù)Jasni 1968 9 16.2 8.7 9 38.1 12.8 Dick 1970 24 17.6 8.0 24 40.7 13.0 Lee 1973 21 30.5 16.5 21 41.4 19.8 Berry 1974 12 13.0 11.0 12 23.7 11.1 Lee 1974 18 14.6 12.4 18 26.4 15.1 Stenber 1992 21 6.3 1.7 21 11.1 2.1 Geital 1995 20 10.8 4.7 20 16.3 7.7 1 在 stata 數(shù)據(jù)編輯器中建

7、立數(shù)據(jù)文件2 變量名稱study 納入的討論year 討論的年份n1 試驗(yàn)組的樣本含量6 / 10 mean1 試驗(yàn)組的均數(shù)sd1 試驗(yàn)組的標(biāo)準(zhǔn)差n2 對(duì)比組的樣本含量mean2 對(duì)比組的均數(shù)sd2 對(duì)比組的標(biāo)準(zhǔn)差3 進(jìn)行 meta 分析metan n1 mean1 sd1 n2 mean2 sd2, labelnamevar=study, yearvar=year輸出的結(jié)果判定異質(zhì)性的方法同上面定性資料的meta,有上圖可以看到,不同的組間存在異質(zhì)性,所以應(yīng)當(dāng)應(yīng)用隨機(jī)效應(yīng)模型 應(yīng)用隨機(jī)效用模型以后輸出的結(jié)果如下7 / 10 stata 同時(shí)輸出森林圖如下4 進(jìn)行發(fā)表偏倚的檢驗(yàn) Stata 命令： metabias _ES _seES,graphbegg 檢驗(yàn)結(jié)果 結(jié)果的說(shuō)明同分類資料的 meta 分析8 / 10 漏斗圖如