lines.predMexhaz: Lines method for a predMexhaz object

In mexhaz: Mixed Effect Excess Hazard Models

Lines 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'
lines(x, which = c("surv", "hazard"), conf.int =
TRUE, lty.pe = "solid", lty.ci = "dashed", ...)

Arguments

x	an object of class predMexhaz, corresponding to (excess) hazard and (net) survival predictions based on a survival model fitted with the mexhaz function. Predictions can be obtained for multiple times for one vector of covariables ("multitime") or for several vectors of covariables at one time point ("multiobs"). The lines() function only applies to the "multitime" type of predictions.
which	type of curve to be plotted. Selection can be made between "surv" (default value) for the (net) survival curve and "hazard" for the (excess) hazard.
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.
lty.ci	type of line used for drawing the confidence limits.
...	additional parameters that are directly passed to the lines function. These parameters, if not already used internally by the lines.predMexhaz function, will apply simultaneously to the point estimate and confidence limit curves.

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")
lines(Pred_Modbs2_2A, which="hazard", conf.int=FALSE)

mexhaz documentation built on Oct. 31, 2022, 5:08 p.m.