After a brief chat with @fonnesbeck in Southampton, UK a few months back, I thought I would try using the
pm.Potential function to help teach MCMC in physics and it seems to work excellently.
I’d be keen to write up a notebook for the documentation if you think it would fit?
Here’s a notebook I produced for our group link. It would need adapting a little but I think an example of the Heisenberg model would be a nice example of how pm.Potential can be used.
I think the important concepts for those looking to use PyMC3 for physics problems would be:
- Showing how energy and the probability distribution link
- Explaining the necessity of the prior distributions
- Plotting the energy and probability surfaces
Any other suggestions?