Next Article in Journal
Adjustments in the Range of Angular Motion during Walking after Amputation of the Toes: A Case Report
Next Article in Special Issue
Convolutions for Bernoulli and Euler–Genocchi Polynomials of Order (r,m) and Their Probabilistic Interpretation
Previous Article in Journal
A Novel (2, 3)-Threshold Reversible Secret Image Sharing Scheme Based on Optimized Crystal-Lattice Matrix
Previous Article in Special Issue
Inventory Models for Non-Instantaneous Deteriorating Items with Expiration Dates and Imperfect Quality under Hybrid Payment Policy in the Three-Level Supply Chain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Alternate Generalized Odd Generalized Exponential Family with Applications to Premium Data

by
Sadaf Khan
1,*,
Oluwafemi Samson Balogun
2,*,
Muhammad Hussain Tahir
1,
Waleed Almutiry
3 and
Amani Abdullah Alahmadi
4
1
Department of Statistics, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
2
School of Computing, University of Eastern Finland, 70211 Kuopio, Northern Europe, Finland
3
Department of Mathematics, College of Science and Arts in Ar Rass, Qassim University, Buryadah 52571, Saudi Arabia
4
College of Science and Humanities, Shaqra University, Shaqra 15572, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(11), 2064; https://doi.org/10.3390/sym13112064
Submission received: 19 September 2021 / Revised: 14 October 2021 / Accepted: 14 October 2021 / Published: 1 November 2021
(This article belongs to the Special Issue Symmetric Distributions, Moments and Applications)

Abstract

:
In this article, we use Lehmann alternative-II to extend the odd generalized exponential family. The uniqueness of this family lies in the fact that this transformation has resulted in a multitude of inverted distribution families with important applications in actuarial field. We can characterize the density of the new family as a linear combination of generalised exponential distributions, which is useful for studying some of the family’s properties. Among the structural characteristics of this family that are being identified are explicit expressions for numerous types of moments, the quantile function, stress-strength reliability, generating function, Rényi entropy, stochastic ordering, and order statistics. The maximum likelihood methodology is often used to compute the new family’s parameters. To confirm that our results are converging with reduced mean square error and biases, we perform a simulation analysis of one of the special model, namely OGE2-Fréchet. Furthermore, its application using two actuarial data sets is achieved, favoring its superiority over other competitive models, especially in risk theory.

1. Introduction

In recent years, there has been a dramatic growth in the number of generalisations of well-known probability distributions. Most notable generalizations are achieved by (i) inducting power parameters in well established parent distributions, (ii) extending the classical distribution by modification in their functions, (iii) introducing special functions such as W[K(x)] as generators and (iv) by compounding of distributions. This heaped surge of generalized families is due to the flexibility in modelling phenomenons related to the changing scenarios of contemporary scientific field including demography, actuarial, survival, biological, ecological, communication theory, epidemiology and environmental sciences. However, a clear understanding of the applicability of these models in most applied areas is necessary if one is to gain insights into systems that can be modeled as random processes. The model, thus obtained, acquires improved empirical results to the real data that is collected adaptively.
Although there exist many functions which act as generators to produce flexible classes of distributions, in this project, we will emphasize generalizations in which a ratio of survival function (sf) has been used in some form, commonly known as the odd ratio. In the reference [1], a proportional odd family viz. a viz. the Marshall Olkin-G (MO-G) was generalized by sf K ¯ ( x ) = 1 K ( x ) , where K ( x ) is the distribution function (cdf) of parent distribution, with the induction of a tilt parameter. Gleaton and Lynch, in the reference [2], used the odd function as generator when they defined a log odd family (OLL-G). In the reference [3], defined the odd Weibull family, as an asymptotically equivalent log-logistic model for larger values of θ , the scale parameter. The reference [4] used the Transformed -Transformer (TX) family, due to the reference [5], to define odd Weibul-G families of distribution. Since then, a myriad of distributions has been generalized using odd function. Some of the important families include [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31], among others.
Focussing on the origins and motivations of our proposed scheme, the authors in [32] proposed the odd generalized exponential family (which we refer to OGE-G) as a better alternative to generalized exponential (GE) family using Lehmann Alternative-I (LA-I). The cdf of the two parameter OGE-G family is mentioned below:
F O G E - G ( x ; α , λ , ψ ) = 1 e λ K ( x ; ψ ) K ¯ ( x ; ψ ) α , x > 0 , α > 0 , λ > 0 .
In the reference [33], an odd family of GE was proposed so-called generalized odd generalized exponential family (which we refer to OGE1-G). The cdf of OGE1-G family is presented as:
F O G E 1 - G ( x ; α , β , ψ ) = 1 e K α ( x ; ψ ) 1 K α ( x ; ψ ) β , x > 0 , α > 0 , β > 0 .
Because of its capacity to simulate variable hazard rate function (hrf) forms of all traditional types in lifetime data analysis, we believe OGE-G offers a sensible combination of simplicity and flexibility. However, the relevance of OGE1-G to lifespan modelling in domains such as reliability, actuarial sciences, informatics, telecommunications, and computational social sciences (just to highlight a few) is still debatable. According to the reference [34], the Lehmann Alternative-II (LA2) approach has received less attention. This motivated us to use LA2 approach to develop the exponentiated odd generalized exponential (OGE2-G), in the same vein as OGE-G and OGE1-G. Adhering to the framework defined in the reference [5], if T follows GE random variable (rv), then the cdf of OGE2-G family is mentioned below:
F O G E 2 - G ( x ; α , β , ψ ) = 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α β , x > 0 , α , β > 0
where α and β are shape parameter and ψ is the vector of baseline parameter.
Consider the following points to emphasise the model’s distinctiveness; (i) In the literature, the proposed model in its current form has not been studied to the best of our knowledge, (ii) From an analytical standpoint, the OGE2-G family has a significantly better configuration and practicality than OGE-G and OGE1-G for inverted models with minimal chance to counter non-identifiability issues, (iii) The OGE2-G has several curious connections to other families. When α approaches 0, F ( x ; α , β , ψ ) tends to GE with λ = 1 , when α = 1 F ( x ; α , β , ψ ) tends to OGE-G, if α 0 and β = 1 then F ( x ; α , β , ψ ) tends to odd exponential (OE) (iv) This new dimension allowed us to explore models which are naturally constituted by LA2. The generalizations, thus attained, produced skewed distributions with much heavier tails enabling its practicality in risk evaluation theory with far better results, (v) The successful application of OGE2-G family motivates future research, as it outperforms nine well-established existing models, (vi) We present a physical explanation for X when α and β are integers. Consider there be a parallel system consisting of β identically independent components. Suppose that the lifetime of a rv Y with a specific K ( x ; ψ ) with α components in a series system such that the risk of failing at time x is represented by the odd function as 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α . Consider that the randomness of this risk is represented by the rv X, then we can assume the following relation holds
P r ( Y x ) = P r ( X 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α ) = F ( x ; α , β , ψ ) ,
explicitly given in Equation (1). The OGE2-G family is offered and explored in this research, emphasising its diversity and scope for application to real life phenomenons. The major features of the OGE2-G family, including the pdf, hrf, qf, and ten unique models from OGE2-G family presented in Table 1, are provided in the first half. Then, certain mathematical properties of the OGE2-G such as series expansion of the exponentiated pdf, moments, parameter estimation, order statistics, Rényi entropy, stress-strength analysis and stochastic dominance results are investigated. Furthermore, Fréchet is specified as baseline model termed as OGE2-Fréchet (denoted as OGE2Fr) and the maximum likelihood (ML) technique is then used to construct statistical applications of the special model. We choose to study OGE2Fr specifically as its nested model include inverse-Rayleigh (IR) and inverse exponential (IE), favoring its suitability over sub-models as well. It is applied to fit two sets of premium data from actuarial field. Using key performance indicators, we reveal that OGE2Fr outperforms nine competing models. A portion pertinent to specific risk measures, with an emphasis on the value at risk (VaR) and the expected shortfall (ES), is presented. Eventually, the estimation of risk measures for the examined data sets is then discussed, with the proposed methodology yielding a rather satisfying result. Equation (1) can be useful in modelling real life survival data with different shapes of hrf. Table 1 lists 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α and the corresponding parameters for some special distributions which are considered to be the potential sub-models of OGE2-G family.
The following is a breakdown of how the paper is constructed. In Section 2, we acquaint the readers to the new family with basic properties and ten potential baseline models which can become members of OGE2-G family. Section 3 is comprised of the mathematical properties of the OGE2-G family. Section 4 progresses by taking Fréchet (Fr) as sub-model to propose OGE2Fr and related statistical and inferential properties. Section 5 specifies two applications of actuarial data sets with emphasis on risk evaluation (premium returns) and the proposed model’s veracity is established. Furthermore, the model is applied to compute some actuarial measures. Section 6 is the final section, with some annotations and useful insights.

2. The OGE2G Family

In this segment, basic statistical properties of the newly proposed family characterized by the cdf, in Equation (1) are presented. Functional forms of ten sub models are also defined.

2.1. Definition of pdf and hrf

The pdf in agreement with Equation (1) is given as (2).
f O G E 2 G ( x ; α , β , ψ ) = α β k ( x ; ψ ) K ¯ ( x ; ψ ) α 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α β 1
Using the results defined in Equations (1) and (2), the hrf is defined as
τ ( x ; α , β , ψ ) = f O G E 2 G ( x ; α , β , ψ ) 1 F O G E 2 G ( x ; α , β , ψ ) = α β k ( x ; ψ ) K ¯ ( x ; ψ ) α 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α β 1 × 1 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α β 1
The hazard rate is just a calculation of the change in survivor rate per unit of time. Hence, its importance in reliability and survival analysis is crucial. The hrf has some characteristic shapes which include monotonic (increasing, decreasing), non-monotonic (bathtub or upside down bathtub) or constant. Standard statistical distribution yield maximum three shapes, but OGE2-G family can yield a diverse range of shapes (including increasing-decreasing-increasing) depending upon the choice of special model. For further details on hrf, see [35].

2.2. Quantile Function and Potential Sub-Models

The OGE2-G family may be readily approximated by reversing Equation (1) as shown below: If indeed the distribution of u is uniform u ( 0 , 1 ) , therefore
x = Q K 1 1 log 1 u 1 / β 1 / α .
Equation (4) can be useful to define statistical measures such as median, skewness, and kurtosis based on quartiles, deciles, or percentiles. These measures facilitates to concisely define the skewness and kurtosis measures which are significant tool to comprehend the shape(s) of the distribution.
Theorem 1 shows how the OGE2 family is related to other distributions.
Theorem 1.
Let X O G E 2 G ( α , β ; ψ ) , then
(a) If Y = 1 K ¯ ( x ; ψ ) α , then F Y ( y ) = 1 e Y 1 Y β , 0 < y < 1 , and
(b) If Y = 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α , then Y G E ( 1 , β ) .

3. Mathematical Properties of OGE2-G Family

To capture the family’s modelling capacity, numerous mathematical features of the OGE2-G are examined in this section. Some of the key results established in this section are then applied in Section 5.

3.1. Linear Expansion of cdf

We provide a useful expansion for (1) in terms of linear combinations of exp-G density functions using the following series expansion as
( 1 z ) η 1 = i = 0 ( 1 ) i Γ ( η ) i ! Γ ( η i ) z i ,
whereas the expansion holds for all | z | < 1 and η > 0 a non-integer value. Then, the cdf of OGE2-G class in (1) can indeed be phrased with
F ( x ) = i = 0 ( 1 ) i Γ ( β + 1 ) i ! Γ ( β + 1 i ) e i 1 K ¯ ff K ¯ ff .
Using series expansion and power series expansion in Equation (5), will yield the following cdf
F ( x ) = F ( x : α , β , ψ ) = = 0 ξ H ( x ) ,
where H ( x ) = K ( x ; ψ ) (for 1 ) denotes the cdf of exp-G distribution with power parameter and
ξ = ( 1 ) i , j = 1 k = 0 j ( 1 ) i + j + k Γ ( β + 1 ) i ! j ! Γ ( β + 1 i ) j k α ( k i ) .
Through differentiating Equation (6) the OGE2-G family density, we may express it as a combination of exp-G densities.
f ( x : α , β , ψ ) = i , j = 1 ξ h ( x ) ,
where h ( x ) = K 1 ( x ; ψ ) k ( x ; ψ ) is the exp-G pdf with power parameter . As a result, numerous features of the proposed model may be deduced from the exp-G distribution’s attributes. Most modern computation frameworks, such as MathCad, Maple, Mathematica, and Matlab, can efficiently handle the formulas derived throughout the article, which can currently operate using the use of analytic formulations of enormous size and complexity.

3.2. Numerous Types of Moments

The fundamental formula for the pth moment of X is supplied by (7) as
μ p = i , j = 0 ξ E ( X p ) .
where E ( X p ) = 0 x p h ( x ) d x . Setting p = 1 in (8) can provide explicit expression for the mean of several parent distributions.
A another expression for μ p is taken from (8) as far as the baseline qf is concerned
μ p = i , j = 0 ξ τ ( p , 1 ) .
where τ ( p , 1 ) = 0 1 Q G ( u ) p u d u .
The central moments ( μ p ) and cumulants ( κ p ) of X can follow from Equation (8) as μ p = k = 0 p p k ( 1 ) k μ 1 p μ p k and κ s = μ s k = 1 s 1 s 1 k 1 κ k μ s k , respectively, where κ 1 = μ 1 .
The rth lower incomplete moment of X can be determined from Equation (7) as
m r ( y ) = i , j = 0 ξ 0 G ( y ) Q G ( u ) r u l 1 d u .
For most G distributions, the final integral may be calculated.

3.3. Inference Related to OGE2 Family

The strategy of maximum likelihood (MLL) approach is used to estimate the unknown parameters of the new class. Let x 1 , , x n be n observations from the OGE2-G density class (2) with parameter vector Θ = ( α , β , ψ ) . Then the likelihood function L ( α , β , δ ) on the domain Θ is defined as
L = n log ( α ) + n log ( β ) + i = 1 n log k ( x i ; ψ ) ( α + 1 ) log K ¯ ( x i ; ψ ) V ( x i ; α , ψ ) + β 1 i = 1 n log 1 e V ( x i ; α , ψ ) ,
where V ( x i ; α , ψ ) = 1 K ¯ ( x i ; ψ ) α K ¯ ( x i ; ψ ) α .
The elements of the score vector U ( Θ ) are as described in the following:
U α = n α i = 1 n log K ¯ ( x i ; ψ ) V ( α ) ( x i ; α , ψ ) + i = 1 n β 1 e V ( x i ; α , ψ ) V ( α ) ( x i ; α , ψ ) 1 e V ( x i ; α , ψ ) , U β = n β + i = 1 n log 1 e V ( x i ; α , ψ ) , U ψ = i = 1 n k ( x i ; ψ ) k ( x i ; ψ ) + ( α + 1 ) i = 1 n k ( x i ; ψ ) k ( x i ; ψ ) K ¯ ( x i ; ψ ) V ( ξ k ) ( x i ; α , ψ ) + ( β 1 ) i = 1 n e V ( x i ; α , ψ ) V ( ξ k ) ( x i ; α , ψ ) 1 e V ( x i ; α , ψ ) ,
where V ( α ) ( . ) and V ( ψ k ) ( . ) means the derivative of the function V with respect to α and ψ , respectively.
The next elements are produced by the components of the score vector J ( Θ ) .
J α α = n α 2 + i = 1 n k ( x i ; α , β , ψ k ) V ( ψ k ) K ¯ ( x i ; α , β , ψ k ) V ( α α ) + ( β 1 ) × i = 1 n e V ( x i ; α , ψ ) V ( x i ; α , α ) + e V ( x i ; α , ψ ) ( V ( x i ; α ) ) 2 ( 1 e V ( x i ; α , ψ ) ) 2 , J α β = i = 1 n e V ( x i ; α , ψ ) V α ( x i ; α , ψ ) [ 1 e V ( x i ; α , ψ ) ] 2 , J α ψ = i = 1 n k ( x i ; α , β , ψ ) k ( x i ; α , β , ψ ) K ¯ ( x i ; α , β , ψ ) V ( α ψ ) ( x i ; α , ψ ) + ( β 1 ) × i = 1 n e V ( x i ; α , ψ ) 2 V ( α ) ( x i ; α , ψ ) V ( ψ ) ( x i ; α , ψ ) e V ( x i ; α , ψ ) V ( α ψ ) ( x i ; α , ψ ) ( 1 e V ( x i ; α , ψ ) ) 2 ( β 1 ) i = 1 n e V ( x i ; α , β , ψ ) V ( α ) ( x i ; α , ψ ) V ( ψ ) ( x i ; α , ψ ) ( 1 e V ( x i ; α , ψ ) ) 2 + ( β 1 ) i = 1 n e V ( x i ; α , ψ ) 2 V ( α ) ( x i ; α , ψ ) V ( ψ ) ( x i ; α , ψ ) ( 1 e V ( x i ; α , ψ ) ) 2 , J β β = n β 2 , J β ψ = i = 1 n e V ( x i ; α , ψ ) V ( ψ ) ( x i ; α , ψ ) 1 e V ( x i ; α , ψ ) ,
J ψ ψ = i = 1 n k ( x i ; ψ ) k ( x i ; ψ ) k ( x i ; ψ ) 2 k ( x i ; ψ ) 2 + ( α + 1 ) i = 1 n k ( x i ; ψ ) k ( x i ; ψ ) { k ( x i ; ψ ) } 2 K ¯ 2 ( ( x i ; ψ ) ) V ( ψ ψ ) ( x i ; ψ ) ( β 1 ) i = 1 n e V ( x i ; α , ψ ) V ( ψ ) ( x i ; α , ψ ) 2 ( 1 e V ( x i ; α , ψ ) ) 2 ( β 1 ) i = 1 n e V ( x i ; α , ψ ) V ( ψ ψ ) ( x i ; α , ψ ) + e V ( x i ; α , ψ ) [ V ( ψ ) ( x i ; α , ψ ) ] 2 ( 1 e V ( x i ; α , ψ ) ) 2 + ( β 1 ) i = 1 n e V ( x i ; α , ψ ) 2 V ( ψ ψ ) ( x i ; α , ψ ) + e V ( x i ; α , ψ ) V ( ψ ) ( x i ; α , ψ ) 2 ( 1 e V ( x i ; α , ψ ) ) 2 ,
where V ( α α ) ( . ) is the derivative of V ( α ) ( . ) with respect to α , V ( α ψ ) ( . ) is the derivative of V α ( . ) with respect to ψ k and V ψ ψ ( . ) is the derivative of V ( ψ ) ( . ) with respect to ψ .

3.4. Entropy

The Rényi entropy due to [36], is characterized as
I R ( γ ) = 1 1 γ log 0 f γ ( x ) d x ,
Let us consider
f γ ( x ) = ( α β ) γ k γ ( x ; ψ ) K ¯ ( x ; ψ ) γ ( α 1 ) e γ 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α × 1 e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α γ ( β 1 ) .
Expanding Equation (12) as in Section 3.1, the Rényi entropy reduces to
I R ( γ ) = 1 1 γ log i = 0 l = 0 i ζ i , l * 0 k γ ( x ) K ¯ { α ( γ + i l ) γ } d x ,
where ζ i , l * = ( 1 ) i + l i ! ( α β ) γ i l j = 0 i ( 1 ) j ( j + i ) i γ ( β 1 ) j .

3.5. Order Statistics

Assume that X 1 , , X n is a random sample (RS) from the OGE2-G. Furthermore, assume that X i : n denote the ith order statistic (OS). Consequently, pdf of X i : n may be interpreted as
f i : n ( x ) = 1 β ( i , n i + 1 ) f ( x ) F ( x ) i 1 1 F ( x ) n i . = 1 β ( i , n i + 1 ) j = 0 n i ( 1 ) j n i j f ( x ) F ( x ) j + i 1 .
Inserting Equations (1) and (2) in the last equation, and expanding it as in Section 3.1, we get
f i : n ( x ) = j = 0 n i η j h m ( x ) ,
where
η j = ( 1 ) j β ( i , n i + 1 ) n i j m = 0 ξ m *
and
ξ m * = ( 1 ) m i , k = 0 l = 0 i ( 1 ) i + k + l ( k + 1 ) i Γ { β ( i + j ) } i ! k ! Γ { β ( i + j ) k } i l α ( i + l + 1 ) 1 m .

3.6. Stress-Strength Reliability

Supp X 1 OGE2 - G ( α , β 1 ; ψ ) and X 2 OGE2 - G ( α , β 2 ; ψ ) are two continuous rvs with pdfs f 1 ( x ) and f 2 ( x ) and cdfs F 1 ( x ) and F 2 ( x ) , therefore the reliability R is supplied via
R = P ( X 1 > X 2 ) = 0 f 1 ( x ) F 2 ( x ) d x .
Theorem 2.
Assume that X 1 and X 2 are two independent rvs established previously with constant parameters β 1 and β 2 . Eventually,
R = β 1 i = 0 β 1 1 i ( 1 ) i Γ ( i + 1 ) β 1 i = 0 β 1 1 i β 2 j ( 1 ) i + j Γ ( i + j + 1 ) .
Proof. 
Using Equations (1) and (2) in Equation (15), we have
0 f 1 ( x ) F 2 ( x ) d x = 0 α β 1 k ( x ; ψ ) K ¯ ( x ; ψ ) ( α + 1 ) e 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α × i = 0 β 1 1 i ( 1 ) i 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α i × 1 j = 0 β 2 j ( 1 ) j 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α j d x .
Equation (16) follows immediately after solving the integral with any mathematical software. □

3.7. Stochastic Ordering

Stochastic ordering has indeed been acknowledged as an essential tool for assessing comparative behavior in reliability theory and other disciplines. Assume X and Y be two rvs via cdfs, sfs and pdfs F 1 ( x ) and F 2 ( x ) , F ¯ 1 ( x ) = 1 F 1 ( x ) and F ¯ 2 ( x ) = 1 F 2 ( x ) , and f 1 ( x ) and f 2 ( x ) , respectively. In the specific planning, the rv X 1 is considered to be lower than X 2 :
  • Stochastic order (symbolized via X 1 s t X 2 ) if F ¯ 1 ( x ) F ¯ 2 ( x ) for all x;
  • LL ratio order (symbolized via X 1 l r X 1 ) if f 1 ( x ) / f 2 ( x ) is decreasing in x 0 ;
  • Hazard rate order (symbolized via X 1 h r X 2 ) if F ¯ 1 ( x ) / F ¯ 2 ( x ) is decreasing in x 0 ;
  • Reversed hazard rate order (symbolized via X 1 r h r X 2 ) if F 1 ( x ) / F 2 ( x ) is decreasing in x 0 .
All these four stochastic orders studied in (1)–(4) are connected to one another as a result of [37] and the accompanying ramifications apply:
( X 1 r h r X 2 ) ( X 1 l r X 2 ) ( X 1 h r X 2 ) ( X 1 s t X 2 ) .
when sufficient conditions are met, the OGE2-G distributions are ordered with regard to the strongest LL ratio ordering, as shown by the next theorem.
Theorem 3.
Assume X 1 O G E 2 ( α 1 , β ; ψ ) and X 2 O G E 2 ( α 2 , β ; ψ ) . If α 1 < α 2 , then X 1 l r X 2 .
Proof. 
First, we have the ratio
f 1 ( x ) f 2 ( x ) = α 1 β k ( x ) K ¯ α 1 1 e 1 K ¯ α 1 K ¯ α 1 1 e 1 K ¯ α 1 K ¯ α 1 β 1 α 2 β k ( x ) K ¯ α 2 1 e 1 K ¯ α 2 K ¯ α 2 1 e 1 K ¯ α 2 K ¯ α 2 β 1 .
After simplification, we obtain
f 1 ( x ) f 2 ( x ) = α 1 K ¯ ( α 2 α 1 ) e 1 K ¯ α 1 K ¯ α 1 + 1 K ¯ α 2 K ¯ α 2 1 e 1 K ¯ α 1 K ¯ α 1 β 1 α 2 1 e 1 K ¯ α 1 K ¯ α 1 β 1 .
Next,
log f 1 ( x ) f 2 ( x ) = log α 2 α 1 + α 2 α 1 log ( K ¯ ) 1 K ¯ α 1 K ¯ α 1 + 1 K ¯ α 2 K ¯ α 2 + β 1 log 1 e 1 K ¯ α 1 K ¯ α 1 β 1 log 1 e 1 K ¯ α 2 K ¯ α 2 .
If a 1 < a 2 , we obtain
d d x log f 1 ( x ) f 2 ( x ) = α 2 α 1 k ( x ) K ¯ ( α 1 + α 2 ) + β 1 k ( x ) K ¯ 1 × α 1 k ¯ α 1 e 1 K ¯ α 1 K ¯ α 1 1 α 2 K ¯ α 2 e 1 K ¯ α 2 K ¯ α 2 1 < 0 .
Thus, f 1 ( x ) / f 2 ( x ) is decreasing in x and hence X 1 l r X 2 . □

4. OGE2-Fréchet Distribution

In this section, we study the first special model defined in Section 2.2, the OGE2-Fréchet (OGE2Fr), in view of its practical application.
The OGE2Fr model can be defined from (1) by taking K ( x ; ψ ) = e ( a / x ) b and k ( x ; ψ ) = b a b x ( b + 1 ) e ( a / x ) b , as cdf and pdf of the baseline Fréchet distribution with a , b > 0 , respectively. The cdf and pdf of OGE2Fr distribution are, respectively, given by
F ( x ; α , β , a , b ) = 1 e 1 1 e a x b α β , x > 0 α , β , a , b > 0 ,
and
f ( x ; α , β , a , b ) = a b α b β x b 1 e 1 a x b 1 e a x b α 1 e a x b α 1 × 1 e 1 1 e a x b α β 1 ,
where α , β and b are shape parameters while a is scale parameter.
The hrf and qf of the OGE2Fr distribution are obtained as
h ( x ) = a b α b β x b 1 e 1 a x b 1 e a x b α 1 e a x b α 1 × 1 e 1 1 e a x b α 1 β 1 1 e 1 1 e a x b α β
Q ( u ) = a log 1 1 log 1 u 1 / β 1 / α 1 / b .
Figure 1 depicts a visualisation of the pdf and hrf functions, exhibiting the range of shapes that all these functions can take at random input parametric values. The OGE2Fr distribution’s pdf can be gradually decreasing, unimodal, and right-skewed, with different curves, tail, and asymmetric aspects, as shown in Figure 1. The hrf, on the other hand, offers an extensive range of increasing, decreasing, unimodal, and increasing-decreasing-increasing (IDI) forms. Given a wide variety of hrf shapes being offered, the OGE2Fr distribution can in fact be a useful tool to model unpredictable time-to-event phenomena.

4.1. Linear Representation and Related Properties

The cdf of the OGE2Fr distribution is quite straightforward and is achieved by using the result defined in Equation (6) as
F ( x ) = F ( x : α , β , a , b ) = = 0 ξ e ( a / x ) b 1 ,
where is the power parameter and noting that = 0 ξ is unity.
For simplicity, we can rewrite the above result as
F ( x ) = F ( x : α , β , a , b ) = = 0 ξ e ( a / x ) b .
By differentiating the last term, we can express the density of OGE2Fr model as follows
f ( x : α , β , a , b ) = i , j = 1 ξ Π ( x ; , a , b ) ,
where Π ( x ; , a , b ) = b a b x b 1 e ( a x ) b represents the Fréchet density function with power parameter . Equation (22) enforces the fact that OGE2Fr density is a linear combination of Fréchet densities. Thus, we can derive various mathematical properties using Fréchet distribution.
Moments are the heart and soul of any statistical analysis. Moments can be used to evaluate the most essential characteristics such as mean, variance, skewness, and kurtosis of a distribution. We now directly present the mathematical expressions for the moments of OGE2Fr model as follows. Let Y be a random variable with density Π ( x ; , a , b ) . Then, core properties of X can follow from those of Y . First, the pth ordinary moment of X can be written as
μ p = i , j = 1 ξ a p 1 p / b Γ ( 1 p / b ) .
Second, the cumulants ( κ p ) of X can be determined recursively from (23) as κ s = μ s k = 1 s 1 s 1 k 1 κ k μ s k , respectively, where κ 1 = μ 1 . The skewness γ 1 = κ 3 / κ 2 3 / 2 and kurtosis γ 2 = κ 4 / κ 2 2 of X can be calculated from the third and fourth standardized cumulants. Plots of mean, variance, skewness and kurtosis of the OGE2Fr distribution are displayed in Figure 2. These plots signifies the significant role of the parameters α and β in modeling the behaviors of X.
Third, the pth incomplete moment of X, denoted by m p ( y ) = E ( X p X y ) = 0 y x p f O G E 2 F r ( x ) d x , is easily found changing variables from the lower incomplete gamma function γ ( v , u x ) = 0 x v 1 e u x dx when calculating the corresponding moment of Y . Then, we obtain
m p ( y ) = i , j = 1 ξ a p 1 + p / b Γ 1 p / b , ( a / x ) b .
Fourth, the first incomplete moment m 1 ( z ) is used to construct the Bonferroni and Lorenz curves as discussed in Section 3.2. Figure 3 provides the income inequality curves (Bonferroni & Lorenz) of the proposed distribution which can easily be derived from (24), respectively, where q = Q ( π ) is the qf of X derived from Equation (20).

4.2. Parameter Estimation

Let x 1 , , x n be a sample of size n from the OGE2Fr distribution given in Equation (18). The log-likelihood function = ( Θ ) for the vector of parameters Θ = ( α , β , a , b ) is -4.6cm0cm
= n b log ( a ) + n log ( α β b ) ( b + 1 ) i = 1 n log ( x i ) i = 1 n a x i b + i = 1 n 1 1 e ( a / x i ) b α α + 1 i = 1 n log 1 e ( a x i ) b + ( β 1 ) i = 1 n log 1 e 1 1 e ( a x i ) b α .
The components of score vectors U( Θ ) are
U α = n α i = 1 n log 1 e a x i b + i = 1 n 1 e a x i b α log 1 e a x i b + ( β 1 ) i = 1 n { 1 e a x i b } α log { 1 e a x i b } 1 e 1 + { 1 e a x i b } α , U β = n β + i = 1 n log 1 e 1 1 e a x i b α , U a = n b a a b b a i = 1 n x i b a b b ( α + 1 ) a i = 1 n x i b e a x i b 1 a b b α a i = 1 n x i b 1 e a x i b α 1 e a x i b + α ( β 1 ) a b b a i = 1 n x i b 1 e a x i b α 1 e a x i b [ 1 e 1 + 1 e a x i b α ] ,
U b = n b + n log ( a ) i = 1 n log x i a b i = 1 n x i b log a x i ( α + 1 ) a b i = 1 n x i b log a x i e a x i b 1 + α a b i = 1 n x i b log a x i 1 e a x i b α e a x i b 1 + α a b ( β 1 ) i = 1 n x i b log a x i 1 e a x i b α e a x i b 1 [ 1 e 1 + 1 e a x i b α ] .

4.3. Order Statistics

For a random sample of X 1 , , X n taken from the OGE2Fr distribution with X i : n as the ith order statistic. For i = 1 , 2 , 3 , , n , the pdf corresponding to X i : n can be expressed as
f i : n ( x ) = 1 β ( i , n i + 1 ) f ( x ) F ( x ) i 1 1 F ( x ) n i . = 1 β ( i , n i + 1 ) j = 0 n i ( 1 ) j n i j f ( x ) F ( x ) j + i 1 ,
where f ( x ) and F ( x ) are the pdf and cdf of OGE2Fr distribution, respectively. Inserting Equations (17) and (18), and using the result defined in Section 3, we have
f i : n ( x ; m a , b ) = j = 0 n i η j m b a b x b 1 e m ( a / x ) b ,
where
η j = ( 1 ) j β ( i , n i + 1 ) n i j m = 0 ξ m *
and
ξ m * = ( 1 ) m i , k = 0 l = 0 i ( 1 ) i + k + l ( k + 1 ) i Γ { β ( i + j ) } i ! k ! Γ { β ( i + j ) k } i l α ( i + l + 1 ) 1 m .
f i : n ( x ; m a , b ) is the probability density function of the OGE2Fr distribution with parameters m a and b.

4.4. Stochastic Ordering

In several areas of probability and statistics, stochastic ordering and disparities are being adhered to at an accelerating rate. For example, in analyzing the contrast of investment returns to random cash flows; two manufacturers may use distinct technologies to make gadgets with the same function, resulting in non-identical life distributions or comparing the strength of dependent structures. Here, we use the term stochastic ordering to refer to any ordering relation on a space of probability measures in a wide sense. Let X and Y be two rvs from OGE2Fr distributions, with assumptions previously mentioned in Section 3. Given that a 1 < a 2 , and for X 1 l r X 2 , f 1 ( x ) / f 2 ( x ) shall be decreasing in x if and only if the following result holds Let X and Y be two rvs from OGE2Fr distributions, with assumptions previously mentioned in Section 3. Given that a 1 < a 2 , and for X 1 l r X 2 , f 1 ( x ) / f 2 ( x ) shall be decreasing in x if the following result holds
d d x log f 1 ( x ) f 2 ( x ) = α 2 α 1 b a b x ( b + 1 ) e ( a / x ) b 1 e ( a / x ) b ( α 1 + α 2 ) + β 1 b a b x ( b + 1 ) e ( a / x ) b 1 e ( a / x ) b 1 × α 1 1 e ( a / x ) b α 1 e 1 e ( a / x ) b a 1 1 1 + 1 β b a b x ( b + 1 ) e ( a / x ) b { 1 e ( a / x ) b } 1 × α 1 1 e ( a / x ) b α 2 e 1 e ( a / x ) b a 2 1 1 < 0 .

4.5. Simulation Study

By using the result defined in Equation (25), we evaluate the sensitivity of the method of estimations using the MLEs of OGE2Fr distribution parameters by Monte Carlo simulation technique. The simulation study is conducted for sample sizes n = 50 , 100 , 200 , 300 , 500 , 600 and parameter combinations, denoted by ω ( . ) , are:
  • ω 1 : α = 0.3 , β = 1.7 , a = 0.8 and b = 0.5 ;
  • ω 2 : α = 0.5 , β = 1.5 , a = 0.8 and b = 0.7 ;
  • ω 3 : α = 1.5 , β = 2.4 , a = 1.5 and b = 0.5 .
We use Equation (20) to generate the random observations. For each ω ( . ) , the empirical bias and MSE values are the average of the values from N = 1000 simulated samples for given sample size n. The formula to evaluate the mean squared error (MSEs) and the average bias (Bias) of each parameter, is given below
M S E ( Θ ^ ) = i = 1 N ( Θ i ^ Θ ) 2 N and B i a s ( Θ ^ ) = i = 1 N Θ i ^ N Θ .
We report the results of the AE, Bias and MSE for the parameters α , β , a and b in Table 2. The MSE of the estimators increases when the assumed model deviates from the genuine model, as anticipated. When the sample size grows larger and the symmetry degrades, the MSE shrinks. Generally speaking, the MSE decreases when the kurtosis grows. Similarly, when the asymmetry rises, the bias grows, and vice versa. As the kurtosis grows, the bias becomes smaller. In conclusion, it is apparent that the MSEs and Biases decrease when the sample size n increases. Thus, we can say that the MLEs perform satisfactorily well in estimating the parameters of the OGE2Fr distribution.

5. Application of OGE2Fr to Premium Data

Most skewed distributions are suitable to measure risk measures associated with actuarial data. The risks involve credit, portfolio, capital, premiums losses, and stocks prices among others. We focus our attention on the stakes based on premiums. Premiums are the payments for insurance that the customer pay to the company to which they are insured. In this section, we apply the OGE2Fr lifetime model for the statistical analysis of two real life data sets both of which include premium losses. Our aim is to compare the fits of the OGE2Fr model with other well-known generalizations of the Fréchet (Fr) models given in Table 3.
The first premium data set, designated as PD1, is derived from complaints upheld against vehicle insurance firms as a proportion of their overall business over a two-year period. The study was conducted by DFR (Darla Fry Ross) insurance and investment company (2009–2016), registered in New York state. The most common complaints are over delays in the settlement of no-fault claims and non-renewal of insurance. Top of the list are insurers with the fewest upheld complaints per million USD of premiums. The companies with the greatest complaint ratios are at the bottom of the list. The data understudy is from the year 2016. The second premium data, denoted by PD2, signifies the net premiums written (in billions of USD) to insurers which, under Article 41 of the New York Insurance Law, are required to meet minimum financial security requirements. Table 4: Descriptives statistics of PD1 and PD2.
The OGE2Fr model is validated through the discriminatory criterions (DCs) we considered for each data set. It includes the negative log-likelihood ( ^ ) of the model taken at the corresponding MLEs, the Akaike Information Criterion (AICs), Bayesian Information Criterion (BICs), Anderson-Darling (AD), Cramér–von Mises (CvM), and Kolmogrov-Smirnov (KS) as well as the p-value (P-KS) of the related KS test. We use the method of maximum likelihood estimation to estimate the unknown parameters as presented in Section 4.2. For each criterion (except p-value (KS)) with highest value), the smallest values is gained by the OGE2Fr model, indicating the best fit among its competitive models.
Some descriptive statistics related to these data are given in Table 4. The skewness and kurtosis are indicative of exponentially tailed data (reversed-J shape). The TTT plots for the both data sets are given in Figure 4. In particular, the TTT plots show largely decreasing hrf, permitting to fit OGE2Fr model on these data sets. The estimated hrf in Figure 5 matches Figure 4. In Table 5, we present the estimates (MLEs) along with their respective standard errors(SEs) while the DCs are listed in Table 6 for PD1 & PD2, respectively. For a more visual view, the estimated pdf, cdf, sf and Q-Q plots of the OGE2Fr model for two data sets are displayed in Figure 6 and Figure 7. Furthermore, the PP-plots of OGE2Fr and its three other competitive 4-parameter models for PD1 and PD2 are displayed in Figure 8 and Figure 9. The log-likelihood function profiles for PD1 and PD2, respectively, are provided in Figure 10 and Figure 11 to highlight the universality of the MLEs of Θ vector. The graphical visualizations are indicative of nice fits for the OGE2Fr model.

5.1. Actuarial Measures

One of the most important duties of actuarial sciences organizations is to assess market risk in a portfolio of instruments, which originates from changes in underlying factors such as equities prices, interest rates, or currency rates.One of the most important duties of actuarial sciences organizations is to assess market risk in a portfolio of instruments, which originates from changes in underlying factors such as equities prices, interest rates, or currency rates. We compute several important risk measures for the suggested distribution in this section, such as Value at Risk (VaR) and Expected Shortfall (ES), which are important in portfolio optimization under uncertainty.

5.1.1. Value at Risk

The quantile premium principal of the distribution of aggregate losses, commonly known as Value at risk (VaR), is the most widely used measure to evaluate exposure to risk in finance. VaR of a rv is the pth quantile of its cdf. If X OGE2Fr denotes a random variable with cdf (17), then its VaR is
V a R p = a log 1 1 log 1 p 1 / β 1 / α 1 / b .

5.1.2. Expected Shortfall

Artzner et al. [47,48] recommended the use of conditional VaR instead of VaR, famously called Expected Shortfall (ES). The ES is a metric that quantifies the average loss in situations where the VaR level is exceeded. It is defined by the following expression
E S p = 1 p 0 p V a R x d x , where 0 < p < 1 .
The ES of OGE2Fr is given by
E S p = 1 p 0 p a log 1 1 log 1 p 1 / β 1 / α 1 / b d x .
Figure 12 illustrates VaR and ES for some random parameter combinations of OGE2Fr.

5.1.3. Numerical Calculation of VaR and ES

The results of OGE2Fr presented in Section 4 allowed us to further explore its application to these risk measures. From Table 5, we take the values of MLEs of PD1 and PD2, respectively, to measure the volatility associated with these measures. Higher values of these risk measures signify heavier tails while lower values indicate a much lighter tail behavior of the model. It is worth mentioning that the OGE2Fr model produced substantially more significant results than its counterparts, indicating that the model has a heavier tail. In Table 7, we show the numerical results of VaRs and ESs of PD1 and PD2, respectively, of the proposed model. For the convenience of the reader, Figure 13 show the results graphically.

6. Discussion

The OGE2-G class of distribution is proposed and studied with some mathematical properties such as ordinary and incomplete moments, mean deviations and generating functions. The maximum likelihood approach is used to estimate the model parameters. Then, we focussed our attention to one of the special member of the family defined with the Fréchet distribution, called the OGE2Fr distribution. We established the optimized maximum likelihood methodology in particular, with the goal of effectively estimating model parameters and validated their convergence by a simulation study, ensuring that the projections have asymptotic properties. To demonstrate the potentiality of the proposed model, two applications to real data sets are provided. The creation of various regression models, Bayesian parameter estimates, and studies of new data sets will all be part of a future effort. We feel that the OGE2-G family can be useful for professionals in statistical analyses beyond the scope of this research because of its several other features.

Author Contributions

Conceptualization: S.K. and M.H.T.; Methodology: W.A.; Software: O.S.B.; Validation: M.H.T., W.A. and A.A.A.; Formal analysis: S.K.; Investigation: S.K.; Resources: M.H.T.; Data curation: O.S.B. and W.A.; Writing—original draft preparation:, S.K.; Writing—review and editing: S.K.; Visualization: S.K.; Supervision: A.A.A.; Project administration: O.S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This manuscript is supported by Digiteknologian TKI-ymparisto project A74338(ERDF, 357 Regional Council of Pohjois-Savo.

Data Availability Statement

https://data.world/datasets/insurance/ (accesed on 28 June 2021).

Acknowledgments

The authors would like to thank the reviewers for their thoughtful remarks and recommendations, which considerably enhanced the paper’s presentation.

Conflicts of Interest

The authors state that they have no conflicting interests to declare in this work.

Abbreviations

The following abbreviations are used in this manuscript:
GEgeneralized exponential
LA-ILehmann alternatives I
LA-IILehmann alternatives II
OGE-GOdd generalized exponential
OGE1-Ggeneralized odd generalized exponentia distributionl
OGE2-Ggeneralized Odd generalized exponential distribution
sfsurvival function
hrfhazard rate function
pdfprobability density function
cdfcumulative distribution function
rvrandom variable
OGE2Frgeneralized Odd generalized exponential distribution
FrFréchet

References

  1. Marshall, A.W.; Olkin, I. A new method for adding parameters to a family of distributions with application to the exponential and Weibull families. Biometrika 1997, 10, 641–652. [Google Scholar] [CrossRef]
  2. Gleaton, J.U.; Lynch, J.D. Properties of generalized log-logistic families of lifetime distributions. J. Probab. Stat. Sci. 2006, 4, 51–64. [Google Scholar]
  3. Cooray, K. Generalization of the Weibull distribution: The odd Weibull family. Stat. Model. Int. J. 2006, 6, 265–277. [Google Scholar] [CrossRef]
  4. Bourguignon, M.; Silva, R.B.; Cordeiro, G.M. The Weibull–G family of probability distributions. J. Data Sci. 2014, 12, 53–68. [Google Scholar] [CrossRef]
  5. Alzatreh, A.; Famoye, F.; Lee, C. A new method for generating families of continuous distributions. Metron 2013, 71, 63–79. [Google Scholar] [CrossRef] [Green Version]
  6. Ahsan ul Haq, M.; Elgarhy, M. The Odd Fréchet-G family of Probability distributions. J. Stat. Appl. Probab. 2017, 7, 185–201. [Google Scholar] [CrossRef]
  7. Alizadeh, M.; Cordeiro, G.M.; Nascimento, A.D.C. Odd-Burr generalized family of distributions with some applications. J. Stat. Comput. Simul. 2016, 87, 367–389. [Google Scholar] [CrossRef]
  8. Afify, A.Z.; Altun, E.; Alizadeh, M.; Ozel, G.; Hamedani, G.G. The Odd Exponentiated Half-Logistic-G Family: Properties, Characterizations and Applications. Chil. J. Stat. 2017, 8, 65–91. [Google Scholar]
  9. Bakouch, H.; Chesneau, C.; Khan, M.N. The Extended Odd Family of Probability Distributions with Practice to a Submodel. Filomat 2018, 33, 3855–3867. [Google Scholar] [CrossRef]
  10. Cordeiro, G.M.; Yousof, H.M.; Ramires, T.G.; Ortega, E.M.M. The Burr XII System of densities: Properties, regression model and applications. J. Stat. Comput. Simul. 2017, 87, 1–25. [Google Scholar] [CrossRef]
  11. Chesneau, C.; Toufik, E.A. Modified Odd Weibull Family of Distributions: Properties and Applications. J. Indian Soc. Probab. Stat. 2019, 21, 259–286. [Google Scholar] [CrossRef] [Green Version]
  12. El-Morshedy, M.; Eliwa, M.S. The odd flexible Weibull-H family of distributions: Properties and estimation with applications to complete and upper record data. Filomat 2019, 33, 2635–2652. [Google Scholar] [CrossRef]
  13. Haghbin, H.; Ozel, G.; Alizadeh, M.; Hamedani, G.G. A new generalized odd log-logistic family of distributions. Commun. Stat. Theory Methods 2016, 46, 9897–9920. [Google Scholar] [CrossRef]
  14. Jamal, F.; Nasir, M.; Tahir, M.H.; Montazeri, N. The odd Burr-III family of distributions. J. Stat. Appl. Probab. 2017, 6, 105–122. [Google Scholar] [CrossRef]
  15. Torabi, H.; Montazari, N.H. The gamma-uniform distribution and its application. Kybernetika 2012, 48, 16–30. [Google Scholar]
  16. Elmorshedy, M.; Eliwa, M.S.; Afify, A.Z. The Odd Chen Generator of Distributions: Properties and Estimation Methods with Applications in Medicine and Engineering. J. Natl. Sci. Found. Sri Lanka 2020, 48, 113–130. [Google Scholar]
  17. Almamy, J.A.; Ibrahim, M.; Eliwa, M.S.; Al-mualim, S.; Yousof, H.M. The Two-Parameter Odd Lindley Weibull Lifetime Model with Properties and Applications. Int. J. Stat. Probab. 2018, 7, 57. [Google Scholar] [CrossRef]
  18. Nasiru, S. Extended Odd Fréchet-G Family of Distributions. J. Probab. Stat. 2018, 2018, 2931326. [Google Scholar] [CrossRef]
  19. Afify, A.Z.; Cordiero, G.M.; Mead, M.E.; Alizadeh, M.; Al-Mofleh, H.; Nofal, Z.M. The Generalized Odd Lindley-G Family: Properties and Applications. Ann. Braz. Acad. Sci. 2019, 91, e20180040. [Google Scholar] [CrossRef] [Green Version]
  20. Korkamaz, M.C. A new family of the continuous distributions: The extended Weibull-G family. Commun. Fac. Sci. Univ. Ank. Ser. Math. Stat. 2018, 68, 248–270. [Google Scholar] [CrossRef]
  21. Alexander, C.; Cordiero, G.M.; Ortega, E.M.M.; Sarabia, J.M. Generalized beta-generated distributions. Comput. Stat. Data Anal. 2012, 56, 1880–1897. [Google Scholar] [CrossRef]
  22. Korkamaz, M.C.; Genc, A.I. A new generalized gwo-sided class of distributions with an emphasis on two-sided generalized normal distribution. Commun. Stat.-Simul. Comput. 2015, 46, 1441–1460. [Google Scholar] [CrossRef]
  23. Korkamaz, M.C.; Cordiero, G.M.; Yousof, H.M.; Pescim, R.R.; Afify, A.Z.; Nadarajah, S. The Weibull Marshall-Olkin family: Regression model and application to censored data. Commun. Stat.-Theory Methods 2018, 48, 4171–4194. [Google Scholar] [CrossRef]
  24. Korkamaz, M.C.; Yousof, H.M.; Hamedani, G.G. The Exponential Lindley Odd Log-Logistic-G Family: Properties, Characterizations and Applications. J. Stat. Theory Appl. 2018, 17, 554–571. [Google Scholar] [CrossRef] [Green Version]
  25. Korkamaz, M.C.; Altun, E.; Yousof, H.M.; Afify, A.Z.; Nadarajah, S. The Burr-X Pareto Distribution: Properties, Applications and VaR Estimation. J. Risk Financ. Manag. 2017, 11, 1. [Google Scholar] [CrossRef] [Green Version]
  26. Korkamaz, M.C.; Alizadeh, M.; Yousof, H.M.; Butt, N.S. The Generalized Odd Weibull Generated Family of Distributions: Statistical Properties and Applications. Pak. J. Stat. Oper. Res. 2018, 542, 541–556. [Google Scholar] [CrossRef] [Green Version]
  27. Korkamaz, M.C.; Yousof, H.M.; Alizadeh, M.; Hamedani, G.G. The Topp-Leone Generalized Odd Log-Logistic Family of Distributions: Properties, Characterizations and Applications. Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. 2018, 68, 1506–1527. [Google Scholar] [CrossRef]
  28. Alizadeh, M.; Korkamaz, M.C.; Almanay, J.A.; Afify, A.Z. Another odd log-logistic logarithmic class of continuous distributions. J. Stat. Stat. Actuar. Sci. 2018, 2, 55–72. [Google Scholar]
  29. Silva, F.G.; Ramos, M.W.; Percontini, A.; Venâncio, R.; De Brito, E.; Cordeiro, G.M. The Odd Lindley-G Family of Distributions. Austrian J. Stat. 2017, 35, 281–308. [Google Scholar] [CrossRef] [Green Version]
  30. Korkamaz, M.C.; Yousof, H.M.; Ali, M.M. Some Theoretical and Computational Aspects of the Odd Lindley Fréchet Distribution. J. Stat. Stat. Actuar. Sci. 2017, 2, 129–140. [Google Scholar]
  31. Korkamaz, M.C.; Altun, E.; Yousof, H.M.; Hamedani, G.G. The Hjorth’s IDB generator of distributions: Properties, characterizations, regression modelling and applications. J. Stat. Theory Appl. 2020, 19, 59–74. [Google Scholar] [CrossRef] [Green Version]
  32. Tahir, M.H.; Cordeiro, G.M.; Alizadeh, M.; Mansoor, M.; Zubair, M.; Hamedani, G.G. The odd generalized exponential family of distributions with applications. J. Stat. Distrib. Appl. 2015, 2, 1. [Google Scholar] [CrossRef] [Green Version]
  33. Morad, A.; Ghosh, I.; Yusuf, H.M.; Rasekhi, M.; Hamedani, G.G. The generalized odd generalized exponential family of distributions. J. Data Sci. 2017, 15, 443–465. [Google Scholar] [CrossRef]
  34. Tahir, M.H.; Nadarajah, S. Parameter induction in continuous univariate distributions: Well-established G families. AABC 2015, 87, 539–568. [Google Scholar] [CrossRef]
  35. Aarset, M.V. How to identify bathtub hazard rate. IEEE 1987, 36, 106–108. [Google Scholar] [CrossRef]
  36. Rényi, A. On measures of entropy and information. In Proceedings of the 4th Berkeley Symposium on Mathematical Statistics and Probability, Berkeley, CA, USA, 20–30 July 1960; Volume 1, pp. 547–561. [Google Scholar]
  37. Shaked, M.; Shantikumar, J.G. Stochastic Orders and Their Applications; Springer: New York, NY, USA, 1994. [Google Scholar]
  38. Nadarajah, S.; Gupta, A.K. The Beta Fréchet distribution. Far East J. Theor. Stat. 2004, 14, 15–24. [Google Scholar]
  39. Souza, W.B.; Cordeiro, G.M.; Simas, A.B. Some Results for Beta Fréchet Distribution. Commun. Stat.-Theory Methods 2011, 40, 798–811. [Google Scholar] [CrossRef]
  40. Cordeiro, G.M.; De-Castro, M. A new family of generalized distributions. J. Stat. Comput. Simul. 2011, 81, 883–898. [Google Scholar] [CrossRef]
  41. Cordeiro, G.M.; Ortega, E.M.M.; Cunha, D.C.C. The exponentiated generalized class of distributions. J. Data Sci. 2013, 11, 1–27. [Google Scholar] [CrossRef]
  42. Krishna, E.; Jose, K.K.; Ristic, M.M. Application of Marshall-olkin Fréchet distribution. Commun. Stat.-Simul. Comput. 2013, 42, 76–89. [Google Scholar] [CrossRef]
  43. Nadarajah, S.; Kotz, S. The exponentiated Fréchet distribution. Interstat Electron. J. 2003, 92, 97–111. [Google Scholar]
  44. Da Silva, R.V.; De Andrade, T.A.N.; Maciel, D.B.M.; Campos, R.P.S.; Cordeiro, G.M. A New Lifetime Model: The Gamma Extended Frechet Distribution. J. Stat. Theory Appl. 2013, 12, 39–54. [Google Scholar] [CrossRef] [Green Version]
  45. Abbas, S.; Taqi, S.A.; Mustafa, F.; Murtaza, M.; Shahbaz, M.Q. Topp-Leone Inverse Weibull Distribution: Theory and Application. Eur. J. Pure Appl. Math. 2017, 10, 1005–1022. [Google Scholar]
  46. Fréchet, M. Sur la loi de probabilite de lecart maximum. APM 1927, 6, 110–116. [Google Scholar]
  47. Artzner, P.; Delbaen, F.; Eber, J.M.; Heath, D. Thinking Coherently. Risk 1997, 10, 68–71. [Google Scholar]
  48. Artzner, P. Application of coherent risk measures to capital requirements in insurance. NAAJ 1999, 3, 11–25. [Google Scholar] [CrossRef]
Figure 1. (a) Plots of density and (b) hazard rate of the OGE2Fr model for some random parameter values.
Figure 1. (a) Plots of density and (b) hazard rate of the OGE2Fr model for some random parameter values.
Symmetry 13 02064 g001
Figure 2. Plots of (a) mean, (b) variance, (c) skewness and (d) kurtosis of the OGE2Fr distribution.
Figure 2. Plots of (a) mean, (b) variance, (c) skewness and (d) kurtosis of the OGE2Fr distribution.
Symmetry 13 02064 g002
Figure 3. Plots of (a) Bonferroni curve and (b) Lorenz curve of OGE2Fr model.
Figure 3. Plots of (a) Bonferroni curve and (b) Lorenz curve of OGE2Fr model.
Symmetry 13 02064 g003
Figure 4. TTT plots of (a) PD1 and (b) PD2.
Figure 4. TTT plots of (a) PD1 and (b) PD2.
Symmetry 13 02064 g004
Figure 5. Estimated hazard rate plots of (a) PD1 and (b) PD2 of OGE2Fr.
Figure 5. Estimated hazard rate plots of (a) PD1 and (b) PD2 of OGE2Fr.
Symmetry 13 02064 g005
Figure 6. Estimated (a) density, (b) cdf, (c) sf, and (d) QQ-plot for PD1.
Figure 6. Estimated (a) density, (b) cdf, (c) sf, and (d) QQ-plot for PD1.
Symmetry 13 02064 g006aSymmetry 13 02064 g006b
Figure 7. Estimated (a) density (b) cdf (c) sf, and (d) QQ-plot for PD2.
Figure 7. Estimated (a) density (b) cdf (c) sf, and (d) QQ-plot for PD2.
Symmetry 13 02064 g007
Figure 8. PP-plots of (a) OGE2Fr alongside competitive (b) BFr, (c) KwFr and (d) EGFr (4-parameter models) for PD1.
Figure 8. PP-plots of (a) OGE2Fr alongside competitive (b) BFr, (c) KwFr and (d) EGFr (4-parameter models) for PD1.
Symmetry 13 02064 g008
Figure 9. PP-plots of (a) OGE2Fr alongside competitive (b) BFr, (c) KwFr and (d) EGFr (4-parameter models) for PD2.
Figure 9. PP-plots of (a) OGE2Fr alongside competitive (b) BFr, (c) KwFr and (d) EGFr (4-parameter models) for PD2.
Symmetry 13 02064 g009aSymmetry 13 02064 g009b
Figure 10. Profiles of the log-likelihood function for the parameters α , β , a and b, respectively, of the OGE2Fr for the PD1.
Figure 10. Profiles of the log-likelihood function for the parameters α , β , a and b, respectively, of the OGE2Fr for the PD1.
Symmetry 13 02064 g010
Figure 11. Profiles of the log-likelihood function for the parameters α , β , a and b, respectively, of the OGE2Fr for the PD2.
Figure 11. Profiles of the log-likelihood function for the parameters α , β , a and b, respectively, of the OGE2Fr for the PD2.
Symmetry 13 02064 g011
Figure 12. Plots of (a) VaR (b) ES for some parameter values.
Figure 12. Plots of (a) VaR (b) ES for some parameter values.
Symmetry 13 02064 g012
Figure 13. Estimated VaR (a,c) with ES (b,d) for PD1 & PD2.
Figure 13. Estimated VaR (a,c) with ES (b,d) for PD1 & PD2.
Symmetry 13 02064 g013
Table 1. Distributions and corresponding 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α functions.
Table 1. Distributions and corresponding 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α functions.
Distribution 1 K ¯ ( x ; ψ ) α K ¯ ( x ; ψ ) α ψ
Fréchet ( x > 0 ) 1 e ( a / x ) b α 1 ( a , b )
Generalized exponential ( x > 0 ) 1 ( 1 e λ x ) θ α 1 ( λ , θ )
Power function ( 0 < x < a ) 1 ( x / a ) b α 1 ( a , b )
Burr III ( x > 0 ) 1 ( 1 + x c ) k α 1 ( c , k )
Half-logistic ( x > 0 ) 2 ( e x + 1 ) 1 α 1 ( α )
Log-logistic ( x > 0 ) 1 + ( x / a ) b α 1 ( a , b )
Inverse Rayleigh ( x > 0 ) 1 e γ / x 2 α 1 ( γ )
Inverse Exponential ( x > 0 ) 1 e γ / x α 1 ( γ )
Normal ( < x < ) 1 ϕ ( x μ σ ) α 1 ( μ , σ )
Gumbel ( < x < ) 1 e e x μ σ α 1 ( μ , σ )
Table 2. AEs, MSEs and Biases for ω 1 , ω 2 & ω 3 .
Table 2. AEs, MSEs and Biases for ω 1 , ω 2 & ω 3 .
ω 1         ω 2 ω 3
θ AEMSEBias AEMSEBias AEMSEBias
n = 50 α 1.142.2670.982 1.970.9971.787 2.142.8920.447
β 2.771.3670.531 2.141.4590.508 3.041.6670.771
a1.701.1190.539 1.451.5170.621 2.121.4860.541
b1.230.9970.408 1.670.8990.447 1.191.0350.546
n = 100 α 0.651.8920.877 1.880.9011.205 2.022.7930.420
β 2.611.2130.509 1.991.3120.472 2.911.5220.656
a1.311.0920.511 1.111.3490.555 1.931.3970.502
b0.820.9010.388 1.400.8350.420 1.080.9230.522
n = 300 α 0.521.4070.853 1.560.8660.773 1.722.7710.407
β 2.331.8090.477 1.831.2920.443 2.671.4870.629
a1.190.9340.483 0.921.2930.507 1.861.2030.488
b0.780.8660.363 1.320.8010.417 0.870.8470.497
n = 400 α 0.390.7750.639 0.890.5190.651 1.551.6360.374
β 1.890.6860.393 1.590.8970.370 2.530.8830.409
a0.930.4730.411 0.850.9450.477 1.590.8600.359
b0.550.5190.287 0.810.5690.374 0.570.5800.338
n = 500 α 0.320.7010.497 0.540.2980.455 1.481.4750.325
β 1.670.5290.376 0.510.7740.323 2.460.6610.373
a0.810.3550.337 0.830.6700.435 1.520.6630.323
b0.490.3930.264 0.750.4110.325 0.490.4440.317
n = 600 α 0.310.4720.316 0.510.2520.424 1.451.4370.311
β 1.650.3330.288 0.490.5540.311 2.450.5160.361
a0.800.2710.298 0.800.6550.409 1.500.5770.305
b0.500.3130.224 0.710.4050.317 0.500.4400.313
Table 3. The comparative fitted models.
Table 3. The comparative fitted models.
DistributionAuthor(s) Θ
BFrNadarajah and Gupta [38]; Barreto and Souza [39]( α , β , a, b)
KwFrCordiero et al. [40]( α , β , a, b)
EGFrCordiero et al. [41]( α , β , a, b)
MOFrKrishna et al. [42]( θ , a, b)
EFrNadarajah and Kotz [43]( θ , a, b)
GaFrDa Silva et al. [44]( γ , a, b)
TLFrAbbas et al. [45]( θ , a, b)
OLiFrSilva et al. [29]( λ , a, b)
FrFréchet [46](a, b)
Table 4. The descriptive statistics related to PD1 & PD2.
Table 4. The descriptive statistics related to PD1 & PD2.
DataSample Size μ σ LowestHighestSkewnessKurtosis
PD18914.08638.381.048204.175.3134.97
PD2343629.4059,372,1617.56736,502.533.069.18
Table 5. Estimates and standard errors for PD1 & PD2.
Table 5. Estimates and standard errors for PD1 & PD2.
Distribution Θ MLE(PD1)SE(PD1)MLE(PD2)SE(PD2)
OGE2Fr α 0.0650.0042.2020.649
β 1.1770.1450.2980.143
a2.3360.09924.06419.485
b5.6850.1090.3270.124
BFr α 5.6234.69230.67113.743
β 1.6991.32434.62216.797
a0.8550.64611.3919.267
b0.7590.4530.1750.126
KwFr α 6.5814.7796.2823.189
β 1.4091.16115.2354.679
a0.6060.4965.7353.726
b0.8580.3900.1570.124
EGFr α 12.9279.1996.7585.059
β 62.07636.09953.63539.338
a21.39917.98222.62916.084
b0.1580.1190.1060.091
MOFr θ 105.26767.873105.56688.076
a0.0600.0280.3020.224
b14.8660.1770.7230.107
EFr θ 1.4091.1611.0970.554
a5.4424.611167.306163.093
b0.8580.3910.4520.128
GaFr γ 5.6423.88914.7155.016
a0.1860.4600.0080.002
b2.0600.7691.6030.395
TLFr θ 0.0750.0090.1220.098
a53.72446.32978.31754.389
b0.9530.1140.4250.089
OLiFr λ 23.26014.66635.7888.969
a6.0360.6538.5831.458
b0.2740.0490.1240.050
Fra4.1780.56410.6303.207
b1.0420.0870.4690.061
Table 6. The statistics ^ , AIC, CAIC, BIC, HQIC, AD, CvM, KS and p-value (KS) for the PD1 & PD2.
Table 6. The statistics ^ , AIC, CAIC, BIC, HQIC, AD, CvM, KS and p-value (KS) for the PD1 & PD2.
Distribution ^ AICCAICBICHQICADCvMKSp-Value (KS)
Premium data 1
OGE2Fr302.22612.44612.92622.39616.450.380.0770.0740.69
BFr306.63621.25621.72631.23625.270.860.1580.0930.39
KwFr306.70621.40621.85631.35625.420.880.1640.0920.41
EGFr306.76621.51621.99631.47625.520.840.1520.0980.34
MOFr311.09628.17628.46635.64631.181.120.1780.1100.22
EFr306.70619.40619.69626.86622.410.880.1640.0920.41
GaFr305.72617.44617.72624.90620.440.770.1460.0900.45
TLFr305.65617.29617.58624.76620.300.710.1330.0900.43
OLiFr308.403622.81623.09630.24625.800.9810.1400.1160.17
Fr306.79617.58617.72622.56619.590.950.1800.0940.39
Premium data 2
OGE2Fr286.09580.18581.56586.29582.260.180.0270.0830.96
BFr287.43582.87584.25588.97584.950.290.0420.1010.84
KwFr287.45582.89584.27588.99584.970.300.0440.1000.86
EGFr288.21584.41585.79590.52586.500.390.0530.1020.84
MOFr288.66583.32584.11587.89584.860.400.0600.1030.83
EFr288.70583.40584.20587.98584.960.470.0630.1190.68
GaFr287.62581.24582.04585.82582.810.330.0450.1010.85
TLFr287.70581.39582.19585.97582.960.350.0470.1020.84
OLiFr287.52581.54581.84586.62582.600.3100.0770.0940.894
Fr288.79581.58581.97584.63582.620.490.0650.1040.82
Table 7. Numerical measures of VaRs and ESs of PD1 & PD2.
Table 7. Numerical measures of VaRs and ESs of PD1 & PD2.
Level of SignificanceVaRs(PD1)ESs(PD1)VaRs(PD2)ESs(PD2)
0.556.7823.245698.117184.508
0.608.0233.589965.385237.819
0.659.5813.9881343.206307.461
0.7011.5884.4561887.655399.670
0.7514.2635.0162693.364523.906
0.8018.0055.7043932.723695.305
0.8523.6386.5815957.8381309.704
0.9033.2387.7689636.2391933.041
0.9554.6109.57818,069.5123659.730
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, S.; Balogun, O.S.; Tahir, M.H.; Almutiry, W.; Alahmadi, A.A. An Alternate Generalized Odd Generalized Exponential Family with Applications to Premium Data. Symmetry 2021, 13, 2064. https://doi.org/10.3390/sym13112064

AMA Style

Khan S, Balogun OS, Tahir MH, Almutiry W, Alahmadi AA. An Alternate Generalized Odd Generalized Exponential Family with Applications to Premium Data. Symmetry. 2021; 13(11):2064. https://doi.org/10.3390/sym13112064

Chicago/Turabian Style

Khan, Sadaf, Oluwafemi Samson Balogun, Muhammad Hussain Tahir, Waleed Almutiry, and Amani Abdullah Alahmadi. 2021. "An Alternate Generalized Odd Generalized Exponential Family with Applications to Premium Data" Symmetry 13, no. 11: 2064. https://doi.org/10.3390/sym13112064

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