如何修改R语言metaprop森林图的Common effect model标题

问题描述 #

我使用R语言的metaprop函数进行比例的Meta分析，在优化最终森林图时遇到两个问题：

• 无法移除“Common effect model”标题
• 无法将该标题重命名为常规的“Fixed effects model”

我的代码如下：

pes.summary=metaprop(cases,total,authoryear,data=dat,sm="PLO",method.ci="NAsm",method.tau="DL",incr=0.5,allincr=FALSE,addincr=FALSE,title="")
png("forestplotMassDXA.png",width=800,height=1200)
forest(pes.summary,
       xlim = c(0,1), pscale =1, 
       rightcols=c("effect", "ci","w.random"),
       rightlabs=c("Proportion", "95% C.I.","Weights"),
       leftcols = c("studlab", "event", "n"), 
       leftlabs = c("Study", "Successful", "Total"),
       xlab = "Proportion", 
       fs.xlab=12,
       fs.study=12,
       fs.study.lables=12,
       fs.heading=12,
       squaresize = 0.5, col.square="navy", col.square.lines="navy",
       col.diamond="navy", col.diamond.lines="navy",
       comb.fixed=FALSE,
       lty.fixed=0,
       lty.random=2, 
       Modality.study="square",
       type.random="diamond",
       ff.fixed="bold.italic",
       ff.random="bold.italic",
       hetlab = "Heterogeneity:",
       fs.hetstat=8,
       smlab="",
       print.Q=TRUE,
       print.pval.Q=TRUE,
       print.I2=TRUE,
       print.tau2=FALSE,
       col.by="black")

解决方法 #

1. 重命名为"Fixed effects model" #

直接在forest()函数中添加fixedlab参数，指定自定义标题：

png("forestplotMassDXA.png",width=800,height=1200)
forest(pes.summary,
       # 新增参数指定标题
       fixedlab = "Fixed effects model",
       # 若需显示固定效应模型结果则设为TRUE，仅改标题但隐藏结果则保持FALSE
       comb.fixed = FALSE,
       # 以下保留原有参数
       xlim = c(0,1), pscale =1, 
       rightcols=c("effect", "ci","w.random"),
       rightlabs=c("Proportion", "95% C.I.","Weights"),
       leftcols = c("studlab", "event", "n"), 
       leftlabs = c("Study", "Successful", "Total"),
       xlab = "Proportion", 
       fs.xlab=12,
       fs.study=12,
       fs.study.lables=12,
       fs.heading=12,
       squaresize = 0.5, col.square="navy", col.square.lines="navy",
       col.diamond="navy", col.diamond.lines="navy",
       lty.fixed=0,
       lty.random=2, 
       Modality.study="square",
       type.random="diamond",
       ff.fixed="bold.italic",
       ff.random="bold.italic",
       hetlab = "Heterogeneity:",
       fs.hetstat=8,
       smlab="",
       print.Q=TRUE,
       print.pval.Q=TRUE,
       print.I2=TRUE,
       print.tau2=FALSE,
       col.by="black")
dev.off()

2. 彻底移除固定效应模型行 #

若不需要显示固定效应模型的任何内容，需在metaprop()阶段就关闭固定效应计算，避免对象中包含相关结果：

# 生成metaprop对象时添加comb.fixed=FALSE
pes.summary=metaprop(cases,total,authoryear,data=dat,sm="PLO",method.ci="NAsm",method.tau="DL",incr=0.5,allincr=FALSE,addincr=FALSE,title="", comb.fixed=FALSE)

# 绘制森林图（其余参数不变）
png("forestplotMassDXA.png",width=800,height=1200)
forest(pes.summary,
       xlim = c(0,1), pscale =1, 
       rightcols=c("effect", "ci","w.random"),
       rightlabs=c("Proportion", "95% C.I.","Weights"),
       leftcols = c("studlab", "event", "n"), 
       leftlabs = c("Study", "Successful", "Total"),
       xlab = "Proportion", 
       fs.xlab=12,
       fs.study=12,
       fs.study.lables=12,
       fs.heading=12,
       squaresize = 0.5, col.square="navy", col.square.lines="navy",
       col.diamond="navy", col.diamond.lines="navy",
       comb.fixed=FALSE,
       lty.fixed=0,
       lty.random=2, 
       Modality.study="square",
       type.random="diamond",
       ff.fixed="bold.italic",
       ff.random="bold.italic",
       hetlab = "Heterogeneity:",
       fs.hetstat=8,
       smlab="",
       print.Q=TRUE,
       print.pval.Q=TRUE,
       print.I2=TRUE,
       print.tau2=FALSE,
       col.by="black")
dev.off()

3. 直接修改metaprop对象属性 #

如果前两种方法无效，可直接修改结果对象的固定效应模型名称属性：

# 修改标题为"Fixed effects model"
pes.summary$fixed.name <- "Fixed effects model"

# 若要移除标题，可将其设为空字符串
# pes.summary$fixed.name <- ""

# 绘制森林图（参数不变）
png("forestplotMassDXA.png",width=800,height=1200)
forest(pes.summary, ...) # 替换为你的原有参数
dev.off()

内容的提问来源于stack exchange，提问作者SirHugo