Abstract

In this paper, a family of statistical models, namely, a new exponential-X family is proposed. A subcase of the introduced family, called the new exponential-Weibull (NE-Weibull) model, is studied. The NE-Weibull model is very competent and possesses heavy-tailed properties. The maximum likelihood estimators of its parameters are derived. The consistency and efficiency of these estimators are assessed in a brief simulation study. Finally, the effectiveness of the NE-Weibull distribution is illustrated by modeling real insurance claims data. The practical analysis shows that the NE-Weibull distribution outclassed other distributions and it can be a better choice for modeling data in the finance sector.

1. Introduction

The extreme value phenomena such as financial returns and other related events can be modeled effectively by extreme value methods. The heavy-tailed (HT) distributions have been proven to be substantial in modeling HT and extreme value data. Researchers have shown a deep concern in the financial sector to study new HT distributions. Among the applicability of the statistical distributions in the applied area, the HT distributions have received much attention for modeling financial phenomena. Most of the data sets in the financial sector possesses HT behavior with long right tail (see Lane [1]; Cooray and Ananda [2]; Wang et al. [3]; Ahn et al. [4]; Jelenković and Tan [5]; Forbes and Wraith [6]; Guo [7]; Punzo et al. [8]; Bhati and Ravi [9]; Punzo [10]; Ke et al. [11]; Dos Reis et al. [12]; and Niu et al. [13]).

Recently, different copulas and new approaches of introducing HT distribution have been studied due to the importance of HT distributions in the financial sector (for example, Bladt et al. [14]; Tikhomirov [15]; Lugosi and Mendelson [16]; and Yousri et al. [17]). For more information about the usefulness of statistical distributions, one can refer to studies by Ramos et al. [18, 19]; He et al. [20]; Alfaer et al. [21]; Afify et al. [22]; and Al Mutairi et al. [23].

A statistical distribution with SF (survival function), say , is said to be a HT model, if its SF verifiesfor all . For more details, see the study by Resnick [24].

An interesting characteristic of the HT distributions is the regularly varying behavior. A statistical distribution is said to possess the regularly varying behavior, if it satisfieswhere , and it is also known as an index of regular variation.

The statistical models possessing such property are very prominent models for modeling HT phenomena in the financial sector [25]. Furthermore, actuaries are very interested in looking for new flexible HT models (see Nadarajah and Bakar [26]; and Ahmad et al. [27]).

Alzaatreh et al. [28] proposed a prominent approach called the T-X family method. Let represent the PDF (probability density function) of T, where T is a RV (random variable) belonging to and . Suppose that is a function of of a RV-X and it has the following conditions:(i)(ii) is a differentiable function as well as monotonically increasing(iii) for and for

The CDF (cumulative distribution function) of the T-X distributions iswhere fulfills the aforementioned conditions. The corresponding PDF of equation (3), say , reduces to

More details about the T-X approach can be found in the work of Ahmad et al. [29]. By implementing the T-X method, the survival distribution family [30] can be obtained via the CDF:where , representing the SF of X.

In this study, a new exponential-X (NE-X) family is proposed based on equation (3). Suppose ; then, its respective CDF is

Linking this to equation (6), the PDF isBy using and in equation (3), the CDF of the NE-X family is as follows:where is a baseline CDF with parametric space . Next, in Propositions 1 and 2, we are going to prove that the expression provided in equation (8) is a CDF.

Proposition 1. For the defined in equation (8), and .

Proof.

Proposition 2. The CDF is right continuous and differentiable.

Proof. From the results proved in Propositions 1 and 2, we arrive at the conclusion that the function defined in equation (8) is differentiable, right continuous, and a compact CDF. Moreover, taking into account the differentiability property of for all , we have the following theorem.

Theorem 1. Let and . With the condition , if and are differentiable, then the quotient is differentiable for all and

Proof. See “Proof of Quotient Rule.”
Taking into account the right continuity of , we have the following theorem.

Theorem 2. If is a right continuous function, so is .

Proof. Assume that for all . Hence, we haveSince ,for each . Therefore, is right continuous.
As we mentioned earlier that the regularly varying tail behavior (RVTB) is a very crucial property to characterize the HT distributions, now we provide the mathematical treatment of RVTB of the NE-X family.

