@article \{st0275,
author = "Crowther, M. J.",
author = "Lambert, P. C.",
title = "Simulating complex survival data",
journal = "Stata Journal",
publisher = "Stata Press",
address = "College Station, TX",
volume = "12",
number = "4",
year = "2012",
abstract = "Simulation studies are essential for understanding and evaluating both
current and new statistical models. When simulating survival times, one often
assumes an exponential or Weibull distribution for the baseline hazard function,
with survival times generated using the method of Bender, Augustin, and Blettner
(2005, Statistics in Medicine 24: 1713\endash1723). Assuming a constant or monotonic
hazard can be considered too simplistic and can lack biological plausibility in
many situations. We describe a new user\endashwritten command, survsim, which allows
the user to simulate survival times from two\endashcomponent parametric mixture
models, providing much more flexibility in the underlying hazard. Standard parametric
distributions can also be used, including the exponential, Weibull, and
Gompertz. Furthermore, survival times can be simulated from the all\endashcause distribution
of cause\endashspecific hazards for competing risks by using the method of
Beyersmann et al. (2009, Statistics in Medicine 24: 956\endash971). A multinomial distribution
is used to create the event indicator, whereby the probability of experiencing
each event at a simulated time t is the cause\endashspecific hazard divided by
the all\endashcause hazard evaluated at time t. Baseline covariates can be included in all
scenarios. We also describe the extension to incorporate nonproportional hazards
in standard parametric and competing\endashrisks scenarios.",
pages = "674-687(14)",
url = "http://www.stata-journal.com/article.html?article=st0275"
}