plotModel now returns the computed fitted model and its components
phi is now not allowed to become less than PHI_MIN in negbinom.rPrior
Global beta and phi now appear as red dots in the hyper-parameters plot