pcp function

Fit a Poisson cluster process