typically, dispersion = 1 / size

Typically dispersion of a negative binomial is 1/size for the size arg of the rnbinom function.

This way, for fixed mu, larger dispersion corresponds to larger variance.

If you accept this change you'd have to fix the man pages for NB and simulate_experiment.