plot median(k) instead of mean(k) in pcnbd.mcmc.plotRegularityRateHeterogeneity