A novel extended inverse Weibull distribution: Statistical analysis and application

Abstract

This paper proposes a new type of exponential-type Weibull distribution based on the inverse Weibull distribution --- the transformed inverse Weibull distribution. This distribution constructs a more flexible parameter structure through mathematical transformation and has a better fitting effect on actual data. We deeply analyzed the key statistical properties of this distribution, including the probability density function, survival function, quantile function, as well as Shannon entropy, Rényi entropy, Tsallis entropy, and Mathai-Haubold entropy, etc. In terms of parameter estimation, various parameter estimation methods such as maximum likelihood estimation and Bayesian estimation were adopted to estimate the parameters of the transformed inverse Weibull distribution, and the performance of various parameter estimation methods was evaluated through Monte Carlo simulation. Finally, two sets of real data were applied to verify the applicability and effectiveness of the model in practical applications. The results show that the transformed inverse Weibull distribution exhibits a superior fitting performance in the goodness-of-fit test compared to the Weibull distribution, weighted exponential distribution, exponential Pareto distribution, flexible Weibull distribution, generalized exponential distribution, and generalized inverse exponential distribution.

1. Introduction

Probability distribution models hold a central position in the fields of statistics and probability theory. They not only provide a theoretical foundation for data modeling, prediction, parameter estimation, and statistical inference but also play a crucial role in describing random processes, constructing complex statistical models, solving optimization problems, and applications in machine learning and artificial intelligence. Although the research on probability distribution models in the existing literature has reached a relatively mature stage, offering robust theoretical support for addressing problems across various domains, traditional probability distribution models do not always achieve optimal fitting in the process of actual data fitting. In light of this, researchers continue to explore more flexible distribution models by extending and transforming classical models, aiming to enhance their fitting performance to meet the evolving needs of data analysis.

In recent years, some progress has been made in converting classical progressive Weibull distribution (WD) models to improve their flexibility. For instance, Alshanbari et al. [1] proposed a new flexible Weibull extension distribution based on the characteristics of extreme data, which can predict and model extreme observations more effectively. In cases where data sets exhibit mixed state faults, traditional probabilistic models are no longer suitable. Therefore, Khan et al. [2] introduced a beta power WD model that is sufficiently flexible to handle multiple failure modes. Al-Marzouki et al. [3] artificially enhanced the accuracy of modeling and prediction by extending the flexible WD (FWD) to develop the modified FWD, and demonstrated through real data verification that its fitting effectiveness surpassed that of the original FWD. Liu et al. [4] proposed a new power FWD by amalgamating the FWD with a novel power transformation method. Simulation studies have shown that this new distribution boasts higher flexibility and superior fitting capabilities compared to the flexible WD. Shi et al. [5] introduced the exponential flexible WD, which integrates the flexible Weibull expansion with the exponential T-X strategy, and it was shown to offer greater flexibility and improved fitting performance over the traditional WD. Gemeay et al. [6] extended the II Laplace semi-logarithmic distribution through power transformation techniques and proposed the power-type II Laplace semi-logarithmic distribution. This new distribution enhances the flexibility and applicability of the original distribution, enabling it to adapt to more complex real-world scenarios. Chaisee et al. [7] proposed a new Gamma-Exponential Weibull Poisson distribution, an extension of the exponential-Weibull Poisson distribution family. It combines the advantages of the exponential, Weibull, and exponential WDs, offering greater flexibility in fitting data and broader applicability, thereby enabling the analysis of more complex datasets. Tu et al. [8] introduced a Weighted Sine-generalized IWD by integrating the generalized IWD with a sine-generated probability framework. This distribution exhibits enhanced flexibility and provides a better fit compared to competing models. Zhu et al. [9] proposed a novel Sine FWD by integrating the traditional WD with a sine function. Compared to the standard Weibull, generalized Weibull, and other competing models, the Sine FWD exhibits superior fitting performance in reliability engineering and related applications. These studies have shown that by modifying the classical probability distribution model, new models can be obtained that are more in line with the characteristics of the actual data, leading to more accurate and reliable results in the field of statistical modeling and data analysis.

The inverse WD (IWD) is a probability distribution extended by WD. Compared with WD, IWD has higher flexibility and is especially suitable for describing the early failure of products. Its failure rate function can present various forms and adapt to different failure modes. This distribution can be applied to a wider set of data, especially atypical failure modes, and can effectively model the reliability of complex systems. In addition, parameter estimates for IWD may be more robust when dealing with extreme values or outliers. Therefore, it has been widely used in survival analysis, reliability analysis and life testing [10–13]. However, with the advancement of technology, we have noticed that IWD also has limitations in practical applications, especially when it comes to fitting real data, it cannot achieve the best fitting effect. Based on this, this study proposes a transformed exponential-type WD, which is a transformation based on IWD. Therefore, we name it the transformed IWD (TIWD). TIWD has a more flexible parameter structure and has a better fitting effect on real data, thereby improving the accuracy of statistical inference. In addition, the proposal of the TIWD model not only expands the existing model theory but also enriches the diversity of probability distributions in statistics. More importantly, the TIWD model also provides new ideas and methods for theoretical researchers, opening a new chapter for the research on probability distribution theory and its applications in various fields.

The rest of this article is as follows. In Section 2, we derived the mathematical expressions for the probability density function (PDF), cumulative distribution function (CDF), survival function (SF), and hazard function (HF) of the TIWD model, discussed the heavy tailed characteristics of the TIWD model, and presented relevant images of the model. In Section 3, we further analyzed the mathematical properties of the TIWD model, such as mixed representation, moments, and quantile functions. In Section 4, we analyzed various entropy measures under the TIWD model. In Section 5, we introduced several parameter estimation methods for the TIWD model, including maximum likelihood (ML) estimation, Bayesian estimation, Anderson Darling (AD) estimation, Cramer-von-Mises (CVM) estimation, and ordinary least squares (OLS) estimation. In Section 6, we compared the performance of ML estimation and three other estimation methods in parameter estimation, mean square error (MSE), and coefficient of variation (CV) through Monte Carlo simulations. In addition, this section also separately analyzed Bayesian estimation to explore the impact of different prior distributions on the estimation results. In Section 7, we applied the TIWD model to two sets of real data for analysis, in order to verify the feasibility of the model in practice. Section 8 provides relevant conclusions, limitations of the model, and future research directions.

2. TIWD model

Let X be a random variable that follows IWD, then the PDF and CDF of IWD are:

(1)(2)

Where is the shape parameter and is a known constant. As shown in Fig 1, when the parameter is fixed, the IWD images under different values show significant differences. For the PDF, when , the image is in an inverted bathtub shape, with a clear peak and a rapid decline after reaching the peak. When , the peak of the PDF significantly decreases, and the downward trend becomes slower. When , the PDF curve becomes more gentle. Thus, as the value increases, the PDF peak gradually decreases and shifts to the right, and the curve becomes smoother. For the CDF, all curves show a trend, and the larger the value, the more gradual the CDF increasing trend.

[Figure omitted. See PDF.]

In this paper, we assume that . Given the transformation , where represents the CDF of the probability distribution model, the PDF and CDF of the TIWD are shown as follows:

(3)(4)

