各位大佬,有没有知道怎么分三类,high,mid,low,然后做生存曲线??
data(myeloma)
head(myeloma)
res.cut <- surv_cutpoint(myeloma, time = "time", event = "event",
variables = c("DEPDC1", "WHSC1", "CRIM1"),
progressbar =F
)
plot(res.cut, "WHSC1", palette = "npg")
res.cat <- surv_categorize(res.cut,labels = c('high','low'))
head(res.cat)
library("survival")
fit <- survfit(Surv(time, event) ~DEPDC1, data = res.cat)
ggsurvplot(fit, data = res.cat, risk.table = TRUE, conf.int = TRUE)