Tips:
- Test your posterior sampling code. See this reference.
- Use debugger. For example, you can
import pdband set break point using commandpdb.set_trace(). - Fix random seed before running your sampling code to be able to reproduce any bugs: i.e.,
np.random.seed(some_seed)