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