powered by
Shrinks towards the median value.
gen.gamma(shape, rate = 1, scale = 1/rate)
same as shape in rgamma
same as rate in rgamma
same as scale in rgamma