Theorem 3. If is a regularly varying distribution, so is .

Proof. Suppose is finite for all . Then, using equation (8), we haveSince , we havefor ; thus, is a regularly varying distribution.
Linking this to equation (8), the PDF represented by isThe RV with PDF (16) is represented by -.
The SF, , and HRF (hazard rate function), , of X arerespectively. In the next section, we discuss a subcase of the NE-X family, called the new exponential-Weibull (NE-Weibull) distribution.

2. NE-Weibull Distribution

This section deals with the NE-Weibull distribution as a subcase of the NE-X family. Suppose that with CDF, , and PDF, , where . By inserting the CDF of the Weibull model in equation (8), the CDF and PDF of the NE-Weibull distribution take the forms

Some possible PDF behaviors for the NE-Weibull distribution are shown in Figure 1. From Figure 1, it is obvious that as the values of and increase, the NE-Weibull model becomes a HT distribution.

As mentioned above, different PDF shapes of the NE-Weibull model for various values of , , and are sketched in Figure 1. When , the NE-Weibull distribution behaves similar to the exponential distribution . For , the NE-Weibull distribution captures the characteristics of . The NE-Weibull model, however, has certain benefits over the Weibull distribution. For instance, it has heavier tails than the Weibull distribution and provides closer fit to data in the finance sector.

Some possible HRF behaviors for the NE-Weibull distribution are shown in Figure 2. The NE-Weibull distribution provides increasing, unimodal, decreasing, and bathtub HRF shapes.

3. Properties

Here, we provide a concise treatment of the mathematical properties of NE-X distribution. These properties include QF (quantile function), moments, SK (skewness), and KUR (kurtosis). Furthermore, different plots for the SK and KUR are also provided.

3.1. Quantile Function

Suppose X denotes the NE-X family with CDF presented in equation (8); then, the QF of NE-X distributions, , iswhere is the solution of .

3.2. Moments

Now, we introduce the moments of the NE-X distributions which can further be used to obtain other characteristics.

The moment of the NE-X family reduces to

By putting equation (16) in (20), we obtainwhere . For , we obtain the first four moments of the NE-X family.

Using the expressions of the moments, we obtain the mathematical form of the SK and KUR measures. The SK and KUR of the NE-Weibull distribution can be calculated, respectively, via the expressions

The effects of , , and on the SK, KUR, variance, and mean of the NE-Weibull distribution are displayed in Figures 35.

4. Estimation

Within this section, we derive the maximum likelihood estimators (MLEs) of the NE-Weibull parameters. Consider as the values of a sample from the NE-Weibull distribution with parameters and . The log-likelihood (LL) function of the NE-Weibull distribution takes the form

The partial derivatives of the LL function are

Equating and to zero and simultaneously solving them yield the MLEs of .

5. Simulation Study

Here, we implement the Monte Carlo simulation approach to address the MLE behavior in estimating the NE-Weibull parameters. The NE-Weibull distribution can be simulated by using equation (8). Let follow the standard uniform distribution; hence, the quantile function reduces to

The simulation is done for (i) , (ii) , and (iii) .

The simulation results are obtained by utilizing the software with the algorithms and “” with . The results obtained are based on replications for samples of size , where . Statistical tools such as biases and mean square errors (MSEs) are obtained as assessing tools. These tools are calculated as follows:where .

The summary measures (SMs) of the three simulated sets of data are provided in Table 1, whereas the histograms, box plots, kernel density estimator, and estimated CDF of the NE-Weibull model are provided in Figures 68. The simulation results are visually displayed in Figures 911.

6. Data Modeling in the Finance Sector

Here, we illustrate the importance of the NE-Weibull distribution in modeling insurance claims data from the financial sector. We also calculated the risk measures such as VaR (value at risk) and TVaR (tail value at risk) for this data.

6.1. Insurance Claims Data

The insurance claims data represents the initial claims of the unemployment insurances per month from 1971 to 2018 [20]. The data set can also be retrieved from https://data.worlddatany-govns8z-xewg.

The summary statistics for the insurance claims data are reported in Table 2. Corresponding to this data, the histogram, box plot, and total time test (TTT) plots are provided in Figure 12. Figure 12 shows that the insurance claims data is unimodal, right-skewed, and HT.

The NE-Weibull distribution is applied to fit the financial data, and it is compared with the Weibull, W-Loss (Weibull-Loss), NHT-Weibull (new heavy-tailed Weibull), Ku-Weibull (Kumaraswamy Weibull), and B-Weibull (beta Weibull) distributions. The CDFs of these models are as follows:(i)Weibull model:(ii)W-Loss model:(iii)NHT-Weibull model:(iv)Ku-Weibull model:(v)B-Weibull model:

To decide about the best fitting of the applied distributions, certain criteria are taken into account. These criteria are as follows:(i)The AIC is(ii)The BIC is(iii)The HQIC is(iv)The CAIC is

Here, denotes the LL function, is a parametric space, represents the model parameters, and is the number of selected samples. In addition to the criteria measures, three goodness-of-fit tests with its corresponding values are also considered. These tests are given by the following:(i)The Anderson–Darling (AD) test statistic iswhere is the observation in samples, calculated after sorting the data in the ascending order.(ii)The Cramer–von Mises (CM) test statistic is(iii)The Kolmogorov–Smirnov (KS) test statistic iswhere represents the supremum of the set of distances and indicates the empirical CDF.

The LL function is optimized, and the goodness-of-fit measures are calculated using the algorithm method “BFGS” via optim() R-function [31]. The MLEs of the fitted models with standard errors are listed in Table 3. Tables 4 and 5 show the goodness-of-fit values for the considered models.

The fitted PDF and CDF of the NE-Weibull model are provided in Figure 13, showing that the right-skewed HT data can be adequately fitted by the NE-Weibull distribution. The Kaplan–Meier survival (KMS) and PP (probability-probability) plots of the NE-Weibull distribution are sketched in Figure 14.

6.2. VaR and TVaR for Insurance Claims Data

Now, we calculate the VaR and TVaR measures for the Weibull and proposed NE-Weibull distributions based on the MLEs and the data provided in Section 6.1. The results of the VaR and TVaR of the NE-Weibull and Weibull distributions are presented in Table 6. We also display these results graphically in Figure 15.

A model having higher values for the TVaR and VaR is considered a HT distribution. The numerical results in Table 6 and their graph in Figure 15 show that the NE-Weibull is a HT model.

7. Final Comments

During the last couple of years, the study and modeling of extreme value data has gained increased interest in numerous areas, particularly in the hydrology and finance sector. Recent studies have shown the importance and potential of the HT distributions in various areas for analyzing real data sets, especially in the finance sector. In this article, a new interesting HT extension of the Weibull distribution, called the NE-Weibull distribution, is proposed using the NE-X approach. The MLEs of the NE-Weibull parameters and are obtained. To show the applicability of the NE-Weibull model, a real-life insurance claims data from the financial sector is considered. Based on certain tests, it is showed that the NE-Weibull model offers an adequate fit to the insurance claims data than the Weibull, W-Loss, NHT-Weibull, Ku-Weibull, and B-Weibull. On the other hand, the numerical results for the VaR and TVaR measures showed that the NE-Weibull model is a HT distribution as compared to the classical Weibull distribution.

Data Availability

To replicate the findings of this paper, the data is available from the corresponding author upon request.

Disclosure

This article was drafted from the PhD work of the first author (Zubair Ahmad).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this work.

Authors’ Contributions

Conceptualization and methodology were carried out by Zubair Ahmad and Eisa Mahmoudi. The original draft was prepared by Zubair Ahmad and Morad Alizadeh. The original draft was edited by Eisa Mahmoudi, Rasool Roozegar, and Ahmed Z. Afify. Formal analysis was performed by Zubair Ahmad. Software was collected by Zubair Ahmad and Morad Alizadeh. The study was supervised by Eisa Mahmoudi, Morad Alizadeh, and Rasool Roozegar. Investigation was conducted by Zubair Ahmad. The manuscript was reviewed and edited by Zubair Ahmad, Eisa Mahmoudi, Morad Alizadeh, Rasool Roozegar, and Ahmed Z. Afify.

Acknowledgments

This work was supported by the Department of Statistics, Yazd University, Yazd, Iran.