Changed parameterisation of Negative Binomial: using phi = 1 / alpha, instead of alpha
Working with rate = 1 / scale in gamma_shape_scale function