用。R语言写生存分析代码

以下是一个用R语言进行生存分析的示例代码：

# 导入生存分析包
library(survival)

# 导入数据
data(lung)

# 查看数据
head(lung)

# 构建生存对象
# status为事件状态，1为死亡，0为存活
# time为事件时间，即患者死亡或最后随访时间
surv_obj <- with(lung, Surv(time, status))

# 基本生存分析
fit <- survfit(surv_obj ~ 1)
summary(fit)

# 绘制生存曲线
plot(fit, main="Survival Curve")

# Cox比例风险回归模型
cox_model <- coxph(surv_obj ~ age + sex + ph.ecog + ph.karno + pat.karno + meal.cal + wt.loss, data=lung)
summary(cox_model)

# 绘制Cox回归系数的森林图
plot(cox_model)

以上代码使用的数据集是R自带的lung数据集，该数据集包含了肺癌患者的基本信息、治疗信息以及生存情况。生存分析基于患者的事件时间和事件状态，使用Kaplan-Meier方法得出生存曲线，使用Cox比例风险回归模型分析影响生存的因素。