提取生存分析预测时返回数字(0)

时间:2017-05-18 21:57:42

标签: r prediction survival-analysis

我使用survreg在R中运行了Weibull AFT生存分析,现在我正在尝试预测直到给定百分比存活率的时间。我已经使用1:98th四分位数创建了累积失败率的向量,并且可以调用任何四分位数的给定失败率的时间。但是,当我尝试调用相应的存活率时,对于低于0.35 R的四分位数,返回数字(0)。任何人都可以解释原因吗?

下面是我的一些代码,虽然不是一个可重现的例子 - 抱歉。我做了一个小的虚拟数据集,但没有遇到同样的问题。

> pct <- (1:98)/100 
> survprobs<-predict(BtCatmodel1, newdata=data.frame(group="B"), type="quantile", p=pct)
> survprobs[(1-pct)==0.4]
[1] 12.39598   # 40% survival at 12.4 days
>survprobs[(1-pct)==0.3]
numeric(0)    # ???

当我为失败率运行相应的百分位时,两者都可以正常工作。

> survprobs[pct==0.6]
[1] 12.39598    # 60% deaths at 12.4 days
> survprobs[pct==0.7]
[1] 18.01726    # 70% deaths at 18 days 

这也有效:

> (1-pct)[70]
[1] 0.3
> survprobs[(1-pct)==((1-pct)[70])]
[1] 18.01726

发生了什么事?

0 个答案:

没有答案