A Novel Modified Alpha Power Transformed Weibull Distribution and Its Engineering Applications

2023-02-27 10:40RefahAlotaibiHassanOkashaMazenNassarandAhmedElshahhat

Refah Alotaibi,Hassan Okasha,Mazen Nassar and Ahmed Elshahhat

1Department of Mathematical Sciences,College of Science,Princess Nourah bint Abdulrahman University,P.O.Box 84428,Riyadh,11671,Saudi Arabia

2Department of Statistics,Faculty of Science,King Abdulaziz University,Jeddah,80200,Saudi Arabia

3Department of Mathematics,Faculty of Science,Al-Azhar University,Nasr City,Cairo,11884,Egypt

4Department of Statistics,Faculty of Commerce,Zagazig University,Zagazig,44519,Egypt

5Faculty of Technology and Development,Zagazig University,Zagazig,44519,Egypt

ABSTRACT This paper suggests a new modified version of the traditional Weibull distribution by adding a new shape parameter utilising the modified alpha power transformed technique.We refer to the new model as modified alpha power transformed Weibull distribution.The attractiveness and significance of the new distribution lie in its power to model monotone and non-monotone failure rate functions,which are quite familiar in environmental investigations.Its hazard rate function can be decreasing, increasing, bathtub and upside-down then bathtub shaped.Diverse structural properties of the proposed model are acquired including quantile function,moments,entropies,order statistics,residual life and reversed failure rate function.The parameters of the distribution were estimated using the maximum likelihood function.The maximum likelihood method is employed to estimate the model parameters and the approximate confidence intervals are also computed.Via a simulation study, the performance of the point and interval estimates are compared using different criteria.Employing real lifetime data sets, we verify that the offered model furnishes a better fit than some other lifetime models including Weibull,gamma and alpha power Weibull models.

KEYWORDS Weibull distribution;modified alpha power transformation method;maximum likelihood;entropy;order statistics

1 Introduction

The Weibull distribution is a wildly favoured lifetime distribution in reliability studies.It is naturally employed for studying biological,hydrological and medical data sets.The Weibull distribution is frequently employed as a suitable alternative to well-known distributions such as exponential,gamma and inverse Weibull distributions.The random variableXis expressed to have a Weibull distribution if its probability density function(PDF)is given by

whereλ >0 andθ >0 are the scale and shape parameters respectively.Also,its cumulative distribution function(CDF)takes the form

One of the major weaknesses of the Weibull distribution is that it does not deliver an adequate fit for some applications, particularly, when the hazard rates are upside-down bathtub or bathtub shapes.To overcome this disadvantage, several investigators have developed various generalizations and modifications of the Weibull distribution to model different types of data in recent years.The generalized Weibull distribution was introduced in reference[1,2]by adding a shape parameter to the Weibull distribution.Likewise,Xie et al.[3]proposed the additive Weibull distribution,the generalized modified Weibull distribution by[4],KumaraswamyWeibull distribution by[5],beta Sarhan-Zaindin modified Weibull distribution by[6],Weibull-Weibull distribution by[7],alpha power Weibull(APW)distribution by [8], log-normal modified Weibull distribution and its reliability implications by [9],generalized extended exponential-Weibull distribution by[10],Poisson modified weibull distribution by [11], alpha logarithmic transformed Weibull distribution by [12] and logarithmic transformed Weibull by[13].

Recently, Mahdavi et al.[14] suggested a new approach to present an additional parameter to a class of distributions for more additional flexibility.The offered technique is named alpha power transformation (APT) and it is worthwhile to incorporate skewness into a family of distributions.They studied the main properties of the APT method and introduced an extension to the exponential distribution using the new approach.Many authors used the same technique to introduce new generalizations of some well-known distribution.For example, alpha power Weibull distribution by[8],alpha power generalized exponential by[15],alpha power transformed inverse Lindley distribution by[16]and alpha power Gompertz distribution by[17].To add more flexibility to the APT method,Alotaibi et al.[18]proposed a new form of the APT method which is called the modified APT(MAPT)method.The CDF and the PDF of the MAPT method are,respectively,given by

whereαis a shape parameter,G(x) is a baseline distribution andg(x) =(x).Reference [18]investigated the major properties of the new generator.They also used the new method to develop a new version of the exponential distribution and studied some of its properties.

