A paper came out recently with a novel accelerated lifetime (AFT) model with cubic splines. This should pique your interest for a few reasons:

1. It helps dethrone the Proportional Hazard (PH) model as the default survival model. People like the PH model because it doesn't make any distributional assumptions. However, like a Trojan horse, there are very strong implicit assumptions that are inherited, often which are too restricting. Suffice to say, I am not a big proponent of the PH model. 

2. A spline-based AFT model weakens the CPH's throne because the model can fit to a larger space of potential models. For example, the Weibull AFT model is a special case of this new AFT model. The authors of the paper also carefully demonstrate that it often has lower bias, standard error, or AIC than other popular AFT models (Generalized Gamma, Generalized F). 

3. AFT models are just simpler to explain. Coefficients of AFT models have a much nicer interpretation than PH or Proportional Odds (PO) models. Simply: a positive (negative) coefficient multiplicatively accelerates (decelerates) a subject's time-to-event. So, a coefficient of 2 means that a subject experiences the event twice as fast as a baseline subject, on average. 

I'm too lazy to give the mathematical details of the model (just drank a strong beer), but what I do want to mention is that the model is implementable in lifelines using our custom model syntax. Here's the code

Update: it's now part of lifelines as `lifelines.CRCSplineFitter`. Happy coding!