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
predictEventProb
prodlim
# 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.