## An Extreme Value Theory Approach for Understanding Queue Length Dynamics in Adaptive Corridors

Shakib Mustavee, 1 Pushkin Kachroo (Senior Member IEEE) 2 and Shaurya Agarwal (Senior Member IEEE) 3

Abstract -This paper introduces a novel approach employing extreme value theory to analyze queue lengths within a corridor controlled by adaptive controllers. We consider the maximum queue lengths of a signalized corridor consisting of nine intersections every two minutes, roughly equivalent to the cycle length. Our research shows that maximum queue lengths at all the intersections follow the extreme value distributions. To the best knowledge of the authors, this is the first attempt to characterize queue length time series using extreme value analysis. These findings are significant as they offer a mechanism to assess the extremity of queue lengths, thereby aiding in evaluating the effectiveness of the adaptive signal controllers and corridor management. Given that extreme queue lengths often precipitate spillover effects, this insight can be instrumental in preempting such scenarios.

## I. INTRODUCTION

Adaptive signal control technologies (ATSC) are critical for managing arterials as they can adjust signal timings according to highly variable traffic demands. ATSC adjusts the timing of red, yellow, and green lights to achieve specific operational goals, such as smooth traffic flow, maximizing throughput, ensuring access equity, and queue management. The effectiveness of adaptive strategies is determined by various measures of effectiveness (MOEs), including throughput capacity, traffic delay, and queue lengths [1], [2]. Among these MOEs, queue lengths are crucial in evaluating signalized intersections' performance. Traffic signal controllers such as InSync utilize artificial intelligence algorithms to optimize traffic flow along signalized arterials. These algorithms are based on cost functions that consider the wait times for each vehicle and queue lengths for each lane. By continuously gathering real-time information from sensors and cameras, the controllers adjust the traffic light's green duration to minimize queue lengths and improve traffic efficiency. It is important to understand the evolution of queue length, particularly cycle-wise maxima, which helps identify spillover conditions, i.e., when the queue length exceeds the capacity. In queue spillover scenarios, upstream traffic conditions are affected by downstream traffic due to the mutual interactions of neighboring intersections. Therefore, estimating the cycle-based maximum queue length is crucial

*This work was not supported by any organization

1 Shakib Mustavee is a Ph.D. candidate with the Department of Civil, Environmental &amp; Construction Engineering, University of Central Florida, Orlando, Florida. shakib.mustavee@ucf.edu

2 Pushkin Kachroo is a Professor with the Department of Electrical &amp; Computer Engineering, University of Nevada, Las Vegas. pushkin@unlv.edu.

3 Shaurya Agarwal is an Assistant Professor with the Department of Civil, Environmental &amp; Construction Engineering, University of Central Florida, Orlando, Florida. shaurya.agarwal@ucf.edu in optimizing signal control and managing congestion in signalized arterials. Long-term analysis of queue length provides insights into the performance and reliability of signal controllers. However, it is difficult to accurately estimate cycle-based queue length due to the lack of high-resolution data [3]. There is a significant gap in the literature regarding the long-term behavior of cycle-to-cycle queue length. This study examines queue length dynamics in a signalized corridor and aims to uncover the long-term behavior of maximum queue length by probabilistic analysis. We draw the maximum queue length of each cycle and characterize its distribution in the light of the extreme value theory (EVT), which is a statistical theory that deals with extreme values and events that occur at the tails of probability distributions. According to EVT, the largest value drawn from a sample of a given size can follow one of three possible distributions: Weibull, Fr´ echet, or Gumbel distributions. Since uncommon and extreme conditions can cause significantly greater consequences, understanding and quantifying rare occurrences is essential for reliability analysis and impact assessment. EVT aims to assess the likelihood of events surpassing previous extremes based on ordered samples of random variables, and it focuses on studying extreme deviations from distribution medians. Extreme value analysis is widely applicable across fields such as structural engineering, finance, economics, geology, telecommunication, and traffic safety. It is crucial for comprehending and preparing for extreme events across various domains. In transportation, the primary application of EVT remains in surrogate safety measures based on traffic safety analysis. Studying the queue length dynamics in the light of extreme value theory is intriguing as it provides insight into extreme incidents when queue lengths exceed capacity.

