Abstract

In this study, a three-parameter modification of the Burr XII distribution has been developed through the integration of the weighted version of the alpha power transformation family of distributions. This newly introduced model, termed the modified alpha power-transformed Burr XII distribution, exhibits the unique ability to effectively model decreasing, right-skewed, or unimodal densities. The paper systematically elucidates various statistical properties of the proposed distribution. The estimation of parameters was obtained using maximum likelihood estimation. The estimator has been evaluated for consistency through simulation studies. To gauge the practical applicability of the proposed distribution, two distinct datasets have been employed. Comparative analyses involving six alternative distributions unequivocally demonstrate that the modified alpha power-transformed Burr XII distribution provides a better fit. Additionally, a noteworthy extension is introduced in the form of a location-scale regression model known as the log-modified alpha power-transformed Burr XII model. This model is subsequently applied to a dataset related to stock market liquidity. The findings underscore the enhanced fitting capabilities of the proposed model in comparison to existing distributions, providing valuable insights for applications in financial modelling and analysis.

1. Introduction

The development of novel or new statistical models is a key area of study in the application of distribution theory. These distributions’ usefulness has led to much research into their theory and the development of new distributions. The idea of generating new continuous distributions by modifying the existing distributions with one or more shape or scale parameters has gained attention in recent years. This parameter introduction has been shown to improve the ability of the developed distributions to fit varied real-life datasets with high degrees of skewness and kurtosis. Some of these newly developed distributions include the modified alpha power transformed Weibull [1], general two-parameter [2], truncated inverse power Ailamujia [3], half-logistic modified Kies exponential [4], truncated inverse power Lindley [5], Marshall–Olkin–Weibull–Burr XII [6], generalised unit half-logistic geometric [7], Chen Burr–Hatke exponential [8], modified XLindley [9], arctan power [10], harmonic mixture Fréchet [11], sine-Weibull geometric [12], bounded odd inverse Pareto exponential [13], new extended Chen [14], power XLindley [15], extended Poisson–Fréchet [16], exponentiated Fréchet loss [17], Gompertz–Makeham [18], and logistic exponential [19] distributions.

The authors of [20] introduced a new method by adding an additional parameter called the alpha power transformation (APT) family. The APT family has been used to develop several modified distributions, including the APT Fréchet [21], APT extended exponential distribution [22], APT inverse Lomax distribution [23], APT log-logistic distribution [24], APT inverse Lindley distribution [25], and APT Pareto distribution [26], among others. With the aim of improving the flexibility of the APT family of distributions, Alotaibi et al. [27] modified the APT family of distributions and obtained a new family of distributions called the modified alpha power transformed method (MAPT). According to Alotaibi et al. [27], the cumulative distribution function (CDF) and probability density function (PDF) of the MAPT are given as follows:and

This study proposes a modification of the Burr XII distribution using the MAPT proposed byAlotaibi et al. [27]. The Burr XII distribution introduced by Burr [28] is widely used in reliability analysis, actuarial studies, medicine, and agriculture. The PDF and CDF of the Burr XII distribution can be expressed, respectively, as follows:and

We are motivated to contribute to the ongoing efforts to enhance the versatility of statistical distributions, thus providing researchers with a powerful tool to analyse and model diverse data scenarios effectively as no single distribution is omnibus. Specifically, our motivations for developing the modified alpha power transformed Burr XII (MAPTBXII) distribution are as follows:(i)Develop an extension of the Burr XII distribution that provides a good parametric fit to data with complex traits(ii)Propose a new Burr XII distribution with closed form CDF and tractable quantile function that facilitates easy generation of random observations for simulation experiments(iii)Formulate a location-scale regression model using the proposed distribution

The subsequent sections of the paper are organised as follows. In Section 2, we develop the MAPTBXII distribution. Section 3 is devoted to deriving various statistical properties of the MAPTBXII distribution. The parameters of the MAPTBXII distribution are estimated through the maximum likelihood estimation method, as detailed in Section 4. Section 5 delves into the discussion of the MAPTBXII regression model. A comprehensive simulation study is presented in Section 6. Real-world applications of the MAPTBXII distribution on two datasets are presented in Section 7. Finally, the paper is concluded in Section 8.

