Contrary to standard survival analysis models, which rely on the assumption that the entire population will eventually experience the event of interest, mixture cure models allow to split the population in susceptible and non-susceptible (cured) individuals.
In this R package, we implement the semi-parametric proportional-hazards (PH) cure model of Sy and Taylor (2000) extended to time-varying covariates. In particular, the penPHcure function allows to:
estimate the regression coefficients and the baseline hazard function (conditional on being susceptible);
compute confidence intervals for the estimated regression coefficients using the basic/percentile bootstrap method;
perform variable selection based on the SCAD-penalized likelihood, as in Beretta and Heuchenne (2019).
Moreover, the penPHcure.simulate function allows to simulate data from a PH cure model, where the event-times are generated on a continuous scale from a piecewise exponential distribution conditional on time-varying covariates, using a method similar to the one described in Hendry (2014).
To install the latest release from CRAN:
install.packages(penPHcure)
To install the latest devel version from GitHub:
library(devtools)
install_github("a-beretta/penPHcure")
If you discover a bug or you have a suggestion to improve the package: https://github.com/a-beretta/penPHcure/issues
For help or more information contact: a.beretta@uliege.be
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.