Next Article in Journal
Credit Valuation Adjustment Compression by Genetic Optimization
Next Article in Special Issue
Loss Reserving Estimation With Correlated Run-Off Triangles in a Quantile Longitudinal Model
Previous Article in Journal
Coherent-Price Systems and Uncertainty-Neutral Valuation
Previous Article in Special Issue
Ruin Probability Functions and Severity of Ruin as a Statistical Decision Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Heavy Tailed Class of Distributions Which Includes the Pareto

by
Deepesh Bhati
1,
Enrique Calderín-Ojeda
2,* and
Mareeswaran Meenakshi
1
1
Department of Statistics, Central University of Rajasthan, Kishangarh 305817, India
2
Department of Economics, Centre for Actuarial Studies, University of Melbourne, Melbourne, VIC 3010, Australia
*
Author to whom correspondence should be addressed.
Risks 2019, 7(4), 99; https://doi.org/10.3390/risks7040099
Submission received: 23 June 2019 / Revised: 8 September 2019 / Accepted: 10 September 2019 / Published: 20 September 2019
(This article belongs to the Special Issue Loss Models: From Theory to Applications)

Abstract

:
In this paper, a new heavy-tailed distribution, the mixture Pareto-loggamma distribution, derived through an exponential transformation of the generalized Lindley distribution is introduced. The resulting model is expressed as a convex sum of the classical Pareto and a special case of the loggamma distribution. A comprehensive exploration of its statistical properties and theoretical results related to insurance are provided. Estimation is performed by using the method of log-moments and maximum likelihood. Also, as the modal value of this distribution is expressed in closed-form, composite parametric models are easily obtained by a mode matching procedure. The performance of both the mixture Pareto-loggamma distribution and composite models are tested by employing different claims datasets.

1. Introduction

In general insurance, pricing is one of the most complex processes and is no easy task but a long-drawn exercise involving the crucial step of modeling past claim data. For modeling insurance claims data, finding a suitable distribution to unearth the information from the data is a crucial work to help the practitioners to calculate fair premiums. In actuarial statistics and finance, the classical Pareto distribution has been deemed better than other models because it provides a good description of the random behaviour of large claims. Loss data mostly have several characteristics such as unimodality, right skewness, and a thick right tail. To accommodate these features in a single model, many probability models have been proposed in the literature.
Over the last decades, different approaches to derive new classes of probability distributions that could provide more flexibility when modelling large losses has been added to the literature. This includes transformation method, the composition of two or more distributions, compounding of distributions or a finite mixture of distributions among other methodologies. In particular for generalizations of the classical Pareto distribution, the reader is referred to the Stoppa distribution (see Stoppa 1990), the Pareto positive stable distribution (see Sarabia and Prieto 2009), the Pareto ArcTan distribution (see Gómez-Déniz and Calderín-Ojeda 2015) and the generalized Pareto distribution proposed in Ghitany et al. (2018). Obviously, adding a parameter to the parent model complicates the parameter estimation.
In this paper, a probabilistic family, the mixture Pareto-loggamma distribution, which belongs to the heavy-tailed class of probabilistic models is introduced. This family is derived by using of an exponential transformation of the generalized Lindley distribution. It is expressed as a convex sum of the classical Pareto and a special case of the loggamma distribution similar to the one presented in Gómez-Déniz and Calderín-Ojeda (2014). The former model is obtained as a particular case whereas the latter one is a limiting case. We further present expressions of statistical and actuarial measures such as moments, variance, cumulative distribution function, hazard rate function, VaR, TVaR and limited expected values. For many of these distributional characteristics, closed-form expressions are obtained. Estimation of the parameters of this distribution can be easily calculated by the maximum likelihood method by using the numerical search of the maximum and also by the method of log-moments.
In many instances, composite distributions give a reasonably good fit as compared to classical distributions since the former models have the advantage of accommodating both low magnitude values with a high frequency as well as large magnitude figures with low frequency. Cooray and Ananda (2005) introduced a composite lognormal-Pareto distribution with restricted mixing weights. Scollnik (2007) improved the composite lognormal-Pareto model by allowing flexible mixing weights. Bakar et al. (2015) considered several probabilistic models in place of lognormal and Pareto distribution using the approach discussed in Scollnik (2007). Recently, Calderín-Ojeda and Kwok (2016) introduced a new class of composite model using mode-matching procedure. They derived composite models using lognormal and Weibull distributions with a Stoppa model, a generalization of the Pareto distribution. As there exists a closed-form expression for the modal value of this new model, in this article, we use the mode-matching procedure to derive new composite models based on this distribution.
The structure of the paper is organized as follows. In Section 2, we present the genesis of the new distribution, and discuss its relationship with other distributions. Besides, the most relevant distributional properties are studied in the same section. Also, different estimation methods are examined. Finally, composite models based on the proposed distribution are derived. Next, in Section 3, some results related to insurance are displayed. Numerical applications are displayed in Section 4 together with the derivation of some income indices. The last section concludes the paper.

2. The MPLG Distribution and Its Properties