2. Modified Alpha Power Transformed Burr XII Distribution

If a random variable X follows the MAPTBXII, then the PDF can be obtained by substituting equation (3) into equation (2), while the CDF can be obtained by substituting equation (4) into equation (1). The CDF and PDF of the MAPTBXII can then be expressed, respectively, as follows:where , , , and .

The hazard function of the MAPTBXII is obtained by finding the ratio of the PDF and complement of the CDF. The hazard function is given by the following equation:

The various shapes (decreasing, right-skewed, left-skewed, or unimodal) of the densities of the MAPTBXII distribution are shown in Figure 1.

The various shapes (decreasing or upside down bathtub) of the hazard function are displayed in Figure 2.

The PDF of the MAPTBXII distribution in linear form is given as follows:where , , , and .

The asymptotic nature of the CDF, PDF, and hazard function is as follows:, as and , as , as and , as , as and , as

The quantile function of MAPTBXII distribution is given as follows:where and is the quantile function.

3. Statistical Properties

3.1. Moments

Moments are a crucial component of statistical theory, and they allow for the examination of many essential characteristics of any distribution.

Mathematically, the moments of are given by

We substitute equation (7) into equation (9), in order to obtain

Letting , which implies and , we obtain the following expression:

Utilising the identity used in [29],then gives

It follows that the moment of the MAPTBXII distribution can be expressed as follows:where is given as the beta function and r = 1, 2, ….

Metrics such as variance , coefficient of variation (CVAR), skewness (CSK), and kurtosis (CKUR) can be derived by utilizing moment-based computations. , , CVAR, CSK, and CKUR, respectively, are obtained using the following expression:

Table 1 displays the , CVAR, CSK, and CKUR of the MAPTBXII distribution, calculated using noncentral moments for specific parameter values. The results suggest that the MAPTBXII distribution exhibits a varying degree of skewness, ranging from highly skewed to moderately skewed. It is noteworthy that the skewness of the distribution may differ depending on the chosen parameter values. Some combinations lead to a positively skewed distribution, while others result in a negatively skewed distribution. This variability underscores the distribution’s flexibility in capturing different skewness patterns based on the specified parameters. The MAPTBXII distribution is leptokurtic (CKUR > 3) and thus has heavier tails and a sharper peak.

3.2. Incomplete Moments

The incomplete moments can be used to calculate the Lorenz curve, the Bonferroni curve, the mean deviation, and the median deviation, among others.

Mathematically, the incomplete moment is given as follows:

We substitute equation (7) into equation (16), in order to obtain

Letting , which implies and , we obtain the following expression:

Utilising the identity used in [29],helps to obtain

It follows that the incomplete moment of the MAPTBXII distribution is given by the following expression:where is given as the incomplete beta function and .

3.3. Moment Generating Function

The moment generating function (MGF) is used to determine the distribution’s moments, if any. Mathematically,

When we substitute equation (14) into equation (22), we obtain the MGF of the MAPTBXII distribution. Hence, the MGF of the MAPTBXII distribution is given by the following expression:

4. Estimation of Parameters

The maximum likelihood estimation (MLE) method is used as an estimator. The log-likelihood function is therefore given as

By maximising equation (24), we obtain the parameter estimates.

5. Modified Alpha Power Transformed Burr XII Regression Model

Suppose follows the MAPTBXII distribution and . The PDF of can be derived by substituting and .where , is the scale parameter, is the shape parameter, , and is the location parameter.

Equation (25) is referred to as the log-MAPT Burr XII (LMAPTBXII) distribution. If , then .

Figure 3 shows various shapes of the densities of the LMAPTBXII distribution.

The survival function of the LMAPTBXII distribution can be expressed as follows:

The proposed location-scale regression model is defined with the dependent variable and predictor variables , where is the intercept and expressed as follows:where , are the regression parameters, and denotes the random error.

The log-likelihood function of the LMAPTBXII model is given by the following expression:

The parameter estimates of the model are obtained by maximising the log-likelihood function. An assessment of the Cox–Snell residuals to ascertain if they behave as a standard exponential distribution would help determine the adequacy of the model. We then diagnose the model using the goodness-of-fit measures (Cramér–von Mises, Anderson–Darling, and Kolmogorov–Smirnov) of the Cox–Snell residuals.

