View source: R/plotPredictEventProb.R
| plotPredictEventProb | R Documentation |
Ploting time-dependent event risk predictions.
plotPredictEventProb(
x,
newdata,
times,
cause = 1,
xlim,
ylim,
xlab,
ylab,
axes = TRUE,
col,
density,
lty,
lwd,
add = FALSE,
legend = TRUE,
percent = FALSE,
...
)
x |
Object specifying an event risk prediction model. |
newdata |
A data frame with the same variable names as those that were
used to fit the model |
times |
Vector of times at which to return the estimated probabilities. |
cause |
Show predicted risk of events of this cause |
xlim |
Plotting range on the x-axis. |
ylim |
Plotting range on the y-axis. |
xlab |
Label given to the x-axis. |
ylab |
Label given to the y-axis. |
axes |
Logical. If |
col |
Vector of colors given to the survival curve. |
density |
Densitiy of the color – useful for showing many (overlapping) curves. |
lty |
Vector of lty's given to the survival curve. |
lwd |
Vector of lwd's given to the survival curve. |
add |
Logical. If |
legend |
Logical. If TRUE a legend is plotted by calling the function
legend. Optional arguments of the function |
percent |
Logical. If |
... |
Parameters that are filtered by |
Arguments for the invoked functions legend and axis are simply
specified as legend.lty=2. The specification is not case sensitive,
thus Legend.lty=2 or LEGEND.lty=2 will have the same effect.
The function axis is called twice, and arguments of the form
axis1.labels, axis1.at are used for the time axis whereas
axis2.pos, axis1.labels, etc. are used for the y-axis.
These arguments are processed via ...{} of
plotPredictEventProb and inside by using the function
SmartControl.
The (invisible) object.
Ulla B. Mogensen ulmo@biostat.ku.dk, Thomas A. Gerds tag@biostat.ku.dk
Ulla B. Mogensen, Hemant Ishwaran, Thomas A. Gerds (2012). Evaluating Random Forests for Survival Analysis Using Prediction Error Curves. Journal of Statistical Software, 50(11), 1-23. DOI 10.18637/jss.v050.i11
predictEventProbprodlim
# competing risk data
library(riskRegression)
library(pec)
set.seed(9)
d=sampleData(80)
csc1 <- CSC(Hist(time,event)~X1+X8, data=d)
nd=sampleData(3)
plotPredictEventProb(csc1, newdata=nd, cause=1, col=1:3)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.