为RE元分析制作轮廓增强漏斗图

时间:2013-04-29 16:13:20

标签: r

我正在尝试为随机效应元分析模型制作轮廓增强漏斗图。我正在使用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)

0 个答案:

没有答案