To our best knowledge, this is the first research paper that uses EVT to analyze queue length dynamics. The contributions of the paper are as follows:

- · The paper introduces extreme value analysis of queue length time series.
- · It provides mathematical justifications for the emergence of extreme value distribution from queueing theory.
- · It provides empirical evidence for the emergence of extreme value distribution using real data.

The structure of the paper is as follows. Section II discusses the relevant literature on EVT analysis in transportation engineering and general queueing systems. Section III provides background on mathematical tools, while

Section IV provides a rich discussion on data and obtained results.

## II. LITERATURE REVIEW

EVT has been used in traffic safety research for over two decades, particularly to analyze the association between traffic conflicts and crashes. A thorough analysis of the studies that have applied EVT in traffic safety can be found in [4]. The review paper offers elaborate discussions on assessing alternate modeling methodologies and associated issues. There are only a few examples of EVT being applied in transportation literature, apart from safety applications. One example is the use of an EVT-based approach to measure extremely prolonged travel time and assess potentially influential factors of emergency vehicles, as proposed in [5]. Another example is the use of EVT in traffic incident detection, as described in [6]. However, to the best of our knowledge, EVT-based approaches have not been used in queue length analysis. While transportation literature largely ignores extreme value analysis of queue length analysis, extensive research exists regarding the extreme value distributions of queue lengths and waiting time in various other queuing systems related to operations research and telecommunication [7], [8], [9]. These works primarily focus on deriving explicit equations and limiting cases for various queuing conditions.

## III. MATHEMATICAL BACKGROUND

## A. Generalized Extreme Value Distribution

EVT focuses on modeling the maximum of a sample of identical and independent random variables (i.i.d r.v). The theory says that after proper renormalization, their distribution can only converge to one of only 3 possible distribution families: type-I (Gumbel), type-II (Fr´ echet), and type-III (Weibull). The Fisher-Tippett-Gnedenko theorem (extreme value theorem) provides the criteria for the maximum value of i.i.d r.v) to converge to a Generalized Extreme Value (GEV) distribution characterized by three parameters: µ, σ, and ξ which combines the Gumbel, Fr´ echet, and Weibull families under one general distribution. It means that under different conditions, these distributions emerge from the GEV distribution.

We consider the set of i.i.d r.v:

$$X _ { 1 }, X _ { 2 }, \dots, X _ { n }$$

Here, we sample each X i from any function F . The distribution of the function is not known. Let Γ n = max( X ,X ,...,X 1 2 n ) .

Let there exists a set of real numbers α , β n n , where α n ≥ 0 , such that Γ n -β n α n has a limiting distribution in the limit of large n and it is non-degenerate. Then extreme value theorem says:

$$\Gamma _ { n } \sim G E V ( \mu, \sigma, \xi )$$

This means the maximum (and minimum) of the i.i.d r.v follows generalized extreme value distribution (GEV).

Generalized Extreme Value (GEV) Distribution: First, we discuss the cumulative distribution function (CDF) of the Generalized Extreme Value (GEV) distribution and then we see how it is related to the three different types of extreme value distributions. The CDF of GEV is:

$$\mathbb { Z } _ { 1 } ^ { \mathbb { Z } } \quad F ( x ; \mu, \sigma, \xi ) = \exp \left \{ - \left [ 1 + \xi \left ( \frac { x - \mu } { \sigma } \right ) \right ] ^ { - 1 / \xi } \right \},$$

where, µ is the location parameter, σ is the scale parameter, ξ is the shape parameter.

Type-I or Gumbel Extreme Value Distribution: The Gumbel distribution emerges when ξ = 0 in the GEV distribution:

$$F _ { \text{Gumbel} } ( x ; \mu, \sigma ) = \exp \left \{ - \exp \left [ - \left ( \frac { x - \mu } { \sigma } \right ) \right ] \right \},$$

where µ and σ are the location and scale parameters, respectively.

