遇到的問(wèn)題:
生存分析進(jìn)行過(guò)程中,一直發(fā)生報(bào)錯(cuò)瑞躺,代碼及 錯(cuò)誤如下:
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
解決辦法:
原來(lái),OS這個(gè)變量我設(shè)置成了因子兴想,而在生存曲線繪制時(shí)幢哨,OS需要數(shù)值型
load("surv_data.Rdata")
library(survival)
library(survminer)
dd1$OS <- as.numeric(dd1$OS) ######轉(zhuǎn)換為數(shù)值型
fit <- survfit(Surv(OS.time,OS)~class_FEV1,dd1)
ggsurvplot(fit)