Modelling real data employing generalized distributions stays vital nowadays.Multiple generalized distributions have been considered and used in different domains.Nevertheless,there still remain considerable necessary issues applying real data,which are not handled by available models.The main purpose of this article is to offer a new unexplored version of the traditional Weibull distribution.To achieve this goal, we utilize the MAPT method to add a new shape parameter to the baseline CDF given by (2).We refer to the new model as modified alpha power transformed Weibull (MAPTW)distribution which contains one scale and two shape parameters.The main properties of the new model including, quantile, mixture expansion, moments, entropies and order statistics are derived.The unknown parameters are estimated using the maximum likelihood estimation method and the approximate confidence intervals (ACIs) of the unknown parameters are also constructed.We are encouraged to introduce the MAPTW distribution because

1.It contains some well-known distributions as special cases including,exponential,Weibull and MAPT exponential (MAPTE) distributions.So, at least, the MAPTW distribution contains the main properties of these important distributions which are commonly used in modelling lifetime data.Also,the main characteristics of these distributions can be obtained directly from the MAPTW distribution.

2.It is appropriate to model positively skewed,negatively skewed and approximately symmetric data which may not be adequately modelled by other competitive models.

3.It can model decreasing,increasing,bathtub and upside-down then bathtub hazard rates which are often faced in real-life problems.

4.It is indicated in Section 2 that the MAPTW distribution can be shown as a combination of the Weibull distribution.This property is especially helpful to derive the properties of the MAPTW distribution directly from the Weibull distribution.

5.The analysis of two real data sets demonstrate that the MAPTW distribution compares satisfactorily with different competing lifetime distributions in modelling engineering data sets.As a results,the MAPTW distribution can be able to model engineering applications rather that some well known and recently proposed generalized models.This is due to the flexibility of its PDF and hazard rate function(HRF).

One of the main advantages of this distribution over many other generalized distributions is that it gives the main forms of the HRF with two shape parameters only.The same advantage can be found in its PDF which can be used to model positively and negatively skewed as well as approximately symmetric data which may not be adequately fitted by other distributions.This makes the new distribution more flexible, especially for modelling data sets when studying reliability experiments,electronics,materials,automotive industries and many engineering applications.Another motivation for this study is examining the behaviour of the maximum likelihood estimators as well as the ACIs by considering different sample sizes and different true parameter values.The rest of the paper is classified as follows: We present the MAPTW distribution in Section 2.The main properties of the MAPTW distribution are derived in Section 3.In Section 4, we discuss the maximum likelihood estimates (MLEs) as well as the ACIs of the model parameter.In Section 5, a simulation study is conducted to compare the performance of point and interval estimates.The effectiveness of the MAPTW distribution is depicted by studying two real data sets from the engineering field in Section 6.Finally,in Section 7 we conclude the paper.

2 The Modified Alpha Power Transformed Weibull Distribution

The MAPTW distribution is discussed in this section.IfXhas a Weibull distribution with PDF and CDF, respectively, provided by (1) and (2), then the CDF of the MAPTW distribution is given by(3).

and the corresponding PDF is given by

whereαandθare the shape parameters whileλis the scale parameter.It is observed thatF(x;α,λ,1)reduces to the MAPTE by[18],F(x;→1,λ,1)tends to the exponential distribution,F(x;→1,λ,2)tends to the Rayleigh distribution,F(x;α,λ,2)tends to the MAPT Rayleigh distribution andF(x;→1,λ,θ)tends to the Weibull distribution.It is to mentioned hare that the MAPTW distribution can be considered as a weighted version of the APW distribution with weight functionThe reliability function(RF)and HRF of the MAPTW distribution are given,respectively,by

Henceforward,Xis used to denote the random variable that has the PDF in(6).Fig.1 shows the various forms of the PDF of the MAPTW distribution given by(6)by assuming different values of the shape parameters and by considering the scale parameterλto be one in all the issues.Fig.1 indicates that the MAPTW distribution can model the left-skewed,right-skewed and approximated symmetric data.Also,Fig.2 presents the different shapes of the HRF given by(8)by considering various values of the shape parameters and by assuming the scale parameter to be one in all the cases.Fig.2 reveals that the HRF of the MAPTW distribution can be decreasing, increasing, bathtub and upside-down then bathtub shaped.

Figure 1:Plots for the PDF of the MAPTW distribution

Figure 2:Plots for the HRF of the MAPTW distribution

To derive the main properties of the MAPTW distribution,we provide a valuable expansion for the PDF in(6).Consider the following power series:

and the generalized binomial expansion in the form

Using the series in(9)and(10),one can write the PDF in(6)as a mixture representation of the Weibull densities as follows:

whereg(x;(j+1)λ,θ)is the PDF of the Weibull distribution with scale parameter(j+1)λand shape parameterθ,where

Moreover,we can write the linear representation of the CDF of the MAPTW distribution after integrating(11),as follows:

whereG(x;(j+1)λ)is the CDF of the Weibull distribution with scale and shape parameters(j+1)λandθ,respectively.

3 Structure Properties

We provide some essential properties for the MAPTW distribution in this section,such as quantile,moments,entropies,order statistics,residual life,and reversed failure rate function.

3.1 Quantile Function

The quantile function(QF)of the MAPTW distribution,sayx=Q(p),can be obtained from the CDF in(5)as

The QF in(13)can be used to generate a random sample from the MAPTW distribution by takingpto beU(0,1).

3.2 Moments

Moments play an essential feature in statistical theory and numerous necessary characteristics of any distribution can be examined via moments.For the MAPTW distribution,thenthmoment can be derived as

whereΓ(.)is the gamma function.Similarly,for the MAPTW distribution we can obtain thenthinverse moment by using(6)as follows:

Practically,the mean of the MAPTW distribution can be obtained from(14)by settingn= 1 as follows:

Similarly,the second moment can be obtained by puttingn=2 in(14)as

Then,once can obtain the variance of the MAPTW distribution using the following relation:

Therthcentral momentμrofXis derived as Moreover, The skewness (Sk) and kurtosis (Ku)measures can be calculated employing the next formulas

The following Lemma explain various kinds of moments of the MAPTW distribution such as incomplete moments(IMs),moment generating function(MGF),characteristic generating function(CGF)and conditional moments(CMs).

Lemma 3.1.IfX~MAPTIW,then

1.The IMs ofXis

whereγ(.,.)is the incomplete gamma function.

2.The MGF ofXis

Figure 3:Mean,variance,skewness and kurtosis of the MAPTW distribution

3.The CGF ofXis

4.The CMs ofXis

whereψn(t)is the IMs andF(x)is the CDF given by(5).

3.3 Entropies

Entropy is one of the numerous popular techniques used to calculate the uncertainty corresponding to a random variable and it was initially employed in the field of physics.Estimating entropy is a significant matter in many fields including statistics and biological phenomenon.High entropy is guided to minor information found in data.For the MAPTW distribution and from (6), the Renyi entropy(RE)can be acquired as follows:

Also,Shannon’s entropy can be derived by Limitingρ↑1 in(15)as

3.4 Order Statistics

Order statistics contain a broad range of applications in statistics, including nonparametric statistics,life testing,quality control monitoring and reliability analysis.LetX1,...,Xnbe a random sample from the MAPTW distribution with CDF and PDF given by(5)and(6),respectively.Suppose also that the order statistics ofX1,...,Xnare denoted byX(1),...,X(n).Then,the PDF of theithorder statistic,X(i),is expressed as

whereB(.,.)is the beta function.Using(5),(6)and a power series expansion,we obtain

Substituting(17)in Eq.(16),the PDF ofX(i)reduces to

whereg(x;(d+1)λ,θ)is the PDF of the Weibull distribution with scale parameter(d+1)λand shape parameterθ,where

3.5 Probability Weighted Moments

The(u,v)thprobability weighted moments can be defined as

For the of the MAPTW distribution,it follows from(5)and(6)that

Substituting (20) in Eq.(19), the probability weighted moments of the MAPTW distribution is given by

3.6 Residual Life and Reversed Failure Rate Function

Therthmoment of the residual life(RL)of the random variableXis defined as follows:

Using(6)and applying the binomial expansion of(x-t)r,we have

whereg(x;(j+1)λ,θ)is the PDF of the Weibull distribution,then

whereΓ(a,y)=dzdenotes the incomplete gamma function.Therthmoment of the reversed RL can be derived using the general formula

Using(5),(6)and applying the binomial expansion of(t-x)r,we can write

3.7 Stress-Strength Model

LetY1~MAPTW(α1,θ,λ)andY2~MAPTW(α2,θ,λ).IfY1represents stress andY2represents strength, then the stress-strength parameter, denoted byR, for the MAPTW distribution.Here, we consider two cases as follows:

Case one:Whenα1≠α2

Using the series expansion in(10),the last equation takes the form

Employing the series expansion in(9),we obtain

Case two:Whenα1≠α2,λ1≠λ2andθ1≠θ2

LetY1~ MAPTW(α1,θ1,λ1) andY2~ MAPTW(α2,θ1,λ2).IfY1represents stress andY2represents strength,then the stress-strength parameter,denoted byR,for the MAPTW distribution is given by

Using the series expansion in(10),we have

Utilizing the series expansion in(9),it comes

4 Estimation of the Parameters