Type-II or Fr´ echet Extreme Value Distribution: The Fr´ echet distribution emerges when ξ &gt; 0 in the GEV distribution:

$$\text{cus} \quad F _ { \text{Fréchet} } ( x ; \mu, \sigma, \xi ) = \exp \left \{ - \left [ 1 + \xi \left ( \frac { x - \mu } { \sigma } \right ) \right ] ^ { - 1 / \xi } \right \},$$

where µ σ , , and ξ are the location, scale, and shape parameters, respectively.

Type-III or Extreme Value (Weibull) Distribution: The Weibull distribution emerges when ξ &lt; 0 in the GEV distribution:

$$\underset { \text{em} } { \text{ition} } \underset { \text{em} } { \text{-III} } \quad F _ { \text{Weibull} } ( x ; \mu, \sigma, \xi ) = \exp \left \{ - \left [ 1 + \xi \left ( \frac { x - \mu } { \sigma } \right ) \right ] ^ { - 1 / \xi } \right \},$$

where µ σ , , and ξ are the location, scale, and shape parameters, respectively.

There are two methods for modeling extreme values within a random process. The first approach, the block minima/maxima, involves dividing the data into non-overlapping, continuous, and uniform intervals. Within each interval, either the maximum or the minimum entry is recorded. These minima/maxima within the intervals follows a generalized extreme value (GEV) distribution, as described by Gumbel in 1958 [10]. The second approach, referred to as threshold exceedances or the peak over the threshold approach, involves selecting observations surpassing/falling short of a threshold value. Using this method, the observations chosen are typically approximated using a generalized Pareto distribution (GPD), as outlined by Pickands III in 1975 [11].

While the block maxima approach is generally easier to apply and interpret, it may not always be directly applicable and often leads to inefficient use of data. Conversely, threshold exceedance can potentially lead to more efficient data utilization, but it is more challenging to implement and ensure that the underlying conditions of the theory are met.

## B. Queuing Process of a Signalized intersection

For several decades, signalized intersections have been modeled as a queueing system. Queueing models can estimate the average queue length and waiting time in a signalized intersection. Different models, such as M/M/1, M/G/1, GI/M/1, and M/D/1 have been proposed to model various attributes of signalized intersections in different scenarios. For example, Pollaczek-Khintchine equation is widely used in transportation literature to model average queue length and delay assuming intersections as an M/G/1 system [12]. A review of various queueing models used in modeling intersections is provided in [13]. Traditional queueing models can not fully capture the features of a signalized intersection due to the alteration of signals. During the red time, signalized intersections suspend service. The queuing process of vehicles at signals can be characterized by queueing models with server vacations [14]. Despite their limitations, these models provide a workable approximation of real queue lengths and delays at signalized intersections.

We use queueing theory to explain the mathematical structure of cycle maxima. Let's assume a signalized intersection with a stable queueing formation and appropriate initial distribution. Here, W i represents the waiting time of the i th vehicle, and L t ( ) denotes the queue length at time t . Here, the objective is to identify approximate distributions for the maximum queue lengths and waiting associated with this process. For sufficiently large values of i and t we assume

$$W _ { i } ^ { * } = \max \{ W _ { k } \colon 0 \leq k \leq i \}, \ \ i \geq 0,$$

and

$$L ^ { * } ( t ) = \max \{ L ( s ) \, \colon 0 \leq s \leq t \}, \ \ t \geq 0,$$

Now, according to extreme-value limit theorems, the maximum queue length process L ∗ ( ) t , that the approximations should follow the mathematical structure [8]:

$$L ^ { * } ( t ) \sim \gamma ( \log t + \log \beta + Z ),$$

where t is presumed to be relatively large, log denotes the natural logarithm, Z follows the Gumbel cumulative distribution function (type-I extreme-value CDF):

$$P ( Z \leq x ) \equiv \Lambda ( x ) = \exp ( - e ^ { - x } ), \ \ - \infty < x < \infty, \quad \stackrel { \alpha } { 2 }$$

and γ and β are positive constants that depend on the particular process being analyzed.

$$I V. \, C A S E \, S T U D Y$$