It can be easily seen that the following expression
f X ( x | θ , λ , x 0 ) = θ 2 x ( θ + λ ) x x 0 θ 1 + λ log x x 0 x x 0 > 0
with λ 0 and θ > 0 is a genuine probability density function (pdf). Here θ , λ are shape parameter and x 0 is scale parameter (see the Appendix A). Note that the pdf (1) can be written as a convex sum of the classical Pareto distribution and a special case of the loggamma distribution. The former model is obtained when λ = 0 and the latter one for λ . Observe that the pdf of the loggamma distribution considered in this work is
f X ( x | θ , γ , x 0 ) = θ γ x 0 Γ ( γ ) x x 0 ( θ + 1 ) log x x 0 γ 1 ,
with θ , γ > 0 and x x 0 . In our model it is assumed that γ = 2 .
The cumulative distribution function (cdf), survival function (sf) and hazard rate function of a random variable (rv) with pdf (1) are respectively given, for x x 0 , as
F X ( x ) = 1 θ + λ + θ λ log x x 0 θ + λ x x 0 θ ,
F ¯ X ( x ) = θ + λ + θ λ log x x 0 θ + λ x x 0 θ ,
h X ( x ) = θ x 1 λ θ + λ + θ λ log x x 0 .
Henceforth, a continuous rv X that follows (1) will be called as mixture Pareto-loggamma distribution and denoted as MPLG ( θ , λ , x 0 ) . The shape of the pdf given in (1) is shown in Figure 1 for different values of the parameters θ and λ and a fixed x 0 . It is observable that the larger is the value of θ , the thinner is the tail. On the contrary, the thickness of the tail decreases when λ drops.
Undoubtedly, the single parameter Pareto distribution is one of the most attractive distributions in statistics; a power-law probability distribution found in a large number of real-world situations inside and outside the field of economics. Furthermore, it is usually used as a basis for Excess of Loss quotations as it gives a pretty good description of the random behavior of large losses. In this sense, many probability distributions can be used for modelling single loss amounts. Then, if the loss is assumed to follow the pdf (1), we have that θ defines the tail behavior of the distribution.
Below, the hazard rate function has been plotted in Figure 2 for different values of the parameters θ and λ for fixed x 0 . It is observable that the hazard rate function increases when the parameter θ grows and the parameter λ declines.
The next result shows the relationship between the MPLG and the Generalized Lindley distribution, which is extensively used in lifetime analysis and reliability.
Remark 1.
The MPLG ( θ , λ , x 0 ) distribution can be also obtained by taking X = x 0 exp { Y } , where Y follows a generalized Lindley distribution with density
g Y ( y ) = θ 2 ( θ + λ ) ( 1 + λ y ) exp { θ y } , y > 0 , λ , θ > 0 .
From this result, it is straightforward to generate random variates from the MPLG ( θ , λ , x 0 ) distribution by observing that the pdf given in (5), can be rewritten as
w ( y ) = p a 1 ( y ) + ( 1 p ) a 2 ( y ) ,
where p = θ θ + λ , a 1 ( · ) is the pdf of the exponential distribution with mean 1 / θ and a 2 ( · ) is the pdf of the Erlang distribution with shape parameter 2 and rate parameter θ . Then, a random variate from the MPLG distribution, x, can be generated following a modification of the algorithm presented in Ghitany et al. (2008) as shown below.
  • Generate two random numbers u 1 and u 2 from the standard uniform distribution, U ( 0 , 1 ) .
  • Generate random variates a ˜ 1 from the exponential distribution with mean 1 / θ and a ˜ 2 from the Erlang distribution with shape parameter 2 and rate parameter θ by using u 1 .
  • If u 2 p , then set y = a ˜ 1 ; otherwise, set y = a ˜ 2 .
  • Generate x = x 0 exp { y } .
Observe that an analogous algorithm could be implemented by using the fact that (1) is a convex sum of the densities of Pareto and loggamma distributions.

2.1. Moments and Log-Moments

The rth raw moment of MPLG ( θ , λ , x 0 ) is given by
E ( X r ) = θ 2 x 0 r ( θ + λ r ) ( θ + λ ) ( r θ ) 2 θ r .
In particular, we have that
μ = E ( X ) = θ 2 x o ( θ + λ 1 ) ( θ 1 ) 2 ( θ + λ ) , θ 1 and E ( X 2 ) = θ 2 x 0 2 ( θ + λ 2 ) ( θ + λ ) ( 2 θ ) 2 , θ 2 .
The variance is given by
V ( X ) = θ 2 x 0 2 ( θ + λ 2 ) ( θ + λ ) ( 2 θ ) 2 θ 2 x 0 ( θ + λ 1 ) ( θ + λ ) ( 1 θ ) 2 2 , θ 1 , 2 .
Remark 1 can be used to obtain the rth raw log-moment of the density (1)
l s = E log X x 0 s = s ! ( θ + λ + λ s ) θ s ( θ + λ ) .
The moment estimator of θ and λ with known x 0 can be simply derived from (7).

2.2. Unimodality

The modal value of the MPLG distribution can be found by taking the first derivative of the pdf (1),
f ( x ) = θ 2 x 0 θ x θ 2 θ + λ λ ( 1 + θ ) 1 + λ log x x 0 ,
then, by setting (8) equal to zero and solving for x, the mode is x m = x 0 exp λ θ 1 λ ( θ + 1 ) . If λ θ + 1 , then the mode is at x m = x 0 .

2.3. Stochastic Ordering

Many parametric families of distributions can be ordered by some stochastic orders according to the value of their parameters. A rv X 1 is said to be stochastically smaller ( X 1 s t X 2 ) than X 2 if F X 1 ( x ) F X 2 ( x ) for all x. The formal definition of stochastic ordering is
Definition 1.
Let X 1 and X 2 be continuous random variables with densities f 1 and f 2 , respectively, such that f 2 ( x ) f 1 ( x ) is non-decreasing over the union of the supports of X 1 and X 2 . Then X 1 is said to be smaller than X 2 in the likelihood ratio order. ( X 1 l r X 2 ) (see Section 1.C of Shaked and Shanthikumar (2007)).
Similarly, a rv can be also ordered for stochastic (distribution function), hazard rate and mean excess orders when the following results hold:
(i)
Stochastic order ( X 1 s t X 2 ) if F X 1 ( x ) F X 2 ( x ) for all x.
(ii)
Hazard rate order ( X 1 h r X 2 ) if h X 1 ( x ) h X 2 ( x ) for all x.
(iii)
Mean excess order ( X 1 m e X 2 ) if e X 1 ( x ) e X 2 ( x ) for all x, where e X ( · ) is mean excess function given in expression (17).
Using Theorem 1.C.1 and Theorem 2.A.1 of Shaked and Shanthikumar (2007), the above stochastic orders hold following implications result
( X 1 l r X 2 ) ( X 1 h r X 2 ) ( X 1 m e X 2 ) ( X 1 s t X 2 ) .
This gives rise to Proposition 1.
Proposition 1.
Let X 1 and X 2 be two rv’s having MPLG distribution with parameters ( θ i , λ i , x 0 i ) i = 1 , 2 . Then following results hold
i. 
If x 01 = x 02 , θ 1 θ 2 and λ 1 λ 2 , then X 1 l r X 2 , X 1 h r X 2 and X 1 s t X 2 .
ii. 
If x 01 x 02 , θ 1 = θ 2 and λ 1 = λ 2 , then X 1 l r X 2 , X 1 h r X 2 and X 1 s t X 2 .
Proof. 
See the Appendix A. □

2.4. Integrated Tail Distribution and Equilibrium Hazard Rate

