Rate-Based Daily Arrival Process Models with Application to Call Centers

We propose, develop, and compare new stochastic models for the daily arrival rate in a call center. Following standard practice, the day is divided into time periods of equal length (e.g., 15 or 30 minutes), the arrival rate is assumed random but constant in time in each period, and the arrivals are from a Poisson process, conditional on the rate. The random rate for each period is taken as a deterministic base rate (or expected rate) multiplied by a random busyness factor having mean 1. Models in which the busyness factors are independent across periods, or in which a common busyness factor applies to all periods, have been studied previously. But they are not sufficiently realistic. We examine alternative models for which the busyness factors have some form of dependence across periods. Maximum likelihood parameter estimation for these models is not easy, mainly because the arrival rates themselves are never observed. We develop specialized techniques to perform this estimation. We compare the goodness-of-fit of these models on arrival data from three call centers, both in-sample and out-of-sample. Our models can represent arrivals in many other types of systems as well. Estimating a model for the vector of counts (the number of arrivals in each period) is generally easier than for the vector of rates, because the counts can be observed, but a model for the rates is often more convenient and natural, e.g., for simulation. We examine and provide insight on the relationship between these two types of modeling. In particular, we give explicit formulas for the relationship between the correlation between rates and that between counts in two given periods, and for the variance and dispersion index in a given period. These formulas imply that for a given correlation between the rates, the correlation between the counts is much smaller in low traffic than in high traffic.

[1]  D. Pollock Smoothing with Cubic Splines , 1999 .

[2]  Peter Kall,et al.  Stochastic Programming , 1995 .

[3]  Ronald F. Boisvert,et al.  NIST Handbook of Mathematical Functions , 2010 .

[4]  W. Rudin Principles of mathematical analysis , 1964 .

[5]  Pierre L'Ecuyer,et al.  Efficient Correlation Matching for Fitting Discrete Multivariate Distributions with Arbitrary Marginals and Normal-Copula Dependence , 2009, INFORMS J. Comput..

[6]  R. Nelsen An Introduction to Copulas , 1998 .

[7]  David R. Cox,et al.  The statistical analysis of series of events , 1966 .

[8]  Louis Anthony Cox,et al.  Wiley encyclopedia of operations research and management science , 2011 .

[9]  Ger Koole,et al.  Call Center Optimization , 2013 .

[10]  Ward Whitt,et al.  Are Call Center and Hospital Arrivals Well Modeled by Nonhomogeneous Poisson Processes? , 2014, Manuf. Serv. Oper. Manag..

[11]  O. Tanir,et al.  Call center simulation in Bell Canada , 1999, WSC'99. 1999 Winter Simulation Conference Proceedings. 'Simulation - A Bridge to the Future' (Cat. No.99CH37038).

[12]  Pierre L'Ecuyer,et al.  A normal copula model for the arrival process in a call center , 2012, Int. Trans. Oper. Res..

[13]  Nabil Channouf,et al.  Modélisation et optimisation d'un centre d'appels téléphoniques : étude du processus d'arrivée , 2008 .

[14]  Seymour Geisser,et al.  8. Predictive Inference: An Introduction , 1995 .

[15]  Pierre L'Ecuyer,et al.  Staffing Multiskill Call Centers via Linear Programming and Simulation , 2008, Manag. Sci..

[16]  Haipeng Shen,et al.  Interday Forecasting and Intraday Updating of Call Center Arrivals , 2008, Manuf. Serv. Oper. Manag..

[17]  Itay Gurvich,et al.  Staffing Call Centers with Uncertain Demand Forecasts: A Chance-Constrained Optimization Approach , 2010, Manag. Sci..

[18]  Sujin Kim,et al.  The stochastic root-finding problem: Overview, solutions, and open questions , 2011, TOMC.

[19]  Ger Koole,et al.  Managing uncertainty in call centres using Poisson mixtures , 2001 .

[20]  Shane G. Henderson,et al.  FORECAST ERRORS IN SERVICE SYSTEMS , 2009, Probability in the Engineering and Informational Sciences.

[21]  Pierre L'Ecuyer,et al.  On the modeling and forecasting of call center arrivals , 2012, Proceedings Title: Proceedings of the 2012 Winter Simulation Conference (WSC).

