修改元分析森林图

编程问答 2022-05-29

问题描述

我想对我制作的一些森林图进行两个具体的更改，而且我对 ggplot 图形更加熟悉，所以我很挣扎。

这是我的数据：

library(Metafor)
library(Meta)        

dat<-structure(list(no = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,25,26,28),ref = c("Smith 2015","Smith 2015","Johnson 2016","Doe 2018","Allens 2016","Happy 2016","PrezFo 2017","Smithy 2008","Bernie 2016","Saint-Cyr 2012","Bernie 2014","Carp 2017","Mountain 2014","Sour 2017","Weinerschnitzel 2016","Word 2016","Park 2015","John 2012","Fansa 2008","Bunt 2011","Vega 2009","Weinerschnitzel 2004","McDonalds 2012","Allens 2012","Amex 2004","Carp 2020"),year = c(2015,2015,2016,2018,2017,2008,2012,2014,2011,2009,2004,2020),Level_of_evidence = c(2,NA,4),Type = c("GUT","SCHMAP","GUT","SCHMAP"),TypeTwo = c("GUT","SCHMAP"
        ),proced = c(1,1,0),Comparative_GUT_vs_SCHMAP = c("yes","yes","no","No"),no_patients = c(39,30,34,96,111,40,56,36,29,138),no_of_flaps = c(54,37,164,154,44,24,101,49,58,33,32,27,31,265),age = c(48,48,48.6,50,46,43,39.4,45.7,31.6,48.5,50.4,52.1,NA),revision = c(NA,NA)),class = c("spec_tbl_df","tbl_df","tbl","data.frame"),row.names = c(NA,-26L),spec = structure(list(
        cols = list(no = structure(list(),class = c("collector_double","collector")),ref = structure(list(),class = c("collector_character",year = structure(list(),Level_of_evidence = structure(list(),Type = structure(list(),TypeTwo = structure(list(),proced = structure(list(),Comparative_GUT_vs_SCHMAP = structure(list(),no_patients = structure(list(),no_of_flaps = structure(list(),age = structure(list(),revision = structure(list(),"collector"))),default = structure(list(),class = c("collector_guess",skip = 1L),class = "col_spec"))
    
    ies.da=escalc(xi=revision,ni=no_of_flaps,measure="PFT",data=dat,add=0)
    subganal.procedure=rma(yi,vi,data=ies.da,mods=~TypeTwo,method="DL")
    pes.summary=Metaprop(revision,no_of_flaps,ref,byvar=TypeTwo,sm="PFT",method.ci="NAsm",method.tau="DL",incr=0.5,allincr=FALSE,addincr=FALSE,tau.common=TRUE,tau.preset=sqrt(subganal.procedure$tau2),title="")



## The plot
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","Cases","Total"),xlab = "Prevalence",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,comb.random=TRUE,lty.fixed=0,lty.random=2,type.study="square",type.random="diamond",ff.fixed="bold.italic",ff.random="bold.italic",hetlab = "",fs.hetstat=10,resid.hetstat = FALSE,smlab="",print.Q=TRUE,print.pval.Q=TRUE,print.I2=TRUE,print.tau2=FALSE,col.by="grey",digits=6)

哪个应该产生这个图：

enter image description here

我正在尝试对图表进行 2 个非常具体的更改：

  1. 将每项“SCHMAP”研究置于每项“GUT”研究之上
  2. 更改底部流行度栏的比例。在这种情况下，我希望它不是 0-1，而是从 0 到 0.5，以便在视觉上“放大”小结果

理想情况下它会像这样结束：

enter image description here

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）