用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)
现在,我们可以使用ggsurvplot()
函数来绘制生存曲线。这个函数使用ggplot2
库来创建一个美观的图形,并根据生存曲线的不同组进行着色。
# 绘制生存曲线
ggsurvplot(surv_fit, data =
还没有评论,来说两句吧...