## A. Data Description

The dataset comprises queue lengths collected from InSync adaptive traffic signals recorded between December 18, 2017, and February 14, 2018, at the Alafaya Trail (SR-434) corridor in East Orlando, FL, shown in Figure 1. This study focuses on queue lengths at nine signalized intersections along this corridor, primarily investigating northward movements (NT) since it is the major road. Prior research demonstrated the effectiveness of adaptive signals for major roads [2]. Traffic sensors may not always produce accurate data due to various factors, such as detector malfunctions, errors during data storage, and unfavorable weather conditions. The study followed the outlier mitigation and data noise reduction methods suggested in [15] and [16].

## B. Data Processing

The extreme value theory assumes i.i.d r.v (independent and identically distributed random variables), meaning that the random variables are- (i) independent: the samples are independent. (ii) identically distributed: each sample comes from the same distribution. EVT often assumes stationarity of the underlying stochastic process. However, the traffic flow dynamics exhibit periodic behavior imposed by peak hours, circadian and weekly cycles. For example, the hybrid model proposed in [17] for traffic flow prediction decomposes time series into four different partsa periodic component, a trend component, a stationary component, and a volatility component. In [16], the queue length time series was decomposed into quasiperiodic and chaotic components. Therefore, seasonality from the traffic flow time series must be adjusted before applying EVT.

In signal processing, a time series is often decomposed as:

$$X _ { t } = S _ { t } + T _ { t } + C _ { t } + R _ { t }$$

where X t represents the data at time t ,

Here, four key components of time series comprise are:

- -( S t ) signifies the periodicity or the seasonal component
- -( T t ) is the gradual changes or the trend component
- -( C t ) denotes the cyclical component
- -( E t ) indicates the residual part of the time series

Distinguishing seasonal and cyclic patterns: Seasonal patterns have a fixed length, while cyclic patterns are variable. Cyclic patterns occur when data fluctuates unpredictably. The average cycle is longer than seasonality, and the magnitude fluctuations are more variable.

Fig. 1. Location of intersections on Alafaya corridor (left) and illustration of the corridor as a dynamical system (right). Northbound (NB) through queue length at each intersection is measured in distance units.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-5ae57a84-6b10-478a-ade5-9d2e007efa7d.jpg)

We observe the distribution of residual components of the queue length time series. Together, seasonality and trend are often known as autoregressive behavior. To obtain residual components, we remove autoregressive behavior from

Fig. 2. Illustrates the original queue length time series and adjusted time series after removing seasonal and trend components from the original queue length time series.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-64bd2e8f-a1e0-4a83-8794-d32c5d77e0ac.jpg)

the time series using seasonal decompose, a Python-based software [18]. It is important to note that the queue length time series dynamics have multiple driving frequencies [16]. The periodic component of the queue length time series has several frequencies. Therefore, removing the seasonality of queue length is not straightforward. Traffic flow dynamics are driven by a circadian circle. Hence, it is a common practice to assume the seasonality period of traffic flow as 24 hours. However, during late night and early morning offpeak hours (10 PM to 07 AM), the queue length usually remains zero. Therefore, we assume the period of seasonality in queue length time series as 15 hours. Figure 2 illustrates the original queue length time series (shown in blue) along with the seasonal and trend components. The seasonal and trend components are added together (shown in red). Figure 3 shows the residual part of the time series.

## C. Distribution of Residual Component

It is observed that the residual components do not follow normal distribution. We tried to find the best-fit distribution for the residual part and used distfit , a Python package, to find the best-fit distribution [19]. The package can test a series of distributions fitted on the empirical data and calculate the residual sum of squares (RSS) for each distribution. RSS quantifies the variance between predicted and observed values in empirical data. Essentially, it highlights the disparities in estimations. It serves as a metric of how well a model aligns with the actual data. A lower RSS suggests a closer model alignment with the data points. The RSS calculation is represented by the formula:

$$R S S = \sum _ { i = 1 } ^ { n } ( y _ { i } - f ( x _ { i } ) ) ^ { 2 }$$

