通過R語言繪制森林圖(僅用HR/OR等,以及CI繪制)


以HR,HigherCI,LowerCI為例,表格中要有這三列

先導入表格lq_2

正態化數據分布

lnhr <- log(lq_2[,"HR"])
lnuci <- log(lq_2[,"HigherCI"])
lnlci <- log(lq_2[,"LowerCI"])
selnhr <- (lnuci-lnlci)/(2*1.96)

meta分析

MetaHR = metagen(TE = lnhr$HR,
                 seTE = selnhr$HigherCI,
                 sm = "HR",
                 data = lq_2,
                 # studlab = paste(lq_1$Clinic_Pt_Im,lq_1$P.value,sep = " P:"),
                 studlab = lq_2$`Clinic-Pt-Rad`,
                 random = T,
                 #backtransf = T
                 )

圖像風格設置

settings.meta('JAMA')
settings.meta('RevMan5')
settings.meta('meta4')

settings.meta('reset') # 重置風格

繪制森林圖(經測試選擇某種風格后修改森林圖參數會無效)

具體參數可用 ?forest 查詢 按照自己的需求調整

forest(MetaHR,
       # xlab = "hahahaha",
       label.left = "C-index = 0.840                                                                        ",
       #label.right = " i am right",
       colgap.studlab = "2cm",
       #colgap.forest.left = "0.2cm",
       colgap.forest.right = "0.5cm",
       colgap.right = "1cm",
       random = FALSE, 
       leftlabs = c("Clinic_Pt_Rad", "HR"),
       leftcols = c("studlab","HR"),
       rightcols = c("ci","P.value"),
       )

效果圖

參考視頻Meta分析自己用-森林圖 R語言(高分SCI繪圖模仿)03


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM