如何移除metafor森林图中的估计刻度以制作堆叠图?
解决方案
要移除堆叠森林图中高亮的汇总估计刻度(即由addpred=TRUE生成的汇总效应菱形标记),只需修改forest()函数的addpred参数为FALSE即可:
dat.bcg dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, slab=paste(author, ", ", year, sep=""), data=dat.bcg) res <- rma(yi, vi, data=dat) predict(res, transf=exp, digits=2) par(mfrow=c(2,1),mar=c(0,4,0,4)) # 保留表头的森林图,移除汇总刻度 forest(res, addpred=FALSE, header=TRUE, xlim=c(-8,5), atransf=exp, at=log(c(.05, .25, 1, 4)),showweights=TRUE) # 无表头的森林图,移除汇总刻度 forest(res, addpred=FALSE, header=FALSE, xlim=c(-8,5), atransf=exp, at=log(c(.05, .25, 1, 4)),showweights=TRUE)
如果仅需移除其中一个图的汇总刻度,仅修改对应forest()调用的addpred参数即可。
关于你提到的亚组分析权重总和无法达到100%的问题:亚组内的权重是基于该亚组内部研究的方差计算的,各亚组权重总和本身就不等于全局的100%(全局汇总分析的权重总和才为100%)。若需让每个亚组内的权重相对总和为100%,可针对每个亚组单独拟合模型,或手动调整权重显示逻辑,但这不属于当前移除估计刻度需求的范畴,如需进一步处理可单独发起提问。
内容的提问来源于stack exchange,提问作者Kitty Cat
相关产品推荐
相关产品推荐

