PyMC version of fitting the robertson model