When one of the parameters is fixed, the PDF image of TIWD shows different trends as the other parameter continues to change, as shown in Fig 2. From the image, it can be intuitively seen that the PDF of the TIWD model shows a decreasing trend with the change of the random variable X within different parameter ranges, and the decay rate gradually slows down. To prove whether the tail of the TIWD model has heavy tail properties, we need to further investigate. We give theorem 1 to prove the heavy tail property of TIWD model.

[Figure omitted. See PDF.]

Theorem 1. Let X be a random variable following TIWD, and the distribution function of TIWD is heavy tailed, if for any , the tail probability satisfies:

where , .

Proof. Since , , then there is . Since , therefore

Obviously, this is an infinitive of form , and , so we rewrite into the Equation (5):

(5)

The Equation (5) then becomes an undetermined form of , which we can use for the L’Hospital’s rule solution limitvalue. Derive the numerator and denominator in Equation (5) separately to obtain:

Then,

(6)

Continuing the L’Hospital’s rule limit on Equation (6), we obtain:

By doing Lopida over and over again, we get:

Therefore, by combining the intuitive judgment of PDF images with the quantitative analysis of the tail probability of CDF, we can comprehensively evaluate and confirm the heavy-tailed characteristic of the TIWD model. As shown in Fig 2, the right tail decay speed of the TIWD model is slower than that of the IWD model. Moreover, as the value increases, all curves show a trend of becoming more concentrated, which means they all have light-tailed characteristics. In other words, a high value leads to a shorter tail of the distribution, while a low value may result in a longer tail or a heavy-tailed phenomenon.

The SF and HF of TIWD are:

(7)(8)

Fig 3 shows the HF curves of the TIWD model under various parameter settings. When observing the image, it is evident that when one parameter remains constant, the HF image exhibits an inverted bathtub shape as the other parameter changes. This indicates that the stability and reliability of the product or system improve over time. In addition, it also indicates that compared to classical models, the model can be adapted to different risk control requirements by setting different parameters, enhancing its significant advantages in adaptability and flexibility. When further analyzing the HF curve, we observed that as the independent variable x increases, the curve as a whole tends to approach zero. This phenomenon indicates that the HF curve displays distinct heavy tailed features. In other words, the curve shows a slower decay rate in areas with higher X values, reflecting a thicker tail than the light tail distribution. The heavy tail characteristic can effectively avoid risk assessment under light tail distribution, thereby ensuring the predictability of risks in high reliability scenarios.

[Figure omitted. See PDF.]

The TIWD shares certain similarities with the inverse Chen distribution. Both are derived through mathematical transformations of their original distributions, and their PDF and HF curves exhibit similar morphological characteristics, demonstrating heavy-tailed features within specific parameter ranges. However, their key difference lies in their transformation mechanisms. The TIWD is obtained by introducing additional parameters for mathematical transformation while keeping one parameter of the original distribution fixed [14], whereas the inverse Chen distribution relies solely on the parameters of the original distribution for its transformation.

Cumulative HF and inverse HF are two important concepts in reliability analysis, closely related to HF. The cumulative HF represents the integral of HF from the initial moment to time t, which is the cumulative amount of failure risk within that time interval. At the same time, the cumulative HF is closely related to SF, which represents the probability that the product has not failed before time t. The reverse HF is used to analyze the probability of failure after time t, and can be used to evaluate the risk of product failure during the remaining life. Here we give the expressions of these two functions [15]:

3. Mathematical properties

3.1. Moments

The r-th moment of the TIWD model can be computed utilizing the subsequent mathematical expression:

(9)

According to Equation (3), when , the exponential function can be expanded by Taylor to obtain . Therefore

then X is said to have a power-law distribution. For the power-law distribution , the existence of its r-order moment is related to , that is, when , there is

According to the necessary and sufficient condition for integral convergence, there is

Therefore, the r-order moment is finite if and only if .

By substituting Equation (3) into Equation (9), we get:

(10)

Let , then Equation (10) is simplified as:

(11)

Therefore, from Equation (11), we can see that the mean and variance of X are:

(12)(13)

3.2. Incomplete moments

Incomplete moments, as an important statistical measure for describing the partial order moments of random variables, play a crucial role in fields such as risk management, financial mathematics, and extreme value theory, especially in the assessment of extreme event risks [16]. In the TIWD model, indepth analysis of the characteristics of incomplete moments will provide us with valuable insights into the distribution characteristics of random variables and expand their potential applications in multiple fields. The definition of incomplete r-th moment for continuous random variable X and its PDF is as follows [17]:

(14)

Therefore, substitute Equation (3) into Equation (14) to obtain the r-th incomplete moment of the random variable X:

(15)

The Lorenz curve is a tool in economics used to describe the degree of inequality in income or wealth distribution. It visualizes the uniformity of income and wealth distribution through images. In statistics, incomplete moments are an important tool used to describe the characteristics of income and wealth distribution. Although the two differ in terms of imagery and mathematical expressions, they both analyze economic distribution characteristics from different perspectives. The corresponding expression is given by the formula for the Lorenz curve below [18]:

(16)

The Bonferroni curve, as a statistical method, is primarily used for multiple hypothesis testing and is defined as the ratio of the Lorenz curve to the CDF [19]:

(17)

3.3 Quantile function

Quantile functions, serving as a pivotal instrument in the realms of statistics and probability theory, facilitate the comprehension and analytical dissection of data attributes. Furthermore, they are instrumental in diverse statistical inferential and decision-making frameworks. These functions can be derived through the process of inverse transformation applied to the CDF:

Among them, w represents the probability value, ranging from 0 to 1. When w equals 0.25, 0.5, and 0.75, they correspond to the first quartile, median and the third quartile, respectively. In Fig 4, we present 3D plots of the first quartile, median and the third quartile. From these plots, we can observe that as the probability values change, the degree of skewness in the distribution gradually becomes smoother. Skewness and kurtosis are important statistical measures used to characterize the state of data distribution. In Fig 5, we present plots of the skewness and kurtosis of the TIWD model across different parameter ranges. We will provide the expressions for skewness and kurtosis:

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

3.4 Order statistics

Order statistics are a powerful and versatile tool in statistics, with extensive applications in data analysis, optimization problems, decision-making, and theoretical research. Let be random samples from the total sample X, where are observations of the random samples. Arrange samples in ascending order to obtain , then it is called is the order statistic of , where is the minimum order statistic of the sample and is the maximum order statistic of the sample. From Equations (3) and (4), we can obtain the distribution density of the j-th order statistic of TIWD, which is:

(18)

So the minimum order statistic is:

The maximum order statistic is:

3.5 Mean-residual life

The mean-residual life is an important concept in reliability analysis and product life testing. It refers to the average time the system can continue to operate correctly after a specific point in time t, which can be calculated using the probability density of the remaining life. The corresponding calculation formula is as follows [20]:

Among them, is the SF, which can be known from Equation (7).

4. Entropy measure under TIWD model

4.1. Shannon entropy

Shannon entropy, a core concept in information theory, quantifies the degree of uncertainty in information and plays a crucial role in various fields, including cryptography, information coding, and communication systems. In statistics, Shannon entropy evaluates probability models, providing rigorous theoretical support for data analysis and model construction, thereby enhancing the depth and accuracy of data analysis. The expression of Shannon entropy within the TIWD model is delineated as follows [21]:

(19)

4.2. Rényi entropy

