View source: R/predict-timereg.r
pava.pred | R Documentation |
Make predictions of predict functions in rows mononotone using the pool-adjacent-violators-algorithm
pava.pred(pred, increasing = TRUE)
pred |
predictions, either vector or rows of predictions. |
increasing |
increasing or decreasing. |
mononotone predictions.
Thomas Scheike
data(bmt); ## competing risks add<-comp.risk(Event(time,cause)~platelet+age+tcell,data=bmt,cause=1) ndata<-data.frame(platelet=c(1,0,0),age=c(0,1,0),tcell=c(0,0,1)) out<-predict(add,newdata=ndata,uniform=0) par(mfrow=c(1,1)) head(out$P1) matplot(out$time,t(out$P1),type="s") ###P1m <- t(apply(out$P1,1,pava)) P1monotone <- pava.pred(out$P1) head(P1monotone) matlines(out$time,t(P1monotone),type="s")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.