making cosmetic changes accroding to input from Florian:

1. Only outputting the posteriors for PIP > 0.0001
2. Additing count for plotting in parallel MJMCMC and GMJMCMC
3. Adding the default for r in gaussian loglik