As we have already seen previously, the proposed distribution can be viewed as a mixture of two thick-tailed distributions, the Pareto and loggamma models, thus it can be used to describe events that include heavy-tailed behaviour. In this particular, we define another distribution function called Integrated tail distribution (ITD) (also known as equilibrium distribution) F e ( x ) which often appears in insurance fields (see Yang 2004). The ITD has many interesting applications, e.g., approximation of the ruin function (see Yang 2004) or characterisation of the tail of the distribution (see Su and Tang 2003) just to name a few. Hence, for MPLG ( θ , λ , x 0 ) distribution bluewith θ > 1 , ITD is obtained as
F e ( x ) = 1 E ( X ) x 0 x F ¯ ( y ) d y = ( θ 1 ) 2 ( θ + λ ) θ 2 x 0 ( θ + λ 1 ) θ λ x 0 x x x 0 θ 1 ( 1 θ ) log x x 0 ( 1 θ ) 2 ( θ + λ ) + x x x 0 θ x 0 1 θ = x x 0 θ ( θ ( 1 θ 2 λ ) + λ ) x x 0 x x 0 θ + ( 1 + θ ) θ λ x log x x 0 θ 2 x 0 ( θ + λ 1 ) .
The associated equilibrium hazard rate r e ( x ) = log 1 F e ( x ) , for θ > 1 is given by
r e ( x ) = ( θ 1 ) 2 θ + λ + θ λ log x x 0 x θ ( θ + 2 λ 1 ) λ + ( θ 1 ) θ λ log x x 0 .
Moreover, it can be easily verified that lim x r e ( x ) = 0 . Hence by using Theorem 2.1 of Su and Tang (2003) we can conclude that MPLG ( θ , λ , x 0 ) is a heavy-tailed distribution. Heavy-tailed distributions are important in non-life insurance when modeling losses related to motor third-party liability insurance, fire insurance or catastrophe insurance.

2.5. Estimation

Let x 1 , , x n be an independent and identically distributed (iid) random sample of size n drawn from a population which follows a MPLG ( θ , λ , x 0 ) distribution. Estimate of parameter x 0 can easily be obtained, as the support of the MPLG ( θ , λ , x 0 ) distribution is greater than x ^ 0 = min ( x 1 , , x n ) . Therefore for x 0 known, we estimate the parameters θ and λ by (i) method of log-moments and (ii) maximum likelihood method.

2.5.1. Method of Log-Moments

The first and second log-moments of MPLG ( θ , λ , x 0 ) distribution obtained by using (7) are given by l 1 = ( θ + 2 λ ) θ ( θ + λ ) and l 2 = 2 ( θ + 3 λ ) θ 2 ( θ + λ ) . Solving these equations for θ and λ , we have that these estimates are provided by
θ ˜ = l 1 2 l 2 + l 1 l 2 and λ ˜ = 3 l 1 l 2 3 l 1 3 3 l 1 2 l 1 2 l 2 + 2 l 2 l 1 2 l 2 l 2 ( 3 l 1 2 4 l 2 ) .

2.5.2. Maximum Likelihood Estimation

The log-likelihood function given the iid random sample x ̲ : = x 1 , , x n of size n is
n ( θ , λ | x ̲ ) = 2 n log θ n log ( θ + λ ) + i = 1 n log 1 + λ log x i x 0 θ i = 1 n log x i x 0 i = 1 n log x i ,
and the normal equations obtained by differentiating (12) with respect to θ and λ are
n ( θ , λ | x ̲ ) θ = 2 n θ n θ + λ i = 1 n log x i x 0 , n ( θ , λ | x ̲ ) λ = n θ + λ i = 1 n log x i x 0 1 + λ log x i x 0 .
The above equations cannot be solved analytically and hence maximum likelihood estimates of θ and λ can be computed numerically using in-built R -function such as nlm ( ) , maxlik ( ) or optim ( ) . Moreover, in all these functions to initialize the program we use initial values obtained from the method of log-moments. The second partial derivative of log-likelihood function with respect to parameter θ and λ are
2 n ( θ , λ | x ̲ ) θ 2 = 2 n θ 2 + n ( θ + λ ) 2 , 2 n ( θ , λ | x ̲ ) θ λ = n ( θ + λ ) 2 , 2 n ( θ , λ | x ̲ ) λ 2 = n ( θ + λ ) 2 i = 1 n log 2 x i x 0 1 + λ log x i x 0 2 .
Furthermore,
E log 2 X x 0 1 + λ log X x 0 2 = θ 2 ( θ + λ ) x 0 log 2 x x 0 1 + λ log x x 0 2 x x 0 θ 1 x 1 + λ log x x 0 d x ,
substituting t = λ log x x 0 , and using Tricomi confluent hypergeometric function U ( a , b , z ) = 1 Γ ( a ) 0 t a 1 ( 1 + t ) a + b 1 e z t d t , we get
E log 2 X x 0 1 + λ log X x 0 2 = 2 θ 2 ( θ + λ ) λ 3 U ( 3 , 3 , θ / λ ) .
Hence, the expected Fisher’s information matrix associated to the parameters θ and λ is given by I i j ( θ , λ ) with i , j = 1 , 2 where
I 11 ( θ , λ ) = E 2 n ( θ , λ | x ̲ ) θ 2 = 2 n θ 2 n ( θ + λ ) 2 , I 12 ( θ , λ ) = I 21 ( θ , λ ) = E 2 n ( θ , λ | x ̲ ) θ λ = n ( θ + λ ) 2 , I 22 ( θ , λ ) = E 2 n ( θ , λ | x ̲ ) λ 2 = n ( θ + λ ) 2 + 2 θ 2 n ( θ + λ ) λ 3 U ( 3 , 3 , θ / λ ) .
The standard errors of the estimates θ ^ and λ ^ can be obtained by inverting the aforementioned matrix and taking the square root of the diagonal entries.

2.6. Composite Models

Composite parametric models are a useful way to describe data that combined loss data of small and moderate size with a high frequency and large observations with a low frequency. They consist of two distributions, the first model is used up to an unknown single threshold value, estimated from the data, and the second distribution beyond this threshold (see Cooray and Ananda 2005). Another approach to derive composite models is utilizing a mode-matching procedure. In this technique, the two distributions are composed at the common modal value which can be estimated from the claims data. Thus, the composite model uses a truncated version of the first model up to the mode and the rest of the model is based on an appropriate truncation of the second distribution from that modal point onwards. The model after composition is similar in shape to either of the models considered but with a thicker tail. By using this methodology, it is guaranteed that the new density is continuous and smooth. These composite models give a significantly better fit as compared to the standard single models for the same empirical data. Here, we derived composite models of our MPLG distribution with lognormal, Weibull and paralogistic distributions.

