日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

R语言ggsurv生存曲线一页多图的实现

發(fā)布時(shí)間:2023/12/20 编程问答 33 豆豆
生活随笔 收集整理的這篇文章主要介紹了 R语言ggsurv生存曲线一页多图的实现 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

R語(yǔ)言ggsurv生存曲線一頁(yè)多圖的實(shí)現(xiàn)

  • Problem
  • Solution


Problem

You want to put multiple graphs on one page.


Solution

The easy way is to use the multiplot function, defined at the bottom of this page. If it isn’t suitable for your needs, you can copy and modify it.

First, set up the plots and store them, but don’t render them yet. The details of these plots aren’t important; all you need to do is store the plot objects in variables.

multiplot <- function(..., plotlist=NULL, file, cols=1, layout=NULL) {library(grid)# Make a list from the ... arguments and plotlistplots <- c(list(...), plotlist)numPlots = length(plots)# If layout is NULL, then use 'cols' to determine layoutif (is.null(layout)) {# Make the panel# ncol: Number of columns of plots# nrow: Number of rows needed, calculated from # of colslayout <- matrix(seq(1, cols * ceiling(numPlots/cols)),ncol = cols, nrow = ceiling(numPlots/cols))}if (numPlots==1) {print(plots[[1]])} else {# Set up the pagegrid.newpage()pushViewport(viewport(layout = grid.layout(nrow(layout), ncol(layout))))# Make each plot, in the correct locationfor (i in 1:numPlots) {# Get the i,j matrix positions of the regions that contain this subplotmatchidx <- as.data.frame(which(layout == i, arr.ind = TRUE))print(plots[[i]], vp = viewport(layout.pos.row = matchidx$row,layout.pos.col = matchidx$col))}} }####創(chuàng)建生存對(duì)象 fit1 <- survfit(Surv(PFS, status)~FIB1, data=df) ggsurv1 <- ggsurvplot(fit1,risk.table=F,#生存統(tǒng)計(jì)統(tǒng)計(jì)表pval=TRUE,ggtheme = theme_bw(),ylab="PFS",legend=c(0.8,0.9),legend.title="",#log-rank檢驗(yàn)pval.method=TRUE)#添加檢驗(yàn)text ###FIB1 fit2 <- survfit(Surv(OS, status)~FIB1, data=df) ggsurv2 <- ggsurvplot(fit2,risk.table=F,#生存統(tǒng)計(jì)統(tǒng)計(jì)表pval=TRUE,ggtheme = theme_bw(),ylab="OS",legend=c(0.8,0.9),legend.title="",#log-rank檢驗(yàn)pval.method=TRUE) ggsurv2 ## multiplot(ggsurv1, ggsurv2, cols=2)

結(jié)果如下:兩張生存曲線圖片組合在一起了

總結(jié)

以上是生活随笔為你收集整理的R语言ggsurv生存曲线一页多图的实现的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。