survreg函数出错:Invalid survival times for this distribution

项目场景:

使用R进行生存分析,建立回归模型,预测患者的生存时间。


问题描述

Cox比例风险模型,是关于风险的模型,自变量是时间t和影响因素X,因变量是风险(危险率)函数。它的特点是,仅仅研究危险因素对生存的影响,而不关心时间对生存的影响。而我要找的,是一种对时间T的回归模型。

调研发现,基于生存时间的威布尔分布的参数AFT(Accelerated Failure Time)生存回归模型,对时间建模。它常使用在工业领域,如研究零件寿命受温度的影响,因此把生存时间称为失效时间(Failure Time)。

然后又发现survreg方法可以自定义数据分布为威布尔分布,并预测生存时间。所以建立survreg函数并预测时间的代码:

# 构建weibull生存模型
fit <- survreg(Surv(OS_MONTHS, OS_STATUS) ~ AGE_AT_EVIDENCE_OF_METS + SEX + TUMOR_PURITY+SAMPLE_COVERAGE,data = train, dist = 'weibull')
# 模型预测
predict(fit, newdata = test, type = "response")

原因分析:

survreg方法不允许时间 = 0,因为在时间 = 0 发生事件,将导致未定义的估计量。(这我也说不清啥意思....)


解决方案:

将生存时间为0的数据都删掉,就不会报错了~

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值