让我们模拟R中的一些数据:
- n < - 10000
- h < - 0.5
- t < - -log(runif(n))/ h
该代码模拟了危险函数的存活时间,即常数。
R语言生存分析Survival analysis原理与晚期肺癌患者分析案例
- 事件< - 1 *(t <5)
- 时间< - t
- obstime [obstime> = 5] < - 5
现在让我们使用R中的生存包绘制估计的生存函数:
survfit < - survfit(Surv(obstime,event)~1)