simPH function

Generate counts under a variant of a 'phenomenological model'