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–1723). Assuming a constant or monotonic hazard can be considered too simplistic and can lack biological plausibility in many situations. We describe http://www.stata-journal.com/software/a new user-written command, survsim, which al- lows the user to simulate survival times from two-component parametric mixture models, providing much more flexibility in the underlying hazard. Standard para- metric distributions can also be used, including the exponential, Weibull, and Gompertz. Furthermore, survival times can be simulated from the all-cause dis- tribution of cause-specific hazards for competing risks by using the method of Beyersmann et al. (2009, Statistics in Medicine 24: 956–971). A multinomial dis- tribution is used to create the event indicator, whereby the probability of expe- riencing each event at a simulated time t is the cause-specific hazard divided by the all-cause hazard evaluated at time t. Baseline covariates can be included in all scenarios. We also describe http://www.stata-journal.com/software/the extension to incorporate nonproportional hazards in standard parametric and competing-risks scenarios. Copyright 2012 by StataCorp LP.