As an extension of Shannon entropy, Rényi entropy shows its strong flexibility through the introduction of parameters , so it has a wide range of applications in machine learning, image processing, bioinformatics and other fields [22–25]. The expression of Rényi entropy within the TIWD model is delineated as follows:

(20)

4.3. Tsallis entropy

Tsallis entropy, conceptualized as a generalization of entropy, was introduced by Constantino Tsallis in 1988. Through the incorporation of the entropy parameter , Tsallis entropy offers a more robust theoretical framework and pragmatic approach for scaling analysis, optimal decision-making, and machine learning, thereby enhancing the modeling of real-world complexities and diversities. The expression of Tsallis entropy within the TIWD model is delineated as follows [21]:

(21)

4.4. Mathai–Haubold entropy

The Mathai-Haubold entropy, advanced within the domain of statistical mechanics, offers a more nuanced representation of certain system configurations. Its distinct mathematical attributes and utility in the context of complex systems render the Mathai-Haubold entropy a pivotal construct in the disciplines of mathematics and physics. The expression of Mathai–Haubold entropy within the TIWD model is delineated as follows [26]:

(22)

5. Parameter estimation of TIWD model

5.1. ML estimation

Let is the n sample size from the TIWD model, denoted , and the likelihood function (LF) can be defined as follows:

(23)

The log-LF is then:

(24)

Subsequently, we derive the parameters in Equation (24) individually to elucidate the likelihood Equations (25) and (26):

(25)(26)

ML estimates for the parameters and can be derived by resolving the Equations (25) and (26). However, direct resolution of these equations is computationally complex, necessitating the use of numerical computation to ascertain the parameters within them. In this study, we employ the dichotomy method to numerically calculate the parameters and obtain the ML estimates and for the parameters and . To ensure the reliability and effectiveness of the ML estimation, we need to demonstrate the existence and uniqueness of the ML estimate. This process not only ensures the accuracy of statistical inference, but also enhances the application value of the ML estimation in practical data analysis and model parameter estimation.

5.1.1 Existence and uniqueness of ML estimation solutions.

Theorem 2.Let the left side of Equation (25)be, then the solution obtained by exists and is unique on .

Proof. (1) Existence

When , it is obvious that , so

When , we have , so

Due to , , , then , , then .

Since is continuous on , the intermediate value theorem shows that has at least one solution on such that , proving existence.

(2) Uniqueness

Taking the derivative of in gives the following formula:

.

It is obvious that , so we can see that is a monotone decreasing function. If a strictly monotonically decreasing function intersects with the x-axis, there is at most one intersection point. Combining ‘existence’, it can be inferred that there exists at least one solution, thus proving uniqueness.

Theorem 3. Let the left side of Equation(26)be

then the solution obtained by exists and is unique on .

Proof. (1) Existence

When , it is obvious that

, , ,

So .

When , we have .

Also because and are related to the values of x, the following discussion now follows:

When , , , , ,

so .

When , , , , ,

so .

When , , , , ,

so we need to compare the growth rates of and .

Let , , so we have

, .

Thus . Because the growth rate of is faster than any other term and , . Therefore, there exists at least one with .

In summary, when , there exists at least one solution with for , and for all . Therefore, according to the intermediate value theorem, it can be concluded that has at least one solution on such that , proving existence.

(2) Uniqueness

The derivative of in yields the following equation:

(27)

Given that , to verify that , it is necessary to compare the magnitudes of the second and third terms. Given that and , it follows that

If , then

If , the growth rate of is faster than the reduction rate of . At this time, the inequality

still holds. So we have

Which implies

Thus,

Since , it follows that is a monotonically decreasing function on . Combining with the “existence” proof, we know that there is at least one solution on , so we have proved the uniqueness. ☐

5.1.2 Asymptotic confidence interval.

In this section, we construct asymptotic confidence intervals (ACIs) for the parameters to evaluate the performance of ML estimates in terms of precision, stability, and applicability. Since ML estimates are asymptotically normal, we can use and to construct ACIs for and . The asymptotic variance needs to be obtained from the inverse of the Fisher information matrix, which is the negative value of the Hessian matrix of the log-LF, i.e.,

Where

We substitute the obtained ML estimates and into the Fisher information moment to obtain:

Then the covariance matrix of and is:

The ACIs for and is:

, .

Where is the percentile of the standard normal distribution.

5.2 Bayesian estimation

Bayesian estimation is a statistical inference method based on Bayes’ theorem, which is based on the principle of combining a priori knowledge with observed data to obtain estimates of parameters [27]. In this process, prior knowledge plays an important role in Bayesian estimation, which not only provides the basis for data analysis, but also determines the posterior distribution of the parameters in combination with the observed data. Therefore, choosing a reasonable prior distribution is decisive for ensuring the accuracy, reliability and stability of parameter estimation. As the conjugate prior of Poisson and exponential distributions, the gamma distribution can flexibly express different prior information by adjusting the parameter values, thus simplifying the computational process of Bayesian estimation and improving the accuracy and reliability of the estimation results. Therefore, in this paper, we choose the gamma distribution as the prior distribution of and . Assuming that and are independent random variables and follow and , respectively, the joint prior distribution of and is:

(28)

The joint posterior density of and is:

(29)

To evaluate the influence of prior information on the estimation results, we also consider non-informative priors in this study. By comparing the results from these two types of priors, we can verify the stability of the TIWD model. When the prior distribution of and is non-informative prior, then the joint prior distribution of and is:

(30)

The joint posterior density of and is:

(31)

In Bayesian estimation, we usually introduce a loss function to quantify the cost loss or decision errors that occur during the estimation process, thereby transforming statistical inference into an optimization problem. In this section, we introduce the precautionary loss function (PLF) to estimate the parameters of the TIWD model. Next, we will analyze Bayesian estimation under PLF. According to Akhtar [28], the definition of PLF is:

Thus, the Bayesian estimator under PLF is:

(32)

Among them, , , and represent functions about and . From Equation (32), it can be seen that the Bayesian estimator of TIWD under PLF is non explicit, and it is relatively complex to directly calculate Equation (32). Therefore, we consider using MCMC sampling to obtain the Bayesian estimates under PLF.

MCMC sampling, a foundational statistical technique, lies at the heart of its ability to efficiently draw samples from complex probability distributions. This is achieved through the integration of Markov chains and Monte Carlo methods. In this approach, a sample from the posterior distribution is obtained by constructing a Markov chain with the target distribution as its stationary state, upon which parameters of the model are estimated and inferred. Presently, popular MCMC algorithms include the Metropolis-Hastings (MHs) algorithm, Gibbs sampling, and slice sampling, among others. In practice, the selection of an appropriate MCMC algorithm often depends on the nature of the target distribution and the dimensionality of the probability space. For instance, the MH algorithm is typically well-suited for low-dimensional probability distributions, whereas Gibbs sampling is favored for high-dimensional scenarios. In this study, we couple the MH algorithm with Gibbs sampling to harness the synergistic benefits of these two methodologies, thereby reducing sampling complexity.

According to Equation (29), the posterior distributions of and are as follows:

(33)(34)

Equations (33) and (34) demonstrate that the posterior distribution of and is implicit, necessitating the application of mixed Gibbs sampling for the extraction of samples from this distribution. The iterative steps for implementing the mixed Gibbs sampling algorithm are delineated as follows:

1. (1) : Set the initial values of the parameters , .

