用R语言绘制生存曲线及其置信区间

用R语言绘制生存曲线及其置信区间

生存曲线是生存分析中常用的可视化工具，用于描述事件发生的概率随时间的变化情况。在R语言中，我们可以使用survival包中的plot函数来绘制生存曲线，并使用置信区间来表示曲线的不确定性。

下面是一个示例，演示如何使用R语言绘制生存曲线及其置信区间。

首先，我们需要安装并加载survival包，以便使用其中的函数和数据集。

install.packages("survival")  # 安装survival包
library(survival)            # 加载survival包

接下来，我们将使用survival包中的lung数据集作为示例数据。该数据集包含了肺癌患者的生存时间和一些其他相关信息。

data(lung)  # 加载lung数据集

现在，我们可以使用survfit函数计算生存曲线的估计值，并使用conf.int参数来计算置信区间。

# 计算生存曲线的估计值和置信区间
surv_fit <- survfit(Surv(time, status) ~ 1, data = lung)
surv_fit_ci <- survfit(Surv(time, status) ~ 1, data = lung, conf.int = TRUE)

接下来，我们可以使用plot函数来绘制生存曲线及其置信区间。