Metafor森林图:双反正弦转换后转为百分比的实现方法咨询
将逆PFT转换结果转为百分比的实现方法
完全可以把逆PFT转换和乘100操作结合,实现森林图显示百分比的需求。你只需要自定义一个复合转换函数,先执行transf.ipft.hm逆转换,再将结果乘以100即可。修改后的代码如下:
# 计算效应量 mort_dat_es <- escalc(xi = deaths, ni = admissions, data = mort_dat, measure = "PFT", add = 0, slab=auth_date) # 合并效应量 mort_dat_pes <- rma(yi, vi, data = mort_dat_es) # 自定义复合转换函数:先逆PFT转换,再乘100转百分比 transf_pct <- function(x, ni) { 100 * transf.ipft.hm(x, ni) } # 转换合并结果为百分比 mort_pes <- predict(mort_dat_pes, transf = transf_pct, targ = list(ni = mort_dat$admissions)) # 绘制森林图,使用自定义转换函数,同时修改x轴标签为百分比 metafor::forest(mort_dat_pes, transf = transf_pct, targ = list(ni = mort_dat$admissions), addpred=TRUE, predstyle = "bar", xlab = "医院死亡率 (%)")
关键改动说明:
- 自定义
transf_pct函数:把transf.ipft.hm的输出直接乘以100,一步完成逆转换和百分比转换 - 把
predict和forest函数中的transf参数替换为自定义的transf_pct - 添加
xlab参数,把x轴标签明确改为百分比形式,让森林图更直观
内容的提问来源于stack exchange,提问作者Megan C
相关产品推荐
相关产品推荐