Fig. 3. Illustrates the residual part computed after removing seasonal and trend components from the original queue length time series.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-b78cea57-116b-4f7e-9fc0-7a7ae0babbbe.jpg)

Here, y i denotes the i th value of the target variable, x i is the i th value of the data, and f ( x i ) denotes the predicted value of y i . Here, the best-fitting distribution is determined based on the lowest RSS score.

In the analysis, we tested a total of 79 different distributions for each intersection, and Table I summarizes the results. The table shows the ranking of the best-fitted distributions based on their RSS scores. Each column lists the best-fit distributions for each intersection, with the corresponding RSS score shown in parentheses. Note that we used acronyms for distribution names in the table. The full form can be found in [19].

We observed that except for Lokanotosa Tr, the distributions of the rest of the intersections could be described by any of the extreme value distributions such as genextreme (generalized extreme value distribution), gumbel r (Gumbel right), gumbel l (Gumbel left), weibull max (Weibull maximum), invweibull (inverse Weibull). This suggests that the distribution of cycle-wise maximum queue length of the corridor follows extreme value distributions.

Figure 4 presents a visual representation of the RSS score of various commonly used distributions, including generalized extreme value, lognormal, t, loggamma, beta, exponential, Pareto, normal, double Weibull, and uniform. Upon examining the figure, we observe that the generalized extreme value distribution has the best RSS score in five intersections and the second-best RSS score in two intersections when compared to the other popular distributions. In the remaining intersections, the RSS score of the generalized extreme value distribution is comparable to the fitted distributions with the best scores.

The histogram of the residual is displayed in Figure 5. We fitted two distributions to the histogram: Gumbel and Fr´ echet,

|   Rank | Challenger Pkwy         | Lokanotosa Tr              | Research Pkwy            | Central FL Blvd        | University Blvd        | Strategy Blvd                         | Centaurus Dr           | Corporate Blvd         | McCullough Rd            |
|--------|-------------------------|----------------------------|--------------------------|------------------------|------------------------|---------------------------------------|------------------------|------------------------|--------------------------|
|      1 | genextreme (0.002835)   | vonmises line (0.00782)    | Fatigue life (0.003468)  | johnsonsu (0.003118)   | mielke (0.00516)       | beta (0.03804)                        | johnsonsu (0.003733)   | johnsonsu (0.001046)   | lognorm (0.010334)       |
|      2 | invweibull (0.002835)   | truncnorm, norm (0.008892) | recipinvgauss (0.003489) | burr (0.003422)        | alpha (0.005407)       | chi, nakagami (0.038059)              | exponnorm (0.005357)   | exponnorm (0.004448)   | genlogistic (0.010349)   |
|      3 | powerlognorm (0.002977) | loggamma (0.009469)        | lognorm (0.003503)       | mielke (0.003619)      | fisk (0.006609)        | gengamma (0.038083)                   | mielke (0.006891)      | mielke (0.004875)      | gumbel r (0.010366)      |
|      4 | invgamma (0.002985)     | rdist (0.009515)           | johnsonsb (0.003512)     | exponnorm (0.005025)   | genextreme (0.00697)   | pearson3, chi2, gamma (0.038166)      | burr (0.007078)        | t (0.009079)           | betaprime (0.010409)     |
|      5 | weibull max (0.003162)  | nakagami (0.009756)        | betaprime (0.003621)     | genextreme (0.010213)  | exponnorm (0.007005)   | betaprime (0.038412)                  | alpha (0.01051)        | fisk (0.009129)        | johnsonsb (0.010677)     |
|      6 | gumbel r (0.003188)     | chi (0.009756)             | f (0.003626)             | invweibull (0.010213)  | genlogistic (0.007018) | f (0.03888)                           | genextreme (0.010805)  | loglaplace (0.00917)   | recipinvgauss (0.010912) |
|      7 | genlogistic (0.003227)  | foldnorm (0.009906)        | invgamma (0.003634)      | genlogistic (0.010808) | gumbel r (0.007145)    | recipinvgauss, fatiguelife (0.038965) | invweibull (0.010805)  | tukeylambda (0.009183) | fatiguelife (0.010938)   |
|      8 | johnsonsu (0.003352)    | logistic (0.010683)        | genextreme (0.003854)    | gumbel r (0.011343)    | weibull max (0.007481) | invgauss (0.039039)                   | exponweib (0.011539)   | alpha (0.01008)        | invgauss (0.010942)      |
|      9 | lognorm (0.003538)      | t (0.01083)                | invweibull (0.003854)    | betaprime (0.013841)   | burr (0.007718)        | powerlognorm (0.039338)               | genlogistic (0.011597) | invweibull (0.010131)  | invgamma (0.011627)      |
|     10 | johnsonsb (0.003822)    | dweibull (0.010993)        | gumbel r (0.003927)      | invgamma (0.013917)    | invgamma (0.008051)    | gumbel r, weibull max, (0.039361)     | gumbel r (0.011838)    | genextreme (0.010131)  | powerlognorm (0.011706)  |

