Modeling time-dependent arrivals to service systems: A case in using a piecewise-polynomial rate function in a nonhomogeneous Poisson process

The authors consider the use of a nonhomogeneous Poisson process in modeling time-dependent arrivals to service systems. In analyzing a set of actual arrival times corresponding to epochs of calls for on-line analysis of electrocardiograms, they found that approximating the rate function by an exponential polynomial (or exponential trigonometric polynomial) suggested in the literature inadequate. This was caused by the nature of the present data which contained frequent cyclic abrupt changes in arrival rates. In this paper, the authors propose the use of a piecewise polynomial to represent the rate function. They present two maximum likelihood estimators for estimating the parameters of the piecewise polynomial-one based on arrival times, and another based on aggregated counts, and a numerical method for carrying out the computation. The authors use a procedure based on thinning for generating arrival times from such a process. For hypothesis testing, these are combined to produce critical values for the Kolmogorov Smirnov statistic by a Monte Carlo simulation. While the present approach is presented and discussed in the context of a specific case, the results are applicable and observations relevant to many other systems sharing similar patterns of variation.