我想省略一些森林情节中存在巨大标准误差的研究,因为它们难以解释。但我不想改变估计。在玩具示例下面:
### load BCG vaccine data
data(dat.bcg)
### meta-analysis of the log relative risks using a random-effects model
res <- rma(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg,
slab=paste(author, year, sep=", "))
### Let's say I want to omit the first study, the rows argument doesn't work as expected
forest(res, rows = c(2:13))
Error in forest.rma(res, rows = c(2:13)) :
Number of outcomes does not correspond to the length of the 'rows' argument.
有什么想法吗?
答案 0 :(得分:1)
您可以使用NSFontAttributeName
构建森林图,将估计值和相应的采样差异传递给函数。使用forest()
参数,您可以省略不希望包含在绘图中的研究。然后将基于模型的摘要估计(使用完整数据集)添加到带有subset
的图中。使用玩具示例:
addpoly()