[22]  Wolfgang Hörmann,et al.  Automatic Nonuniform Random Variate Generation , 2011 .

[23]  Shane G. Henderson,et al.  Optimizing Call Center Staffing Using Simulation and Analytic Center Cutting-Plane Methods , 2008, Manag. Sci..

[24]  van der,et al.  Proceedings of the 2012 winter simulation conference , 2001, WSC 2008.

[25]  Ward Whitt,et al.  An Introduction to Stochastic-Process Limits and their Application to Queues , 2002 .

[26]  Pierre L'Ecuyer,et al.  Inter-dependent, heterogeneous, and time-varying service-time distributions in call centers , 2016, Eur. J. Oper. Res..

[27]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[28]  W. Whitt,et al.  Improving Service by Informing Customers About Anticipated Delays , 1999 .

[29]  Tevfik Aktekin,et al.  Call center arrival modeling: A Bayesian state‐space approach , 2011 .

[30]  J. Michael Harrison,et al.  A Method for Staffing Large Call Centers Based on Stochastic Fluid Models , 2005, Manuf. Serv. Oper. Manag..

[31]  Pierre L'Ecuyer,et al.  Call-type dependence in multiskill call centers , 2013, Simul..

[32]  W. Whitt,et al.  Choosing arrival process models for service systems: Tests of a nonhomogeneous Poisson process , 2014 .

[33]  G. P. Bhattacharjee,et al.  The Incomplete Gamma Integral , 1970 .

[34]  Haipeng Shen,et al.  FORECASTING TIME SERIES OF INHOMOGENEOUS POISSON PROCESSES WITH APPLICATION TO CALL CENTER WORKFORCE MANAGEMENT , 2008, 0807.4071.

[35]  Pierre L'Ecuyer,et al.  Forecasting Call Center Arrivals: Fixed-Effects, Mixed-Effects, and Bivariate Models , 2012, Manuf. Serv. Oper. Manag..

[36]  Zeynep Akşin,et al.  The Modern Call Center: A Multi‐Disciplinary Perspective on Operations Management Research , 2007 .

[37]  Ward Whitt,et al.  Dynamic staffing in a telephone call center aiming to immediately answer all calls , 1999, Oper. Res. Lett..

[38]  Pierre L'Ecuyer,et al.  Modeling and simulation of call centers , 2005, Proceedings of the Winter Simulation Conference, 2005..

[39]  R. L. Iman,et al.  An iterative algorithm to produce a positive definite correlation matrix from an approximate correlation matrix (with a program user's guide) , 1982 .

[40]  Avishai Mandelbaum,et al.  Workload forecasting for a call center: Methodology and a case study , 2009, 1009.5741.

[41]  Jonathan Weinberg,et al.  Bayesian Forecasting of an Inhomogeneous Poisson Process With Applications to Call Center Data , 2007 .

[42]  Pierre L'Ecuyer,et al.  A Java library for simulating contact centers , 2005, Proceedings of the Winter Simulation Conference, 2005..

[43]  Avishai Mandelbaum,et al.  Telephone Call Centers: Tutorial, Review, and Research Prospects , 2003, Manuf. Serv. Oper. Manag..

[44]  Pierre L'Ecuyer,et al.  Staffing multi-skill call centers via search methods and a performance approximation , 2009 .

[45]  Avishai Mandelbaum,et al.  Statistical Analysis of a Telephone Call Center , 2005 .

[46]  Pierre L'Ecuyer,et al.  Modeling Daily Arrivals to a Telephone Call Center , 2003, Manag. Sci..

[47]  Pierre L'Ecuyer,et al.  Quasi-Monte Carlo methods with applications in finance , 2008, Finance Stochastics.

[48]  Haipeng Shen Statistical Analysis of Call‐Center Operational Data: Forecasting Call Arrivals, and Analyzing Customer Patience and Agent Service , 2010 .

[49]  Michael C. Fu,et al.  Guest editorial , 2003, TOMC.

[50]  C. Reinsch Smoothing by spline functions , 1967 .

[51]  Refik Soyer,et al.  Modeling and Analysis of Call Center Arrival Data: A Bayesian Approach , 2008, Manag. Sci..