ENH: make model.plot more complete, add example

This PR adds an example comparing DisMod-MR fits of simulated data for PD with a range of different parameters, and also includes some additions to the code that complement this example.