2. (2) : Sample from the distribution using the MH algorithm. Choose the normal distribution as the proposal distribution for , where represents the current state of and represents its variance.

1. (i) A sample is selected from the and resampled when . The acceptance probability is then calculated based on :

2. (ii) Sample is taken from the uniform distribution U (0,1), so that:

3. (iii) Let j = j + 1, return (i) to continue the above step.

4. (3): Sample from the distribution using the MH algorithm. Choose the normal distribution as the proposal distribution for , where represents the current state of and represents its variance.

1. (i) A sample is selected from the and resampled when . The acceptance probability is then calculated based on :

2. (ii) Sample is taken from the uniform distribution U (0,1), so that:

3. (iii) Let j = j + 1, return (i) to continue the above step.

Samples and of and are obtained through the above sampling steps, so the parameters of the TIWD model under PLF are Bayes estimator as:

(35)

In the above situation, we consider sampling under the condition of informative prior. When considering the non-informative prior, we only need to change the posterior distribution of and under the condition of informative prior to the posterior distribution under the condition of non-informative prior. It should be noted that during the iteration process of the Markov chain before reaching the target distribution, the sampling samples may be affected by the initial values. In order to ensure that the sampling samples used in the analysis can reliably represent the probability distribution being studied, we need to discard some samples and use to represent them.

5.3. AD estimation

AD estimation is a statistical test usually used to determine whether a set of data conforms to a particular probability distribution model, based on the difference between the cumulative and empirical distribution functions of the data. This test method was proposed by Anderson and Darling [29] in 1952. Compared with the traditional goodness of fit test, AD test has some unique advantages. It analyzes the cumulative function of the distribution, thus avoiding the direct analysis of the PDF. This makes the AD test more flexible in practical applications, adapting to a wide variety of data types and distributions. In the existing literature, many researchers used the AD estimation method to test the goodness of fit, so as to understand and analyze the data. For example, Alsadat et al [30] used AD estimation to estimate the reliability of unit Gompertz distribution. Shafiq et al. [31] used AD test to measure the parametric efficiency of the semilog unit Gamboz type I distribution. Aboraya et al. [32] used AD estimation to estimate small samples in the new compound Lomax model. In this paper, the AD estimates for and can be obtained by minimizing the and functions: In this paper, the AD estimates for and can be obtained by minimizing the and functions:

(36)

Among them, represents the parameters and .

In addition, we can also obtain and by solving the following nonlinear equations:

(37)(38)

Here

(39)(40)

5.4. CVM estimation

The CVM estimation is an important alternative to the Kolmogorov-Smirnov test, proposed by Harald Cramer and Richard von Mises, and it holds a significant place in statistics. This test offers a flexible and effective statistical method to verify whether a sample set conforms to a specific distribution by quantifying the discrepancy between the CDF of the sample and the CDF of the hypothetical distribution, thus providing an indepth comparison and evaluation of the hypothesis that the sample data set follows a particular distribution. As the CVM test is a non-parametric test, it possesses wide applicability and flexibility for small sample sizes, rendering it extensively used across various fields. In current literature, the CVM test is widely applied to validate various statistical models and data distributions. For instance, Wang and Zhu [33] employed CVM estimation to test the fitness of the lognormal distribution for power functions, corroborating the method with real-world data. Kutzker et al. [34] utilized the CVM test to assess the induction class of conditional distribution functions, while Cavaliere et al. [35] utilized the CVM test statistic to verify the correct specification of conditional variance functions within the GARCH model. Anis et al. [36] considered the application of the CVM test for parameter estimation in the Rayleigh distribution. These studies underscore the practical significance and broad applicability of the CVM test across numerous fields of research. In this paper, the CVM estimates for and can be obtained by minimizing the and functions:

(41)

In addition, we can also obtain and by solving the following nonlinear equations:

(42)(43)

Among them, , j = 1, 2 can be obtained from Equations (39) and (40).

5.5. OLS estimation

OLS is an estimation method used to estimate unknown parameters in linear regression models. It does this by minimizing the sum of squares of the residuals observed between the dependent and predictor variables in a given data set. Owing to its simple principle, strong optimality, and high efficiency, the OLS method has become an important tool for data analysis and is widely used in economics, finance, engineering, and other fields [37–40]. In this paper, the OLS estimates for and can be obtained by minimizing the and functions:

(44)

In addition, we can also obtain and by solving the following nonlinear equations:

(45)(46)

Among them, , j = 1, 2 can be obtained from Equations (39) and (40).

Due to the complex forms of Equations (37), (38), (42), (43), (45), and (46), it is difficult to solve them directly. Therefore, we use numerical methods to solve these equations. As an efficient numerical analysis method, Newton Raphson method has significant advantages in solving root problems of nonlinear equations. Therefore, in this paper, we use Newton Raphson method to numerically solve these complex equations in order to obtain accurate solutions.

6. Simulation

To evaluate the uncertainty and complexity of the TIWD model, enhance the accuracy of statistical inference, and improve the efficiency of data processing, Matlab software is employed to compute various entropy measures under the TIWD model. We set the parameters , and entropy parameters , with their specific values presented in Tables 1 and 2, respectively.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Then, this paper deeply analyzes the ML estimation, Bayesian estimation and other three estimation methods in TIWD model parameter estimation through simulation. Specifically, we first set the initial values of different parameters, and respectively, and selected the sample size n = 20, 30, 50, 80, 100 for 1000 cycles of simulation. In order to obtain ML estimation and other three estimation methods, the parameter average estimates (AEs) and the corresponding MSEs are obtained. In order to further evaluate the accuracy of each parameter estimation method, we also calculate the CVs of the parameters under each method. The CV is an important statistical indicator used to measure the relative dispersion of data, which eliminates the influence of different units or means on the dispersion by calculating the ratio of standard deviation to mean. Detailed results are shown in Tables 3 and 4. In the discussion of Bayesian estimation, since its performance is affected by the selection of prior distributions, we discuss it as a separate part. In this part, we consider two cases of informative prior and non-informative prior, and set the initial parameter values as and , and keep the sample size unchanged. Through 1000 cycle simulation, we adopted the hyperparameter under the informative prior, and applied MCMC algorithm to carry out 5000 iterations, among which the first 500 iterations were used as the combustion period to eliminate the influence of the initial value, thus obtaining the AEs and MSEs of parameter Bayesian estimation, refer to Tables 5 and 6. It should be noted that for the sake of simplicity, we refer to INF and Non-INF as Bayesian estimates under informative prior and non-informative prior, respectively. To further analyze the performance of the estimation methods used, we set and , construct ACIs for ML estimates at different confidence levels of , and calculate the coverage probability (CP) and average width (AW) of the parametric confidence intervals. The specific results are shown in Table 7.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

According to the above tables, the following conclusions can be drawn:

1. (1) When and are fixed, the measures of entropy show an upward trend with the increase of . When and are fixed, the entropy measures tend to decrease with the increase of .

2. (2) With the increase of sample size, both MSEs and CVs of TIWD model parameter estimates showed a decreasing trend, indicating that the accuracy of estimation increased with the increase of sample size.

3. (3) Among various parameter estimation methods, ML estimation shows superior performance in overall accuracy compared with the other three estimation methods. At the same time, the OLS estimation is generally superior to other estimation methods in terms of precision.

