写一段生存分析R语言的代码

下面是一段生存分析的 R 语言代码：

library(survival)

# 读取数据
data(lung)

# 构建生存模型
fit <- survfit(Surv(time, status) ~ age + sex, data = lung)

# 生存曲线
plot(fit, lty = 1:2, xlab = "时间(天)", ylab = "生存概率")

# 半衰期
library(durations)
halflife(fit)

这段代码使用了 survival 包来进行生存分析，首先导入了该包，然后读取了 lung 数据集，之后建立了一个生存模型，模型的自变量包括年龄和性别。最后，通过 plot 函数绘制了生存曲线，并通过 durations 包中的 halflife 函数计算了半衰期。