NamastePymc3 is an update to my previous Namaste code using the "exoplanet" models of Dan Foreman-Mackey et al. This uses a PyMC3/Theano back-end, a differentiable Keplerian orbital model, and the Gaussian process modeling of Celerite to fit exoplanetary transits.
Here we simply provide a pre-written "exoplanet" model specific to single transits, and plot/analyse the results.
See Example.ipynb for an example of how to use it. If that doesn't render, try this nbviewer link.
NamastePymc3 natively includes eccentricity and fits for period (rather than scaled transit velocity as before). The period prior can be applied according to your prefered index - P^-8/3 (as suggested by Kipping et al ), or the less-steep P^-5/3 (as used by Osborn et al 2016 and Sandford et al 2019)