plot.predMexhaz: Plot method for a predMexhaz object In mexhaz: Mixed Effect Excess Hazard Models

Description

Function for plotting the predicted (excess) hazard or (net) survival based on a `predMexhaz` object.

Usage

 ```1 2 3``` ```## S3 method for class 'predMexhaz' plot(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 `plot()` 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 value 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 `plot` function. These parameters will apply simultaneously to the point estimate and confidence limit curves.

`predict.mexhaz`, `points.predMexhaz`
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15``` ```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") ```