6. Monte Carlo Simulations

In this section, we ascertain the consistency of the estimators of the MAPTBXII distribution through a simulation study. The results were obtained using sample sizes of 50, 100, 200, 250, 300, 350, 500, and 600 with parameter values , , and , respectively. It can be observed that the average biases (AB) and root mean square error (RMSE) decrease as the sample size increases as shown in Tables 2.

7. Applications

In this section, we provide the applications of the MAPTBXII distribution using two uncensored datasets.

The precipitation (in inches) in the Minneapolis dataset was used in [30, 31]. The second dataset represents runoff amounts at Jug Bridge, Maryland, and was used by Makubate et al. [32]. The performance of the MAPTBXII distribution is compared with the Topp–Leone Burr-XII distribution (TLBXII) [33], Weibull Burr XII distribution (WBXII) [29], inverse Weibull distribution (INW) [34], exponentiated exponential Burr XII distribution (EEBXII) [35], Cauchy [36], and inverse Gompertz distribution (IGD) [37].

Analytical measures such as Kolmogorov–Smirnov (K-S) test, Anderson–Darling test (AD), and Cramér–von Mises test (CVM), Akaike Information Criterion (AIC), Consistent Akaike Information Criterion (AICc), and Bayesian Information Criterion (BIC) are considered in evaluating the goodness of fit of the proposed distribution and the other fitted models.

The maximum likelihood estimates and other goodness-of-fit statistics for the two datasets are presented in Tables 3 and 4.

The MLEs of the MAPTBXII model for datasets 1 and 2 are both unique and represent real maxima, as demonstrated by the profile log-likelihood plots shown in Figures 4 and 5.

The fitted densities are presented in Figures 6 and 7, whereas the fitted CDFs are shown in Figures 8 and 9. It can be observed from these results that the MAPTBXII provides a better fit to the two datasets than the other competing lifetime distributions.

7.1. Application of Modified Alpha Power Transformed Burr XII Regression Model

The LMAPTBXII model was employed in the analysis of a dataset related to stock market liquidity. Dataset can be retrieved from https://instruction.bus.wisc.edu/jfrees/jfreesbooks/Regression20Modeling/BookWebDec2010/data.html (accessed on 8 January 2023). The competing models are the log-harmonic mixture Burr XII (LHMBXII) distribution [31] and the log-Gumbel Burr XII (LGBXII) distribution [38]. The response variable is the total number of shares that were traded on an exchange during a specific period (volume), while the covariate is the number of shares outstanding as of December 31, 1984, expressed in millions of shares (shares) . The fitted model is given by the following expression:

The response variable data are 16.221, 5.693, 11.965, 3.834, 13.235, 0.658, 13.794, 2.009, 27.600, 18.515, 23.466, 18.192, 28.163, 4.428, 4.912, 9.802, 15.513, 10.781, 11.149, 7.561, 14.969, 2.726, 17.275, 14.570, 11.967, 4.155, 3.072, 5.872, 17.744, 4.351, 16.995, 1.902, 19.727, 12.986, 2.823, 27.106, 46.832, 2.584, 18.056, 9.180, 7.999, 27.070, 8.914, 7.141, 28.574, 5.992, 2.118, 37.720, 13.080, 12.249, 19.874, 15.006, 4.503, 6.110, 12.238, 6.268, 5.810, 16.576, 33.565, 64.572, 3.000, 3.842, 16.242, 34.453, 21.718, 25.377, 8.317, 15.255, 4.971, 4.277, 2.473, 6.798, 10.848, 2.908, 3.695, 12.074, 37.923, 10.646, 19.624, 15.165, 4.054, 21.249, 8.726, 8.409, 34.073, 14.464, 8.903, 11.556, 5.067, 14.148, 4.712, 14.855, 6.180, 9.457, 6.170, 7.238, 9.877, 13.001, 12.967, 30.597, 15.134, 1.365, 39.273, 19.387, 15.505, 6.810, 18.001, 7.093, 3.791, 2.638, 18.960, 6.927, 19.436, 5.098, 16.257, 2.805, 10.961, 6.299, 5.971, 13.477, 11.716, 40.585, and 20.430.

