Simulate Bespoke Time-to-Event Models Using ODEs

Simulates time-to-event (survival) datasets for clinical trial design and analysis using ordinary differential equation (ODE) models solved via the 'mrgsolve' backend. Built-in Weibull and flexible M-spline baseline hazard models are provided out of the box, and fully bespoke hazard models can be implemented as custom 'mrgsolve' ODE systems. Event times are generated by inverse transform sampling from the resulting cumulative hazard functions. See Bender et al. (2005) for the inverse transform sampling methodology and Royston and Parmar (2002) for flexible parametric survival models.