2.6.1. Composite Lognormal- MPLG Model

Let us consider, for x > 0 , the two-parameter lognormal distribution with pdf and cdf respectively given by
f 1 ( x ) = 1 2 π x σ exp 1 2 log x μ σ 2 and F 1 ( x ) = Φ log x μ σ ,
with μ R , σ > 0 where Φ ( · ) is the cdf of the standard normal distribution. Now by taking expressions (1) and (2) as the pdf, f 2 ( x ) , and cdf, F 2 ( x ) of the MPLG distribution respectively with x x 0 and setting equal the modes of lognormal and MPLG distributions we obtain
σ = μ log x 0 λ θ 1 λ ( θ + 1 ) .
Note that we must impose the additional constraint μ > log x 0 + λ θ 1 λ ( θ + 1 ) . Now, the unrestricted mixing weight is given by
r = θ 2 x m ( θ + λ ) x m x 0 θ 1 + λ log x m x 0 Φ log x m μ σ × θ 2 x m ( θ + λ ) x m x 0 θ 1 + λ log x m x 0 Φ log x m μ σ + 1 2 π x m σ exp 1 2 log x m μ σ 2 θ + λ + θ λ log x m x 0 θ + λ x m x 0 θ 1 .
Then, the pdf of the four-parameter composite lognormal- MPLG distribution is given by
f ( x ) = r F 1 ( x m ) 1 1 2 π x σ exp 1 2 log x μ σ 2 0 x x m ( 1 r ) ( 1 F 2 ( x m ) ) 1 θ 2 x ( θ + λ ) x x 0 θ 1 + λ log x x 0 x m x .

2.6.2. Composite Weibull- MPLG Model

Let
f 1 ( x ) = τ x x ϕ τ exp x ϕ τ and F 1 ( x ) = 1 exp x ϕ τ
be the pdf and cdf of a two-parameter Weibull distribution with ϕ , τ > 0 and x > 0 . Let us again consider the pdf and cdf of the MPLG distribution given by (1) and (2) respectively. By equating the modes of the two distributions, we have that
ϕ = x 0 exp λ θ 1 λ ( θ 1 ) τ τ 1 1 τ .
where the restriction τ > 1 must be imposed to guarantee that the mode of the Weibull distribution is greater than 0. Now, the pdf of the four-parameter composite Weibull- MPLG distribution is provided by
f ( x ) = r F 1 ( x m ) 1 τ x x ϕ τ exp x ϕ τ 0 x x m ( 1 r ) ( 1 F 2 ( x m ) ) 1 θ 2 x ( θ + λ ) x x 0 θ 1 + λ log x x 0 x m x , .
where the unrestricted mixing weight is given by
r = θ 2 x m ( θ + λ ) x m x 0 θ 1 + λ log x m x 0 1 exp x m ϕ τ × θ 2 x m ( θ + λ ) x m x 0 θ 1 + λ log x m x 0 1 exp x m ϕ τ + τ x m x m ϕ τ exp x m ϕ τ θ + λ + θ log x m x 0 θ + λ x m x 0 θ 1 .

2.6.3. Composite Paralogistic- MPLG Model

Finally, let us assume that
f 1 ( x ) = α 2 ( x τ ) α x ( 1 + ( x τ ) α ) α + 1 and F 1 ( x ) = 1 1 ( 1 + ( x τ ) α ) α
are the pdf and cdf of the two-parameter paralogistic distribution with α, τ > 0, and x > 0 . Once again, we consider the pdf and cdf of the MPLG model given by (1) and (2). By setting equal the modal values of the two distributions, we have that
τ = x 0 exp λ θ 1 λ ( θ 1 ) 1 α 1 α 2 + 1 1 α .
where the restriction α > 1 is again established to ensure that the mode of the paralogistic distribution is larger than 0. The unrestricted mixing weight is now given by
r = θ 2 x m ( θ + λ ) x m x 0 θ 1 + λ log x m x 0 1 1 1 + ( x m τ ) α α × θ 2 x m ( θ + λ ) x m x 0 θ 1 + λ log x m x 0 1 1 1 + ( x m τ ) α α + α 2 ( x m τ ) α x m ( 1 + ( x m τ ) α ) α + 1 θ + λ + θ log x m x 0 θ + λ x m x 0 θ 1 .
The pdf of the four-parameter composite paralogistic- MPLG distribution is provided by
f ( x ) = r F 1 ( x m ) 1 α 2 ( x τ ) α x ( 1 + ( x τ ) α ) α + 1 0 x x m ( 1 r ) ( 1 F 2 ( x m ) ) 1 θ 2 x ( θ + λ ) x x 0 θ 1 + λ log x x 0 x m x .

3. Insurance Results

In the following, several theoretical results related to insurance for the MPLG ( θ , λ , x 0 ) distribution and related composite models are derived.

3.1. Mean Excess Function

The mean excess function measures the expected payment per claim on a policy with a fixed amount deductible of x, ignoring the claims with amounts less than or equal to x. It is defined as
e ( x ) = E ( X x | X > x ) = x F ¯ ( u ) d u F ¯ ( x ) ,
which can also be obtained by inverting the equilibrium hazard rate (Su and Tang 2003), hence for θ > 1 the mean excess function for the model with pdf (1) is given by
e ( x ) = x θ ( θ + 2 λ 1 ) λ + ( θ 1 ) θ λ log x x 0 ( θ 1 ) 2 θ + λ + θ λ log x x 0 .

3.2. Excess of Loss Reinsurance

Let X be a rv denoting the individual claim size taking values greater than d. Let us also assume that X follows (1), and the expected cost per claim to the reinsurance layer when the loss in excess of m subject to a maximum of l is given by
E ( min ( l , max ( 0 , X m ) ) = m m + l ( x m ) f ( x ) d x + l F ¯ ( m + l ) .
Now by replacing f ( · ) and F ¯ ( · ) by expressions (1) and (3) and solving the integral, we have that
E ( min ( l , max ( 0 , X m ) )   = θ 2 x 0 θ θ + λ ( l + m ) θ λ ( θ l + m ) log l + m x 0 ( θ 1 ) θ ( θ + λ 1 ) ( l + m ) ( θ 1 ) 2 + m ( θ + λ ) θ 2   + m 1 θ θ 2 + θ ( 2 λ 1 ) λ + ( θ 1 ) θ λ log m x 0 ( θ 1 ) 2 θ 2   + l x 0 θ ( l + m ) θ θ + λ + θ λ log l + m x 0 θ + λ .

3.3. Value-at-Risk and Tail Value-at-Risk

In this subsection, we firstly discuss the most widely used risk measure, the Value-at-Risk (VaR). It is defined as the minimum value of the distribution such that the probability of the loss larger than this value does not exceed a given probability. In statistical terms, VaR is a quantile of a random variable and the formal definition is as follows.
Definition 2.
Let X be a loss rv with a continuous cdf F X ( · ) , and δ be a probability level such that 0 < δ < 1 , the Value-at-Risk at probability level δ, denoted by VaRδ(X), is the δ-quantile of X. That is
V a R δ ( X ) = inf { x R : F X ( x ) δ } .
Hence for rv having MPLG ( θ , λ , x 0 ) distribution the Value-at-Risk is given in next Proposition.
Proposition 2.
For θ , λ > 0 , the VaRδ(X) of the MPLG ( θ , λ , x 0 ) distribution is
V a R δ ( X ) = x 0 exp θ + λ θ λ 1 θ W 1 ( 1 δ ) ( θ + λ ) exp θ + λ λ λ
where W 1 is the negative branch of Lambert-W function.
Proof. 
See the Appendix A. □
The median of MPLG ( θ , λ , x o ) will be obtained by taking δ = 0.5 in VaRδ(X). As in general, the loss distributions are typically skewed, the VaR is a non-coherent risk measure due to the lack of subadditivity (see Klugman et al. 2012), for that reason, the Tail Value-at-Risk (TVaR) of X (see Acerbi and Tasche 2002) is usually considered as a more informative and more useful risk measure. The TVaR is given by
TVaRd δ ( X ) = 1 1 δ δ 1 VaR z ( X ) d z ,
which is a coherent risk measure. If X is continuous, Pr ( X VaR δ ( X ) ) = δ , and then the TVaR is the conditional tail expectation TVaR δ ( X ) = E ( X | X > VaR δ ( X ) ) .
Proposition 3.
For θ , λ > 0 , the TVaRδ(X) for MPLG ( θ , λ , x 0 ) distribution is given by
T V a R δ ( X ) = 1 ( 1 δ ) λ x 0 θ exp ( 1 θ ) ( θ + λ ) θ λ ( 1 θ ) 2 ( θ + λ ) × ( 1 + ( 1 θ ) K ( δ ; θ , λ ) ) exp ( 1 θ ) K ( δ ; θ , λ ) θ 1
where K ( δ ; θ , λ ) = W 1 ( 1 δ ) exp θ + λ λ ( θ + λ ) λ .
Proof. 
See the Appendix A. □

3.4. Limited Loss Variable for Composite Models

Finally, in this subsection the kth moment of the limited loss variable for composite models based on the MPLG distribution is presented.
Definition 3.
The kth incomplete moment transform of a rv X with cdf F ( x ) and pdf f ( x ) , given that μ k = E ( X k ) < , is a rv with pdf given by f ( k ) ( x ) = x k f ( x ) μ k and cdf given by F ( k ) ( x ) = 0 x t k f ( t ) d t μ k , 0 < x < .
Definition 4.
The limited loss variable is defined as, X u = X , X < u , u , X u . where u is the maximum benefit paid by the insurance policy.
Proposition 4.
The moment of order k of the limited loss variable for the composite models based on the MPLG distribution is given by
E ( ( X u ) k ) = r E ( X 1 k ) F 1 ( k ) ( u ) F 1 ( x m ) + u k 1 r F 1 ( u ) F 1 ( x m ) ,
if 0 < u x m and
E ( ( X u ) k ) = r E ( X 1 k ) F 1 ( k ) ( x m ) F 1 ( x m ) + ( 1 r ) E ( X 2 k ) F 2 ( k ) ( u ) F 2 ( k ) ( x m ) 1 F 2 ( x m ) + u k ( 1 r ) 1 F 2 ( u ) 1 F 2 ( x m ) ,
if x m < u < .
Proof. 
The kth moment of the limited loss variable can be derived as,
E [ ( X u ) k ] = 0 u x k f X ( x ) d x + u k [ 1 F X ( u ) ] .
If 0 < u < x m ,
E [ ( X u ) k ] = 0 u r x k f 1 ( x ) F 1 ( x m ) d x + u k 1 r F 1 ( u ) F 1 ( x m ) ,
and we get (20). Now, if x m < u < ,
E [ ( X u ) k ] = 0 x m r x k f 1 ( x ) F 1 ( x m ) d x + x m u ( 1 r ) x k f 2 ( x ) 1 F 2 ( x m ) d x + u k 1 ( r + ( 1 r ) F 2 ( u ) F 2 ( x m ) 1 F 2 ( x m ) ) ,
and (21) is obtained. □
Observe that in expressions (20) and (21), E ( X i k ) and F i ( k ) ( · ) denotes the kth incomplete moment and the cdf of the kth incomplete moment for lognormal, Weibull or paralogistic distributions ( i = 1 ) and MPLG distribution ( i = 2 ).

4. Case Studies

In this section, we use two claims datasets to assess the performance of MPLG distribution and the MPLG composite models. Finally, some results related to income indices are given.

4.1. An Application to Automobile Claims Data

First, we examine the performance of MPLG distribution as compared to other heavy-tailed distributions available in the literature by employing of a real automobile claims dataset (see De Jong and Heller 2008). This dataset describes a one-year vehicle insurance policies taken out in 2004 or 2005. There are 67,856 policies, of which 4624 (6.8%) had at least one claim. The variable of interest is the size of the claims. The minimum claim amount is $200 and the maximum value is $55,922.13. We have fitted to this dataset the MPLG distribution by setting the value of x 0 equal to the minimum value to explain the claims amount distribution. In Table 1, parameter estimates and standard errors (S.E.) for the MPLG and other two-parameter and three-parameter heavy-tailed probabilistic models with support in R + are illustrated. For a review of these distributions the reader is referred to Hogg and Klugman (2009) or Klugman et al. (2012). For the classical Pareto, PAT the location parameter has been set equal to the minimum value, i.e., x 0 = 200 . For the loggamma distribution the location parameter can be chosen in the neighborhood of this value. For comparison purposes, three measures of model selection has been included in this Table, the negative of the log-likelihood function (NLL), Akaike’s information criterion (AIC) and Bayesian information criterion (BIC). It is observable that the MPLG provides the best fit to data in terms of these measures of model selection. The figures of the measures of model selection reveal that the Pareto distribution is preferable to the loggamma distribution for this dataset.
We have also fitted the three-parameter Burr distribution to this dataset. For this model the value of the NLL is 38,028.13. However, the algorithm used to search for the maximum of the log-likelihood surface to find the estimates stopped at the boundary of the parameter space, and thus these estimates are no longer displayed in Table 1.

4.2. An Application to Fire Insurance Claims

In our second application, the composite models derived from the MPLG distribution are compared to the composite Pareto and composite Lomax families. For that reason, we consider the well-known Danish fire insurance dataset that consists of 2492 fire insurance losses in millions of Danish kroner (DKr) from the years 1980 to 1990 (both inclusive), adjusted to reflect 1985 values. This dataset may be found in the ‘SMPrcacticals’ add-on package for R, available from the CRAN website http://cran.r-project.org/. Parameter estimation for all the models considered has been completed by the method of maximum likelihood (which is implemented using the function ‘mle’/‘mle2’ in R). In Table 2, parameter estimates and standard errors (S.E.) for all the composite distributions and the three measures of model validation considered in our first example are exhibited. Among the composite models, we can see that the Weibull- MPLG composite model gives the best fit overall for this dataset in terms of NLL, AIC and BIC values.
Finally, to select composite models that provide an acceptable description of the loss process, we must verify that the first-order moment of the limited loss variable and the empirical counterpart, given by E n ( u ) = 1 n i = 1 n min ( x i , u ) are essentially in agreement. Obviously when u tends to infinity, the former quantity and E n ( u ) converge to E ( X ) and the sample mean, respectively. In Table 3 empirical and fitted limited expected value for the seven composite models considered are displayed. As it can be observed, the composite lognormal-Pareto and Weibull-Pareto distributions tend to overestimate the empirical limited expected value when the policy limit u increases. Although, the remainder models stay closer to the empirical limited expected value for different values of u, the composite Weibull-Lomax displays the best behaviour.

4.3. Theil’s Income Indices

Among several economic inequality measures the Gini Index is commonly used. However, it does not allow us to vary the sensitivity of the index to redistributional movements at specific income ranges, which can point out different evolutions of disparities when no Lorenz dominance is found (see Sarabia and Castillo 2005). Yitzhaki (1983) further generalized this index which is sensitive to changes in the right tail of the distribution as the parameter increases. The two limiting cases of this family of inequality measures, known as Theils’s Indices ( T 1 and T 0 ), corresponding to the Theil entropy index (TEI) and the mean log deviation (MLD), are defined by
T 1 ( X ) = E log X μ and T 2 ( X ) = E X μ log X μ .
Proposition 5.
The TEI and MLD indices for a rv X that follows the MPLG ( θ , λ , x 0 ) distribution are provided by
T 1 ( X ) = θ + 2 λ θ ( θ + λ ) + log θ 2 ( θ + λ 1 ) ( θ 1 ) 2 ( θ + λ ) , θ > 1
and
T 2 ( X ) = 2 θ 1 1 θ + λ 1 log θ 2 ( θ + λ 1 ) ( θ 1 ) 2 ( θ + λ ) , θ > 1 .
Proof. 
The proof follows directly after some computation. □

5. Conclusions

In this paper, we have proposed a new heavy-tailed class of distribution, that is obtained by using an exponential transformation of the generalized Lindley distribution. This class generalizes both the classical Pareto model and a special case of the loggamma distribution. Some of its most relevant statistical properties were examined. The model is very flexible since it allows for closed-form expressions for many results related to insurance. Besides, as the mode of this new model can be written in closed-form, composite models based on this distribution can be easily derived. The numerical illustrations reveal that the mixture Pareto-loggamma distribution provides a good fit to loss data, and it is a competitive model with other existing heavy-tailed distributions in the literature.

Author Contributions

The authors contributed equally to this work.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Expression (1) is a genuine pdf:
f X ( x | θ , λ , x 0 ) = θ θ + λ θ x x x 0 θ + λ θ + λ θ 2 x x x 0 θ log x x 0 = θ θ + λ θ x 0 θ x θ + 1 + λ θ + λ θ 2 x 0 x x 0 ( θ + 1 ) log x x 0 = θ θ + λ g 1 ( x ) + λ θ + λ g 2 ( x ) ,
where g 1 ( x ) and g 2 ( x ) are the densities of the Pareto and loggamma distribution respectively.
Proof of Proposition 1.
Let X i be a MPLG rv’s with parameter ( θ i , λ i , x 0 i ) , i = 1 , 2 . Then
d d x log f X 1 ( x ) f X 2 ( x ) = λ 1 x 1 + λ 1 log x x 01 λ 2 x 1 + λ 2 log x x 02 θ 1 x + θ 2 x
it is easy to see that, if x 01 = x 02 , (A1) is negative when θ 1 θ 2 and λ 1 λ 2 which implies that X 1 l r X 2 . Further note that, when θ 1 = θ 2 and λ 1 = λ 2 , (A1) is negative when x 01 x 02 and therefore X 1 l r X 2 . The other results follows from (9). □
Proof of Proposition 2.
By assuming p = log x x 0 , the cdf can be written as
F X ( x ) = 1 θ + λ + θ λ p θ + λ exp { θ p } ,
for fixed θ , λ > 0 and δ ( 0 , 1 ) , the δ th quantile function is obtained by solving F X ( x ) = δ . By re-arranging the above, we obtain
( 1 δ ) ( θ + λ ) = ( θ + λ + θ λ p ) exp { θ p }
Now, by multiplying both sides of (A2) by exp ( θ + λ ) / λ , we obtain
θ + λ λ + θ p exp θ + λ λ + θ p = ( 1 δ ) ( θ + λ ) λ exp θ + λ λ .
From expression (A3), we see that θ + λ λ + θ p is the Lambert-W function of real argument ( 1 δ ) ( θ + λ ) λ exp θ + λ λ . Thus, we have
W ( 1 δ ) ( θ + λ ) λ exp θ + λ λ = θ + λ λ + θ p .
Moreover, for any θ , λ > 0 , it is immediate that θ + λ λ + θ p > 1 , and it can also be checked that ( 1 δ ) ( θ + λ ) λ exp θ + λ λ ( 1 / e , 0 ) since δ ( 0 , 1 ) . Therefore, by taking into account the properties of the negative branch of the Lambert-W function, we deduce the following
W 1 ( 1 δ ) ( θ + λ ) λ exp θ + λ λ = θ + λ λ + θ p .
Again, substituting p = log x x 0 , and solving for x, we obtain
VaR δ ( X ) = x 0 exp θ + λ θ λ 1 θ W 1 ( 1 δ ) ( θ + λ ) . exp θ + λ λ λ .
This completes the proof of Proposition 2. □
Proof of Proposition 3.
Let X follows MPLG ( θ , λ , x 0 ) . Then, the T V a R δ ( X ) is defined as
TVaR δ ( X ) = 1 1 δ δ 1 VaR z ( X ) d z ,
substituting the value VaR z ( X ) from (18), we get
TVaR δ ( X ) = 1 1 δ δ 1 x 0 exp θ + λ θ λ 1 θ W 1 ( 1 z ) exp θ + λ λ ( θ + λ ) λ d z = x 0 1 δ exp ( θ + λ ) ( 1 θ ) θ λ δ 1 exp 1 θ W 1 ( 1 z ) exp θ + λ λ ( θ + λ ) λ d z .
Now, by denoting ( 1 z ) e θ λ 1 ( θ + λ ) λ = u , we get
TVaR δ ( X ) = x 0 λ 1 δ exp ( θ + λ ) ( 1 θ ) θ λ 0 δ exp 1 θ W 1 ( u ) d u ,
where δ = ( 1 δ ) e θ λ 1 ( θ + λ ) λ . By letting W 1 ( u ) = t , it gives t exp { t } = u , and d u = exp { t } t + 1 d t , then we have
TVaR δ ( X ) = x 0 λ 1 δ exp ( θ + λ ) ( 1 θ ) θ λ 0 W 1 ( δ ) exp 1 θ 1 ( 1 + t ) d t = 1 ( 1 δ ) λ x 0 θ exp ( 1 θ ) ( θ + λ ) θ λ ( 1 θ ) 2 ( θ + λ ) × ( 1 + ( 1 θ ) K ( δ ; θ , λ ) ) exp ( 1 θ ) K ( δ ; θ , λ ) θ 1 ,
where K ( δ ; θ , λ ) = W 1 ( 1 δ ) exp θ + λ λ ( θ + λ ) λ . Hence the proposition. □

References

  1. Acerbi, Carlo, and Dirk Tasche. 2002. On the coherence of expected shortfall. Journal of Banking & Finance 26: 1487–503. [Google Scholar]
  2. Bakar, Shaiful A., Nor A. Hamzah, Mastoureh Maghsoudi, and Saralees Nadarajah. 2015. Modeling loss data using composite models. Insurance: Mathematics and Economics 61: 146–54. [Google Scholar]
  3. Calderín-Ojeda, Enrique, and Chun F. Kwok. 2016. Modeling Claims Data with Composite Stoppa Models. Scandinavian Actuarial Journal 9: 817–36. [Google Scholar] [CrossRef]
  4. Cooray, Kahadawala, and Malwane M. A. Ananda. 2005. Modeling actuarial data with a composite lognormal–Pareto model. Scandinavian Actuarial Journal 5: 321–34. [Google Scholar] [CrossRef]
  5. De Jong, Piet, and Heller Gillian H. 2008. Generalized Linear Models for Insurance Data. International Series on Actuarial Science; Cambridge: Cambridge University Press. [Google Scholar]
  6. Ghitany, Mohamed E., Barbra Atieh, and Saralees Nadarajah. 2008. Lindley distribution and its application. Mathematics and Computers in Simulation 78: 493–506. [Google Scholar] [CrossRef]
  7. Ghitany, Mohamed E., Emilio Gómez-Déniz, and Saralees Nadarajah. 2018. A new generalization of the Pareto distribution and its application to insurance data. Journal of Risk and Financial Management 11: 10. [Google Scholar] [CrossRef]
  8. Gómez-Déniz, Emilio, and Enrique Calderín-Ojeda. 2014. A Suitable Alternative to the Pareto Distribution. Hacettepe Journal of Mathematics and Statistics 43: 843–60. [Google Scholar]
  9. Gómez-Déniz, Emilio, and Enrique Calderín-Ojeda. 2015. Modeling insurance data with the Pareto ArcTan distribution. Astin Bulletin 45: 639–60. [Google Scholar] [CrossRef]
  10. Hogg, Robert V., and Stuart A. Klugman. 2009. Loss Distributions. Vol. 249 of Wiley Series in Probability and Statistics; New York: Wiley. [Google Scholar]
  11. Klugman, Stuart A., Harry H. Panjer, and Gordon E. Willmot. 2012. Loss Models: From Data to Decisions, 4th ed. Wiley Series in Probability and Statistics; Hoboken: John Wiley and Sons, Inc. [Google Scholar]
  12. Sarabia, José M., and Enrique Castillo. 2005. About a class of max–stable families with applications to income distributions. Metron LXIII: 505–27. [Google Scholar]
  13. Sarabia, José M., and Faustino Prieto. 2009. The Pareto-positive stable distribution: A new descriptive model for city size data. Physica A 388: 4179–91. [Google Scholar] [CrossRef]
  14. Scollnik, David P. M. 2007. On Composite Lognormal-Pareto Models. Scandinavian Actuarial Journal 1: 20–33. [Google Scholar] [CrossRef]
  15. Shaked, Moshe, and J. George Shanthikumar. 2007. Stochastic Orders. New York: Springer Science and Business Media. [Google Scholar]
  16. Stoppa, Gabriele. 1990. Proprietà campionarie di un nuovo modello Pareto generalizzato. In Atti XXXV Riunione Scientifica della Società Italiana di Statistica. Padova: Cedam, pp. 137–44. [Google Scholar]
  17. Su, Chun, and Qihe H. Tang. 2003. Characterizations on heavy-tailed distributions by means of hazard rate. Acta Mathematicae Applicatae Sinica 19: 135–42. [Google Scholar] [CrossRef]
  18. Yang, Hailiang. 2004. Cramér-Lundberg asymptotics. In Encyclopedia of Actuarial Science. New York: Wiley. [Google Scholar]
  19. Yitzhaki, Shlomo. 1983. On an Extension of the Gini Inequality Index. International Economic Review 24: 617–28. [Google Scholar]
Figure 1. Graphs of the pdf (1) for different values of parameter θ , λ and fixed value of x 0 .
Figure 1. Graphs of the pdf (1) for different values of parameter θ , λ and fixed value of x 0 .
Risks 07 00099 g001
Figure 2. Graphs of the hazard rate function (4) for different values of parameter θ , λ and for fixed value of x 0 .
Figure 2. Graphs of the hazard rate function (4) for different values of parameter θ , λ and for fixed value of x 0 .
Risks 07 00099 g002
Table 1. Estimated values of different heavy-tailed models, corresponding standard errors (S.E.), negative of the log-likelihood function (NLL), Akaike’s information criterion (AIC) and Bayesian information criterion (BIC) for automobile claims dataset (De Jong and Heller 2008).
Table 1. Estimated values of different heavy-tailed models, corresponding standard errors (S.E.), negative of the log-likelihood function (NLL), Akaike’s information criterion (AIC) and Bayesian information criterion (BIC) for automobile claims dataset (De Jong and Heller 2008).
DistributionParameter Estimates (S.E.)NLLAICBIC
Pareto θ ^ = 0.661 ( 0.010 ) 38,024.8076,051.6176,058.05
lognormal μ ^ = 6.810 ( 0.017 ) , σ ^ = 1.189 ( 0.012 ) 38,852.1577,708.3177,721.19
loggamma α ^ = 1.115 ( 0.021 ) , β ^ = 5.109 ( 0.118 ) 38,998.1878,000.3678,013.23
Frećhet α ^ = 0.659 ( 0.019 ) , s ^ = 254.0 ( 11.87 ) 38,408.3076,822.6076,841.92
b ^ = 157.1 ( 4.279 )
Weibull μ ^ = 0.786 ( 0.008 ) , σ ^ = 1690.8 ( 33.64 ) 39,491.6078,987.1979,000.07
Lomax α ^ = 2.047 ( 0.088 ) , λ ^ = 2205.1 ( 133.1 ) 39,169.8578,343.7078,356.58
PAT α ^ = 0.895 ( 0.095 ) , θ ^ = 0.740 ( 0.017 ) 38,006.3076,016.6176,029.49
Inverse Weibull α ^ = 1.053 ( 0.012 ) , σ ^ = 518.8 ( 7.636 ) 38,595.6177,195.2277,208.09
MPLG θ ^ = 0.943 ( 0.018 ) , λ ^ = 0.698 ( 0.073 ) 37,965.9975,935.9875,948.86
Table 2. Estimated values of different heavy-tailed models, corresponding standard errors (S.E.), NLL, AIC, and BIC of different composite models for the Danish fire losses dataset.
Table 2. Estimated values of different heavy-tailed models, corresponding standard errors (S.E.), NLL, AIC, and BIC of different composite models for the Danish fire losses dataset.
DistributionParameter Estimates (S.E.)NLLAICBIC
lognormal-Lomax μ ^ = 0.104 ( 0.020 ) σ ^ = 0.182 ( 0.011 ) 3860.477728.947752.223
λ ^ = 0.365 ( 0.123 ) θ ^ = 1.144 ( 0.029 )
lognormal-Pareto μ ^ = 0.137 ( 0.019 ) σ ^ = 0.197 ( 0.012 ) 3865.867737.727755.182
θ ^ = 1.208 ( 0.030 )
lognormal- MPLG μ ^ = 0.045 ( 0.017 ) θ ^ = 2.060 ( 0.039 ) 3872.407752.817776.09
x ^ 0 = 0.745 ( 0.070 ) λ ^ = 65.804 ( 377.425 )
Weibull-Lomax τ ^ = 15.345 ( 0.671 ) ϕ ^ = 0.969 ( 0.007 ) 3823.707655.407678.68
λ ^ = 0.561 ( 0.127 ) θ ^ = 0.971 ( 0.007 )
Weibull-Pareto τ ^ = 14.048 ( 0.502 ) ψ ^ = 0.997 ( 0.008 ) 3840.387686.767704.22
θ ^ = 1.003 ( 0.008 )
Weibull- MPLG τ ^ = 18.763 ( 0.986 ) x ^ 0 = 0.787 ( 0.006 ) 3823.307654.607677.88
θ ^ = 1.938 ( 0.034 ) λ ^ = 4.614 ( 0.159 )
paralogistic- MPLG α ^ = 16.719 ( 1.311 ) x ^ 0 = 0.901 ( 0.020 ) 3824.4827656.967680.25
θ ^ = 1.989 ( 0.035 ) λ ^ = 3.379 ( 0.266 )
Table 3. Limited expected value for the distributions considered and different values of the policy limit u.
Table 3. Limited expected value for the distributions considered and different values of the policy limit u.
Composite Models
PolicyEmpiricallognormallognormalWeibullWeibulllognormalWeibullParalogistic
Limit u ParetoLomaxParetoLomax MPLG MPLG MPLG
10.9890.9980.9850.9890.9890.9870.9860.989
21.5651.7261.5251.5601.5481.3331.2661.467
52.1382.4892.1222.1552.1432.0151.9432.093
102.4472.9322.3622.4822.4642.3412.2922.398
202.7073.2852.5662.7812.6742.5422.5222.591
1002.9583.8512.8263.2992.9232.7402.7782.793
2003.0374.0172.8843.4632.9722.7702.8232.826
3.0634.6653.0054.2893.0602.8042.8852.868

Share and Cite

MDPI and ACS Style

Bhati, D.; Calderín-Ojeda, E.; Meenakshi, M. A New Heavy Tailed Class of Distributions Which Includes the Pareto. Risks 2019, 7, 99. https://doi.org/10.3390/risks7040099

AMA Style

Bhati D, Calderín-Ojeda E, Meenakshi M. A New Heavy Tailed Class of Distributions Which Includes the Pareto. Risks. 2019; 7(4):99. https://doi.org/10.3390/risks7040099

Chicago/Turabian Style

Bhati, Deepesh, Enrique Calderín-Ojeda, and Mareeswaran Meenakshi. 2019. "A New Heavy Tailed Class of Distributions Which Includes the Pareto" Risks 7, no. 4: 99. https://doi.org/10.3390/risks7040099

APA Style

Bhati, D., Calderín-Ojeda, E., & Meenakshi, M. (2019). A New Heavy Tailed Class of Distributions Which Includes the Pareto. Risks, 7(4), 99. https://doi.org/10.3390/risks7040099

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