Shrinks towards the median value.
gen.beta(shape1, shape2, ncp = 0)
shape1
same as shape1 in rbeta
shape2
same as shape2 in rbeta
ncp
same as ncp in rbeta