單組率meta分析，發(fā)表偏倚檢驗、敏感性分析、亞組分析怎么做？

先說重點

1 、對于可以用metan完成森林圖的數(shù)據(jù)（也就是無極端值或偏態(tài)數(shù)據(jù)），可以實現(xiàn)發(fā)表偏倚檢驗、敏感性分析、亞組分析的操作；

2 、對于只能用metaprop完成森林圖的數(shù)據(jù)（數(shù)據(jù)存在極端值或非正態(tài)分布），發(fā)表偏倚檢驗、敏感性分析都做不了（做了結(jié)果也不對），但可以做亞組分析；

3 、其實單組率的meta分析，沒必要做發(fā)表偏倚檢驗，敏感性分析的意義也不大。

接著聊操作

1、 發(fā)表偏倚檢驗

單組率meta分析，不需要做發(fā)表偏倚檢驗，為什么？因為單組率只是一個描述性結(jié)果，而不是差異比較結(jié)果，不存在所謂的“陽性”結(jié)果或有統(tǒng)計學(xué)意義的結(jié)果。
然而，一定要做發(fā)表偏倚檢驗的話，操作上是可行的。以下圖的這些數(shù)據(jù)為例，我們已經(jīng)知道它們可以使用“metan”完成效應(yīng)量合并的森林圖。

在這個過程中，我們生成了p和se，利用它們就可以實現(xiàn)發(fā)表偏倚檢驗。

運行代碼：metabias? p se, egger graph得到egger檢驗的結(jié)果，P = 0.010 < 0.05，提示有顯著的發(fā)表偏倚。

這個結(jié)果沒有出乎意料，事實上，單組率的meta分析，“存在”發(fā)表偏倚的可能性，我認為是在95%以上，因為它們肯定都是“陽性”結(jié)果。這也是為什么我不認同單組率的meta分析需要做發(fā)表偏倚檢驗。

2 、敏感性分析

單組率的meta分析，如果數(shù)據(jù)可以通過“metan”實現(xiàn)的，敏感性分析跟常規(guī)meta分析一樣。
按照截圖的設(shè)置，選擇主菜單（Main）和二級菜單（Effect Opts）的相關(guān)參數(shù)。

以上菜單操作，就是以率(p)和標(biāo)準(zhǔn)誤(se)為數(shù)據(jù)，通過隨機效應(yīng)模型完成敏感性分析，點擊“OK”得到結(jié)果。

同樣的，結(jié)果毫無疑問是“穩(wěn)定”。不管剔除哪個研究，剩余研究的合并結(jié)果都不會變成沒有統(tǒng)計學(xué)意義。
用“metaprop”完成單組率meta分析的數(shù)據(jù)，為什么不能按上面的方法做發(fā)表偏倚檢驗或敏感性分析呢？我們用實例驗證。

數(shù)據(jù)中，“p, se”是通過metan完成森林圖時，通過公式計算得到的率和標(biāo)準(zhǔn)誤；而“_es, _sees”是通過metaprop對應(yīng)的率和標(biāo)準(zhǔn)誤。
既然森林圖的合并結(jié)果是用“_es, _sees”計算得到的，用“p, se”做敏感性分析就不合適。那能不能用前者做敏感性分析呢？答案也是不行的，我們可以操作一下。

此時，我們要看的是，總合并結(jié)果：P (95%CI) = 0.22 (0.05, 0.39)。事實上，用metaprop得到的合并結(jié)果是P (95%CI) = 0.20 (0.04, 0.44)。
結(jié)果的差異在于這是兩種不同的算法，因此，metaprop完成的單組率meta分析，不能用metaninf或metainf做敏感性分析。

3 亞組分析

因為要做亞組分析，因此需要整理一個分組變量，如下圖的“group”，分為A、B兩個亞組。

操作很簡單，我們已經(jīng)知道了metaprop實現(xiàn)單組率meta分析森林圖的代碼：metaprop case n, random cimethod(exact) lcols( study) xlab(0.2, 0.4, 0.6, 0.8) dp(3)
只需要在代碼的最后加上：by(xx)，也就是：metaprop case n, random cimethod(exact) lcols( study) xlab(0.2, 0.4, 0.6, 0.8) dp(3) by(group)

Ok，亞組分析完成了。
要提醒大家的一點是：metaprop的亞組分析，分組后，亞組的研究數(shù)量在3篇或以下，是不顯示異質(zhì)性檢驗結(jié)果的。如果審稿人有疑問，解釋一下即可。

今天的分享就到這了，meta分析學(xué)習(xí)過程，我們會碰到各式各樣的問題，如果你想更快、更高效地掌握meta分析，早日發(fā)表SCI，請不要錯過爾云間的meta分析培訓(xùn)班！