用R语言绘制优美的生存曲线

用R语言绘制优美的生存曲线

生存曲线是生物学和医学研究中常用的可视化工具，用于描述个体或群体在一段时间内存活的概率。在R语言中，我们可以使用一些库来绘制漂亮的生存曲线，例如survival和ggplot2。

首先，我们需要加载所需的库并准备数据。假设我们有一个包含个体生存时间和事件状态（是否发生事件）的数据集。以下是一个示例数据集，包含了10个个体的生存时间和事件状态：

# 加载所需的库
library(survival)
library(ggplot2)

# 创建一个示例数据集
time <- c(12, 15, 18, 24, 30, 36, 42, 48, 54, 60)
status <- c(1, 1, 0, 1, 1, 0, 0, 1, 0, 1)

# 创建生存时间和事件状态的数据框
data <- data.frame(time, status)

接下来，我们可以使用survfit()函数计算生存曲线的估计值。该函数可以接受一个生存对象作为参数，我们可以使用Surv()函数创建一个生存对象。

# 计算生存曲线的估计值
surv_fit <- survfit(Surv(time, status) ~ 1, data = data)