First, specify the two parameters,
shape1 and shape2, and the
desired population mean and standard
deviation. The random numbers, drawn
from a beta distribution by
stats::rbeta() will then be
rescaled with the desired population
mean and standard deviation.