2023年2月6日 星期一

[R 語言] continuity correction 細解

最近跟合作夥伴仔細討論了R 語言裡面, 診斷性統合分析的一些技術細節, 在這邊做個筆記, 診斷性統合分析主要會遇到的問題就是因為使用odds ratio去做計算, 所以可能遇到值為零而導致不能取log的問題, 在這邊就會使用 continuity correction 的技術去處理, 其實就是把零的值去加上0.5

但是在診斷性統合分析裡面常用的兩個package, 彼此之間的預設校正方法會有些不同, 因此在生成圖表的時候需要特別注意

首先是拿來生成forest plot 跟 sROC 的 mada package, 裡面常用的madad 函式是只要有一個值為零, 就會把所有研究的值都加上0.5

再來是生成funnel plot 的 metabin 函式, 其預設是只會針對單筆研究裏面有零的, 去把單筆研究裏面的值都加0.5

差別雖然細微, 但是使用的時候還是要盡量上下的圖形都是一致的, 這裡寫個小筆記避免忘記

ref.

1. https://cran.r-project.org/web/packages/mada/mada.pdf
2. https://search.r-project.org/CRAN/refmans/meta/html/metabin.html