2022年Stata在meta分析中的應用

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

2、勸慰劑組總例數死亡例數總例數死亡例數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 把數據輸入 stata 軟件1 / 10 2 變量的說明Study 納入的討論Year 年份Death1 Aspirin 組的死亡人數Live1 Aspirin 組的存活人數Death

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

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

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

6、者關節(jié)壓痛指數試驗發(fā)表試驗組均數標準差對比組均數標準差名稱年份例數例數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 數據編輯器中建

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