How to update PyMC model with counterfactuals?

Yes, does that make sense in the context of your model?