This study proposes an efficient approach to detect one or more change points for gamma distribution. We plug a closed-form estimator into the gamma log-likelihood function to obtain a sharp approximation to the maximum of log-likelihood. We further derive a closed form calibration of approximate likelihood which is asymptotically equivalent to the exact log-likelihood. This circumvents iterative optimization procedures to find maximum likelihood estimates which can be a burden in detecting multiple change points. The simulation study shows that the approximation is accurate and the change points can be detected much faster. Two case studies on the time between events arising from industrial accidents are presented and extensively investigated.

Reliability data obtained from life tests and degradation tests have been extensively used for purposes such as estimating product reliability and predicting warranty costs. When there is more than one candidate model, an important task is to discriminate between the models. In the literature, the model discrimination was often treated as a hypothesis test and a pairwise model discrimination procedure was carried out. Because the null distribution of the test statistic is unavailable in most cases, the large sample approximation and the bootstrap were frequently used to find the acceptance region of the test. Although these two methods are asymptotically accurate, their performance in terms of size and power is not satisfactory in small sample size. To enhance the small-sample performance, we propose a new method to approximate the null distribution, which builds on the idea of generalized pivots. Conventionally, the generalized pivots were often used for interval estimation of a certain parameter or function of parameters in presence of nuisance parameters. In this study, we further extend the idea of generalized pivots to find the acceptance region of the model discrimination test. Through extensive simulations, we show that the proposed method performs better than the existing methods in discriminating between two lifetime distributions or two degradation models over a wide range of sample sizes. Two real examples are used to illustrate the proposed methods.