遇到的问题:
生存分析进行过程中,一直发生报错,代码及 错误如下:
load("surv_data.Rdata")
library(survival)
library(survminer)
fit <- survfit(Surv(OS.time,OS)~class_FEV1,dd1)
ggsurvplot(fit)
ggsurvplot(fit)
Error in data.frame(..., check.names = FALSE) : arguments imply differing number of rows: 59, 0, 118
解决办法:
原来,OS这个变量我设置成了因子,而在生存曲线绘制时,OS需要数值型
load("surv_data.Rdata")
library(survival)
library(survminer)
dd1$OS <- as.numeric(dd1$OS) ######转换为数值型
fit <- survfit(Surv(OS.time,OS)~class_FEV1,dd1)
ggsurvplot(fit)
网友评论