4. (4) In Bayesian estimation, Bayesian estimation based on informative prior is superior to Bayesian estimation based on non-informative prior in accuracy.

5. (5) With the increase of confidence level, CP and AW of confidence interval both show an increasing trend.

7. Real data analysis

In this section, we evaluate the applicability of the TIWD model in practical applications by empirically analyzing two sets of real data. The first set of data represents the fatigue life of metal components (in cycles) [41]. The second set of data reflects the survival time of patients with head and neck cancer [42]. The specific details of these two sets of data are detailed in Table 8.

[Figure omitted. See PDF.]

Before conducting in-depth analysis on the two sets of real data we provided, it is necessary to verify whether the real data conforms to the characteristics of algebraic decay. Whether the real data satisfies algebraic decay directly affects the accuracy and reliability of the TIWD model in statistical inference, which is the key to our subsequent research. Through Equation (7), we provide the following analysis:

When , there is , and is expanded by Taylor expansion to obtain:

.

So there is . That is to say, when , there is . It indicates that exhibits algebraic decay characteristics. Next, we will further verify that the real data satisfies the algebraic decay characteristic by plotting the log-log SF graph. Please refer to Fig 6. From Fig 6, it can be observed that the tail data is closer to the theoretical asymptotic line. The blue data points gradually align with the red dotted line, indicating that the tail conforms to the asymptotic characteristic of algebraic decay, which is consistent with the thick-tailed characteristic of TIWD.

[Figure omitted. See PDF.]

In order to further evaluate the performance of the TIWD model in fitting two sets of real data, we used Kolmogorov Smirnov (KS) test, Akaike information criterion (AIC), bias correction AIC (AICc) [43], and Bayesian information criterion (BIC) [44] to quantitatively evaluate the model. The calculation formulas for AIC, AICc, and BIC are as follows:47

Here, represents the log-likelihood function of the parameters given the data, k indicates the number of parameters in the model, and n represents the sample size. During this process, we use the p-value of the KS test as the key indicator to evaluate the goodness of fit of the model. To further highlight the outstanding performance of the TIWD model, we compared it with the WD, weighted exponential distribution (WED) [45], exponentiated Pareto distribution (EPD) [46], FWD [47], generalized exponential distribution (GED) [48], and generalized inverse exponential distribution (GIED) [49]. The specific test results are detailed in Table 9. To visually demonstrate the advantages of the TIWD model in fitting the two sets of data, we plotted the empirical distribution function of the real data and the CDFs of each model, as shown in Fig 7. Additionally, we created PP plots and QQ plots of the TIWD model for the two sets of real data to more clearly reveal the fitting effect. The relevant graphical displays are presented in Fig 8.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Through these specialized analysis charts, we can clearly observe that TIWD has a good fitting effect, which enables us to more accurately evaluate the performance of the TIWD model in practical applications. We used all the estimation methods in this study to obtain the parameter values and corresponding KS test values and p-values for these two sets of real data, as shown in Table 10. According to Table 10, the ML estimation performance is the best in real Data I, and the Bayesian estimation performance under informative prior is the best in real Data II.

[Figure omitted. See PDF.]

In addition to the two sets of real data mentioned above, the durability test data of deep groove ball bearings under working conditions and the dataset of infection occurrence times of hemodialysis patients within several months can also be considered for analysis. These two sets of data have good fitting characteristics with the TIWD model, so similar research can be conducted on them to further expand the application scope of the TIWD model.

8. Conclusions

In this study, we conducted a series of transformations on the IWD, proposed a novel TIWD, and conducted in-depth analysis of the statistical and mathematical properties of this distribution model. Firstly, we discussed the mathematical properties and statistical characteristics of TIWD. Subsequently, we used five estimation methods to estimate the parameters of the model and conducted a comprehensive evaluation of these estimation methods through simulation. The research results indicate that Bayesian estimation performs outstandingly in parameter estimation accuracy with the introduction of prior information, effectively reducing the impact of prior on the results. To verify the practical application performance of the TIWD model, we selected two sets of real data for empirical testing. Through goodness of fit testing, we observed that the TIWD model performed better in terms of data fitting accuracy and adaptability compared to other distribution models. Previous studies have shown that traditional IWD models have certain limitations in processing complex data in certain situations, and cannot accurately describe the distribution characteristics of specific data. The TIWD model obtained through mathematical transformation introduces new parameters, improves the flexibility and adaptability of the model, and can effectively improve the fitting effect of the model on data. From a theoretical perspective, the TIWD model is a further extension and generalization of the IWD model. It not only enriches the probability distribution theory system and improves the theoretical tools for describing complex data distributions, but also provides a deeper understanding of the properties and characteristics of the IWD model and its extended models. At the same time, the proposal of the TIWD model has strengthened the connection between statistics and other disciplines, promoted the development of interdisciplinary research, and provided new methodological support for solving complex problems in practice.

Although the TIWD model has certain advantages in expanding the applicability and fitting accuracy of traditional IWD models, it still has certain limitations. Firstly, due to the introduction of new parameters in the TIWD model, parameter estimation is more complex compared to traditional IWD models. When analyzing ML estimation, there may be difficulties in convergence, often requiring consideration of more numerical methods for computation. Secondly, the TIWD model requires specific data for calibration and validation when fitting real data, otherwise it will affect the performance of the model. Therefore, collecting and organizing data that meets the requirements of the model may require a significant amount of time and manpower, and may face difficulties in obtaining data in practical applications. Finally, due to the fixed parameter structure of the TIWD model, it may exhibit poor adaptability to new data in some practical application scenarios, making it difficult to flexibly adapt to different conditions and requirements. Continuous adjustments are needed to adapt to new scenarios.

Given the limitations of the TIWD model, future research needs to develop more efficient numerical calculation methods to obtain parameter estimates. Before performing data fitting, it is necessary to preprocess the data to remove outliers and improve data quality. At the same time, clarify the data standards that meet the requirements of the TIWD model, providing guidance for collecting and organizing data to reduce time and manpower. Finally, we should actively explore the research of TIWD models in emerging fields such as artificial intelligence and big data analysis. By expanding the research field and discovering the advantages and disadvantages of TIWD in different fields, the model can be improved and perfected, providing new ideas and methods for data analysis in various fields.

Supporting information

S1 File. Support information explanation.

https://doi.org/10.1371/journal.pone.0335555.s001

(TIF)

References

1. 1. Alshanbari HM, Odhah OH, Al-Mofleh H, Ahmad Z, Khosa SK, El-Bagoury AA-AH. A new flexible Weibull extension model: different estimation methods and modeling an extreme value data. Heliyon. 2023;9(11):e21704. pmid:38027837

* View Article

* PubMed/NCBI

* Google Scholar

2. 2. Khan F, Ahmad Z, Khosa SK, Alomair MA, Alomair AM, Alsharidi AK. A new modification of the flexible Weibull distribution based on power transformation: Monte Carlo simulation and applications. Heliyon. 2023;9(6):e17238. pmid:37426796

* View Article

* PubMed/NCBI

* Google Scholar

3. 3. Al-Marzouki S, Alrashidi A, Chesneau C, Elgarhy M, Khashab RH, Nasiru S. On improved fitting using a new probability distribution and artificial neural network: application. AIP Advances. 2023;13(11).

* View Article

* Google Scholar

