How to update PyMC model with counterfactuals?

So something more like this?