Function for adding to an already existing graphical window
the predicted (excess) hazard or (net) survival based on a predMexhaz
object.
1 2 3 
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. 
lty.ci 
type of line used for drawing the confidence limits. 
... 
additional parameters that are directly passed to the

predict.mexhaz
, plot.predMexhaz
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16  data(simdatn1)
## Fit of a fixedeffect hazard model, with the baseline hazard
## described by a linear Bspline 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.mexhaz(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", conf.int=FALSE)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.