4. 4. Liu X, Ji J, Alrashidi A, Almulhim FA, Alshawarbeh E, Seong J-T. A new probabilistic model with mixed-state failure rates: modeling time-to-event scenarios in reliability and music engineering. Alexandria Eng J. 2024;96:99–111.

* View Article

* Google Scholar

5. 5. Shi X, Hu J, Gao R. Integrated modeling of the sports and reliability data: implications of the probabilistic model and deep learning approaches. Alexandria Eng J. 2024;94:274–86.

* View Article

* Google Scholar

6. 6. Gemeay AM, Bashiru SO, Sapkota LP, Kayid M, Dutta S, Mohammad S. A new power transformed distribution with applications to radiotherapy and environmental datasets. J Rad Res Appl Sci. 2025;18(2):101339.

* View Article

* Google Scholar

7. 7. Chaisee K, Khamkong M, Paksaranuwat P. A new extension of the exponentiated Weibull–Poisson family using the gamma-exponentiated weibull distribution: development and applications. Symmetry. 2024;16(7):780.

* View Article

* Google Scholar

8. 8. Tu X, Kong J, Fu Q, Chang S, Zhang K, Alballa T, et al. A new trigonometric-inspired probability distribution: a simulation study and applications in reliability and hydrology. Alexandria Eng J. 2025;113:181–94.

* View Article

* Google Scholar

9. 9. Zhu JQ, El-Din MMM, Seong JT, Tashkandy YA, Bakr ME, Kumar A. A new sine-arisen probabilistic model and artificial neural network methods for statistical modeling of the music engineering and reliability data. Alex Eng J. 2024;101:108–17.

* View Article

* Google Scholar

10. 10. Almazah MMA, Ullah K, Hussam E, Hossain MdM, Aldallal R, Riad FH. New statistical approaches for modeling the COVID‐19 data set: a case study in the medical sector. Complexity. 2022;2022(1).

* View Article

* Google Scholar

11. 11. AbaOud M, Almuqrin MA. The weighted inverse Weibull distribution: Heavy-tailed characteristics, Monte Carlo simulation with medical application. Alexandria Eng J. 2024;102:99–107.

* View Article

* Google Scholar

12. 12. Afify AZ, Shawky AI, Nassar M. A new inverse Weibull distribution: properties, classical and Bayesian estimation with applications. KJS publishes peer-review articles in Mathematics, Computer Science, Physics, Statistics, Biology, Chemistry, and Earth and Environmental Sciences. 2021;48(3).

* View Article

* Google Scholar

13. 13. Hassan MKH, Aslam M. DUS-neutrosophic multivariate inverse Weibull distribution: properties and applications. Complex Intell Syst. 2023;9(5):5679–91.

* View Article

* Google Scholar

14. 14. Shakhatreh MK, Lemonte AJ, Cordeiro GM. On the generalized extended exponential-Weibull distribution: properties and different methods of estimation. Inter J Comp Math. 2019;97(5):1029–57.

* View Article

* Google Scholar

15. 15. Aslam M, Tahir M, Hussain Z. Reliability analysis of 3-component mixture of distributions. Scientia Iranica. 2017;0(0):0–0.

* View Article

* Google Scholar

16. 16. Hassan AS, Alsadat N, Chesneau C, Elgarhy M, Kayid M, Nasiru S, et al. Inverse power XLindley distribution with statistical inference and applications to engineering data. Sci Rep. 2025;15(1):4385. pmid:39910109

* View Article

* PubMed/NCBI

* Google Scholar

17. 17. Ogunde AA, Dutta S, Almetawally EM. Half logistic generalized rayleigh distribution for modeling hydrological data. Ann Data Sci. 2024;12(2):667–94.

* View Article

* Google Scholar

18. 18. Gemeay AM, Karakaya K, Bakr ME, Balogun OS, Atchadé MN, Hussam E. Power Lambert uniform distribution: statistical properties, actuarial measures, regression analysis, and applications. AIP Adv. 2023;13(9).

* View Article

* Google Scholar

19. 19. Irshad MR, Aswathy S, Maya R, Al-Omari AI, Alomani G. A flexible model for bounded data with bathtub shaped hazard rate function and applications. MATH. 2024;9(9):24810–31.

* View Article

* Google Scholar

20. 20. Kayid M. Further results on the proportional vitalities model. Entropy (Basel). 2021;23(9):1201. pmid:34573826

* View Article

* PubMed/NCBI

* Google Scholar

21. 21. Estañón CR, Montgomery HE Jr, Angulo JC, Aquino N. The confined helium atom: an information–theoretic approach. Int J Quant Chem. 2024;124(4).

* View Article

* Google Scholar

22. 22. Shrahili M. Some new results involving residual Renyi’s information measure for $ k $-record values. MATH. 2024;9(5):13313–35.

* View Article

* Google Scholar

23. 23. Chennaf S, Ben Amor J. Renyi entropy of uncertain random variables and its application to portfolio selection. Soft Comput. 2023;27(16):11569–85.

* View Article

* Google Scholar

24. 24. Bialas P, Burda Z, Johnston DA. Rényi entropy of zeta urns. Phys Rev E. 2023;108(6–1):064108. pmid:38243505

* View Article

* PubMed/NCBI

* Google Scholar

25. 25. Fazlollahi HR. Rényi entropy correction to expanding universe. Eur Phys J C. 2023;83(1).

* View Article

* Google Scholar

26. 26. El-Raouf MMA, AbaOud M. A novel extension of generalized Rayleigh model with engineering applications. Alexandria Eng J. 2023;73:269–83.

* View Article

* Google Scholar

27. 27. Dutta S, Alqifari HN, Almohaimeed A. Bayesian and non-bayesian inference for logistic-exponential distribution using improved adaptive type-II progressively censored data. PLoS One. 2024;19(5):e0298638. pmid:38753595

* View Article

* PubMed/NCBI

* Google Scholar

28. 28. Akhtar N, Khan SA, Ismail EAA, Awwad FA, Gul T, Khan AA. Bayesian estimation of 3-component mixture geometric distribution under different loss functions. AIP Adv. 2024;14(1).

* View Article

* Google Scholar

29. 29. Anderson TW, Darling DA. Asymptotic theory of certain “Goodness of Fit” criteria based on stochastic processes. Ann Math Statist. 1952;23(2):193–212.

* View Article

* Google Scholar

30. 30. Alsadat N, Hassan AS, Elgarhy M, Chesneau C, Mohamed RE. An efficient stress–strength reliability estimate of the unit gompertz distribution using ranked set sampling. Symmetry. 2023;15(5):1121.

* View Article

* Google Scholar

31. 31. Shafiq A, Sindhu TN, Dey S, Lone SA, Abushal TA. Statistical features and estimation methods for half-logistic unit-gompertz type-I model. Mathematics. 2023;11(4):1007.

* View Article

* Google Scholar

32. 32. Aboraya M, Ali MM, Yousof HM, Ibrahim M. A novel lomax extension with statistical properties, copulas, different estimation methods and applications. Bull Malays Math Sci Soc. 2022;45(S1):85–120.

* View Article

* Google Scholar

33. 33. Wang C, Zhu H. Tests of fit for the power function lognormal distribution. PLoS One. 2024;19(2):e0298309. pmid:38386634

* View Article

* PubMed/NCBI