Here, the method of maximum likelihood is employed to estimate the unknown parameters of the MAPTW distribution.Moreover, the approximate confidence intervals (ACIs) of the unknown parameters are obtained.Assume that we have a random sample of sizentaken from the MAPTW distribution with PDF given by (6), then we can express the log-likelihood function in this case as follows:

To get the maximum likelihood estimates(MLEs)ofα,λandθdenoted by ˆα, ˆλand ˆθ,one can maximize the objective function in(21)with respect toα,λandθ.An alternative approach to compute the required estimates is to solve the following normal equations simultaneously:

It is to note that there are no closed-form solutions for Eqs.(22)-(24).To get the MLEs in this case,one can direct to any numerical technique for this purpose.Using the asymptotic properties of the MLEs,we can obtain the ACIs ofα,λandθ.It is known thatwhereI0-1(α,λ,θ)is the asymptotic variance-covariance(AVC)of the MLEs.Practically,it is not easy to computeI0-1(α,λ,θ),thus,the approximate AVC of the MLEs denoted bycan be used as follows:

whereIαα,Iαλ=Iλα,Iαθ=Iθα,Iλθ=Iθλ,IλλandIθθare the second derivatives of(21)and given by

Thus,the(1-γ)%ACIs ofα,λandθcan be obtained as follows:

wherezγ/2is the upper(γ/2)thpercentile point of a standard normal distribution.

5 Simulation Study

To evaluate the performance of the MLEs of the MAPTW distribution,a Monte Carlo simulation is conducted based on a sufficiently large 5,000 independent from the MAPTW distribution using different sample sizesn(= 50,100,150,200) and different choices for true values of the model parametersα,λandθasα(= 0.1,0.3,0.5,0.7),λ(= 0.5,0.8,1.2,1.6)andθ(= 0.5,0.8,1.5,2.5).The simulation procedure is carried out based on following the steps:

Step 1:Determine the sample size and the starting values for the parameters.

Step 2:Generate a random sample of sizenfrom the MAPTW distribution from(9).

Step 3:Compute the average estimates with their root mean squared-errors(RMSEs)and relative absolute biases(RABs)ofα,λandθ.

Step 4:Obtain the(1-γ)%CLs of the parametersα,λandθ.

Step 5:Repeat steps 2-5 5000 times.

Step 6:Calculate the average values(AV)of MLEs,RMSEs,RABs and CLs of any function ofα,λandθ(sayφ)are given,respectively,by

whereφiLandφiUare the lower and upper ACI bounds,respectively.

All required numerical results are obtained via R software version 4.0.4 utilizing‘maxLik’package proposed by[19]and later recommended by[20,21].The simulation results are shown in Tables 1 and 2.From Table 1,we observe that the maximum likelihood estimates of the parametersα,λandθare quite satisfactory in terms of minimum RMSEs and RABs.It can also be shown that,asnincreases,the AVMLEs of the parametersα,λandθclose to the true parameter values.This result implies the fact that the MLEs behaved as asymptotically unbiased estimators.For each set of parameter values,asnincreases,the AV-RMSEs and AV-RABs decrease which implying that the MLEs are also consistent.To be more specific,it is noted that the AV-RMSE and AV-RAB for the parameterθdecrease rapidly after the sample size exceeds 50.On the other hand,the AV-RMSEs and AV-RABs for the parametersαandλneed a larger sample size to get a small value for these measures.After revising the computation results, we have observed that some samples from the 1000 generated samples give high estimated values forαandλand this is the reason for the high values of AV-RMSEs and AV-RABs in this case,but the majority of the generated samples get satisfactory results.Further,Table 2 shows that the AV-CLs decrease are narrowed down with the increase in the sample sizenas expected.Graphically,the AV-MSEs, AV-RABs and AV-CLs of the parametersα,λandθare displayed in Fig.4 (when(α,λ,θ)=(0.1,0.5,0.5)as an example).This conclusion refers to the fact that whennincreases more additional information is gathered.

Table 1:The AV-MLEs(first-line),AV-RMSEs(second-line)and AV-RABs(third-line)of α,λ and θ

Table 1 (continued)

Table 1 (continued)

Table 2: The AV-CLs for 95%ACIs of α,λ and θ

6 Engineering Applications

To show how our proposed model works in practice,in this section we present the analysis of two real data sets from engineering science for illustrative purposes.The first data set(denoted by Data-I)consists of 40 observations of the active repair times for airborne communication transceiver,see[22].The second data set(denoted by Data-II)represents the failure times of 84 Aircraft Windshields and been obtained from[23].The values of both data sets I and II are arranged with ascending order and reported in Table 3.

Figure 4:Plot of the RMSEs,RABs(left-panel)and CLs(right-panel)of α,λ and θ

