I know this is a really old thread but can someone please help me figure out how I can use PyMC3 to do exactly as quoted?