Hello,

I am just wondering your help which R function or code should I use to model a three stage hierarchical model for a count data (a Poisson Process With Change Point)?

Is MCMC the best choice to estimate the time when the rate of an event changed? What should be specified in the MCMC model?

It will be great if I can get some worked examples too.

Regards