Table 3: The failure times of communication transceiver and aircraft windshields

First, to identify the shape of the HRF of MAPTW distribution, we shall consider a graphical method,given by[24],based on the total time on test(TTT)plot.However,the scaled TTT transform based on both observed data sets I and II is given bywherez(i)is theith order statistic of the observed sample.Graphically,the scaled TTT transform is displayed by plotting(m/n,Gn(m/n)).Using both data sets I and II in Table 3,plots of the empirical and estimated scaled TTT-Transform of the MAPTW distribution are provided in Fig.5.It shows that the TTT plot for the data sets I and II is convex and concave,respectively.So,the failure rate shape for the Data-I set decreases while it for the Data-II set increases.

Figure 5:Empirical and estimated scaled TTT-transform plot of the MAPTW distribution

To monitored the applicability and flexibility of the MAPTW distribution, we compare fit the MAPTW distribution to both data sets I and II with several other competitive models exhibiting various hazard rates, namely: alpha power transformed generalized-exponential (APTGE), alpha power exponential (APE), APW, generalized-exponential (GE), Nadarajah-Haghighi (NH), gamma(G)and Weibull(W)distributions.The corresponding PDFs of the competing models(forx >0)are listed in Table 4.

Table 4: Some competing lifetime models of the MAPTW distribution

To show the feasibility and validity of the MAPTW model along with other seven competing models, several goodness-of-fit measures namely; estimated negative log-likelihood (NL), Akaike information criterion(AIC),consistent Akaike information criterion(CAIC),Bayesian information criterion(BIC)and Hannan-Quinn information criterion(HQIC)are considered.Based on both data sets I and II,the selection criteria as well as the MLEs along with their standard errors(SEs)of the model parameters are calculated and provided in Table 5.It shows that the MAPTW distribution has the smallest values of NCL, AIC, CAIC, BIC, HQIC, K-S and the highestP-value.Moreover,Fig.6 shows graphically the quantile-quantile (QQ) plots of all competitive distributions for given data sets.Furthermore, the relative histograms of both data sets and the fitted densities, as well as the fitted/empirical survival functions,are plotted in Fig.7.It can be seen from Figs.6 and 7,that the graphical presentations support the same numerical findings presented in Table 5.Thus,the proposed MAPTW distribution provides the best fit, for the given data sets, than the APTGE, APTW and other popular lifetime distributions.Since the empirical HRFs of the two data sets are decreasing and increasing,respectively,and most of the competitive models such as Weibull and APW distributions can model such HRFs,the results show that the MAPTW distribution provides a better fit for the two engineering applications in terms of goodness of fit statistics.This is because the MAPTW distribution has a very flexible HRF with different shape behaviour as well as a flexible PDF,especially for skewed data.To sum up,we recommend the use of MAPTW for real practical purposes.

Table 5: The MLEs(SEs) and selection criteria of the MAPTW distribution and other competing models

Figure 6:The QQ plots of MAPTW distribution and its competing models for Data-I and Data-II

Figure 7: Histogram and fitted densities (left panel), empirical and fitted survival functions (right panel)of MAPTW distribution and some competing models

7 Conclusion

In this paper,we have offered a fresh version of Weibull distribution named the modified alpha power transformed Weibull distribution.The proposed distribution is acquired by taking the Weibull distribution as the baseline distribution in the modified alpha power transformed method.Some properties of the new distribution are derived.The hazard rate function of the new distribution can take different shapes possessing decreasing, increasing, bathtub and upside-down then bathtub shaped.Accordingly,it can be viewed perfectly effectively in modelling lifetime data.Via the maximum likelihood method, the point and interval estimates of the model parameters are evaluated.The performance of the point and interval estimates is assessed through a simulation study.The simulation yields verified that the estimates are asymptotically unbiased and consistent.In addition,two real data sets to active repair times for airborne communication transceiver and the failure times of aircraft windshields are investigated.From the empirical results, we can conclude that the new distribution provides a more adequate fit than some known distributions including gamma, Weibull and alpha power Weibull distributions.For future work,it is of interest to study the proposed distribution under different censoring schemes including Type-II and progressive Type-II censoring.Another future work is to investigate the Bayesian estimation of the unknown parameters of the suggested model.

Acknowledgement:The authors would like to thank the Editor, Assistant Editor and Reviewers for their constructive recommendations which assisted us to improve the manuscript.The authors extend their appreciations to Princess Nourah bint Abdulrahman University Researchers Supporting Project No.(PNURSP2022R50),Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia.

Funding Statement:This research was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project No.(PNURSP2022R50), Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.