* Google Scholar

34. 34. Kutzker T, Klein N, Wied D. Flexible specification testing in quantile regression models. Scandinavian J Stat. 2023;51(1):355–83.

* View Article

* Google Scholar

35. 35. Cavaliere G, Perera I, Rahbek A. Specification tests for GARCH processes with nuisance parameters on the boundary. J Busin Econ Stat. 2023;42(1):197–214.

* View Article

* Google Scholar

36. 36. Anis MZ, Okorie IE, Ahsanullah M. A review of the Rayleigh distribution: properties, estimation and application to COVID-19 data. Bull Malays Math Sci Soc. 2023;47(1).

* View Article

* Google Scholar

37. 37. Lim JS, Kim YD, Woo J-C. Approximated uncertainty propagation of correlated independent variables using the ordinary least squares estimator. Molecules. 2024;29(6):1248. pmid:38542885

* View Article

* PubMed/NCBI

* Google Scholar

38. 38. Mardones C, Orellana A. Estimating production functions to determine the marginal economic value and price elasticity of water in the Chilean industrial sector. Utilities Policy. 2023;85:101677.

* View Article

* Google Scholar

39. 39. Ahsan-ul-Haq M, Farooq MU, Nagy M, Mansi AH, Habineza A, Marzouk W. A new flexible distribution: statistical inference with application. AIP Adv. 2024;14(3).

* View Article

* Google Scholar

40. 40. Dong Y, Soale A-N, Power MD. A selective review of sufficient dimension reduction for multivariate response regression. J Stat Plan Inference. 2023;226:63–70.

* View Article

* Google Scholar

41. 41. Jia J-M, Yan Z-Z, Peng X-Y. A new discrete extended Weibull distribution. IEEE Access. 2019;7:175474–86.

* View Article

* Google Scholar

42. 42. Unal C, Cakmakyapan S, Ozel G. Alpha power inverted exponential distribution: properties and application. Gazi Univ J Sci. 2018;31:954–65.

* View Article

* Google Scholar

43. 43. Kudo S, Fujimoto M, Sato T, Nagano A. Determination of the optimal number of linked rigid-bodies of the trunk during walking and running based on Akaike’s information criterion. Gait Posture. 2020;77:264–8. pmid:32087596

* View Article

* PubMed/NCBI

* Google Scholar

44. 44. Ndeba BS, El Alani O, Ghennioui A, Benzaazoua M. Comparative analysis of seasonal wind power using Weibull, Rayleigh and Champernowne distributions. Sci Rep. 2025;15(1):2533. pmid:39833214

* View Article

* PubMed/NCBI

* Google Scholar

45. 45. Niu R, Tian W, Zhang Y. Discriminating among generalized exponential, weighted exponential and weibull distributions. Mathematics. 2023;11(18):3847.

* View Article

* Google Scholar

46. 46. Ali S, Shah I, Wang L, Yue Z. A comparison of shewhart-type time-between-events control charts based on the renewal process. IEEE Access. 2020;8:113683–701.

* View Article

* Google Scholar

47. 47. Alshanbari HM, Al-Mofleh H, Seong J-T, Khosa SK. A new tangent-generated probabilistic approach with symmetrical and asymmetrical natures: monte carlo simulation with reliability applications. Symmetry. 2023;15(11):2066.

* View Article

* Google Scholar

48. 48. Raza A, Ali S, Shah I, Al-Rezami AY, Almazah MMA. A comparative analysis of mean charts assuming Weibull and generalized exponential distributions. Heliyon. 2024;10(21):e40001. pmid:39568830

* View Article

* PubMed/NCBI

* Google Scholar

49. 49. Al-Saiary ZA, Bakoban RA. The topp-leone generalized inverted exponential distribution with real data applications. Entropy. 2020;22(10):1144.

* View Article

* Google Scholar

Citation: Gong Q, Zhang Z, Zeng L, Ren H (2025) A novel extended inverse Weibull distribution: Statistical analysis and application. PLoS One 20(10): e0335555. https://doi.org/10.1371/journal.pone.0335555

About the Authors:

Qin Gong

Roles: Conceptualization, Data curation, Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing

Affiliation: College of Science, Jiangxi University of Science and Technology, Ganzhou, China

Ziwen Zhang

Roles: Data curation, Investigation, Methodology, Supervision, Visualization, Writing – original draft

Affiliation: College of Science, Jiangxi University of Science and Technology, Ganzhou, China

Lihua Zeng

Roles: Formal analysis, Funding acquisition, Project administration, Validation, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Teaching Department of Basic Subjects, Jiangxi University of Science and Technology, Nanchang, China

ORICD: https://orcid.org/0009-0007-6666-5813

Haiping Ren

Roles: Investigation, Project administration, Resources, Validation, Writing – review & editing

Affiliation: Teaching Department of Basic Subjects, Jiangxi University of Science and Technology, Nanchang, China

References

1. Alshanbari HM, Odhah OH, Al-Mofleh H, Ahmad Z, Khosa SK, El-Bagoury AA-AH. A new flexible Weibull extension model: different estimation methods and modeling an extreme value data. Heliyon. 2023;9(11):e21704. pmid:38027837

2. Khan F, Ahmad Z, Khosa SK, Alomair MA, Alomair AM, Alsharidi AK. A new modification of the flexible Weibull distribution based on power transformation: Monte Carlo simulation and applications. Heliyon. 2023;9(6):e17238. pmid:37426796

3. Al-Marzouki S, Alrashidi A, Chesneau C, Elgarhy M, Khashab RH, Nasiru S. On improved fitting using a new probability distribution and artificial neural network: application. AIP Advances. 2023;13(11).

4. Liu X, Ji J, Alrashidi A, Almulhim FA, Alshawarbeh E, Seong J-T. A new probabilistic model with mixed-state failure rates: modeling time-to-event scenarios in reliability and music engineering. Alexandria Eng J. 2024;96:99–111.

5. Shi X, Hu J, Gao R. Integrated modeling of the sports and reliability data: implications of the probabilistic model and deep learning approaches. Alexandria Eng J. 2024;94:274–86.

6. Gemeay AM, Bashiru SO, Sapkota LP, Kayid M, Dutta S, Mohammad S. A new power transformed distribution with applications to radiotherapy and environmental datasets. J Rad Res Appl Sci. 2025;18(2):101339.

7. Chaisee K, Khamkong M, Paksaranuwat P. A new extension of the exponentiated Weibull–Poisson family using the gamma-exponentiated weibull distribution: development and applications. Symmetry. 2024;16(7):780.

8. Tu X, Kong J, Fu Q, Chang S, Zhang K, Alballa T, et al. A new trigonometric-inspired probability distribution: a simulation study and applications in reliability and hydrology. Alexandria Eng J. 2025;113:181–94.

9. Zhu JQ, El-Din MMM, Seong JT, Tashkandy YA, Bakr ME, Kumar A. A new sine-arisen probabilistic model and artificial neural network methods for statistical modeling of the music engineering and reliability data. Alex Eng J. 2024;101:108–17.

10. Almazah MMA, Ullah K, Hussam E, Hossain MdM, Aldallal R, Riad FH. New statistical approaches for modeling the COVID‐19 data set: a case study in the medical sector. Complexity. 2022;2022(1).