## TABLE I

## RANKING OF BEST-FIT DISTRIBUTIONS BASED ON RSS ERROR

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-21303778-f4be-40d5-b404-5d2eaa945147.jpg)

Fig. 4. Illustrates the RSS score of some popular distributions. The general extreme value distribution is the best fit for most intersections.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-4552a53d-f609-48ed-a391-eb2a25d2d64c.jpg)

which are special cases of the generalized extreme value distribution. Visually, both distributions fit the histogram well. However, to further verify this, we applied a few graphical tests.

The Probability-Probability (P-P) plots and QuantileQuantile (Q-Q) plots are two graphical tools used in data analysis to assess the similarity between two distributions. The P-P plot is utilized to evaluate quality of match between a given data set and a particular probability distribution. It compares the empirical CDF of the data with the hypothesized actual CDF. If the scatter plot of the two distributions exhibit a roughly linear pattern, it suggests that the hypothesized true distribution is a good fit for the data set being analyzed. On the contrary, the Q-Q plot is a graph used to

Fig. 5. Illustrates the probability distribution of the residual part of the queue length time series. Gumbel and Fr´ echet distributions are fitted into the histogram

identify the true probability distribution of a given set of data in a similar way to P-P plots. However, this graph compares the quantiles of the empirical distribution of the data with the quantiles of the assumed true probability density function. If the resulting graph of these two distributions follows a linear pattern, it suggests that the assumed probability density function reasonably matches the given data.

Figure 6 and Figure 7 present the P-P and Q-Q plots for all intersections respectively. These plots assume that the residual follows the Gumbel distribution. As per the figures, it is evident that all intersections follow a linear pattern in their respective P-P and Q-Q plots. This indicates that using the Gumbel distribution is a suitable way to approximate the

Fig. 6. P-P plots of the data assuming it follows Gumbel distribution

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-b9f8120e-9d09-4e24-b879-1a1baa453112.jpg)

queue length time series of the corridor.

## V. CONCLUSION

The paper discusses the relationship between the cycleto-cycle maximum queue length and extreme value theory. Characterizing long-term queue length patterns as GEV could lead to potential applications in arterial management. However, this paper only serves as a first step in utilizing EVT for queue length prediction and corridor management. In future studies, we plan to demonstrate the practical applications of this idea. Additionally, the paper only used the block maxima approach to characterize the extreme behavior of queue lengths. In future work, we plan to show the threshold exceedance analysis.

## REFERENCES

