I did an update of this model for the latest PyMC here but never shared it. Now’s as good a time as any to do so!