Rstanarm

[ermeel]

In the meantime you can try to use a workaround with rstanarm or brms in case you are interested in semi-parametric piece-wise constant proportional hazard models. These you can cast into the form of a poisson regression, following the ideas outlined here

Using stan_jm for parametric proportional hazards regression only

This is an awesome idea! Just to make sure I understand correctly (and for future readers) let me try and formulate the idea in my own words (or mathematical notation): W.l.o.g. set the origin of time to 0. Suppose the event (either 0 or 1) for individual i has occurred at time t_i. Suppose there are n_i unique event times before t_i and denote them by t_1 < t_2 < \dots <t_{n_i}. For individual i we construct n_i+1 observations of the form (t_1, 0), ( t_2 - t_1, 0), (t_3-t_2, 0),…, (t_{n_{i}…