3. MLE change point estimator
Here, it is assumed that the underlying process initially
operates in a state of statistical control, with observations
coming from a Poisson distribution with the known parameters
β = β0 (β is a p-dimensional vector); so, the mass probability
function is: