Next Article in Journal
Three-Way Decisions Based on Q-Rung Orthopair Fuzzy 2-Tuple Linguistic Sets with Generalized Maclaurin Symmetric Mean Operators
Next Article in Special Issue
On the Condition of Independence of Linear Forms with a Random Number of Summands
Previous Article in Journal
Acceleration of Boltzmann Collision Integral Calculation Using Machine Learning
Previous Article in Special Issue
Characterization of Probability Distributions via Functional Equations of Power-Mixture Type
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Extended Log-Logistic Distribution: Inference and Actuarial Applications

1
Department of Mathematics & Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
2
Department of Mathematics, Faculty of Science, Tanta University, Tanta 31527, Egypt
3
Department of Statistics, Mathematics and Insurance, Benha University, Benha 13511, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(12), 1386; https://doi.org/10.3390/math9121386
Submission received: 5 May 2021 / Revised: 2 June 2021 / Accepted: 8 June 2021 / Published: 15 June 2021
(This article belongs to the Special Issue Characterization of Probability Distributions)

Abstract

:
Actuaries are interested in modeling actuarial data using loss models that can be adopted to describe risk exposure. This paper introduces a new flexible extension of the log-logistic distribution, called the extended log-logistic (Ex-LL) distribution, to model heavy-tailed insurance losses data. The Ex-LL hazard function exhibits an upside-down bathtub shape, an increasing shape, a J shape, a decreasing shape, and a reversed-J shape. We derived five important risk measures based on the Ex-LL distribution. The Ex-LL parameters were estimated using different estimation methods, and their performances were assessed using simulation results. Finally, the performance of the Ex-LL distribution was explored using two types of real data from the engineering and insurance sciences. The analyzed data illustrated that the Ex-LL distribution provided an adequate fit compared to other competing distributions such as the log-logistic, alpha-power log-logistic, transmuted log-logistic, generalized log-logistic, Marshall–Olkin log-logistic, inverse log-logistic, and Weibull generalized log-logistic distributions.

1. Introduction

Modeling insurance losses data has received significant interest from actuaries and risk managers who often evaluate and study the unlikely outcomes that the value-at-risk may express by chance. The insurance data are usually unimodal [1], right-skewed [2], positive [3], have a heavy-tailed density [4], and have a unimodal hump shape [5].
There is a clear need to develop and propose more flexible distributions by extending the well-known classical distributions or by introducing a new family to model several insurance datasets such as financial returns, unemployment insurance data, insurance losses data, and risk management data, among others.
The log-logistic (LL) distribution is also known as the Fisk distribution in the income distribution literature [6]. Some authors, such as [7,8], have referred to the Fisk distribution as the LL distribution, whereas Arnold [9] referred to it as the Pareto Type III distribution and included an additional location parameter. Further details about the LL model can be found in [10]. Several authors have studied different generalized forms of the LL distribution to improve its capability and flexibility. Some notable examples are the following: Kumaraswamy-LL [11], beta-LL [12], Marshall–Olkin LL [13], McDonald LL [14], Zografos-Balakrishnan LL [15], and odd Lomax LL distributions [16].
This article suggests a new version of the LL distribution called the extended log-logistic (Ex-LL) distribution, which can provide more flexibility in modeling insurance data than other competing models. Hence, the aim of the paper was three-fold. The first was devoted to proposing the Ex-LL model as a new form of the LL distribution via the Type-I heavy-tailed G family [17]. The Ex-LL model has some desirable characteristics as follows: It is a more flexible extension with only an extra shape parameter for the LL distribution, and it improves goodness-of-fit to real-life data; it produces negative and positive skewness, and also, it causes the kurtosis to be more flexible as compared to the baseline LL model. The Ex-LL model has a heavier tail than the LL and exponentiated log-logistic (ELL) models (see Section 6). The Ex-LL hazard rate function (HRF) exhibits an increasing shape, an upside-down bathtub shape, a decreasing-J shape, and a reversed-J shape. Its density exhibits a reversed-J shape, a symmetrical shape, an asymmetrical shape (right-skewed or left-skewed), a unimodal shape, and a J shape. Furthermore, it can be adopted to analyze various data in the engineering and actuarial sciences. Two sets of real data from the insurance and engineering fields were fitted using the Ex-LL model, showing its flexibility over other competing distributions.
Second, we estimated the Ex-LL parameters using some estimation approaches: the least-squares estimators (LSEs), the Anderson–Darling estimators (ADEs), the maximum likelihood estimators (MLEs), the weighted least-squares estimators (WLSEs), and the Cramér–von Mises estimators (CVMEs). The proposed estimators were compared through extensive simulations to determine the best estimation approach to estimate the Ex-LL parameters.
Estimating the exposure to market risk, due to changes in the prices of equity exchange rates and interest rates, is considered an important subject in the actuarial sciences. Therefore, our third objective was devoted to deriving some important actuarial or risk measures based on the Ex-LL distribution called the value-at-risk (VaR), tail value-at-risk (TVaR), tail variance premium (TVP), tail variance (TV), and expected shortfall (ES). These measures are very important in portfolio optimization under uncertainty.
The article is outlined as follows. The new Ex-LL distribution is defined in Section 2. In Section 3, its distributional properties were obtained. Five classical estimation methods were adopted to estimate the Ex-LL parameters, as shown in Section 4. We provide, in Section 5, detailed numerical results for the introduced methods. In Section 6, we explored five important risk measures based on the Ex-LL model, and provided numerical simulation results for these risk measures. In Section 7, two real-life data from the engineering and actuarial sciences were analyzed to explore the importance and flexibility of the Ex-LL distribution. Some final concluding remarks were provided in Section 8.

2. The Ex-LL Distribution

The cumulative distribution function (CDF) of the two-parameter LL model is given (for x > 0 ) by G ( x ) = 1 + λ x β 1 , β , λ > 0 , and its probability density function (PDF) reduces to g ( x ) = β λ x β 1 1 + λ x β 2 .
We define the CDF of the Ex-LL model based on the extended family [17] which is specified (for x R ) by the CDF
F x = 1 1 G x ; φ 1 1 α G x ; φ α , α > 0 ,
and the PDF
f x = α 2 g x ; φ 1 G x ; φ α 1 1 1 α G x ; φ α + 1 , α > 0 ,
where G x ; φ is any baseline CDF that depends on φ R . The CDF (1) reduces to the baseline CDF for α = 1 .
The CDF of the Ex-LL distribution follows, by inserting the LL CDF in Equation (1), as
F ( x ) = 1 λ λ + α x β α , α , β , λ > 0 , x > 0 .
The corresponding PDF of the Ex-LL distribution takes the form
f ( x ) = α 2 β λ α x β 1 1 λ + α x β α + 1 , α , β , λ > 0 , x > 0 .
The LL model follows as a special case from (4) with α = 1 . The EX-LL HRF reduces to
h ( x ) = α 2 β x β 1 λ + α x β , α , β , λ > 0 , x > 0 .
Plots of the PDF and HRF of the Ex-LL distribution are depicted respectively in Figure 1 and Figure 2. These plots show that the Ex-LL density exhibits a symmetrical shape, an asymmetrical shape, a J shape, a unimodal shape, and a reversed-J shape. Furthermore, its HRF exhibits an increasing shape, an upside-down bathtub shape, a decreasing-J shape, and a reversed-J shape.

3. Mathematical Properties

3.1. Mode and Quantile Function

On differentiating the logarithm of (4) with respect to x and equating to zero, the unique mode of the Ex-LL distribution follows as
x 0 = ( β 1 ) λ α ( α β + 1 ) 1 / β , β > 1 .
For β 1 , the Ex-LL distribution has no mode.
The quantile function (QF) of the Ex-LL distribution follows as
Q ( p ) = λ ( 1 p ) 1 / α 1 α 1 / β , 0 < p < 1 .
Let p Uniform ( 0 , 1 ) , then the QF of the Ex-LL distribution can be adopted to generate its random data by the formula
x i = λ ( 1 p i ) 1 / α 1 α 1 / β , i = 1 , 2 , , n .

3.2. Moments and Moment Generating Function

The rth ordinary moment of the Ex-LL distribution is given by
μ r = E ( X r ) = 0 x r f ( x ) d x = Γ r + β β α λ r β Γ α r β Γ ( α ) , r β < α .
The first four ordinary moments of the Ex-LL distribution can be obtained directly by setting r = 1 , 2 , 3 , and 4 in the last equation.
The mean ( μ 1 ), variance ( V a r ( X ) ), skewness ( γ 1 ( X ) ), and kurtosis ( γ 2 ( X ) ) of the Ex-LL distribution are obtained numerically for some choices of its parameters using the Wolfram Mathematica program version 12.0. These results are reported in Table 1. It is noted that the Ex-LL skewness varies in the interval ( 0.84559 , 45.48282 ) , whereas the LL skewness varies only within the interval ( 0.35216 , 1.81999 ) . Furthermore, the kurtosis of the Ex-LL distribution spreads much larger in the interval ( 4.10572 , 2573.493 ) , whereas the kurtosis of the LL distribution spreads only in the interval ( 4.50847 , 14.76564 ) for same parameter values.
The rth incomplete moment of the Ex-LL distribution has the form
I r ( t ) = 0 t x r f ( x ) d x = α 2 β λ α 1 t β + r 1 λ + α t β α 2 F 1 1 , r β α + 1 ; r β + 2 ; t β α λ β + r ,
where 2 F 1 1 , r β α + 1 ; r β + 2 ; t β α λ is the hyper geometric function.
The moment generating function of the Ex-LL distribution takes the form
M ( t ) = k = 0 t k Γ k + β β α λ k β Γ α k β k ! Γ ( α ) , k β < α .

3.3. Mean Residual Life, Mean Inactivity Time and Inequality Curves

The mean residual life of the Ex-LL distribution at age t has the form
M R L = 1 I 1 ( t ) S ( t ) t = ( β + 1 ) α 2 β λ α 1 t β + 1 1 λ + α t β α 2 F 1 1 , 1 α + 1 β ; 2 + 1 β ; t β α λ ( β + 1 ) λ λ + α t β α t ,
where I 1 ( t ) is the first incomplete moments.
The mean inactivity time of the Ex-LL distribution takes the form
M I T = t I 1 ( t ) F ( t ) = α 2 β λ α 1 t β + 1 1 λ + α t β α 2 F 1 1 , 1 α + 1 β ; 2 + 1 β ; t β α λ ( β + 1 ) λ λ + α t β α 1 + t ,
The Lorenz, Bonferroni and Zenga curves (see, Lorenz [18] and Bonferroni [19]) are considered the most important inequality curves and their useful applications are common in insurance, reliability, medicine, and economics.
The Lorenz curve is defined for the Ex-LL distribution as follows:
L ( p ) = I 1 ( x p ) μ = α β λ α Γ ( α ) α λ 1 β + 1 x p β + 1 1 λ + α x p β α 2 F 1 1 , 1 α + 1 β ; 2 + 1 β ; α x p β λ ( β + 1 ) Γ 1 + 1 β Γ α 1 β ,
where F ( x p ) = p , x p is the QF, and I 1 ( t ) refers to the first incomplete moments.
The Bonferroni and Zenga inequality curves can be determined, through their relationship with the Lorenz curve, by the following formulae ([20])
B ( p ) = L ( p ) p and Z ( p ) = L ( p ) p p [ 1 L ( p ) ] .

3.4. Some Entropies

The entropy of the random variable X has important applications in many applied fields including statistics for testing hypotheses [21] and engineering, physics, and information theory to describe dynamical systems or nonlinear chaotic [22]. Furthermore, Song [23] developed the log-likelihood-based distribution measure based on the Rényi information. Song’s measure is exist and can be defined for all distributions. Song’s measure provides meaningful comparisons between distributions as compared with traditional kurtosis measure.
The Rényi P X ( k ) , Tsallis L X ( k ) , and Shannon H X ( 1 ) entropies of the Ex-LL distribution can be derived for the Ex-LL distribution by the following formulae.
P X ( k ) = 1 1 k log x = 0 f k ( x ) d x , k > 0 , k 1 = 1 k 1 k log ( α β ) ( k 1 ) log α λ β + log β Γ ( α k + k ) Γ k ( β 1 ) + 1 β Γ α β k + k 1 β
and
L X ( k ) = 1 1 k x = 0 f k ( x ) d x 1 = ( α β ) k Γ k ( β 1 ) + 1 β α λ k 1 β Γ α β k + k 1 β β Γ ( α k + k ) β ( 1 k ) Γ ( α k + k ) .
The Shannon entropy, say H X ( 1 ) , follows from P X ( k ) as r 1 . Then, H X ( 1 ) follows for the Ex-LL distribution as
H X ( 1 ) = lim r 1 P X ( k ) = β α β log ( α β ) + ( β 1 ) Φ ( α ) log ( α ) + Υ ( β 1 ) + β + log ( λ ) β ,
where Φ ( z ) = d d z log [ Γ ( z ) ] and Υ refers to the Euler Mascheroni constant.

3.5. Order Statistics

The PDF of the k th order statistic, X k : n , for the Ex-LL distribution is defined by
f k : n ( x ) = n ! ( n k ) ! ( k 1 ) ! f ( x ) [ F ( x ) ] k 1 [ 1 F ( x ) ] n k = α 2 β n ! x β 1 λ α ( n + 1 k ) 1 λ + α x β α i α + α n + 1 1 λ λ + α x β α k 1 Γ ( k ) Γ ( n + 1 k ) .
The CDF of X k : n for the Ex-LL distribution takes the form
F k : n ( x ) = l = k n ( l n ) [ 1 F ( x ) ] n l [ F ( x ) ] l = n k 1 λ λ + α x β α k λ λ + α x β α ( n k ) × 2 F 1 1 , k n ; k + 1 ; 1 λ α x β + λ α ,
where 2 F 1 1 , k n ; k + 1 ; 1 λ α x β + λ α denotes the hyper geometric function.
The PDF and CDF of the minimum, say ( W n ) , and the maximum order statistics, say ( Z n ) , follows by setting k = 1 and k = n , respectively. The limiting distributions (Theorem 2.1.5 [24]) of W n and Z n are given by
lim n + P ( W n < d n x ) = 1 exp ( x β ) , d n = F 1 1 n
and
lim n + P ( Z n < b n x ) = exp ( x α β ) , b n = F 1 1 1 n .

4. Estimation Methods

Here, the Ex-LL parameters are estimated using some estimation approaches.

4.1. Maximum Likelihood Estimation

Let x 1 , x 2 , , x n be a random sample from the PDF (4), then the log-likelihood function takes the form
L = n log α 2 β λ α ( α + 1 ) i = 1 n log α x i β + λ + ( β 1 ) i = 1 n log x i .
The first derivatives with respect to α , β and λ follows by differentiating Equation (8) and equating them to zero. We obtain
L α = n λ α α 2 β λ α log ( λ ) + 2 α β λ α α 2 β ( 1 + α ) i = 1 n x i β α x i β + λ i = 1 n log α x i β + λ = 0 , L λ = α n λ ( α + 1 ) i = 1 n 1 α x i β + λ = 0 , L β = n β ( 1 + α ) i = 1 n α x i β log x i α x i β + λ + i = 1 n log x i = 0 .
Solving the previous equations, we obtain the MLEs of the Ex-LL parameters α , β and λ . However, these equations cannot be solved analytically, hence statistical software including Maple, R or Mathematica are adopted to solve them numerically using iterative methods such as Newton-Raphson algorithm.

4.2. Least-Squares and Weighted Least-Squares Estimation

Let x 1 : n , x 2 : n , , x 2 : n be the order statistics of the Ex-LL distribution. Hence, the LSEs of the Ex-LL parameters are provided by minimizing:
O = i = 1 n F ( x i : n ) i n + 1 2 = i = 1 n 1 λ α x i : n β + λ α i n + 1 2 ,
Besides, by solving the following equations we obtain the LSEs of the Ex-LL parameters:
i = 1 n 1 λ α x i : n β + λ α i n + 1 Δ s ( x i : n ) = 0 , s = 1 , 2 , 3 ,
where
Δ 1 ( x i : n ) = α F ( x i : n ) = λ α x i : n β + λ α log λ α x i : n β + λ + α x i : n β α x i : n β + λ ,
Δ 2 ( x i : n ) = λ F ( x i : n ) = α λ α x i : n β + λ α 1 λ α x i : n β + λ 2 1 α x i : n β + λ ,
Δ 3 ( x i : n ) = β F ( x i : n ) = α 2 λ α x i : n β log x i : n α x i : n β + λ 1 α α x i : n β + λ 2 .
The WLSEs of the Ex-LL parameters are calculated by minimizing:
W = i = 1 n ( n + 1 ) 2 ( n + 2 ) i ( n i + 1 ) F ( x i : n ) i n + 1 2 = i = 1 n ( n + 1 ) 2 ( n + 2 ) i ( n i + 1 ) 1 λ α x i : n β + λ α i n + 1 2 .
Besides, the WLSEs of α , β and λ can be obtained by solving:
i = 1 n ( n + 1 ) 2 ( n + 2 ) i ( n i + 1 ) F ( x i : n ) i n + 1 Δ s ( x i : n ) = 0 ,
where Δ s ( x i : n ) , s = 1 , 2 , 3 are given in (9)–(11), respectively.

4.3. Anderson—Darling Estimation

The ADEs of the Ex-LL parameters are calculated by minimizing:
A = n 1 n i = 1 n ( 2 i 1 ) [ log F ( x i : n ) + log S ( x i : n ) ] .
The ADEs are also obtained by solving:
i = 1 n ( 2 i 1 ) Δ s ( x i : n ) F ( x i : n ) Δ s ( x n + 1 i : n ) S ( x n + 1 i : n ) = 0 ,
where Δ s ( x i : n ) , s = 1 , 2 , 3 are given in (9)–(11), respectively.

4.4. Cramér—Von Mises Estimation

The CVMEs of the Ex-LL parameters are calculated by minimizing:
C V = 1 12 n + i = 1 n F ( x i : n ) 2 i 1 2 n 2 = 1 12 n + i = 1 n 1 λ α x i : n β + λ α 2 i 1 2 n 2 ,
These estimators can be also obtained by solving the nonlinear equations:
i = 1 n 1 λ α x i : n β + λ α 2 i 1 2 n Δ s ( x i : n ) = 0 ,
where Δ s ( x i : n ) , s = 1 , 2 , 3 are given in (9)–(11), respectively.

5. Numerical Simulations for the Estimation Methods

Now, we provided detailed simulation results to explore the performances of the introduced estimation methods in estimating the parameters of the Ex-LL model. We considered several sample sizes and different values of the parameters, that is, n = { 20 , 60 , 100 , 200 , 400 } and α = { 0.25 , 0.50 , 0.75 , 1.50 , 2.0 , 2.50 } , β = { 0.25 , 0.50 , 0.75 , 1.50 , 2.0 , 2.50 , 3.0 } , and λ = { 0.25 , 0.50 , 0.75 , 1.50 , 2.0 , 3.0 } . We generated N = 5000 random samples using Equation (7). The behavior of the different estimates is compared with respect to their: average absolute biases ( | B I A S | ), | B I A S | = 1 N i = 1 N | ϑ ^ ϑ | , mean square errors (MSEs), M S E s = 1 N i = 1 N ( ϑ ^ ϑ ) 2 , and mean relative errors (MREs), M R E s = 1 N i = 1 N | ϑ ^ ϑ | / ϑ .
The Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9 show the simulation results, average estimates of the parameters (AVEs), | B I A S | , MSEs, and MREs, of the Ex-LL parameters using different approaches. These results showed that estimates are very close to their true values and have small biases, MSEs and MREs. The results illustrated that the biases, MSEs, and MREs decrease as n increases, showing that the introduced estimators are consistent. We can conclude that the introduced estimation methods performed very well in estimating the Ex-LL parameters. Generally, based on our study, the ordering performance of these estimators, in terms of their MSEs, is the MLEs, WLSEs, ADEs, CVMEs, and LSEs. Hence, the maximum likelihood (ML) method is the best estimation approach to estimate the Ex-LL parameters and it will be adopted in the subsequent section to estimate the parameters of the Ex-LL model from real datasets.

6. Actuarial Measures

In this section, we discuss the mathematical and computational aspects of five important actuarial measures in portfolio optimization under uncertainty namely VaR, TVaR, TV, TVP, and ES for the Ex-LL distribution.

6.1. VaR Measure

The VaR has some other names such as the quantile premium principle or quantile risk measure, and it can be specified with a typical degree of confidence q, where q = 90 % , 95 % or 99 % ).
The VaR, say V a R q , of a random variable X is the q th QF of its CDF, i.e., V a R q = Q ( q ) (Artzner [25]).
The VaR of the Ex-LL model is derived as
V a R q = λ ( 1 q ) 1 / α 1 α 1 / β .

6.2. TVaR and TV Measures

The TVaR quantifies the expected value of loss given that an event outside a certain probability level has occurred. The TVaR is defined by the formula
T V a R q = 1 ( 1 q ) V a R q x f ( x ) d x .
The TVaR of the Ex-LL distribution has the form
T V a R q = α 1 α β λ α V a R q 1 α β 2 F 1 α + 1 , α 1 β ; α 1 β + 1 ; V a R q β λ α ( α β 1 ) ( 1 q ) ,
where 2 F 1 α + 1 , α 1 β ; α 1 β + 1 ; V a R q β λ α is the hyper geometric function.
The TV measure (Landsman [26]) represents the variance of a loss distribution beyond a particular critical value, and it pays attention to the TV beyond VaR. The TV of the Ex-LL distribution is defined as
T V q ( X ) = E ( X 2 | X > x q ) ( T V a R q ) 2 = 1 ( 1 q ) V a R q x 2 f ( x ) d x ( T V a R q ) 2 ,
where
E ( X 2 | X > x q ) = α 1 α β λ α V a R q 2 α β 2 F 1 α + 1 , α 2 β ; α 2 β + 1 ; V a R q β λ α ( α β 2 ) ( 1 q ) .
Using Equation (12)–(14), we obtain the TV of the Ex-LL distribution.

6.3. TVP and ES Measures

The TVP plays an important role in insurance field. The TVP of the Ex-LL distribution is defined by
T V P q ( x ) = T V a R q + λ T V q ,
where 0 < λ < 1 . The TVP of the Ex-LL distribution follows by substituting the expressions (12) and (13) in Equation (15).
Artzner [25] introduced another important measure of financial risk called the ES. The ES of the Ex-LL distribution takes the form
E S q ( x ) = 1 q 0 q V a R t d t = ( λ α ) 1 / β q 0 q ( 1 t ) 1 / α 1 1 / β d t .

6.4. Numerical Computations for Actuarial Measures

In this section, we presented numerical simulations for the studied risk measures, VaR, TVaR, TV, TVP, and ES, of the Ex-LL, ELL, and LL distributions. We generated a random sample, n = 100 , from the Ex-LL, ELL, and LL distributions, and the parameters were estimated by the ML method. The results were obtained after 1000 repetitions to calculate the five risk measures. The numerical results of these measures were reported in Table 10 and Table 11 for the Ex-LL, ELL, and LL distributions. For visual comparisons, we displayed the results, in Table 10 and Table 11, graphically as shown in Figure 3 and Figure 4.
The values and plots in Table 10 and Table 11 and Figure 3 and Figure 4 reveal that the introduced Ex-LL model has a heavier tail than the tails of the ELL and LL distributions. Hence, it may be adopted to model heavy-tailed datasets.

7. Modeling Real Data from the Engineering and Insurance Fields

In this section, we analyzed two real datasets from the engineering and insurance fields to explore the usefulness of the Ex-LL distribution. The first data was studied by [27]. It contains 74 observations and it refers to gauge lengths of 20 mm. This data was analyzed by [28,29].
The second data set refers to losses from private passenger automobile insurance policies in United Kingdom. It consists of 32 observations and 4 variables. We particularly analyzed the variable number 4 which represents number of claims. These data is available on R©software library.
The two data sets are displayed in Table 12 and Table 13.
The Ex-LL distribution is compared with some competing distributions including the alpha power log-logistic (APLL) [30], transmuted log-logistic (TLL) [31], generalized log-logistic (GLL) [32], Marshall–Olkin log-logistic (MOLL) [13], Poisson Burr-X log-logistic (PBXLL) [33], transmuted inverse log-logistic (TILL) [34], inverse log-logistic (ILL) [34], Weibull generalized log-logistic (WGLL) [35], and LL distributions.
The competing distributions are checked using some goodness-of-fit measures including Anderson–Darling (AD), Cramér–von Mises (CM), and Kolmogorov–Smirnov (KS) with its p-value (KS-p-value).
The parameters of the competing models are estimated via the ML method. The estimates and analytical measures are obtained using the Mathematica program version 12.0. Table 14 and Table 15 provide the analytical measures along with the ML estimates and their standard errors (SEs) in parenthesis. The fitted PDF, CDF, SF, and P-P plots of the Ex-LL model for the two data sets are shown in Figure 5 and Figure 6, respectively. The results in Table 14 and Table 15 indicate that the Ex-LL distribution provides adequate fit than other competing models for the two datasets.

8. Conclusions

In this article, we introduce a flexible extension of the log-logistic distribution called the extended log-logistic (Ex-LL) distribution. The EX-LL distribution can be adopted to model heavy-tailed actuarial data. The EX-LL distribution exhibits increasing, reversed-J, decreasing, upside-down bathtub, and J-shapes hazard rate functions. We derive some of its basic mathematical properties. Some risk measures are obtained for the Ex-LL distribution along with their detailed simulation results which illustrate that the tail of the Ex-LL distribution is heavier than the tails of the log-logistic and exponentiated log-logistic distributions. The parameters of the Ex-LL model are estimated using five classical estimation approaches and simulation results are conducted to explore their behavior for different sample sizes. The simulation results show that the maximum likelihood is the best estimation method for the Ex-LL parameters. The practical importance of the Ex-LL distribution is illustrated by two real data sets from the engineering and insurance sciences, showing its superiority fit as compared by nine competing distributions. We hope that the Ex-LL model will attract wider applications in other applied fields such as medicine, economics, reliability, life testing, and survival analyses.

Author Contributions

Investigation, N.M.A., A.M.G. and H.M.A.; Methodology, N.M.A., A.M.G. and H.M.A.; Project administration, A.Z.A.; Resources, A.Z.A.; Software, A.M.G. and A.Z.A.; Writing—original draft, A.M.G. and A.Z.A.; Writing—review & editing, H.M.A. and A.Z.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Taif University Researchers Supporting Project number (TURSP-2020/279), Taif University, Taif, Saudi Arabia.

Acknowledgments

The authors would like to thank the Editorial Board, and two reviewers for their constructive suggestions which greatly improved the final version of this manuscript. The authors extend their thanks to Taif University Researchers Supporting Project number (TURSP-2020/279), Taif University, Taif, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cooray, K.; Ananda, M.M.A. Modeling actuarial data with a composite lognormal-Pareto model. Scand. Actuar. J. 2005, 2005, 321–334. [Google Scholar] [CrossRef]
  2. Lane, M.N. Pricing risk transfer transactions 1. ASTIN Bull. J. IAA 2000, 30, 259–293. [Google Scholar] [CrossRef] [Green Version]
  3. Klugman, S.A.; Panjer, H.H.; Willmot, G.E. Loss Models: From Data to Decisions; John Wiley & Sons: Hoboken, NJ, USA, 2012; Volume 715. [Google Scholar]
  4. Ibragimov, R.; Prokhorov, A. Heavy Tails and Copulas: Topics in Dependence Modeling in Economics and Finance; World Scientific: Singapore, 2017. [Google Scholar]
  5. Punzo, A.; Mazza, A.; Maruotti, A. Fitting insurance and economic data with outliers: A flexible approach based on finite mixtures of contaminated gamma distributions. J. Appl. Stat. 2018, 45, 2563–2584. [Google Scholar] [CrossRef]
  6. Fisk, P.R. The graduation of income distributions. Econom. J. Econom. Soc. 1961, 29, 171–185. [Google Scholar] [CrossRef]
  7. Dagum, C. A model of income distribution and the conditions of existence of moments of finite order. Bull. Int. Stat. Inst. 1975, 46, 199–205. [Google Scholar]
  8. Shoukri, M.M.; Mian, I.U.H.; Tracy, D.S. Sampling properties of estimators of the log-logistic distribution with application to Canadian precipitation data. Can. J. Stat. 1988, 16, 223–236. [Google Scholar] [CrossRef]
  9. Arnold, B.C. Pareto Distributions Fairland; International Cooperative Publishing House: Silver Spring, MD, USA, 1983. [Google Scholar]
  10. Kleiber, C.; Kotz, S. Statistical Size Distributions in Economics and Actuarial Sciences; John Wiley & Sons: Hoboken, NJ, USA, 2003; Volume 470. [Google Scholar]
  11. De Santana, T.V.F.; Ortega, E.M.M.; Cordeiro, G.M.; Silva, G.O. The Kumaraswamy-log-logistic distribution. J. Stat. Theory Appl. 2012, 11, 265–291. [Google Scholar]
  12. Lemonte, A.J. The beta log-logistic distribution. Braz. J. Probab. Stat. 2014, 28, 313–332. [Google Scholar] [CrossRef]
  13. Gui, W. Marshall–Olkin extended log-logistic distribution and its application in minification processes. Appl. Math. Sci. 2013, 7, 3947–3961. [Google Scholar] [CrossRef]
  14. Tahir, M.H.; Mansoor, M.; Zubair, M.; Hamedani, G. McDonald log-logistic distribution with an application to breast cancer data. J. Stat. Theory Appl. 2014, 13, 65–82. [Google Scholar] [CrossRef] [Green Version]
  15. Hamedani, G. The Zografos-Balakrishnan log-logistic distribution: Properties and applications. J. Stat. Theory Appl. 2013, 12, 225–244. [Google Scholar]
  16. Cordeiro, G.M.; Afify, A.Z.; Ortega, E.M.; Suzuki, A.K.; Mead, M.E. The odd Lomax generator of distributions: Properties, estimation and applications. J. Comput. Appl. Math. 2019, 347, 222–237. [Google Scholar] [CrossRef]
  17. Zhao, W.; Khosa, S.K.; Ahmad, Z.; Aslam, M.; Afify, A.Z. Type-I heavy tailed family with applications in medicine, engineering and insurance. PLoS ONE 2020, 15, e0237462. [Google Scholar]
  18. Lorenz, M.O. Methods of measuring the concentration of wealth. Publ. Am. Stat. Assoc. 1905, 9, 209–219. [Google Scholar] [CrossRef]
  19. Bonferroni, C.E. Elementi di Statistica Generale; Seeber: Firenze, Italy, 1930. [Google Scholar]
  20. Arcagni, A.; Porro, F. The graphical representation of inequality. Rev. Colomb. De Estad. 2014, 37, 419–437. [Google Scholar] [CrossRef] [Green Version]
  21. Morales, D.; Pardo, L.; Vajda, I. Some new statistics for testing hypotheses in parametric models. J. Multivar. Anal. 1997, 62, 137–168. [Google Scholar] [CrossRef] [Green Version]
  22. Kurths, J.; Voss, A.; Saparin, P.; Witt, A.; Kleiner, H.J.; Wessel, N. Quantitative analysis of heart rate variability. Chaos Interdiscip. J. Nonlinear Sci. 1995, 5, 88–94. [Google Scholar] [CrossRef]
  23. Song, K.S. Rényi information, loglikelihood and an intrinsic distribution measure. J. Stat. Plan. Inference 2001, 93, 51–69. [Google Scholar] [CrossRef]
  24. Galambos, J. The Asymptotic Theory of Extreme Order Statistics; R.E. Krieger Pub. Co.: Malabar, FL, USA, 1987. [Google Scholar]
  25. Artzner, P. Application of coherent risk measures to capital requirements in insurance. N. Am. Actuar. J. 1999, 3, 11–25. [Google Scholar] [CrossRef]
  26. Landsman, Z. On the tail mean—variance optimal portfolio selection. Insur. Math. Econ. 2010, 46, 547–553. [Google Scholar] [CrossRef]
  27. Kundu, D.; Raqab, M.Z. Estimation of R = P(Y < X) for three parameter Weibull distribution. Stat. Probab. Lett. 2009, 79, 1839–1846. [Google Scholar]
  28. Afify, A.; Cordeiro, G.; Butt, N.; Ortega, E.; Suzuki, A. A new lifetime model with variable shapes for the hazard rate. Braz. J. Probab. Stat. 2017, 31, 516–541. [Google Scholar] [CrossRef]
  29. Afify, A.; Mohamed, O. A new three-parameter exponential distribution with variable shapes for the hazard rate: Estimation and applications. Mathematics 2020, 8, 135. [Google Scholar] [CrossRef] [Green Version]
  30. Aldahlan, M.A. Alpha power transformed log-logistic distribution with application to breaking stress data. Adv. Math. Phys. 2020, 2020, 2193787. [Google Scholar] [CrossRef]
  31. Granzotto, D.C.T.; Louzada, F. The transmuted log-logistic distribution: Modeling, inference, and an application to a polled tabapua race time up to first calving data. Commun. Stat. Theory Methods 2015, 44, 3387–3402. [Google Scholar] [CrossRef]
  32. Adeyinka, F.S.; Olapade, A.K. On transmuted four parameters generalized log-logistic distribution. Int. J. Stat. Distrib. Appl. 2019, 5, 32–37. [Google Scholar]
  33. Almamy, J.A. Extended Poisson log-logistic distribution. Int. J. Stat. Probab. 2019, 8, 56–69. [Google Scholar] [CrossRef]
  34. Para, B.A.; Jan, T.R. Transmuted inverse log logistic model: Properties and application in medical sciences and engineering. Math. Theory Model. 2017, 7, 157–181. [Google Scholar]
  35. Abouelmagd, T.H.M.; Hamed, M.S.; Almamy, J.A.; Ali, M.M.; Yousof, H.M.; Korkmaz, M.C. Extended Weibull log-logistic distribution. J. Nonlinear Sci. Appl. 2019, 12, 523–534. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Possible shapes of the Ex-LL density function for various values of α , β and λ .
Figure 1. Possible shapes of the Ex-LL density function for various values of α , β and λ .
Mathematics 09 01386 g001
Figure 2. Possible shapes of the Ex-LL hazard function for various values of α , β and λ .
Figure 2. Possible shapes of the Ex-LL hazard function for various values of α , β and λ .
Mathematics 09 01386 g002aMathematics 09 01386 g002b
Figure 3. Shapes of the VaR, TVaR, TV, TVP, and ES using the numerical values in Table 10.
Figure 3. Shapes of the VaR, TVaR, TV, TVP, and ES using the numerical values in Table 10.
Mathematics 09 01386 g003
Figure 4. Shapes of the VaR, TVaR, TV, TVP, and ES using the numerical values in Table 11.
Figure 4. Shapes of the VaR, TVaR, TV, TVP, and ES using the numerical values in Table 11.
Mathematics 09 01386 g004
Figure 5. Histogram of data set I along with the fitted functions of the Ex-LL model.
Figure 5. Histogram of data set I along with the fitted functions of the Ex-LL model.
Mathematics 09 01386 g005
Figure 6. Histogram of data set II along with the fitted functions of the Ex-LL model.
Figure 6. Histogram of data set II along with the fitted functions of the Ex-LL model.
Mathematics 09 01386 g006
Table 1. The numerical values for μ 1 , V a r ( X ) , γ 1 ( X ) , and γ 2 ( X ) of the Ex-LL distribution for several choices of α , β , and λ .
Table 1. The numerical values for μ 1 , V a r ( X ) , γ 1 ( X ) , and γ 2 ( X ) of the Ex-LL distribution for several choices of α , β , and λ .
Parameters μ 1 Var ( X ) γ 1 ( X ) γ 2 ( X )
( α = 0.5 , β = 15 , λ = 1.5 ) 1.198520.051222.1318115.34159
( α = 0.75 , β = 7.5 , λ = 0.75 ) 1.112440.113142.3446222.42470
( α = 1.5 , β = 3 , λ = 3 ) 1.177750.438112.8218349.06847
( α = 1.5 , β = 3 , λ = 0.75 ) 0.741940.173862.8218349.06840
( α = 1.75 , β = 15 , λ = 2.5 ) 0.974400.009770 0.05100 3.75522
( α = 2 , β = 9 , λ = 3 ) 0.949000.023780.109113.67376
( α = 2 , β = 6 , λ = 3 ) 0.933670.051040.434494.10571
( α = 3 , β = 1.5 , λ = 2 ) 0.410120.144783.5818263.19551
( α = 3 , β = 1.5 , λ = 5 ) 0.755460.491243.581863.19552
( α = 3.75 , β = 15 , λ = 2.5 ) 0.869250.00610 0.46932 3.66150
( α = 5 , β = 1.5 , λ = 1.5 ) 0.156100.015601.9938711.18786
( α = 10 , β = 10 , λ = 2.5 ) 0.825170.01221 0.27530 3.39749
( α = 15 , β = 0.5 , λ = 0.75 ) 1.1  × 10 6 2.4  × 10 9 45.482822573.493
( α = 15 , β = 25 , λ = 5 ) 0.841420.00184 0.84557 4.31163
Table 2. Numerical values of the Ex-LL distribution for the parameters α = 0.25 , β = 0.75 , λ = 0.5 .
Table 2. Numerical values of the Ex-LL distribution for the parameters α = 0.25 , β = 0.75 , λ = 0.5 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
200.455440.823220.428820.200920.578170.207490.058070.386380.433970.401840.770900.82998
600.497670.854030.278820.135190.467440.056230.025820.281130.010140.270390.623260.22491
MLEs1000.504700.844750.263610.109010.383990.036130.017090.207150.002500.218010.511980.14453
2000.505990.823690.255480.081500.302280.024380.009840.140930.001020.163000.403040.09753
4000.504110.795500.252950.062340.224060.018470.006010.083000.000560.124690.298750.07387
200.470060.850470.365500.179610.552060.153210.046720.363120.320440.359230.736090.61282
600.497280.857960.272380.132830.459310.052220.024690.272270.005360.265650.612410.20888
ADEs1000.497940.832690.265120.112800.398680.039580.017860.220220.003060.225610.531570.15832
2000.506830.832780.255930.086940.324440.026280.011330.159560.001160.173880.432580.10510
4000.511330.817060.251030.067480.246310.019210.006920.096110.000590.134960.328410.07682
200.437190.832500.454110.203640.591800.234530.059340.399880.367340.407290.789070.93813
600.478160.842290.296550.148410.500600.074140.030610.307160.015610.296810.667470.29658
CVMEs1000.495900.840210.269770.123050.435320.046920.020870.252720.004450.246110.580430.18766
2000.497690.809710.261280.100320.357000.032280.014070.181180.001850.200630.476010.12913
4000.507190.808580.252990.077330.280060.022370.008810.122010.000810.154650.373410.08949
200.458470.825780.372920.206370.589050.171740.058440.395510.161330.412740.785400.68696
600.486780.816180.280390.150450.480060.064980.031250.291380.010550.300910.640080.25991
LSEs1000.502860.848740.264540.122200.432610.044490.020820.249670.004030.244400.576810.17796
2000.505440.826220.257900.097670.347950.031390.013830.174680.001740.195340.463940.12557
4000.505570.802990.252800.073080.265560.021060.008280.113980.000750.146160.354070.08425
200.421280.629190.373460.166180.367690.161540.042170.176870.204310.332360.490260.64616
600.460210.685790.280470.101230.300330.053670.016680.117340.009090.202460.400440.21467
WLSEs1000.473830.698040.265720.087360.271350.037020.011950.094990.002760.174720.361800.14808
2000.481500.717480.259400.067350.223100.025590.006910.064950.001180.134690.297470.10237
4000.489880.731400.255190.054570.185470.017900.004280.045740.000540.109150.247290.07159
Table 3. Numerical values of the Ex-LL distribution for the parameters α = 1.5 , β = 0.5 , λ = 0.75 .
Table 3. Numerical values of the Ex-LL distribution for the parameters α = 1.5 , β = 0.5 , λ = 0.75 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
201.411530.527910.927370.529990.371350.239620.376760.161230.133720.353330.742710.31950
601.481450.557510.811530.406450.344240.113640.211490.141250.028640.270970.688480.15152
MLEs1001.501030.563980.785530.359530.316800.079820.167580.125070.012600.239690.633600.10643
2001.504310.548740.770720.294130.271350.057850.115290.099630.005960.196090.542700.07713
4001.546520.580120.754780.246130.239910.041020.083720.082600.002760.164090.479810.05469
201.311100.473180.906980.514270.373890.230580.366180.162950.113410.342850.747780.30744
601.444580.544310.806160.411580.342720.111760.225920.141160.027060.274390.685440.14901
ADEs1001.467570.547770.785390.372950.328570.085490.178880.130980.013740.248640.657140.11399
2001.490990.547530.767320.311290.283370.059730.128040.105960.006050.207520.566740.07964
4001.513910.551800.758760.252250.239140.043590.089400.082300.003080.168170.478270.05812
201.272580.470910.974730.569180.397460.294620.443070.177600.190270.379450.794920.39282
601.442400.558260.817070.443890.367430.125050.258990.156890.033410.295930.734870.16674
CVMEs1001.458900.556260.800780.416120.352770.100530.222590.146690.021850.277410.705540.13404
2001.477050.550320.777500.360310.316800.068890.166150.124630.008860.240210.633600.09186
4001.495390.548160.767660.304430.273810.050740.123640.100770.004630.202950.547630.06766
201.258550.489200.908730.550360.404740.258490.428800.182980.143530.366910.809480.34465
601.402720.544300.810880.454380.374630.131350.274030.161060.038220.302920.749250.17513
LSEs1001.425170.538750.793240.415600.352210.099480.224500.146370.020110.277070.704420.13265
2001.467600.552550.774570.360100.319550.071630.168700.126820.009790.240070.639090.09550
4001.501270.555800.762570.303600.277330.048430.121200.102840.004140.202400.554660.06457
201.281690.486530.882150.542490.392460.222830.403960.173760.117420.361660.784920.29711
601.375710.509670.804920.438800.355960.115270.254210.148860.029460.292540.711920.15369
WLSEs1001.468300.553190.779850.374770.329930.083290.180070.132560.013140.249850.659860.11105
2001.493210.554310.769600.318530.291480.060660.133400.109880.006630.212350.582950.08088
4001.516730.557440.759320.266040.249560.044600.097540.087780.003250.177360.499120.05946
Table 4. Numerical values of the Ex-LL distribution for the parameters α = 2 , β = 1.5 , λ = 0.25 .
Table 4. Numerical values of the Ex-LL distribution for the parameters α = 2 , β = 1.5 , λ = 0.25 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
201.853141.489990.311580.737851.024710.080480.746101.158260.017360.368930.683140.32193
601.900101.511240.273750.570380.911810.037310.436090.951680.002940.285190.607870.14923
MLEs1001.929901.529210.264070.493150.841420.027200.322310.823520.001390.246580.560950.10879
2002.006221.617200.256490.397650.720070.017270.211500.644930.000560.198830.480040.06909
4001.994471.567950.254600.339970.631420.013380.153060.513630.000300.169980.420950.05354
201.722001.386960.297400.698041.022080.073230.706561.165390.012840.349020.681380.29294
601.828691.446910.268580.562240.908900.036730.436590.947010.002800.281120.605940.14693
ADEs1001.870631.454460.262880.505960.837930.026890.347330.829770.001410.252980.558620.10754
2001.919451.495140.257410.429400.762740.018870.242650.697130.000630.214700.508500.07549
4001.956981.513780.254610.362800.651600.013260.173090.538830.000290.181400.434400.05304
201.616671.272730.336390.813591.071180.107980.938121.290570.028620.406790.714120.43192
601.787051.422090.281170.661401.002750.047120.592671.115120.004940.330700.668500.18848
CVMEs1001.846971.460060.268190.571220.923110.033300.440080.965530.002370.285610.615410.13318
2001.906681.513280.259880.488340.834060.022290.314000.812930.000940.244170.556040.08916
4001.932951.507410.255680.411970.730340.015720.223060.648110.000420.205990.486900.06286
201.600991.321670.309460.769271.054510.092920.892641.255850.021880.384640.703010.37167
601.787121.440560.271760.638290.990050.043110.559171.086960.003950.319150.660030.17243
LSEs1001.838451.466570.265260.560440.922000.031260.426320.961970.002090.280220.614670.12505
2001.871121.473760.260700.501980.847590.023080.333930.840250.001060.250990.565060.09230
4001.956001.549850.254850.412930.733540.015680.222280.655670.000430.206470.489020.06273
201.620261.331650.301140.744601.061070.080860.822601.243710.016520.372300.707380.32345
601.791921.413440.268630.602380.950580.038330.495141.012110.003070.301190.633720.15332
WLSEs1001.869711.478010.261660.514540.865380.028500.362430.873110.001610.257270.576920.11401
2001.937411.540870.256360.443500.787290.019050.255770.730880.000650.221750.524860.07619
4001.978941.565570.253810.359770.656640.013940.171460.550770.000330.179880.437760.05575
Table 5. Numerical values of the Ex-LL distribution for the parameters α = 2.5 , β = 3 , λ = 1.5 .
Table 5. Numerical values of the Ex-LL distribution for the parameters α = 2.5 , β = 3 , λ = 1.5 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
202.245602.649031.861980.772881.472820.468440.980412.898010.570680.309150.490940.31230
602.257292.631471.618630.597371.369520.191910.562212.389950.079310.238950.456510.12794
MLEs1002.338102.746891.578060.526031.263380.144490.419492.009360.040520.210410.421130.09632
2002.424852.935351.539020.414831.084890.093200.253311.469110.015630.165930.361630.06213
4002.425922.899771.527150.364390.973870.070370.187361.200440.008670.145760.324620.04691
201.966432.222081.818090.865041.673170.461061.177723.604220.536680.346010.557720.30737
602.159672.493011.634610.644381.453520.222800.678192.704680.105430.257750.484510.14853
ADEs1002.232782.601751.583870.564461.339840.155380.518672.322780.048370.225780.446610.10358
2002.312632.721851.546400.474811.189990.109900.351071.823510.022180.189920.396660.07327
4002.369092.795131.525760.395471.039670.074020.230151.381180.009670.158190.346560.04935
201.919102.211652.015480.939241.626480.622341.444213.736310.951560.375700.542160.41489
602.058562.357791.693670.726071.545650.274750.881493.140120.192110.290430.515220.18317
CVMEs1002.140482.446081.628590.663681.488030.200090.701782.822500.081010.265470.496010.13339
2002.275892.668181.564180.530371.278970.128950.440652.099300.030150.212150.426320.08596
4002.352262.793801.536310.447411.149700.087340.291891.645280.013350.178960.383230.05823
201.851722.189121.869520.914481.605530.567181.423353.730970.780030.365790.535180.37812
602.062102.403491.645210.723441.522990.265790.885583.081530.153340.289380.507660.17719
LSEs1002.152662.498611.594540.630511.450320.189400.648172.699590.080130.252210.483440.12627
2002.218302.563881.566150.561821.331450.131330.499942.288950.033680.224730.443820.08755
4002.311872.700321.533890.457191.152810.084930.314981.701770.013030.182880.384270.05662
201.874512.171241.841220.930471.741850.498801.375063.874230.680650.372190.580620.33253
602.099052.389311.638690.692611.520410.232530.772232.965590.121640.277040.506800.15502
WLSEs1002.195622.520311.571880.590811.381960.155740.543072.405930.044810.236330.460650.10383
2002.286382.650471.551640.499051.241210.113690.376381.942730.024050.199620.413740.07579
4002.372142.793791.527990.399521.050130.075440.231211.395080.009970.159810.350040.05030
Table 6. Numerical values of the Ex-LL distribution for the parameters α = 0.5 , β = 2 , λ = 3 .
Table 6. Numerical values of the Ex-LL distribution for the parameters α = 0.5 , β = 2 , λ = 3 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
200.502432.243413.282240.137480.943060.678560.028811.026690.582930.274960.471530.22619
600.486962.112863.211490.108200.788430.500110.016520.768280.363280.216400.394220.16670
MLEs1000.494482.072433.147150.096610.688490.412550.013300.614900.262790.193220.344240.13752
2000.498512.076413.078940.076510.573910.301770.008490.455470.149120.153020.286950.10059
4000.502122.077113.040370.059950.477000.215110.005470.330870.076030.119900.238500.07170
200.442742.011934.201960.168951.003391.607890.042911.1899915.310040.337910.501690.53596
600.481392.057363.278010.116000.809840.600480.020030.806570.786940.232010.404920.20016
ADEs1000.496312.116803.157420.096980.734660.437980.013990.679130.387490.193960.367330.14599
2000.494082.060023.093910.075800.603460.303840.008600.491720.165490.151600.301730.10128
4000.504342.083533.031780.061620.479800.227240.005700.332260.085750.123240.239900.07575
200.471912.154013.344050.140590.974360.737880.029101.097710.665180.281180.487180.24596
600.484882.135093.247550.114220.838680.560690.018280.843670.427950.228440.419340.18690
CVMEs1000.489092.074223.167340.103940.783300.448530.015360.756470.305840.207880.391650.14951
2000.492442.070163.118850.083380.655960.349370.010040.561220.198390.166760.327980.11646
4000.505622.108753.037110.067240.516060.253950.006680.376350.104480.134470.258030.08465
200.484502.009723.200690.157331.014020.724150.035161.181990.653210.314660.507010.24138
600.486412.004713.153640.120650.862350.540570.020680.894730.409180.241300.431180.18019
LSEs1000.490402.036663.136910.106530.766660.458690.015950.727570.313290.213050.383330.15290
2000.495552.038123.095140.089970.678090.349660.011390.586590.195120.179940.339050.11655
4000.502412.074923.045810.070530.537690.269440.007290.403410.119770.141050.268840.08981
200.489961.972733.142950.153840.975970.705150.034621.116210.620330.307680.487990.23505
600.493692.081733.140920.110130.789460.493560.017500.772940.354240.220260.394730.16452
WLSEs1000.495012.053563.094850.097710.719780.413660.013560.657210.255560.195430.359890.13789
2000.496022.036503.071050.076120.598060.295130.008370.475200.142430.152240.299030.09838
4000.500162.058933.041930.060490.474670.223670.005500.328860.083100.120980.237340.07456
Table 7. Numerical values of the Ex-LL distribution for the parameters α = 0.75 , β = 0.25 , λ = 2 .
Table 7. Numerical values of the Ex-LL distribution for the parameters α = 0.75 , β = 0.25 , λ = 2 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
200.834920.462972.249330.331190.371600.494960.147170.213100.362540.441591.486420.24748
600.818410.405492.089860.244510.286050.324780.089610.147440.169900.326021.144210.16239
MLEs1000.793250.365022.091120.215100.244910.286370.070330.113620.139630.286800.979620.14318
2000.786360.326512.026680.151420.171650.186650.038760.062340.059590.201890.686610.09333
4000.765200.284872.020400.106560.114510.136710.018700.025520.031350.142080.458030.06836
200.798890.435822.212860.319130.361830.508470.135920.200870.371460.425501.447310.25424
600.795110.380072.109990.248570.277580.335940.090530.138520.182400.331431.110320.16797
ADEs1000.792540.366122.077090.211960.245460.273750.069340.115010.125170.282610.981840.13688
2000.786890.325952.023690.151820.170200.184310.038800.060730.056130.202430.680800.09215
4000.766880.289612.015990.112010.121210.143620.020590.028910.033190.149350.484850.07181
200.824720.465332.268760.336540.387870.549580.149410.224900.427980.448721.551480.27479
600.802020.402892.159380.275910.312660.380800.105540.166720.234980.367881.250620.19040
CVMEs1000.808910.392742.099500.250050.281030.317470.089810.139850.169160.333401.124130.15874
2000.792650.353822.058500.194090.219920.240600.059960.097390.098820.258790.879670.12030
4000.774110.305892.021150.136850.149920.160240.029920.045060.041030.182460.599690.08012
200.793700.450122.168790.330670.377390.536750.141760.215630.403260.440891.509570.26837
600.784990.386492.126630.271310.301380.388990.101370.154670.236660.361741.205540.19449
LSEs1000.804860.390982.063590.243470.277410.311880.086360.137760.159020.324631.109640.15594
2000.781740.334502.042870.181780.199880.224080.053140.080860.084150.242370.799510.11204
4000.782160.316062.012070.139030.154620.163890.032440.050260.043810.185370.618480.08195
200.814420.467102.163860.331180.386940.512760.144300.226230.373330.441581.547760.25638
600.806060.408892.081380.264300.304030.350640.099410.158480.200940.352401.216120.17532
WLSEs1000.790000.359232.062130.213680.240010.282310.069980.109460.130100.284910.960030.14115
2000.776230.320282.034880.159350.178060.200800.040680.063830.066550.212470.712250.10040
4000.767260.289262.014690.108640.117950.139610.019830.028920.030800.144850.471790.06980
Table 8. Numerical values of the Ex-LL distribution for the parameters α = 0.5 , β = 2.5 , λ = 0.75 .
Table 8. Numerical values of the Ex-LL distribution for the parameters α = 0.5 , β = 2.5 , λ = 0.75 .
AVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
200.461532.712770.954550.189051.337910.301810.049322.089630.158800.378090.535170.40241
600.483632.640100.834400.127211.092990.165350.024271.496590.054320.254410.437200.22047
MLEs1000.501332.686080.788930.105270.946190.113570.016601.165920.023340.210540.378480.15143
2000.504992.702070.773990.083250.763980.080870.010130.823620.012010.166500.305590.10782
4000.501722.611690.762540.059460.576680.054620.005440.516920.004940.118920.230670.07282
200.472052.765430.897520.171821.297640.260270.042501.964570.119970.343650.519060.34702
600.490222.697520.818840.127081.110960.156570.023151.517540.047350.254160.444380.20877
ADEs1000.495422.642610.791740.107480.973520.114820.016831.231720.024100.214960.389410.15310
2000.502142.632230.767150.084710.797130.081820.010750.874690.011630.169430.318850.10910
4000.503912.604250.759340.062250.579410.056070.005850.510210.005080.124510.231760.07475
200.453662.782920.956770.187051.356600.307730.047642.112190.158950.374100.542640.41030
600.483932.799420.846240.135551.147990.181610.026321.610910.065530.271090.459200.24214
CVMEs1000.489052.678780.812340.115201.011290.137160.019381.303450.036670.230400.404510.18289
2000.495582.641280.785360.093940.851060.099810.013021.000940.018550.187870.340430.13308
4000.504742.635470.763200.074430.682580.068200.008300.686280.007780.148870.273030.09093
200.464852.487930.899690.200451.408150.294100.053782.294460.144250.400890.563260.39213
600.479642.606320.834850.143231.179690.182240.028761.679400.065660.286460.471870.24298
LSEs1000.489332.612500.806030.120781.049000.137040.020531.372010.036650.241570.419600.18272
2000.500762.623960.775560.096390.857700.097880.013460.978960.017490.192780.343080.13050
4000.506052.631790.757870.072620.682150.065070.007860.685520.006760.145230.272860.08676
200.480812.621030.877450.192721.352980.274480.051682.136990.132280.385450.541190.36598
600.487622.612030.812170.128271.103420.151560.024111.518950.044820.256540.441370.20208
WLSEs1000.501432.686940.785200.108230.997370.118200.017091.264560.026370.216460.398950.15760
2000.505812.697350.769130.087090.820570.083840.011090.919420.011600.174190.328230.11178
4000.506392.608100.754960.060650.586350.053950.005710.523270.004620.121310.234540.07194
Table 9. Numerical values of the Ex-LL distribution for the parameters α = 2 , β = 0.75 , λ = 1.5 .
Table 9. Numerical values of the Ex-LL distribution for the parameters α = 2 , β = 0.75 , λ = 1.5 .
MethodnAVEs|BIAS|MSEsMREs
α β λ α β λ α β λ α β λ
202.020560.871231.752680.753610.600400.374410.721760.406770.246700.376800.800530.24961
602.045960.906061.582180.593710.550200.189820.453350.358240.068280.296850.733600.12655
MLEs1002.035510.884571.561590.531180.509360.147180.367240.318730.043480.265590.679140.09812
2002.034330.859371.534110.458600.450670.107240.273910.265480.020190.229300.600890.07150
4002.055910.865501.511650.364120.375040.076500.184310.200530.009720.182060.500060.05100
201.597650.536851.760360.660740.432790.382110.665280.236660.254940.330370.577060.25474
601.734840.614431.612130.488790.369010.203260.390900.174900.082940.244400.492020.13550
ADEs1001.809130.650801.576230.405310.331060.153080.270710.141460.047610.202650.441410.10205
2001.865910.676831.546180.334650.292950.106950.177780.110820.021070.167320.390590.07130
4001.900850.693791.530280.287380.260200.076600.124180.087400.010730.143690.346930.05107
201.600210.539821.840330.711410.434410.457490.768880.244320.336330.355710.579220.30499
601.700610.596761.692930.563190.391830.277270.514220.199220.152180.281600.522430.18485
CVMEs1001.762600.628551.615110.473240.365250.190520.365680.170170.073080.236620.487000.12701
2001.794830.640841.580120.412770.336410.136740.274980.145440.038130.206390.448550.09116
4001.873480.681781.539100.321800.282660.085950.160650.102930.014380.160900.376880.05730
201.582410.584881.695640.651380.415160.376370.683320.223050.244300.325690.553540.25091
601.685500.602711.636530.530910.382730.235330.472550.189630.111020.265460.510310.15689
LSEs1001.746430.626111.594270.470120.365680.182900.366830.169590.066560.235060.487580.12193
2001.800750.645441.560000.396260.328860.126460.254670.138000.031480.198130.438480.08431
4001.885730.697501.532860.326040.287310.086980.162300.104550.013680.163020.383080.05799
201.613200.584531.697650.643980.412600.361270.664790.219890.233510.321990.550130.24085
601.733190.617951.610110.486550.367080.207330.394130.173570.084450.243270.489440.13822
WLSEs1001.752030.613161.594440.440460.348360.166680.318730.158080.053980.220230.464480.11112
2001.868350.679691.550630.348660.302530.108770.187610.115740.021040.174330.403380.07251
4001.924920.713591.524270.273310.252800.072290.111140.081220.009180.136650.337070.04819
Table 10. Simulation results of the five risk measures for the Ex-LL, ELL, and LL distributions.
Table 10. Simulation results of the five risk measures for the Ex-LL, ELL, and LL distributions.
DistributionParametersSignificance LevelVaRTVaRTVTVPES
Ex-LL α = 0.9 , λ = 0.5 , β = 2.5 0.600.991082.0126617.0703612.254870.58926
0.651.082792.1522119.3495714.729420.62361
0.701.190842.3217122.3683417.979550.66015
0.751.323152.5351526.5612322.456070.69977
0.801.493962.8178132.7900429.049840.74382
0.851.732443.2218943.0447639.809930.79451
0.902.112683.8803363.2138660.772800.85620
0.952.920955.30891122.11669121.319770.94018
LL λ = 0.5 , β = 2.5 0.600.891251.700533.513463.808610.53540
0.650.970231.810643.915544.355740.56574
0.701.062441.943244.441145.052040.59783
0.751.174212.108625.159655.978360.63240
0.801.316762.325226.206227.290190.67052
0.851.512872.630687.886119.333870.71395
0.901.819553.1197211.0764113.088490.76609
0.952.452664.1524819.8766723.035310.83561
ELL λ = 0.5 , β = 2.5 , c = 0.75 0.600.766271.499884.520124.211950.43416
0.650.839111.599595.082344.903110.46245
0.700.923661.719495.823195.795720.49229
0.751.025541.868826.845767.003140.52433
0.801.154752.064218.353018.746620.55951
0.851.331602.3396610.8095911.527810.59941
0.901.606922.7809315.5747316.798190.64706
0.952.173693.7154829.1837931.440080.71019
Table 11. Simulation results of the five risk measures for the Ex-LL, ELL, and LL distributions.
Table 11. Simulation results of the five risk measures for the Ex-LL, ELL, and LL distributions.
DistributionParametersSignificance LevelVaRTVaRTVTVPES
Ex-LL α = 0.5 , λ = 1.5 , β = 5 0.601.731162.9438910.753309.395871.23942
0.651.840313.1095412.0610310.949211.28134
0.701.969503.3106513.7755912.953571.32574
0.752.128593.5635816.1283615.659851.37375
0.802.335183.8977719.5723019.555611.42707
0.852.625164.3734225.1369825.739851.48847
0.903.088905.1423435.8059737.367721.56331
0.954.070346.7832465.7172369.214611.66531
LL λ = 1.5 , β = 5 0.601.175751.561930.204701.684750.89085
0.651.226901.613510.212341.751530.91470
0.701.284051.673280.222331.828910.93900
0.751.350081.744700.235661.921450.96415
0.801.429871.833740.254172.037080.99067
0.851.532861.952050.281652.191461.01939
0.901.681282.127270.327732.422221.05173
0.951.952202.455470.429382.863381.09106
ELL λ = 1.5 , β = 5 , c = 1.5 0.601.547622.550353.774614.815121.05488
0.651.652982.686244.158885.389511.09674
0.701.775112.848554.657636.108901.14072
0.751.921913.049065.333407.049111.18774
0.802.107243.308706.306838.354161.23916
0.852.358943.669697.8472010.339811.29719
0.902.745554.2366510.7163513.881361.36600
0.953.520445.3980918.3937622.872161.45596
Table 12. The observations of gauge lengths of 20 mm data.
Table 12. The observations of gauge lengths of 20 mm data.
1.3121.3141.4791.5521.7001.8031.8611.8651.9441.958
1.9661.9972.0062.0212.0272.0552.0632.0982.1402.179
2.2242.2402.2532.2702.2722.2742.3012.3012.3592.382
2.4262.4342.4352.3822.4782.5542.5142.5112.4902.535
2.5662.5702.5862.6292.8002.7732.7702.8093.5852.818
2.6422.7262.6972.6842.6482.6333.1283.0903.0963.233
2.8212.8802.8482.8183.0672.8212.9542.8093.5853.084
3.0122.8802.8483.433
Table 13. The losses from private passenger automobile insurance policies data.
Table 13. The losses from private passenger automobile insurance policies data.
2140235631719244140343
318129123448361169151479381166
245970719304266859504162260578
31296
Table 14. Measures of goodness-of-fit and estimates of the Ex-LL distribution and other distributions for data set I.
Table 14. Measures of goodness-of-fit and estimates of the Ex-LL distribution and other distributions for data set I.
ModelADCMKSKS-p-ValueEstimates (SEs)
Ex-LL0.187080.024540.053070.98524 α ^ = 5.35622 ( 6.59867 )
λ ^ = 14044.4 ( 27477.7 )
β ^ = 6.44473 ( 0.98951 )
LL0.560270.067560.059310.95704 λ ^ = 2407.73 ( 1907.88 )
β ^ = 8.63147 ( 0.83741 )
APLL0.244470.033370.054320.98112 α ^ = 0.0025 ( 0.0054 )
a ^ = 7.00987 ( 0.5545 )
b ^ = 3.32516 ( 0.22732 )
TLL0.281210.039380.058080.96415 α ^ = 2.79405 ( 0.12052 )
β ^ = 7.39915 ( 0.27229 )
λ ^ = 1.0000 ( 0.37615 )
GLL28.25656.11810.523160.0000 α ^ = 0.41192 ( 0.10455 )
β ^ = 1.00000 ( 0.11506 )
θ ^ = 0.51734 ( 0.08067 )
MOLL0.560270.067560.059310.95704 c ^ = 0.68853 ( 747.74 )
λ ^ = 2.57367 ( 323.811 )
β ^ = 8.63147 ( 0.83741 )
PBXLL1.106820.176310.091590.56396 c ^ = 24.5304 ( 4.10689 )
λ ^ = 5.84672 × 10 8 ( 0.00018 )
β ^ = 0.72432 ( 0.02922 )
TILL29.94006.559660.515790.0000 λ ^ = 1.00000 ( 0.29515 )
α ^ = 1.96889 ( 0.16327 )
ILL54.745811.50140.665890.0000 α ^ = 1.71725 ( 0.15506 )
WGLL0.214430.026480.059620.95509 c ^ = 0.28789 ( 0.23415 )
λ ^ = 2.36831 ( 2.22062 )
β ^ = 4.55449 ( 1.61755 )
Table 15. Measures of goodness-of-fit and estimates of the Ex-LL distribution and other distributions for data set II.
Table 15. Measures of goodness-of-fit and estimates of the Ex-LL distribution and other distributions for data set II.
ModelADCMKSKS-p-ValueEstimates (SEs)
Ex-LL0.146860.023320.077210.99108 α ^ = 28.8527 ( 63.2613 )
λ ^ = 558573 ( 2.27773 × 10 6 )
β ^ = 1.15130 ( 0.17685 )
LL0.438280.048240.087130.96832 λ ^ = 4708.8 ( 6204.98 )
β ^ = 1.60132 ( 0.23852 )
APLL0.438120.048160.087040.96861 α ^ = 1.3504 ( 7.91583 )
a ^ = 1.5993 ( 0.25100 )
b ^ = 179.106 ( 326.737 )
TLL0.438160.048180.087060.96853 α ^ = 187.685 ( 198.581 )
β ^ = 1.59980 ( 0.24797 )
λ ^ = 0.07504 ( 1.68415 )
GLL6.058791.230070.347400.00088 α ^ = 12.7009 ( 9.22939 )
β ^ = 1.00000 ( 0.16110 )
θ ^ = 0.36290 ( 0.09022 )
MOLL0.438280.048240.087120.96832 c ^ = 113.207 ( 7.31161 × 10 6 )
λ ^ = 10.258 ( 1.76675 × 10 6 )
β ^ = 1.60132 ( 0.57739 )
PBXLL0.846840.144360.147820.48651 c ^ = 26.6401 ( 6.83793 )
λ ^ = 6.65389 × 10 8 ( 0.00095 )
β ^ = 0.12657 ( 0.00757 )
TILL13.41332.931890.511870.00000 λ ^ = 1.00000 ( 0.44793 )
α ^ = 0.33840 ( 0.04258 )
ILL24.09765.042010.67941 0.00000 α ^ = 0.29516 ( 0.04048 )
WGLL0.304830.044630.089400.96016 c ^ = 0.00178 ( 8.97093 )
λ ^ = 69.1613 ( 347363 )
β ^ = 4.26012 ( 0.62882 )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alfaer, N.M.; Gemeay, A.M.; Aljohani, H.M.; Afify, A.Z. The Extended Log-Logistic Distribution: Inference and Actuarial Applications. Mathematics 2021, 9, 1386. https://doi.org/10.3390/math9121386

AMA Style

Alfaer NM, Gemeay AM, Aljohani HM, Afify AZ. The Extended Log-Logistic Distribution: Inference and Actuarial Applications. Mathematics. 2021; 9(12):1386. https://doi.org/10.3390/math9121386

Chicago/Turabian Style

Alfaer, Nada M., Ahmed M. Gemeay, Hassan M. Aljohani, and Ahmed Z. Afify. 2021. "The Extended Log-Logistic Distribution: Inference and Actuarial Applications" Mathematics 9, no. 12: 1386. https://doi.org/10.3390/math9121386

APA Style

Alfaer, N. M., Gemeay, A. M., Aljohani, H. M., & Afify, A. Z. (2021). The Extended Log-Logistic Distribution: Inference and Actuarial Applications. Mathematics, 9(12), 1386. https://doi.org/10.3390/math9121386

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