在R中是否有一种方法可以显示元对象或森林对象中影响大小的自定义摘要？

Question

我需要提出疫苗有效性研究的元分析。本质上，我需要作为疫苗效力或"VE“的简要统计，即"1-RR”而不是"RR"，或者是"1-OR“而不是" or”。目前，元对象允许这些标准的效果大小大小为RR和OR，但是是否有一种方法来表示自定义输出？

元none包有大量的转换列表(请参阅"transf“选项)，这些转换都不允许自定义输入，"Meta”包也没有自定义输入。

是否有其他程序或方法需要注意？

它可能如下所示(使用元from包中的数据)：

### load metafor package
library(metafor)

### load BCG vaccine dataset
data(dat.bcg)
### Example of what it might look like with a custom summary estimate
dat <- escalc(measure="1-RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)

Answer 1

在元计算的predict()函数中，可以指定自定义转换。因此，您使用(日志转换的)风险比率进行元分析，然后根据您认为合适的情况进行转换。例如：

library(metafor)

# compute log risk ratios and corresponding sampling variances
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)

# fit random-effects model
res <- rma(yi, vi, data=dat)
res

# estimated pooled risk ratio (with 95% CI/PI)
predict(res, transf=exp)

# estimated pooled 1-RR (with 95% CI/PI)
predict(res, transf=function(x) 1-exp(x))