The covariate’s data are 81.141, 27.088, 189.680, 13.492, 72.600, 6.736, 107.743, 20.851, 220.776, 79.964, 97.225, 151.508, 141.309, 21.433, 27.309, 50.338, 137.269, 73.200, 46.544, 29.718, 38.301, 19.683, 92.775, 72.127, 83.009, 16.049, 13.689, 15.161, 176.210, 13.524, 39.761, 8.724, 190.192, 20.232, 26.879, 155.600, 172.200, 12.402, 46.600, 78.946, 46.870, 126.268, 67.662, 43.206, 454.876, 50.366, 20.402, 315.451, 70.039, 35.408, 108.496, 34.787, 280.173, 49.342, 46.713, 40.222, 55.816, 96.849, 256.478, 612.686, 36.720, 26.756, 49.942, 139.747, 182.845, 125.229, 52.777, 130.825, 41.947, 23.261, 27.168, 29.495, 87.213, 16.209, 30.065, 116.298, 407.704, 72.103, 189.167, 78.281, 9.121, 99.637, 46.101, 62.933, 87.766, 115.561, 21.716, 19.896, 29.574, 95.962, 22.525, 93.908, 30.959, 16.399, 44.673, 86.376, 38.180, 103.353, 84.354, 361.610, 59.854, 11.217, 300.702, 271.429, 250.052, 53.415, 54.348, 32.016, 17.147, 19.276, 53.827, 24.616, 34.743, 30.653, 121.370, 13.610, 60.899, 25.436, 34.305, 139.916, 129.852, 783.051, and 95.704.

Table 5 presents the maximum likelihood estimates, standard errors (in parentheses), and values for the fitted model. The model selection criteria clearly indicate that the LMAPTBXII model is the most suitable. Leveraging the parameter estimates derived from the LMAPTBXII model, we formulate the following equation:

We can deduce that the effect shares had on the total number of shares that were traded on an exchange during a specific period was positively significant.

To assess the appropriateness of the LMAPTBXII, LGBXII, and LHMBXII models, Cox–Snell residuals were generated. Upon examination of the probability-probability (P-P) plot illustrated in Figure 10, it is evident that the residuals of the LMAPTBXII model exhibit closer alignment to the diagonal line in comparison to those of the LGBXII and LHMBXII models. This observation signifies that the LMAPTBXII model provides a better fit to the data.

Table 6 summarises the diagnostics results ( values in parentheses), which affirm the conclusions made using the P-P plots.

8. Conclusion

In this study, we thoroughly investigated the tractability, performance, and flexibility of a novel three-parameter modified alpha power transformed Burr XII distribution. One notable feature of this model is the discovery of a closed-form expression for its quantile function, adding to its analytical convenience. The parameter estimation for the proposed model was carried out through the rigorous maximum likelihood estimation method. This approach ensures that the model parameters are optimised to best capture the characteristics of the data under consideration. To assess the performance and flexibility of the new model, a comprehensive simulation study was conducted, and the model’s applicability was demonstrated through its application to two lifetime datasets. The robust results obtained from both the simulation study and the practical applications convincingly validate the enhanced flexibility of the proposed modified alpha power transformed Burr XII distribution. Building upon the original model, we introduced a logarithmic transformation to create a new log-modified alpha power transformed Burr XII model. Through thorough development and rigorous validation, we demonstrated the viability of this new model, offering an alternative perspective for researchers working with diverse data patterns. It is worth noting that the motivation behind introducing modifications to existing distributions, such as the Burr XII, is rooted in the need for increased flexibility. Recognising that no single distribution can perfectly fit all datasets, these modifications empower researchers to account for a broader range of data patterns. In this context, the modified alpha power transformed Burr XII stands out as a valuable addition, capable of accommodating decreasing, skewed, and near-symmetric datasets. In summary, our study makes a significant contribution to the continuous endeavour of improving the versatility of statistical distributions. As we continue to refine and expand our toolkit in statistical modelling, the MAPTBXII distribution stands out as a valuable addition, empowering researchers to address a wide range of data challenges with increased precision and adaptability. The modified alpha power transformed Burr XII distribution, along with its log-transformed counterpart, represents a valuable asset for researchers seeking a flexible and robust framework for data analysis and modelling.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest with regard to this article.