我对 10 项研究进行了荟萃分析,现在正在将结果绘制在森林图上。然而,我遇到了一个恼人的功能,即具有较小 SE(适合重量框)的研究结果变成白色,并且似乎不遵守 color.study 参数。还有其他方法可以绕过这个吗? 如果有必要的话,我不介意完全失去重量盒,因为重量几乎相同。
df<-data.frame(study=c('s1','s2','s3','s4','s5','s6','s7','s8','s9','s10'),
SMD=c(0.016416141,0.299488949,0.024292283,0.110415656 ,0.160754906,-0.018112735,0.027915812,-0.025253501,-0.024464375,0.001146802),
se=c(0.013495131,0.021345836,0.010902911,0.006237354,0.014503043,0.012365650,0.009825055,0.006861558,0.018569984,0.021838104))
## Run meta-analysis
m.gen <- metagen(TE = SMD,
seTE = se,
studlab = study,
data = df,
sm = "SMD",
fixed = FALSE,
random = TRUE,
method.tau = "REML",
hakn = FALSE,
prediction = TRUE,
title = "Test")
## Plot forest plot
forest.meta(m.gen,
sortvar = TE,
prediction = TRUE,
print.tau2 = FALSE,
smlab="Test",
print.I2=FALSE,
print.I2.ci=FALSE,
col.study='black'
)