The change point is estimated as follows:
• Using consecutive event times, disjoint intervals are defined.
• Each point within such an interval can be considered to be a change point defining a piecewise Weibull NHPP model with two phases. So long as the two phases defined by that point each consist of at least two events, the algorithm can compute MLEs for the parameters of that model. The log-likelihood for that model can also be computed.
• Within each of the disjoint intervals, a constrained optimization routine is used to find a local optimum for the log-likelihood function.
• These local optima are compared, and the point corresponding to the largest is chosen as the estimated change point.
Note that this procedure differs from the grid-based approach described in Guo et al. (2010).