R Kaplan Meier曲线按性别和危险因素划分的-divide曲线

install.packages("survival")
install.packages("survminer")
library(survival)
library(survminer)

# 示例数据
df <- data.frame(
  time = c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
  status = c(1, 1, 0, 1, 0, 1, 0, 1, 0, 1),
  gender = c("Male", "Female", "Male", "Female", "Male", "Female", "Male", "Female", "Male", "Female"),
  risk_factor = c("High", "Low", "High", "Low", "High", "Low", "High", "Low", "High", "Low")
)

# 拟合基本的Kaplan-Meier模型
km_model <- survfit(Surv(time, status) ~ 1, data = df)

# 绘制按性别划分的Kaplan-Meier曲线
ggsurvplot(km_model, data = df, risk.table = TRUE, pval = TRUE, break.time.by = 20,
            ggtheme = theme_minimal(), legend.title = "Gender",
            facet.by = "gender", palette = "jco")

# 绘制按危险因素划分的Kaplan-Meier曲线
ggsurvplot(km_model, data = df, risk.table = TRUE, pval = TRUE, break.time.by = 20,
            ggtheme = theme_minimal(), legend.title = "Risk Factor",
            facet.by = "risk_factor", palette = "jco")

# 绘制按性别和危险因素划分的Kaplan-Meier曲线
ggsurvplot(km_model, data = df, risk.table = TRUE, pval = TRUE, break.time.by = 20,
            ggtheme = theme_minimal(), legend.title = "Gender and Risk Factor",
            facet.by = c("gender", "risk_factor"), palette = "jco")