R语言Metafor包:计算并显示比值比而非对数比值比
如何在metafor包的森林图中直接显示比值比(OR)而不是对数比值比(logOR)
当然有超简便的实现方法!metafor的forest()函数自带参数可以直接帮你完成尺度转换,完全不用手动计算OR再重新建模绘图。
核心方法:使用trans=exp参数
不管你已经拟合了logOR的模型,还是刚完成效应量计算,只需要在调用forest()时加上trans=exp,就能把x轴和效应值都转换为原始的OR尺度。
举个完整的示例代码:
# 加载metafor包 library(metafor) # 用内置数据计算logOR(measure="OR"返回的就是对数转换后的效应量) dat <- escalc(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg) # 拟合meta分析模型 model <- rma(yi, vi, data=dat) # 绘制森林图,直接显示OR forest(model, trans=exp, # 关键参数:将logOR转换为OR xlab="Odds Ratio (OR)", # 修改x轴标签为OR digits=c(2,2,2) # 控制效应值和置信区间的小数位数 )
额外说明
trans=exp是通用转换参数,如果你拟合的是其他对数转换的效应量(比如log相对风险logRR),用这个参数也能转回到原始尺度。- 加上
digits参数可以让显示的OR数值更美观,三个数值分别对应单个研究的效应值、汇总效应值、置信区间的小数位数。 - 森林图里的所有数值(包括研究的效应量、汇总效应量)都会自动转换为OR,不需要你手动处理原始数据。
内容的提问来源于stack exchange,提问作者ragzoxaim
相关产品推荐
相关产品推荐

