用R中的竞争风险模型预测每种失效类型的可能性

我希望在历史数据上运行一个竞争风险模型，并在一个新的数据集中预测特定时期（比如一个时期）每种类型死亡的可能性。到目前为止，我已经研究了

timereg

软件包中的

comp.risk

和

cmprsk

中的

crr

，但我很难弄清楚如何使用他们的预测方法返回这些可能性

以

bmt

数据集（来自

timereg

软件包）和

comp.risk

为例，我想做如下操作：

m <- comp.risk(Surv(time, cause>0)~platelet+age+tcell, data=bmt,
               bmt$cause, causeS=1, resample.iid=1)
ndata <- data.frame(platelet=c(1,0,0), age=c(0,1,0), tcell=c(0,0,1), 
                    start.time=c(1, 1, 1), end.time=c(2, 2, 2))
out <- predict(m, newdata=ndata)

m <- comp.risk(Surv(rep(0, nrow(bmt)), time, cause>0)~platelet+age+tcell, data=bmt,
               bmt$cause, causeS=1, resample.iid=1)