R生存分析可视化-CSDN博客

本文链接：https://blog.csdn.net/qq_42424756/article/details/81605405
library(forestplot)
> RFS.forestplot
                                 V1         V2                 V3   V4   V5   V6              V7      V8
1                            Factor     Cohort Numbers of patient   NA   NA   NA       HR(95%CI) P Value
2                                                                   NA   NA   NA                        
3            Age (>60 : ≤ 60 years)                                 NA   NA   NA                        
4                                     Training              79/84 1.23 0.77 1.99 1.23(0.77-1.99)     0.4
5                                   Validation              84/69 1.02 0.63 1.64 1.02(0.63-1.64)     0.9
6                                        Whole            163/153 1.13 0.81 1.58 1.13(0.81-1.58)     0.5
7            Gender (Male : Female)                                 NA   NA   NA                        
8                                     Training             108/55 0.79 0.48 1.28 0.79(0.48-1.28)     0.3
9                                   Validation             109/44 1.27 0.74 2.20 1.27(0.74-2.20)     0.4
10                                       Whole             217/99 1.00 0.70 1.43 1.00(0.70-1.43)       1
11 Family cancer history (Yes : No)                                 NA   NA   NA                        
12                                    Training              47/96 0.73 0.42 1.25 0.73(0.42-1.25)     0.2
13                                  Validation              45/88 0.89 0.51 1.53 0.89(0.51-1.53)     0.7
14                                       Whole             92/184 0.79 0.54 1.16 0.79(0.54-1.16)     0.2
15          AFP  (≥ 20 : < 20 μg/L)                                 NA   NA   NA                        
16                                    Training              55/67 0.77 0.43 1.39 0.77(0.43-1.39)     0.4
17                                  Validation              53/62 1.88 1.07 3.28 1.88(1.07-3.28)    0.03
18                                       Whole            108/129 1.15 0.77 1.71 1.15(0.77-1.71)     0.5
19              Child grade (B : A)                                 NA   NA   NA                        
20                                    Training             10/103 1.56 0.61 3.97 1.56(0.61-3.97)     0.3
21                                  Validation               6/87 1.63 0.50 5.35 1.63(0.50-5.35)     0.4
22                                       Whole             16/190 1.54 0.74 3.20 1.54(0.74-3.20)     0.2
23 Pathologic stage (III-IV : I-II)                                 NA   NA   NA                        
24                                    Training             35/118 3.04 1.80 5.12 3.04(1.80-5.12)  <0.001
25                                  Validation             39/103 2.80 1.67 4.68 2.80(1.67-4.68)  <0.001
26                                       Whole             74/221 2.96 2.05 4.26 2.96(2.05-4.26)  <0.001
27 Histologic grade (G3-G4 : G1-G2)                                 NA   NA   NA                        
28                                    Training             56/105 0.96 0.58 1.60 0.96(0.58-1.60)     0.9
29                                  Validation              52/98 1.15 0.70 1.89 1.15(0.70-1.89)     0.6
30                                       Whole            108/203 1.04 0.73 1.48 1.04(0.73-1.48)     0.8
31     Vascular invasion (Yes : No)                                 NA   NA   NA                        
32                                    Training              39/96 1.59 0.89 2.84 1.59(0.89-2.84)     0.1
33                                  Validation              51/83 1.52 0.89 2.60 1.52(0.89-2.60)     0.1
34                                       Whole             90/179 1.54 1.05 2.28 1.54(1.05-2.28)    0.03
35       MSH (high risk : low risk)                                 NA   NA   NA                        
36                                    Training              77/86 2.64 1.61 4.35 2.64(1.61-4.35)  <0.001
37                                  Validation              57/96 1.74 1.07 2.81 1.74(1.07-2.81)    0.02
38                                       Whole            134/182 2.08 1.48 2.91 2.08(1.48-2.91)  <0.001



RFS.forestplot$V5 <- as.numeric(RFS.forestplot$V5)
RFS.forestplot$V4 <- as.numeric(RFS.forestplot$V4)
RFS.forestplot$V6 <- as.numeric(RFS.forestplot$V6)
View(RFS.forestplot)
RFS.forestplot[21,7] <- "1.63(0.50-5.35)"
RFS.forestplot[38,7] <- "2.08(1.48-2.91)"
RFS.forestplot[29,7] <- "1.15(0.70-1.89)"
RFS.forestplot[9,7] <- "1.27(0.74-2.20)"
RFS.forestplot[25,7] <- "2.80(1.67-4.68)"
RFS.forestplot[10,7] <- "1.00(0.70-1.43)"
RFS.forestplot[22,7] <- "1.54(0.74-3.20)"
View(RFS.forestplot)
d <- c(NA,NA,NA,0.2,0.2,0.4,rep(c(NA,0.2,0.2,0.4),8))

forestplot(as.matrix(RFS.forestplot[,c(1,2,3,7,8)]),mean=RFS.forestplot$V4,lower = RFS.forestplot$V5,
           upper = RFS.forestplot$V6,graph.pos = 4,is.summary = c(T,T,T,F,F,F,rep(c(T,F,F,F),8)),
           txt_gp = fpTxtGp(ticks = gpar(cex=1),
                            summary = gpar(cex = 0.8),cex=0.8,label=gpar(cex=0.9,col="black"),
                            xlab=gpar(cex = 0.8),
                            title=gpar(cex = 0.8)),boxsize = d,zero=1,
           lineheight = unit(3.8,"mm"),
           col=fpColors(box="steelblue", lines="black", zero = "black"),
           lwd.xaxis=1.5,graphwidth = unit(6,"cm"), clip =c(0, 4),lwd.ci=2,colgap = unit(0.5,"cm"))






forestplot.OS[22,7] <- "1.57(0.74-3.30)"
forestplot.OS[33,7] <- "1.26(0.72-2.20)"
forestplot.OS$V5 <- as.numeric(forestplot.OS$V5)
forestplot.OS$V4 <- as.numeric(forestplot.OS$V4)
forestplot.OS$V6 <- as.numeric(forestplot.OS$V6)

d <- c(NA,NA,NA,0.2,0.2,0.4,rep(c(NA,0.2,0.2,0.4),8))
forestplot(as.matrix(forestplot.OS[,c(1,2,3,7,8)]),mean=forestplot.OS$V4,lower = forestplot.OS$V5,
           upper = forestplot.OS$V6,graph.pos = 4,is.summary = c(T,T,T,F,F,F,rep(c(T,F,F,F),8)),
           txt_gp = fpTxtGp(ticks = gpar(cex=1),
                            summary = gpar(cex = 0.8),cex=0.8,label=gpar(cex=0.9,col="black"),
                            xlab=gpar(cex = 0.8),
                            title=gpar(cex = 0.8)),boxsize = d,zero=1,
           lineheight = unit(3.8,"mm"),
           col=fpColors(box="steelblue", lines="black", zero = "black"),
           lwd.xaxis=1.5,graphwidth = unit(6,"cm"), clip =c(0, 4),lwd.ci=2,colgap = unit(0.5,"cm"))
森林图
