Simulates time-to-event (survival) datasets for clinical trial
design and analysis, including datasets driven by a joint
pharmacokinetic/pharmacodynamic (PK/PD) model solved via the
'mrgsolve' backend. Event times can be located directly during ODE
integration ('in-solver' event detection) using a library of
built-in PK/PD hazard models or a user-supplied 'mrgsolve' model
converted for this purpose, or generated by inverse transform
sampling from a closed-form Weibull or flexible M-spline baseline
hazard, or from any externally simulated survival trajectory. See
Bender et al. (2005)