Next Article in Journal
Application of Adaptive MOMEDA with Iterative Autocorrelation to Enhance Weak Features of Hoist Bearings
Next Article in Special Issue
Causal Information Rate
Previous Article in Journal
The Impact of Linear Filter Preprocessing in the Interpretation of Permutation Entropy
Previous Article in Special Issue
Information Geometric Theory in the Prediction of Abrupt Changes in System Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extreme Value Theory in Application to Delivery Delays

by
Marcin Fałdziński
1,
Magdalena Osińska
2,* and
Wojciech Zalewski
3
1
Department of Econometrics and Statistics, Nicolaus Copernicus University, Gagarina 11, 87-100 Toruń, Poland
2
Department of Economics, Nicolaus Copernicus University, Gagarina 11, 87-100 Toruń, Poland
3
Department of Logistics, Nicolaus Copernicus University, Gagarina 11, 87-100 Toruń, Poland
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(7), 788; https://doi.org/10.3390/e23070788
Submission received: 5 May 2021 / Revised: 18 June 2021 / Accepted: 19 June 2021 / Published: 22 June 2021

Abstract

:
This paper uses the Extreme Value Theory (EVT) to model the rare events that appear as delivery delays in road transport. Transport delivery delays occur stochastically. Therefore, modeling such events should be done using appropriate tools due to the economic consequences of these extreme events. Additionally, we provide the estimates of the extremal index and the return level with the confidence interval to describe the clustering behavior of rare events in deliveries. The Generalized Extreme Value Distribution (GEV) parameters are estimated using the maximum likelihood method and the penalized maximum likelihood method for better small-sample properties. The findings demonstrate the advantages of EVT-based prediction and its readiness for application.

1. Introduction

The Extreme Value Theory (EVT) evaluates both the magnitude and frequency of rare events and supports long-term forecasting. Delivery delays in transportation occur stochastically and rarely. However, they always increase costs, decrease consumers’ satisfaction, and lower the confidence of subcontractors in the supply chain. Thus, this paper refers to information-based vehicle planning and route analysis in road freight transportation. Telematics systems track the vehicles and send data to integrated transport management systems (TMS). They help to anticipate extreme time delays when both the transport operator and subcontractors have delivered a fleet.
This study aims to model and forecast rare events that appear as delivery delays in road transport using the tools provided by the EVT. The detailed research questions refer to the following issues:
RQ1: Is the Generalized Extreme Values (GEV) distribution able to model delivery delays data properly?
RQ2: What are the implications of estimated return level and extremal index for delivery planning?
Events that are rare in frequency and massive in magnitude have become important due to their negative impact. The Extreme Value Theory constitutes a framework allowing one to describe and measure such events. The EVT has had a vast array of applications. It covers, among many others, environmental data [1,2,3,4,5,6,7,8,9], risk management [10,11,12,13,14,15,16,17,18], but also earthquake size distribution [19] and athletic records [20].
In the literature on logistics and supply chain management, the Extreme Value Theory is rarely applied. One of the reasons is that logistics means the right organization. Supply chain management provides techniques to provide adequate services at every stage of the supply chain, including delivery. The focus is on optimizing activities to ensure timeliness, the security of goods, and service quality throughout the supply chain. Although quality management is currently highly developed and carefully evaluated, the risk of delivery delays remains an essential factor for the entire supply chain’s instability.
In recent publications on transportation, one can find several contributions to the theory and practice of the EVT. Firstly, the extreme value theory has been widely used in the case of safety on road analysis. Zheng and Sayed [21] proposed using the extreme value theory approach to conduct traffic conflict-based before-after analysis. The capability of providing a confident estimation of extreme events by the EVT approach drives the before-after analysis to shift from normal traffic conflicts to more severe conflicts, which are relatively rare but have more in common with actual crashes. The generalized Pareto distributions of traffic conflicts with a time-to-collision (TTC) of less than 0.5 s were applied. The results show an approximately 34% reduction in total extreme-serious conflicts in the analyzed area. Orsini et al. [22] used the EVT to predict the annual number of crashes in motorway cross-sections. Two traditional approaches, block maxima and peak-over-threshold, were applied to estimate EVT parameters for each cross-section. The authors found that both methods produced reliable predictions of annual rear-end collisions. Secondly, travel delays were the subject of the EVT. Esfeh et al. [23] proposed a new class of extreme value distribution called compound generalized extreme value (CGEV) distribution for investigating the effects of monthly and seasonal variation on extreme travel delays in road networks. They assumed that the frequency and severity of extreme events are highly correlated to the variation in weather conditions. The change in driving behavior varies according to road/weather conditions and contributes to the monthly and seasonal variation in observed extreme travel times.
Moreover, time delays in road transportation are typically considered from the perspective of congestion in the urban area within a temporal and spatial characteristics framework [24,25], intelligent transport systems [26], motorway traffic control [27], and the deviation between the actual and estimated shipment demand of containers [28].
The paper allows filling the gap between the usefulness of the EVT and its application in delivery operations being a part of a supply chain.
The methodology is applied based on real-time delay data obtained from a large brewery company operating in a European trademark network, which contracted a road transport company to provide delivery services. The data comes from the integrated transport management system.
The remaining of the paper is as follows. Section 2 introduces the methodology based on the extreme value theory. Section 3 details the data characteristics and empirical results. Section 4 contains discussion, and Section 5 concludes.

2. Extreme Value Theory—A Methodological Background