11. AbaOud M, Almuqrin MA. The weighted inverse Weibull distribution: Heavy-tailed characteristics, Monte Carlo simulation with medical application. Alexandria Eng J. 2024;102:99–107.

12. Afify AZ, Shawky AI, Nassar M. A new inverse Weibull distribution: properties, classical and Bayesian estimation with applications. KJS publishes peer-review articles in Mathematics, Computer Science, Physics, Statistics, Biology, Chemistry, and Earth and Environmental Sciences. 2021;48(3).

13. Hassan MKH, Aslam M. DUS-neutrosophic multivariate inverse Weibull distribution: properties and applications. Complex Intell Syst. 2023;9(5):5679–91.

14. Shakhatreh MK, Lemonte AJ, Cordeiro GM. On the generalized extended exponential-Weibull distribution: properties and different methods of estimation. Inter J Comp Math. 2019;97(5):1029–57.

15. Aslam M, Tahir M, Hussain Z. Reliability analysis of 3-component mixture of distributions. Scientia Iranica. 2017;0(0):0–0.

16. Hassan AS, Alsadat N, Chesneau C, Elgarhy M, Kayid M, Nasiru S, et al. Inverse power XLindley distribution with statistical inference and applications to engineering data. Sci Rep. 2025;15(1):4385. pmid:39910109

17. Ogunde AA, Dutta S, Almetawally EM. Half logistic generalized rayleigh distribution for modeling hydrological data. Ann Data Sci. 2024;12(2):667–94.

18. Gemeay AM, Karakaya K, Bakr ME, Balogun OS, Atchadé MN, Hussam E. Power Lambert uniform distribution: statistical properties, actuarial measures, regression analysis, and applications. AIP Adv. 2023;13(9).

19. Irshad MR, Aswathy S, Maya R, Al-Omari AI, Alomani G. A flexible model for bounded data with bathtub shaped hazard rate function and applications. MATH. 2024;9(9):24810–31.

20. Kayid M. Further results on the proportional vitalities model. Entropy (Basel). 2021;23(9):1201. pmid:34573826

21. Estañón CR, Montgomery HE Jr, Angulo JC, Aquino N. The confined helium atom: an information–theoretic approach. Int J Quant Chem. 2024;124(4).

22. Shrahili M. Some new results involving residual Renyi’s information measure for $ k $-record values. MATH. 2024;9(5):13313–35.

23. Chennaf S, Ben Amor J. Renyi entropy of uncertain random variables and its application to portfolio selection. Soft Comput. 2023;27(16):11569–85.

24. Bialas P, Burda Z, Johnston DA. Rényi entropy of zeta urns. Phys Rev E. 2023;108(6–1):064108. pmid:38243505

25. Fazlollahi HR. Rényi entropy correction to expanding universe. Eur Phys J C. 2023;83(1).

26. El-Raouf MMA, AbaOud M. A novel extension of generalized Rayleigh model with engineering applications. Alexandria Eng J. 2023;73:269–83.

27. Dutta S, Alqifari HN, Almohaimeed A. Bayesian and non-bayesian inference for logistic-exponential distribution using improved adaptive type-II progressively censored data. PLoS One. 2024;19(5):e0298638. pmid:38753595

28. Akhtar N, Khan SA, Ismail EAA, Awwad FA, Gul T, Khan AA. Bayesian estimation of 3-component mixture geometric distribution under different loss functions. AIP Adv. 2024;14(1).

29. Anderson TW, Darling DA. Asymptotic theory of certain “Goodness of Fit” criteria based on stochastic processes. Ann Math Statist. 1952;23(2):193–212.

30. Alsadat N, Hassan AS, Elgarhy M, Chesneau C, Mohamed RE. An efficient stress–strength reliability estimate of the unit gompertz distribution using ranked set sampling. Symmetry. 2023;15(5):1121.

31. Shafiq A, Sindhu TN, Dey S, Lone SA, Abushal TA. Statistical features and estimation methods for half-logistic unit-gompertz type-I model. Mathematics. 2023;11(4):1007.

32. Aboraya M, Ali MM, Yousof HM, Ibrahim M. A novel lomax extension with statistical properties, copulas, different estimation methods and applications. Bull Malays Math Sci Soc. 2022;45(S1):85–120.

33. Wang C, Zhu H. Tests of fit for the power function lognormal distribution. PLoS One. 2024;19(2):e0298309. pmid:38386634

34. Kutzker T, Klein N, Wied D. Flexible specification testing in quantile regression models. Scandinavian J Stat. 2023;51(1):355–83.

35. Cavaliere G, Perera I, Rahbek A. Specification tests for GARCH processes with nuisance parameters on the boundary. J Busin Econ Stat. 2023;42(1):197–214.

36. Anis MZ, Okorie IE, Ahsanullah M. A review of the Rayleigh distribution: properties, estimation and application to COVID-19 data. Bull Malays Math Sci Soc. 2023;47(1).

37. Lim JS, Kim YD, Woo J-C. Approximated uncertainty propagation of correlated independent variables using the ordinary least squares estimator. Molecules. 2024;29(6):1248. pmid:38542885

38. Mardones C, Orellana A. Estimating production functions to determine the marginal economic value and price elasticity of water in the Chilean industrial sector. Utilities Policy. 2023;85:101677.

39. Ahsan-ul-Haq M, Farooq MU, Nagy M, Mansi AH, Habineza A, Marzouk W. A new flexible distribution: statistical inference with application. AIP Adv. 2024;14(3).

40. Dong Y, Soale A-N, Power MD. A selective review of sufficient dimension reduction for multivariate response regression. J Stat Plan Inference. 2023;226:63–70.

41. Jia J-M, Yan Z-Z, Peng X-Y. A new discrete extended Weibull distribution. IEEE Access. 2019;7:175474–86.

42. Unal C, Cakmakyapan S, Ozel G. Alpha power inverted exponential distribution: properties and application. Gazi Univ J Sci. 2018;31:954–65.

43. Kudo S, Fujimoto M, Sato T, Nagano A. Determination of the optimal number of linked rigid-bodies of the trunk during walking and running based on Akaike’s information criterion. Gait Posture. 2020;77:264–8. pmid:32087596

44. Ndeba BS, El Alani O, Ghennioui A, Benzaazoua M. Comparative analysis of seasonal wind power using Weibull, Rayleigh and Champernowne distributions. Sci Rep. 2025;15(1):2533. pmid:39833214

45. Niu R, Tian W, Zhang Y. Discriminating among generalized exponential, weighted exponential and weibull distributions. Mathematics. 2023;11(18):3847.

46. Ali S, Shah I, Wang L, Yue Z. A comparison of shewhart-type time-between-events control charts based on the renewal process. IEEE Access. 2020;8:113683–701.

47. Alshanbari HM, Al-Mofleh H, Seong J-T, Khosa SK. A new tangent-generated probabilistic approach with symmetrical and asymmetrical natures: monte carlo simulation with reliability applications. Symmetry. 2023;15(11):2066.

48. Raza A, Ali S, Shah I, Al-Rezami AY, Almazah MMA. A comparative analysis of mean charts assuming Weibull and generalized exponential distributions. Heliyon. 2024;10(21):e40001. pmid:39568830

49. Al-Saiary ZA, Bakoban RA. The topp-leone generalized inverted exponential distribution with real data applications. Entropy. 2020;22(10):1144.