Different approaches to the Challenger Example

I think this is usually modelled using a Poisson point process. I have never fit one myself, but using the GP module you can do something similar: http://docs.pymc.io/notebooks/GP-slice-sampling.html#Gaussian-Process-Classification