points.predMexhaz {mexhaz} | R Documentation |
Points method for a predMexhaz object
Description
Function for adding to an already existing graphical window
the predicted (excess) hazard or (net) survival based on a predMexhaz
object.
Usage
## S3 method for class 'predMexhaz'
points(x, which = c("surv", "hazard"), conf.int =
TRUE, lty.pe = "solid", lty.ci = "dashed", ...)
Arguments
x |
an object of class |
which |
type of curve to be plotted. Selection can be made
between |
conf.int |
logical values allowing the user to decide whether to plot the confidence limits of the survival (or hazard). |
lty.pe |
type of line used for drawing the hazard/survival estimate (used when 'type="l"'). |
lty.ci |
type of line used for drawing the confidence limits (used when 'type="l"'). |
... |
additional parameters that are directly passed to the
|
See Also
predict.mexhaz
, plot.predMexhaz
Examples
data(simdatn1)
## Fit of a fixed-effect hazard model, with the baseline hazard
## described by a linear B-spline with two knots at 1 and 5 year and with
## effects of age (agecr), deprivation index (depindex) and sex (IsexH)
Mod_bs2_2 <- mexhaz(formula=Surv(time=timesurv,
event=vstat)~agecr+depindex+IsexH, data=simdatn1, base="exp.bs",
degree=2, knots=c(1,5), verbose=0)
## Prediction at several time points for one vector of covariates
Pred_Modbs2_2A <- predict(Mod_bs2_2, time.pts=seq(0,10,by=0.1),
data.val=data.frame(agecr=0,depindex=0.5,IsexH=1))
plot(Pred_Modbs2_2A, which="hazard", col="red")
points(Pred_Modbs2_2A, which="hazard", type="l", conf.int=FALSE)