perturbate_poisson function

Generate perturbated counts by the Poisson distribution