- [1] M. Essa and T. Sayed, 'Self-learning adaptive traffic signal control for real-time safety optimization,' Accident Analysis &amp; Prevention , vol. 146, p. 105713, 2020.
- [2] M. S. I. Shafik, 'Field evaluation of insync adaptive traffic signal control system in multiple environments using multiple approaches,' 2017.
- [3] J. Yao and K. Tang, 'Cycle-based queue length estimation considering spillover conditions based on low-resolution point detector data,' Transportation research part C: emerging technologies , vol. 109, pp. 1-18, 2019.
- [4] Y. Ali, M. M. Haque, and F. Mannering, 'Assessing traffic conflict/crash relationships with extreme value theory: Recent developments and future directions for connected and autonomous vehicle and highway safety research,' Analytic methods in accident research , p. 100276, 2023.
- [5] Z. Zhang, Q. He, J. Gou, and X. Li, 'Analyzing travel time reliability and its influential factors of emergency vehicles with generalized extreme value theory,' Journal of Intelligent Transportation Systems , vol. 23, no. 1, pp. 1-11, 2019.
- [6] S. Tuli, M. R. Wilkinson, and C. Kettell, 'Radnet: Incident prediction in spatio-temporal road graph networks using traffic forecasting,' arXiv preprint arXiv:2206.05602 , 2022.
- [7] S. Asmussen, 'Extreme value theory for queues via cycle maxima,' Extremes , vol. 1, no. 2, pp. 137-168, 1998.
- [8] A. W. Berger and W. Whitt, 'Maximum values in queueing processes,' Probability in the Engineering and Informational Sciences , vol. 9, no. 3, pp. 375-409, 1995.
- [9] K. Akbash, N. Doronina, and I. Matsak, 'On extreme values of the queue length in some queuing systems,' Georgian Mathematical Journal , vol. 31, no. 1, pp. 1-16, 2024.
- [10] E. J. Gumbel, Statistics of extremes . Columbia university press, 1958.
- [11] J. Pickands III, 'Statistical inference using extreme order statistics,' the Annals of Statistics , pp. 119-131, 1975.
- [12] N. ROUPHAIL15, A. TARKO16, and J. LI17, 'Traffic flow at signalized intersections,' 1992.
- [13] F. Gunes, S. Bayrakli, and A. H. Zaim, 'Smart cities and data analytics for intelligent transportation systems: an analytical model for scheduling phases and traffic lights at signalized intersections,' Applied Sciences , vol. 11, no. 15, p. 6816, 2021.
- [14] Q. Yang and Z. Shi, 'The evolution process of queues at signalized intersections under batch arrivals,' Physica A: Statistical Mechanics and its Applications , vol. 505, pp. 413-425, 2018.
- [15] R. Rahman and S. Hasan, 'Real-time signal queue length prediction using long short-term memory neural network,' Neural Computing and Applications , vol. 33, pp. 3311-3324, 2021.
- [16] S. Das, S. Mustavee, S. Agarwal, and S. Hasan, 'Koopman-theoretic modeling of quasiperiodically driven systems: Example of signalized traffic corridor,' IEEE Transactions on Systems, Man, and Cybernetics: Systems , 2023.
- [17] H. Zhang, X. Wang, J. Cao, M. Tang, and Y. Guo, 'A hybrid shortterm traffic flow forecasting model based on time series multifractal characteristics,' Applied Intelligence , vol. 48, pp. 2429-2440, 2018.
- [18] 'Statsmodels: statistical modeling and econometrics in Python.' [Online]. Available: https://github.com/statsmodels/statsmodels
- [19] E. Taskesen, 'distfit is a python library for probability density fitting.' Jan. 2020. [Online]. Available: https://erdogant.github.io/distfit

Fig. 7. Q-Q plots of the data assuming it follows Gumbel distribution

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00867v1-figure-279331b8-60e8-45a1-afc4-aec186fc4cf7.jpg)
Shakib Mustavee Pushkin Kachroo Shaurya Agarwal
An Extreme Value Theory Approach for Understanding Queue Length Dynamics in Adaptive Corridors
This paper introduces a novel approach employing extreme value theory to
analyze queue lengths within a corridor controlled by adaptive controllers. We
consider the maximum queue lengths of a signalized corridor consisting of nine
intersections every two minutes, roughly equivalent to the cycle length. Our
research shows that maximum queue lengths at all the intersections follow the
extreme value distributions. To the best knowledge of the authors, this is the
first attempt to characterize queue length time series using extreme value
analysis. These findings are significant as they offer a mechanism to assess
the extremity of queue lengths, thereby aiding in evaluating the effectiveness
of the adaptive signal controllers and corridor management. Given that extreme
queue lengths often precipitate spillover effects, this insight can be
instrumental in preempting such scenarios.