我想使用R中的Survival软件包,用威布尔分布假设拟合加速故障时间模型后,计算预期的存活率。
时间t的生存率公式为 S(t)= exp(-lambda * t ^ p)
lambda在哪里 lambda = exp(B * X) B是系数的向量 X是预测变量的向量
t是时间 p是Weibull分布的比例参数。
使用type =“ linear”或等效选项调用predict.survreg时,是否与上面的lambda相同?
在上述公式中,将predict.survreg(data,type =“ linear”)的输出插入lambda不会给出概率。传递分位数的向量并内插生存率的计算量太大。
predict.survreg(data,type =“ quantile”,p = seq(.01,.99,by = .01)
关于从线性预测变量到生存率的数学转换是什么的任何想法?