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

We 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, we 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 our data which contained frequent cyclic abrupt changes in arrival rates. In this paper, we propose the use a piecewise polynomial to represent the rate function. We 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. We 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 our 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.

[1]  Gwilym M. Jenkins,et al.  Time series analysis, forecasting and control , 1972 .

[2]  Herbert P. Galliher,et al.  Nonstationary Queuing Probabilities for Landing Congestion of Aircraft , 1958 .

[3]  Edward P. C. Kao,et al.  Simulating a Computer System for Processing Electrocardiograms , 1980, IEEE Transactions on Systems, Man, and Cybernetics.

[4]  P. A. W. Lewis,et al.  Simulation of Nonhomogeneous Poisson Processes with Degree-Two Exponential Polynomial Rate Function , 1979, Oper. Res..

[5]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[6]  Paul Bratley,et al.  A guide to simulation , 1983 .

[7]  W. Fuller.  GRAFTED POLYNOMIALS AS APPROXIMATING FUNCTIONS , 1969 .

[8]  G. Shedler,et al.  Simulation of Nonhomogeneous Poisson Processes by Thinning , 1979 .

[9]  D. Cox,et al.  The statistical analysis of series of events , 1966 .

[10]  Sheldon M. Ross,et al.  Stochastic Processes , 2018, Gauge Integral Structures for Stochastic Calculus and Quantum Electrodynamics.

[11]  Paul Bratley,et al.  A guide to simulation (2nd ed.) , 1986 .

[12]  P. Lewis.  Recent results in the statistical analysis of univariate point processes , 1971 .

[13]  George A. F. Seber,et al.  Linear regression analysis , 1977 .

[14]  V. Klema.  LINPACK user's guide , 1980 .

[15]  L. Dickson.  First Course in the Theory of Equations , 1922, Nature.

[16]  P. Lewis.  Remarks on the theory, computation and application of the spectral analysis of series of events , 1970 .

[17]  C. Maclean,et al.  Estimation and testing of an exponential polynomial rate function within the nonstationary Poisson process , 1974 .

[18]  H. Lilliefors.  On the Kolmogorov-Smirnov Test for the Exponential Distribution with Mean Unknown , 1969 .

[19]  Peter A. W. Lewis,et al.  Statistical Analysis of Non-Stationary Series of Events in a Data Base System , 1976, IBM J. Res. Dev..

[20]  D. G. Beech,et al.  Handbook of Statistical Tables. , 1962 .

[21]  Linus Schrage,et al.  A guide to simulation , 1983 .

[22]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[23]  Bernard O. Koopman,et al.  Air-Terminal Queues under Time-Dependent Conditions , 1972, Oper. Res..