Chris Jackson

There's no way of doing this using `flexsurvspline`, I think, as that fits a cubic spline basis. At the moment you'd have to use `flexsurvreg` and construct a custom piecewise...

For putting covariates in, I'd consider parameterising as `(bhaz, hr2, hr3, ...)` where `bhaz` is the hazard in period 1, and `hr2`, `hr3` are the hazard ratios between period 2,...

Does this do what you want, or do I misunderstand? This gives a sample of 10 alternative values for the 365-day survival probability, based on 10 random draws from the...

OK - I think it would be a worthwhile enhancement to summary.flexsurvreg to have the option to return a sample of B values for any required summary statistic, rather than...

@t-patton The code for `pmatrix.fs` actually does this sampling if you request a confidence interval. It calls the function `bootci.fmsm`, which draws a sample, calculates the P matrix (or some...

Thanks for this report - yes, one would expect there to be equivalent information in the messy and tidy results. For the moment I will at least fix the documentation...

Interesting examples! In the first one, the MLE of the Weibull scale parameter (using scale in the `dweibull` sense not the `survreg` sense) in this example is `exp(4e+05)`, which is...

The error happens when trying to find initial values for the `flexsurvreg` fit. It does this by calling `survreg(..., dist="lognormal")` on the natural-scale survival times. This results in the `invalid...

I agree this would be a beneficial feature, if it were implemented in a clean modular way. For advanced users who want to might try different optimisation algorithms, it's already...