add a function to simulate from beta distributions along a single gradient with either Poisson or neg binomial random counts