我有一个数字time变量,我希望使用R的Kaplan实现来观察我的数据。
km <- survfit(Surv(time,event)~1)
summary(km)但是我有特定的处理时间范围,比如1-3,4-6等等,但是我的方法返回一个错误。
km <- survfit(Surv(time=c(c(1,3),c(4,6),c(7,9)),event)~1)
summary(km)这可以用另一段代码实现吗?还是我必须更改time变量以满足我的需要(即重新编码为类别)?
发布于 2013-12-19 14:26:04
您提供的time变量应该包括数据集中各个观测的生存时间。
不确定我是否正确地理解了您的问题,但是如果您希望得到特定时间点的KM估计,您可以使用summary函数中的survfit时间参数来获得这些估计值。
summary(km, times=c(1,3,4,6,7,9))当然,您也可以用plot(km)绘制KM曲线来可视化它。
https://stackoverflow.com/questions/20683590
复制相似问题