Would be better to do some more specific unit tests rather than just testing expection See suggestions in https://github.com/UCL/causalprog/pull/42#discussion_r2251222506 - [ ] Experiment with making a new graph with `do()` and calling numpyro's `do` - [ ] #91