rm(list = ls())
library(scales)
library(readxl)
library(survival)
library(survminer)
library(magrittr)
library(ggpubr)
library(ggplot2)
library(lattice)
library(readr)
#建立导出集合
p<-list()
kmfit<-list()
m<-list()
NO3ACT0221 <- read_csv("F:/paper/NO3ACT0221.csv")
GC<-NO3ACT0221
setwd("F:/paper/NO2MachinelearningGCRF/")
colnames(GC)
GC[,c(3:4,6:7)]<-lapply(GC[,c(3:4,6:7)],as.factor)
GC<-GC[,-2]
GC<-data.frame(GC)
m <- split(GC, GC$center)
length(m)
for (i in 1:3) {
print(percent(i/length(m)))
n <- m[[i]] %>%
surv_group_by("HRG")
length(n)
for (a in 1:2) {
D<-n[[2]][[a]]
kmfit[[a]]<-survfit(Surv(failuretime,failure) ~ ACT, data=D)
#save the plot
p[[a]]<-ggsurvplot(kmfit[[a]],
risk.table = T,risk.table.y.text=F,tables.col="BLACK",
pval = T,pval.size=4,
title=" ",legend.labs=c("good","poor"),
xlab="Time(month)",ylab = "Survival Probability",
conf.int=T,conf.int.style="ribbon",linetype=1,
surv.median.line="hv", xlim=c(0,60),
palette = c("#E7B800", "#2E9FDF"),break.y.by=0.25,break.x.by=12,
pval.coord=c(48,0.9),ncensor.plot = FALSE)
ggpubr::ggexport(p,filename = paste0("survival","m=",i,"n=",a,".tiff"), device=cairo_ps,width = 1800,height = 1800,res = 300)
}
}
R语言for循环结果异常,想根据 center 和 HRG对数据集进行亚组分析,但是跑出来的12张图中有好几张空白,只有6张图是想要得到的,请问如何在代码设计中就去掉冗余的图片,或者此代码是否存在什么问题