variance
. This requires
that the original value is always positive. When using this, the simulators
are called separately for each locus, which can dramatically increase the
time needed to simulate models with many loci.par_variation(par, variance)