When modeling the maxima of a random variable, extreme value theory serves the same role as the central limit theorem plays when modeling sums of random variables. In both cases, the theory indicates limiting distributions. The Fisher and Tippett [29] theorem states that rescaled sample maxima converge in distribution to a variable having a distribution of Generalized Extreme Value (GEV) distribution:
H γ , μ , σ ( z ) = { e x p { ( 1 + γ z μ σ ) 1 / γ }   γ 0 e x p { e x p ( ( z μ σ ) ) }   γ = 0
where 1 + γ ( z μ ) / σ > 0 , and < μ < , σ > 0 and < μ < .
The most common approach for describing extreme events for stationary data is the block maximum method, which models the maxima of a set of contiguous blocks of observations using GEV distribution. The blocks method assumes a given sample, X 1 , , X N , divided into m blocks of length n , where the whole sample is N = n m . Then, the j-th block comprises X ( j 1 ) n + 1 , , X j n . For each block, the maximum M n , j = m a x ( X ( j 1 ) n + 1 , , X j n ) is calculated. The next step is to fit GEV, H γ , μ , σ , distribution to the block maxima, M n , 1 , , M n , m , that generates estimates for the parameters ( γ ^ , μ ^ , σ ^ ) is short H ^ [30]. Note that GEV comprises the following three subclasses: for the parameter γ > 0 named Fréchet (the distribution is heavy-tailed), for γ < 0 and for γ = 0 named Weibull. The theoretical return level when γ = 0 for the GEV distribution indicates how the return level would look when the γ parameter equals zero, and all other coefficients remain constant.
The likelihood function for independent maxima Z 1 , , Z m , is given:
L ( γ , μ , σ ) = m log σ ( 1 + 1 γ ) i = 1 m log [ 1 + γ ( Z i μ σ ) ] i = 1 m [ 1 + γ ( Z i μ σ ) ] 1 / γ
provided that 1 + γ ( z μ ) / σ > 0 , for i = 1 , , m .
Coles and Dixon [31] proposed a penalty function for the maximum likelihood method for the GEV distribution. The penalty is in terms of both bias and variance magnitudes. The penalty function is of the form:
P ( γ ) = { 1 if   γ 0 exp { λ ( 1 1 γ 1 ) β } if   0 < γ < 1 0 if   γ 1
for a range of non-negative values of λ and β . Therefore, the penalized likelihood function is as follows:
L p e n ( γ , μ , σ ) = L ( γ , μ , σ ) × P ( γ )
Coles and Dixon [31] advise setting λ = β = 1 , which leads to reasonable performance. The Penalized Maximum Likelihood (PML) estimator appears to be slightly better than, or at least as good as, the PWM estimator. Following Coles and Dixon [31], the two approaches are applied for comparison, i.e., the penalized maximum likelihood (PML) and the maximum likelihood (ML) estimators.
The selection of the blocks corresponds to a period of a specific time length (e.g., one year). In the case of the block method, an alternative type of quantile is of interest. Generally, a return level R n , T with a return period T is the quantile of the maxima distribution M n , which satisfies P ( M n > R n , T ) = 1 / T [32]. The return level R n , T is expected to exceed on average once every 1 / T given period (e.g., for annual data 1 / p years) with probability p . The block in which the exceedance happens is referred to as the stress period. If one assumes that there are m blocks and each of length n , these m blocks follows the GEV distribution, then estimates of extreme quantiles R n , T of the maximum distribution are obtained by inverting the GEV distribution, as follows:
R n , T = { H γ ^ , μ ^ , σ ^ 1 ( 1 1 T ) = μ ^ σ ^ γ ^ ( 1 ( l o g ( 1 1 T ) ) γ ^ )   for   γ 0 H γ ^ , μ ^ , σ ^ 1 ( 1 1 T ) = μ ^ σ ^ ( l o g ( 1 1 T ) ) for   γ = 0
H γ ^ , μ ^ , σ ^ 1 is the quantile function for GEV distribution (in short H ^ 1 ).
Suppose that n = 250 , i.e., approximately 250 working days in a year and T = 10 , then P ( M 250 > R 250 , 10 ) = 1 / 10 . Then, one can infer that the return level is expected to be exceeded on average once every ten years. The variance of the return level R n , T is obtained by using the delta method [33]:
V a r ( R n , T ) = V a r ( R p ) R p T V R p
where R p T = [ R p μ , R p σ , R p γ ] = [ 1 , γ 1 ( 1 y p γ ) , σ γ 2 ( 1 y p γ ) σ γ 1 y p γ l n y p ] , y p = ln ( 1 p ) and V is the variance-covariance matrix of ( γ ^ , μ ^ , σ ^ ) .
The confidence interval for R n , T with a significance level α equals:
P { R n , T z α V a r ( R n , T ) < R ^ n , T < R n , T + z α V a r ( R n , T ) } = 1 α
To assess the estimates, one can use a return level plot, which comprises a graph of R n , T = μ σ γ ( 1 ( l o g ( 1 1 / T ) ) γ ) against y p = log ( 1 p ) using a logarithmic scale.
Confidence intervals and empirical estimates can be added to the plot. If the GEV model is suitable for the data, then empirical and theoretical estimates should be the same. Generally, more accurate confidence intervals are obtained by evaluating the profile likelihood for the return level. The profile likelihood is a method for making inference on a particular parameter of θ i of parameters vector θ . The log-likelihood for θ can be written as l ( θ i ) , where θ i denotes all parameters of θ excluding θ i . The profile log-likelihood for θ i is defined as l p ( θ i ) = max θ i l ( θ i ) . The profile log-likelihood is the maximum log-likelihood with respect to all remaining parameters in θ excluding θ i .
In this approach, one parameter (e.g., γ = γ 0 ) is fixed, and the log-likelihood function is maximized with respect to the remaining parameters ( μ , σ ) . This process is repeated for different values of γ 0 .
To obtain confidence intervals for the return level, it is required to re-parameterize the GEV distribution so that R n , T is one of the parameters of the GEV. It is done by rearranging Equation (5) as follows:
μ ^ = R n , T + σ ^ γ ^ ( 1 ( l o g ( 1 1 / T ) ) γ ^ )
and then replacing μ in Equation (1) with Equation (8). Consequently, it is possible to perform maximization for the return level as one of the parameters of the GEV distribution. The GEV parameter vector ψ = ( γ , R n , T , σ ) is partitioned into two components i.e., ψ 1 = R n , T and ψ 2 = ( γ , σ ) , and the log-likelihood profile is defined as:
L ( ψ 1 ) = m a x ψ 2 ( ψ 1 , ψ 2 )
Under suitable regularity conditions for large samples, the deviance statistics is defined as:
D ( ψ 1 ) = 2 ( L ( ψ ^ ) L ( ψ 1 ) ) ~ χ 1 2
The set of values C L α for which D ( ψ 1 ) C L α gives a ( 1 α ) confidence interval for R n , T , where C L α is the ( 1 α ) quantile of the χ 1 2 distribution.
If the extremes show a tendency to form clusters (dependency), it is reasonable to inquire about their sizes. To do that, following [31], define the extremal index. Let 0 θ 1 and τ > 0 . Then, one can find a sequence u n ( τ ) such that the following hold:
l i m n n ( 1 F ( u n ( τ ) ) ) l i m n P { M n u n ( τ ) } = e x p ( θ τ )
where u n is a non-decreasing sequence of real numbers, and F ( · ) denotes the marginal distribution. A straightforward way to estimate the extremal index θ is to use the blocks method. The asymptotic estimator of θ is:
θ ^ B M = n 1 l o g ( 1 K u / m ) l o g ( 1 N u / m n )
where K u is the number of blocks in which the threshold is exceeded, m is the number of blocks, N u is the number of exceedances over the threshold u and n is the number of observations in the m -th block. Smith and Weissman [34] investigate the statistical properties of such an estimator as the problem of choosing n (or m ) and u . Note that both m and n should be large. Setting n = 20 observations in the m -th block is deemed fine [34].
O’Brien [35] proposes the runs method where a high threshold u is used as well, but the estimator of θ is based on the runs, which means observations that are below the threshold u are utilized:
θ ^ R M = i = 1 n r I A i , n i = 1 n I { X i > u } = i = 1 n r I A i , n N u
where A i , n = { X i > u , X i + 1 u , , X i + r u } and r is the number of observations below the threshold. Computations by [34] suggest that the runs estimator has a lower bias than the blocks method.
A more recent proposal comes from [36]. They propose an interval estimator based on the rescaled inter-exceedance intervals. Suppose N u exceedances over a threshold value u and 1 t 1 < < t N u n is the time of these exceedances; then the proposed estimator takes the following form:
θ ^ 1 = 2 ( i = 1 N u 1 T i ) 2 ( N u 1 ) ( i = 1 N u 1 T i 2 )
where T i = t i + 1 t i for i = 1 , , N u 1 .
Ferro and Segers [36] propose another estimator of θ with the bias correction:
θ ^ 2 = 2 ( i = 1 N u 1 ( T i 1 ) ) 2 ( N u 1 ) ( i = 1 N u 1 ( T i 1 ) ( T i 2 ) )
The estimator of θ ensures the estimate is within the region of 0 and 1 and is defined as:
θ ^ F S = { 1 θ ^ 1     i f     m a x { T i : 1 i < N u 1 } 2 1 θ ^ 2     i f     m a x { T i : 1 i < N u 1 } > 2
The findings of Süveges [37] support the conjecture that the inter-exceedances estimator possesses the best characteristics. The mean cluster size of extreme values is obtained as reciprocal of θ .

3. EVT in Modeling and Forecasting of Delivery Delays

3.1. Data Characteristics

Using the tools described in Section 2, actual delivery data are analyzed. High-frequency data covering 3770 observations over approximately three months were taken into account.
It is crucial how the estimated time of arrival is calculated. The starting point is to set the actual time of the vehicle loading in the distribution center. Then, using a digital map, the distance and route are determined to the place of distribution included in the delivery order. Calculations take into account the restrictions related to the gross vehicle weight and road speed limitations. The average speed of travel on individual routes is determined based on history. This allows selecting the estimated time of arrival (ETA). In addition, due to the provisions of EC Regulation WE561, the rest time is set for routes that cannot be completed in the primary dimension of driving a vehicle by a driver, i.e., during 9 h of driving a day. If this happens, we will add 11 h of driver rest time to the ETA.
In the described delivery process, about 70–100 vehicles are engaged in the transportation, depending on the demand for transportation services. Thus, delays are not generated by one driver or one vehicle. They occur in different periods in real-time due to vehicle position registration and loading/unloading using GPS/GPRS technology (via intelligent transport system; ITS). The process is homogenous in that one transport company is offering its service to one supplier (a brewery company). The brewery product (beer) can be thought of as homogenous in transportation. It is packed and loaded in a standardized way. The number of receivers is defined as well. In a short period (three months), the receivers can be assumed constant (possible changes are irrelevant). Additionally, the routes are established in advance.
One can consider the causes of delays in transportation from three perspectives:
  • The management perspective, i.e., cargo management, determining optimal routes, analysis of the delivery time from the standpoint of planning, the technical and operational readiness of vehicles, and proper structure of contracts with subcontractors.
  • The road conditions perspective related to the possibility of dynamic changes on the road such as the occurrence of random events and hazards on the road (accidents, collisions), congestion, and changes in weather conditions.
  • The driver’s behavior perspective, including defective implementation of the plan, received from transport managers, psychophysical condition of the driver (exhaustion, stress), and irrational behavior on the road.
Delays appear stochastically, and they are treated as a realization of a stochastic process. Unless the transport management system remains unchanged, delays are time-dependent, although not assigned to an equally spaced time set.
The time series comprises the delays, i.e., the difference between the expected and the actual delivery time, expressed in hours (Figure 1). Due to many transport operations and a relatively small number of delays (51 delays in total), the ratio of delays equals 1.35% of the total.
Table 1 compared the characteristics of the entire set of the data (3770 observations) and delay data (51 observations). The minimum for all deliveries equals zero hours, i.e., delivery was before the expected time, so there is no delay.
The comparison illustrates that when the total data set is considered, delays are not significant. The median is equal to zero, and the mean value equals 8 min 48 s. However, standard deviation and shape parameters (skewness and kurtosis) indicate that the rare data distort the distribution. Looking at delay data (51 observations), it becomes clear that the observed delays are from 15 min to 42 h, which corresponds to the scale of possible consequences for a given carrier and the entire supply chain. The average delay is around 11 h (10:59:45), taking the mean, and 12 h taking the median. The standard deviation of almost 8 h shows considerable variability in delays, while shape parameters indicate departures from ideal values in Gaussian distribution. The latter is supported by the Jarque-Bera test [38] results, allowing rejecting the hull hypothesis assuming normality. Our findings are in line with [23]. They described time delays in travel time. State-of-the-art practices often use the mean and variance of travel time to derive buffer time and planning time indices as measures of road reliability. However, empirical travel time analysis showed that travel time distribution is not symmetrical but highly skewed to the right with a heavy tail, especially in the presence of road network disruptions resulting from adverse weather conditions, car collisions, or other incidents.
Figure 1 presents two plots with the delivery delays distributed over time. The upper one depicts all deliveries (namely 3770 observations) where deliveries with no delays are assigned zero hours. We can see that some of the delays happened on the same day, i.e., indicating some link between delays. The lower plot depicts only delays that are presented consecutively for better visibility. Most of the delays are less than 24 h, with one notable exception where the delay reached 42 h.

3.2. The Empirical Results

The analyzed type of time series data is called a randomly distributed zero-inflated time series and has specific properties [39,40,41]. The time series can be viewed as two processes when delays occur or not. In practice, we are interested in modeling the delays only, so the focus is directed at the maxima series rather than the expected value. That is why the time series consisting of maxima, i.e., delays in hours, are used. Then, the EVT is employed. The maximum likelihood (ML) method is applied to estimate the parameters of the GEV distribution. For small samples, the maximum likelihood estimator is biased, and its variance is higher than for the case of other methods, such as the Penalized Maximum Likelihood method (PML).
Table 2 shows estimates of the GEV distribution parameters with standard errors and p-values of significance for both cases, ML and PML. Only the parameter γ is insignificant (ML), while the others are highly significant. The maximized log-likelihood and the Bayesian Information Criterion (BIC) indicate the penalized ML method estimates’ advantage [42].
Figure 2 and Figure 3 depict various plots, i.e., panel (a)—probability, panel (b)—quantile, and panel (c) and (d) return level plots, as diagnostic tools (More about diagnostic plots and their interpretations can be found in [33,43,44,45].). For the probability plot and quantile plot, the points should lie close to the unit diagonal. Significant deviations from linearity mean that the model is not adequate for the data. The probability plot is the empirical distribution function evaluated at z i , which are ordered extreme values. Consequently, the probability plot consists of the points { ( H ^ ( z i ) ,   ( i / ( m + 1 ) ) } ,   i = 1 , , m where m is the sample size. The quantile plot consists of the points { ( H ^ 1 ( i / ( m + 1 ) ) ) ,   z i } ,   i = 1 , , m where H ^ 1 is the quantile function for GEV distribution line.
The plotted empirical points are close to linear, indicating that the GEV model is valid for extrapolation purposes, although there are some deviations from linearity, especially in the probability plot. The quantile plots are quite linear except for one observation. The variability of the observations likely causes the departures from linearity. Nonetheless, all plots seem to favor the PML over the ML methods. The probability and the quantile plots are similar. However, a small distinction exists. A significant difference is observed for the return level plots, where confidence intervals are narrower for penalized ML than for the ML method.
To check the validity of the observed plots, we employed the Anderson-Darling [46], Cramer [47] and von Mises [48] tests to check the goodness-of-fit of the fitted GEV distribution. The Anderson-Darling test is one of the most powerful tests based on the empirical distribution function. Table 3 presents results for both cases, i.e., the maximum likelihood method and the penalized maximum likelihood method. At 5% significance level, we do not reject the null hypothesis, so there is no evidence against assuming that the empirical distribution is from the GEV distribution for both estimation methods. However, the Anderson-Darling test at 10% significance level indicates rejection of the null so that it may be improved upon.
Figure 4 shows the profile log-likelihood estimated for the return level when the return period is 5 and 10 days from the PML. The 95% confidence intervals (orange lines) are rather symmetrical, and the bell-curve shape of the profile log-likelihood indicates that both estimates are fairly accurate.
Figure 5 displays a range of return levels and their confidence intervals obtained from the profile likelihood method. Firstly, it should be noted that Confidence Intervals (CIs) using the Penalized ML method are closer to the return levels than those obtained from the ML method. Secondly, the CIs for the PML are close to being symmetrical, contrary to the CIs from the ML method, where considerable asymmetry is visible. It indicates that estimates of the return levels and CIs from the PML are less uncertain in estimating the bands, meaning they have smaller differences between the upper and the lower bands than those obtained from the ML method. The PML estimates indicate that for every five delays, one expects to have a maximum delay of about 18 h and 16 min with 95% confidence intervals of (16:12; 20:38).
This indicates that estimates of the return levels and CIs from the PML are more appropriate than those obtained from the ML method. For every five delays, one expects to have a maximum delay of about 18 h and 16 min with 95% confidence intervals of (16:12; 20:38). For every ten delays, this is, respectively, about 22 h and 22 min with a 95% confidence interval of (20:02; 25:00).
It is interesting to infer the implied cluster size of the series of extreme delays. It is reasonable to assume that one delay can cause another, which means that delays may form clusters. For the original time series (3770 observations), there are on average 60 deliveries per day, so setting 63 blocks gives nearly 60 observations per block. By contrast, the median of deliveries equals 66, generating 57 blocks. Setting the latter number of observations per block enables us to make interpretations in terms of days. It is common practice when it comes to the extremal index [30,32]. Table 4 shows the extremal index estimates using three different estimators, i.e., blocks method, runs method, and Ferro and Segers method.
The results show that delays tend to form clusters with a mean length between two and three days. Since the [36] estimator generates more accurate values than others, it is expected to be more trustworthy. These findings indicate that it is reasonably likely that one delay causes further delays.

4. Discussion

Divery is inevitably a part of the supply chain. Supply chain management is typically described as a proactive relationship and integration among various tiers in the chain [49]. Delivery represents a link between different elements and subcontractors [50,51,52,53].
The empirical results show that the extreme value theory applied to delivery delays is a promising tool of a supply chain distortion analysis. The presented approach is experience-based. It makes the assessment of delivery dependent on Key Performance Indicators [54]. One of them is timeliness, expressed by the ratio of ‘correct’ deliveries measured by the number of deliveries made on time to the total number of completed transports. The second is the delivery quality measured by the number of damaged cargo units to the total number of transported units in a certain number of transport operations.
Since the EVT application in the present study is relatively novel in the specified area, we recommend including the EVT module in the integrated transport management system, being a part of the supply chain management system. It would be beneficial at the stage of anticipating delays and planning for reserves of vehicles and drivers. As the entire process is repeatable, finding optimal paths and time windows for deliveries within the framework of different constraints is possible. However, everyday practice indicates that rare but highly impactful delays are responsible for the distortion of the supply chain system, including optimizing the transportation processes.
Using the EVT as a tool of permanent monitoring and planning transportation in the supply chain requires a perfect information system integrated over all parties involved [55]. The optimal solution is using autonomous information collection without human intervention. Such systems are already utilized in just-in-time logistics, particularly in the production sector, e.g., automotive [54]. As concerns the theoretical model, it is straightforward that the delivery delays follow the process:
y t = { z t       if   there   is   a   delay 0             otherwise
It is a randomly distributed zero-inflated time series [39,40,41]. The time series can be viewed as two processes, whether delays occur or not. It is specific in the sense that we focus only on z t i.e., delays because zero is a normal state of no delays. Delays are distributed randomly. Thus, the process dynamics are different from the Ornstein-Uhlenbeck process, as was considered in [56]. It is often called a jump-diffusion process or the Ornstein–Uhlenbeck processes driven by the Lévy process [57,58]. Here, the state of diffusion is stable and equal to zero. The Process (20) can be modeled in different ways. One of the most often representations is the Markov model with Poisson distribution for count data [59].
Future works would include a comparison of the EVT model with the Markov zero-inflated Poisson model. The possible development is related to identifying factors of delivery delays as well. In addition, it would also be interesting to extend the analysis by studying the possibility that deliveries are the nonstationary extremes.

5. Conclusions

The purpose of this paper was to model and forecast rare events that appear as delivery delays in road transport using the Extreme Value Theory (EVT). We evaluated both the magnitude and frequency of delivery delays. The actual data set originated from the intelligent transport system (ITS) from a large European brewery company, and the cooperating road transport company has been studied. The results show that out of 3770 transport operations under study, only 51 delays have occurred, i.e., a 1.35% ratio, and the discrepancy between the plan and its implementation is low. This number justified applying the extreme value approach.
The empirical results positively answer the research questions. The extreme value theory is the appropriate tool for analyzing delivery delays since it gives estimated Generalized Extreme Value (GEV) distribution quantiles that fit the actual distribution. Estimating the extremal index in a considered process shows that the logistics operator should plan additional resources (vehicle and driver) every two to three days to guarantee a smooth delivery and flexible supply chain management.
The advantage of the EVT is that it allows considering any period. If the business environment is stable, the results of the EVT are stable as well. Its application can be easily extended for prediction, planning, and decision-making in any means of transportation as a part of a supply chain.

Author Contributions

Conceptualization, M.F., M.O., and W.Z.; methodology, M.F.; software, M.F.; validation, M.F., M.O., and W.Z.; formal analysis, W.Z.; investigation, M.O., and W.Z.; resources, W.Z.; data curation, W.Z.; writing—original draft preparation, M.O.; writing—review and editing, M.F., M.O., and W.Z.; visualization, M.F., and W.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in RepOD, V1 at https://doi.org/10.18150/FZS9EI (accessed on 21 June 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davison, A.C.; Smith, R.L. Models for exceedances over high thresholds. J. R. Stat. Soc. Ser. B (Methodol.) 1990, 52, 393–442. Available online: https://www.jstor.org/stable/2345667 (accessed on 21 June 2021). [CrossRef]
  2. Tawn, J.A. Estimating probabilities of extreme sea-levels. Appl. Stat. 1992, 41, 77–93. [Google Scholar] [CrossRef]
  3. De Haan, L.; de Ronder, J. Sea and wind: Multivariate extremes at work. Extremes 1998, 7, 7–45. [Google Scholar] [CrossRef]
  4. Barão, M.I.; Tawn, J.A. Extremal analysis of short series with outliers: Sea-levels and athletic records. J. R. Stat. Soc. 1999, 48, 469–487. [Google Scholar] [CrossRef]
  5. Nadarajah, S.; Shiau, J.T. Analysis of extreme flood events for the pachang river Taiwan. Water Resour. Manag. 2005, 9, 363–374. [Google Scholar] [CrossRef]
  6. Katz, R.W.; Parlange, M.B.; Naveau, P. Statistics of extremes in hydrology. Adv. Water Resour. 2002, 25, 1287–1304. [Google Scholar] [CrossRef] [Green Version]
  7. Caires, S.; Swail, V.R.; Wang, X.L. Projection and analysis of extreme wave climate. J. Clim. 2006, 19, 5581–5605. [Google Scholar] [CrossRef]
  8. Panagoulia, D.; Economou, P.; Caroni, C. Stationary and nonstationary generalized extreme value modelling of extreme precipitation over a mountainous area under climate change. Environmetrics 2014, 25, 29–43. [Google Scholar] [CrossRef]
  9. Masingi, V.N.; Maposa, D. Modelling long-term monthly rainfall variability in selected provinces of South Africa: Trend and extreme value analysis approaches. Hydrology 2021, 8, 70. [Google Scholar] [CrossRef]
  10. McNeil, A.J. Extreme value theory for risk managers. In Internal Modelling and CAD II; RISK Books: London, UK, 1999; pp. 93–113. [Google Scholar]
  11. Danielsson, J.; DE Vries, C.G. Value-at-risk and extreme returns. Ann. D’economie Stat. 2000, 60, 239–270. [Google Scholar] [CrossRef] [Green Version]
  12. McNeil, J.A.; Frey, F. Estimation of tail-related risk measures for heteroscedastic financial time series: An extreme value approach. J. Empir. Financ. 2000, 7, 271–300. [Google Scholar] [CrossRef]
  13. Fałdziński, M.; Osińska, M.; Zdanowicz, T. Detecting risk transfer in financial markets using different risk measures. Cent. Eur. J. Econ. Model. Econom. 2012, 1, 45–64. [Google Scholar]
  14. Bień-Barkowska, K. Looking at extremes without going to extremes: A new self-exciting probability model for extreme losses in financial markets. Entropy 2020, 22, 789. [Google Scholar] [CrossRef]
  15. Božović, M. Portfolio tail risk: A multivariate extreme value theory approach. Entropy 2020, 22, 1425. [Google Scholar] [CrossRef]
  16. Fałdziński, M.; Osińska, M. The use of range-based volatility estimators in testing for Granger causality in risk on international capital markets. J. Risk Model Valid. 2020, 14. [Google Scholar] [CrossRef]
  17. Gijbels, I.; Kika, V.; Omelka, M. Multivariate tail coefficients: Properties and estimation. Entropy 2020, 22, 728. [Google Scholar] [CrossRef]
  18. Trucíos, C.; Tiwari, A.V.; Alqahtani, F. Value-at-risk and expected shortfall in cryptocurrencies’ portfolio: A vine copula–based approach. Appl. Econ. 2020, 52, 2580–2593. [Google Scholar] [CrossRef]
  19. Kagan, Y.Y. Earthquake size distribution and earthquake insurance. Commun. Stat. Stoch. Models 1997, 13, 775–797. [Google Scholar] [CrossRef]
  20. Strand, M.; Boes, D. Modeling road racing times of competitive recreational runners using extreme value theory. Am. Stat. 1998, 52, 205–210. [Google Scholar]
  21. Zheng, L.; Sayed, T. Application of extreme value theory for before-after road safety analysis. Transp. Res. Rec. 2019, 2673, 1001–1010. [Google Scholar] [CrossRef]
  22. Orsini, F.; Gecchele, G.; Gastaldi, M.; Rossi, R. Large-scale road safety evaluation using extreme value theory. IET Intell. Transp. Syst. 2020, 14. [Google Scholar] [CrossRef]
  23. Esfeh, M.A.; Kattan, L.; Lam, W.H.K.; Esfe, R.A.; Salari, M. Compound generalized extreme value distribution for modeling the effects of monthly and seasonal variation on the extreme travel delays for vulnerability analysis of road network. Transp. Res. Part C 2020, 120, 102808. [Google Scholar] [CrossRef]
  24. Abadi, A.; Ioannou, P.A.; Dessouky, M.M. Multimodal dynamic freight load balancing. IEEE Trans. Intell. Transp. Syst. 2016, 2, 356–366. [Google Scholar] [CrossRef]
  25. Chen, D.; Ignatius, J.; Sun, D.; Goh, M.; Zhan, S. Impact of congestion pricing schemes on emissions and temporal shift of freight transport. Transp. Res. Part E 2018, 118, 77–105. [Google Scholar] [CrossRef] [Green Version]
  26. Cavone, G.; Dotoli, M.; Seatzu, C. A survey on petri net models for freight logistics and transportation systems. IEEE Trans. Intell. Transp. Syst. 2017, 6, 1795–1813. [Google Scholar] [CrossRef]
  27. Iordanidou, G.; Papamichail, I.; Roncoli, C.; Papageorgiou, M. Feedback-based integrated motorway traffic flow control with delay balancing. IEEE Trans. Intell. Transp. Syst. 2017, 9, 2319–2329. [Google Scholar] [CrossRef]
  28. Wang, T.; Xing, Z.; Hu, H.; Qu, X. Overbooking and delivery-delay-allowed strategies for container slot allocation. Transp. Res. Part E 2019, 433–447. [Google Scholar] [CrossRef]
  29. Fisher, R.A.; Tippett, L.H.C. Limiting forms of the frequency distribution of the largest or smallest member of a sample. Math. Proc. Camb. Philos. Soc. 1928, 24, 180–290. [Google Scholar] [CrossRef]
  30. Embrechts, P.; Klüppelberg, C.; Mikosch, T. Modelling Extremal Events for Insurance and Finance; Springer: New York, NY, USA, 2003. [Google Scholar]
  31. Coles, S.G.; Dixon, M.J. Likelihood-based inference for extreme value models. Extremes 1999, 1, 5–23. [Google Scholar]
  32. McNeil, A.J. Calculating Quantile Risk Measures for Financial Time Series Using Extreme Value Theory; ETH Zurich: Zurich, Switzerland, 1998. [Google Scholar] [CrossRef]
  33. Coles, S.G. An Introduction to Statistical Modeling of Extreme Values; Springer: New York, NY, USA, 2001. [Google Scholar]
  34. Smith, R.L.; Weissman, I. Estimating the extremal index. J. R. Stat. Soc. Ser. B (Methodol.) 1994, 3, 515–528. [Google Scholar] [CrossRef] [Green Version]
  35. O’Brien, G. Extreme values for stationary and Markov sequences. Ann. Probab. 1987, 15, 281–291. [Google Scholar] [CrossRef]
  36. Ferro, C.A.T.; Segers, J. Inference for clusters of extreme values. J. R. Stat. Soc. Ser. B (Methodol.) 2003, 65, 545–556. [Google Scholar] [CrossRef]
  37. Süveges, M. Likelihood estimation of the extremal index. Extremes 2007, 10, 41–55. [Google Scholar] [CrossRef]
  38. Jarque, C.M.; Bera, A.K. A test for normality of observations and regression residuals. Int. Stat. Rev. 1987, 55, 163–172. [Google Scholar] [CrossRef]
  39. Cameron, A.C.; Trivedi, P.K. Essentials of count data regression. In A Companion to Theoretical Econometrics; Baltagi, P.H., Ed.; Blackwell: Oxford, UK, 2001; pp. 331–348. [Google Scholar]
  40. Winkelmann, R. Econometric Analysis of Count Data; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  41. Hilbe, J.M. Modeling Count Data; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  42. Neath, A.A.; Cavanaugh, J.E. The Bayesian information criterion: Background, derivation, and applications. Wiley Interdiscip. Rev. Comput. Stat. 2012, 4, 199–203. [Google Scholar] [CrossRef]
  43. Fawcett, L.; Walshaw, D. Extreme Values: Statistical Analysis Using R; Wiley: Hoboken, NJ, USA, 2016. [Google Scholar]
  44. Alam, M.A.; Emura, K.; Farnham, C.; Yuan, J. Best-fit probability distributions and return periods for maximum monthly rainfall in Bangladesh. Climate 2018, 6, 9. [Google Scholar] [CrossRef] [Green Version]
  45. Iyamuremye, E.; Mung’atu, J.; Mwita, P. Extreme value modelling of rainfall using poisson-generalized pareto distribution: A case study tanzania. Int. J. Stat. Distrib. Appl. 2019, 5, 67–75. [Google Scholar] [CrossRef] [Green Version]
  46. Anderson, T.W.; Darling, D.A. A Test of goodness-of-fit. J. Am. Stat. Assoc. 1954, 49, 765–769. [Google Scholar] [CrossRef]
  47. Cramér, H. On the composition of elementary errors. Scand. Actuar. J. 1928, 1, 13–74. [Google Scholar] [CrossRef]
  48. Mises, R.E. Wahrscheinlichkeit, Statistik und Wahrheit; Springer: Berlin/Heidelberg, Germany, 1928; Available online: https://www.springer.com/de/book/9783662418635 (accessed on 21 June 2021).
  49. Trkman, P.; Stemberger, M.I.; Jaklic, J.; Groznik, A. Process approach to supply chain integration. Supply Chain. Manag. 2007, 12, 116–128. [Google Scholar] [CrossRef]
  50. Zeng, A.Z. Global sourcing: Process and design for efficient management. Supply Chain Manag. 2003, 8, 367–379. [Google Scholar] [CrossRef]
  51. Christopher, M. Logistics and Supply Chain Management: Creating Value-Adding Networks; Prentice Hall: London, UK, 2005; Available online: https://www.mscceducation.com/wp-content/uploads/2019/06/SCM-Logistics-and-SCM-Creating-value-added-networks-Pg-1-to-17.pdf (accessed on 21 June 2021).
  52. Lambert, D.M.; García-Dastugue, S.J.; Croxton, K.L. An evaluation of process-oriented supply chain management frameworks. J. Bus. Logist. 2005, 26, 25–51. [Google Scholar] [CrossRef]
  53. Gunasekaran, A.; Lai, K.-H.; Cheng, T.C.E. Responsive supply chain: A competitive strategy in a networked economy. Omega 2008, 36, 549–564. [Google Scholar] [CrossRef]
  54. Dörnhöfer, M.; Schröder, F.; Günthner, W.A. Logistics performance measurement system for the automotive industry. Logist. Res. 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  55. Kumar, R.S.; Pugazhendhi, S. Information sharing in supply chains: An overview. Procedia Eng. 2012, 38, 2147–2154. [Google Scholar] [CrossRef] [Green Version]
  56. Guel-Cortez, A.J.; Kim, E.-J. Information length analysis of linear autonomous stochastic processes. Entropy 2020, 22, 1265. [Google Scholar] [CrossRef]
  57. Önalan, O. The Ornstein-uhlenbeck processes driven by lévy process and application to finance. Electron. Eng. Comput. Technol. 2010, 60, 443–453. [Google Scholar] [CrossRef]
  58. Kou, S. A jump-diffusion model for option pricing. Manag. Sci. 2002, 48, 1086–1101. [Google Scholar] [CrossRef] [Green Version]
  59. Wang, P. Markov zero-inflated Poisson regression models for a time series of counts with excess zeros. J. Appl. Stat. 2001, 28, 623–632. [Google Scholar] [CrossRef]
Figure 1. Delivery delays. Panel (a) covers all observations (3770 obs), and panel (b) covers delivery delays (51 obs).
Figure 1. Delivery delays. Panel (a) covers all observations (3770 obs), and panel (b) covers delivery delays (51 obs).
Entropy 23 00788 g001
Figure 2. Diagnostic plots for Generalized Extreme Values distribution (ML method).
Figure 2. Diagnostic plots for Generalized Extreme Values distribution (ML method).
Entropy 23 00788 g002
Figure 3. Diagnostic plots for Generalized Extreme Values distribution (PML method).
Figure 3. Diagnostic plots for Generalized Extreme Values distribution (PML method).
Entropy 23 00788 g003
Figure 4. The profile likelihood for the return level when the return period is 5, panel (a), and 10 days, panel (b).
Figure 4. The profile likelihood for the return level when the return period is 5, panel (a), and 10 days, panel (b).
Entropy 23 00788 g004
Figure 5. Return level estimates with 95% confidence intervals obtained from the profile likelihood for the ML method, panel (a) and the PML, panel (b).
Figure 5. Return level estimates with 95% confidence intervals obtained from the profile likelihood for the ML method, panel (a) and the PML, panel (b).
Entropy 23 00788 g005
Table 1. Descriptive statistics [HH:MM:SS].
Table 1. Descriptive statistics [HH:MM:SS].
No. of Obs.MinMaxMedianMeanStd. dev.KurtSkewJB StatJB
p-Value
377000:00:0042:00:0000:00:0000:08:4801:32:52228.93913.3668,131,0730.000
5100:15:0042:00:0012:00:0010:59:4507:54:335.6650.94222.6340.000
Note: The JB denotes the Jarque-Bera test for normality.
Table 2. Estimates of the Generalized Extreme Values (GEV) distribution.
Table 2. Estimates of the Generalized Extreme Values (GEV) distribution.
MLCoefficientStd. Errorp-Value
γ−0.04520.09810.6449
μ0.31290.04310.0000
σ0.27310.03160.0000
Log likelihood−12.7468BIC0.7312
Penalized MLCoefficientStd. Errorp-Value
γ−0.18020.04340.0000
μ0.34350.04730.0000
σ0.31820.04020.0000
Log likelihood−6.8146BIC0.4985
Table 3. The goodness-of-fit tests.
Table 3. The goodness-of-fit tests.
MLPML
Statisticp-ValueStatisticp-Value
Anderson-Darling 1.87940.10712.01730.089
Cramer-von Mises 0.30960.12690.29230.1420
Table 4. Extremal index estimation results.
Table 4. Extremal index estimation results.
Blocks MethodRuns MethodFerro and Segers
Extremal index θ ^   ( m = 63 ) 0.44030.29410.3193
θ ^ 1   ( m = 63 ) 2.273.403.13
Extremal index θ ^ ( m = 57 ) 0.39310.31370.3193
θ ^ 1 ( m = 57 ) 2.543.193.13
Note: The assumed number of exceedances over a high threshold Nu = 51.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fałdziński, M.; Osińska, M.; Zalewski, W. Extreme Value Theory in Application to Delivery Delays. Entropy 2021, 23, 788. https://doi.org/10.3390/e23070788

AMA Style

Fałdziński M, Osińska M, Zalewski W. Extreme Value Theory in Application to Delivery Delays. Entropy. 2021; 23(7):788. https://doi.org/10.3390/e23070788

Chicago/Turabian Style

Fałdziński, Marcin, Magdalena Osińska, and Wojciech Zalewski. 2021. "Extreme Value Theory in Application to Delivery Delays" Entropy 23, no. 7: 788. https://doi.org/10.3390/e23070788

APA Style

Fałdziński, M., Osińska, M., & Zalewski, W. (2021). Extreme Value Theory in Application to Delivery Delays. Entropy, 23(7), 788. https://doi.org/10.3390/e23070788

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop