我试图使用生存包中的第一个coxph,然后使用survMisc包中的cutp()来找到变量的最佳切割点。我的数据被截断,因此我的coxph估算如下:
cph1 <- coxph(Surv(randtime, endtime, recurrence) ~ con + x + y, data=data)
然而,当我尝试使用cutp()
时:
cutp(cph1, var=con, plot=TRUE) it gives me an error:
1:D:NA / NaN参数出错 另外:有50个或更多警告(使用警告()查看前50个)
警告() 警告信息: 1:在Surv中(s1 [,1],s1 [,2]):无效状态值,转换为NA
如果我处理的数据没有被截断,只使用Surv(endtime, recurrence)..)
cutp()
效果很好。
有没有办法找到带有左截断数据的cutpoint
的{{1}}?
提前致谢!