変量効果メタ分析モデルの等高線強化ファンネル プロットを作成しようとしています。パッケージのfunnel.rma()機能を使用しています。meta私が直面している問題は、要約推定値を表す垂直線から外れた漏斗状の対角線 (信頼区間) が、固定効果推定値に対してのみ表示されていることです。FE 推定値の垂直線を消すことはできますが、RE 推定値の垂直線に対角線を付けることができません。
私が使用しているデータは
私が使用しているコードは次のとおりです。
require(meta)
g <- c(0.88, 0.66, -0.02, 0.34, 0.30, -0.12, 0.70, 0.42, 2.66, 0.10, 0.73, -0.54, 0.72, -0.54)
g_v <- c(0.12, 0.05, 0.06, 0.06, 0.09, 0.08, 0.03, 0.07, 0.10, 0.06, 0.11, 0.07, 0.13, 0.06)
g_se <- sqrt(g_v)
FC<-metagen(g, g_se, method.tau = "REML", hakn = TRUE, comb.fixed = FALSE)
funnel(FC, lty.random = 1, xlim = c(-3, 3), xlab = "Effect size", ylab = "Standard Error", cex = .75, col = 1, bg = 1, contour=c(0.00001, 0.95), col.contour=c("grey", "white"), family = "serif", cex.lab = 1.5)