@bob-carpenter I moved the thread to its own topic.
You can use pm.observe to define observed nodes in pre-existing model thse days. You can also do a lot of model mutation with pm.do, including setting parameters to specific values, taking prior predictive draws and then conditioning on those for calibration workflows.
Our README shows an example of such workflow: pymc/README.rst at main · pymc-devs/pymc · GitHub