Alpha-Power Exponentiated Inverse Rayleigh distribution and its applications to real and simulated data

The main goal of the current paper is to contribute to the existing literature of probability distributions. In this paper, a new probability distribution is generated by using the Alpha Power Family of distributions with the aim to model the data with non-monotonic failure rates and provides a better fit. The proposed distribution is called Alpha Power Exponentiated Inverse Rayleigh or in short APEIR distribution. Various statistical properties have been investigated including they are the order statistics, moments, residual life function, mean waiting time, quantiles, entropy, and stress-strength parameter. To estimate the parameters of the proposed distribution, the maximum likelihood method is employed. It has been proved theoretically that the proposed distribution provides a better fit to the data with monotonic as well as non-monotonic hazard rate shapes. Moreover, two real data sets are used to evaluate the significance and flexibility of the proposed distribution as compared to other probability distributions.


Introduction
In statistical theory, the development of new distributions has become a common practice in recent decades; this is done generally by adding an extra parameter [1] to the baseline distribution, using generators [2,3], or by combining two distributions [4]. Ramadan and Magdy [5] produced a new probability distribution by applying the Inverse Weibull (IW) to the Alpha Power Family of distribution. Alzaatreh et al. [2] introduced T-X family of continuous distributions by interchanging the probability density function of any continuous random variable with the probability density function of Beta distribution. Lee et al. [3] developed a technique of generating single variable continuous distributions. Jones [6] applied the Beta distribution to the family of distribution presented by Eugene et al. [7].
The main purpose of such an amendment to the existing distributions is to model the real data both with monotonically and non-monotonically hazard rate functions. Secondly, to increase the model flexibility of the complex data structures as compared to existing probability distributions. Because the existing distribution has some limitations so model the complex data structures, for example, Exponential and Weibull distributions fail the real data following a non-monotonic failure rate functions. In this aim of presenting the paper is to contribute a new probability distribution that will model the data with both monotonically and non-monotonically hazard rate functions. The proposed model will also increase the model flexibility as compared to other models.

Alpha Power Transformation
In the Recent past, Mahdavi and Kundu [8] suggested a new technique, called Alpha Power Transformation (APT), for including an additional parameter in the life time distributions. The primary purpose of this family was to utilize the non-symmetrical behavior of the parent distribution. The Alpha Power Transformation is defined by Let X is a continuous random variable with F(x) as Cumulative Distribution Function, the Cumulative Distribution Function of Alpha Power Transformation is as follows; The associated Probability density function is given below The Alpha Power transformation has been used by many researchers, for example, Dey et.al [9] explored the new probability distributions by applying the Exponential and Rayleigh distribution to the Alpha Power Family of distributions. By using the same Transformation Nassar et al. [10] produced Alpha Power Weibull distribution, Alpha Power Inverse Weibull distribution was produced by Ramadan and Magdy [5], Alpha Power Transformed Extended Exponential distribution by Hassan et al [11].
The main aim of the paper is to produce a new probability distribution by using the Alpha Power family of distribution. In this paper, we considered the Exponentiated Inverse Rayleigh distribution as a baseline distribution presented by Rehman and Sajjad [12]. The Exponentiated Inverse Rayleigh distribution is the extension of the Inverse Rayleigh distribution presented by Voda [13]. He discussed various statistical properties such as moment generating function, survival function, and order statistics. A random variable X is said to be Inverse Rayleigh if it possesses the following Pdf and Cdf FðyÞ The Exponentiated Inverse Rayleigh (EIR) distribution has the following pdf and cdf; f ðxÞ ¼ 2xy FðxÞ ¼ e À ay The current study is linked with the introduction of a novel distribution which is stated as Alpha Power Exponentiated Inverse Rayleigh (APEIR) distribution. Various statistical properties of the APEIR distribution are studied such as quantile function, median, mode, moment generating function and r th moment, order statistics, mean residual life function, and stress strength parameter are obtained and discussed. Furthermore, an expression for the Renyi entropy and for the Mean Waiting Time has been explored. The estimation of the parameters is done by using the maximum likelihood. In order to prove the flexibility of the model, we considered the application by using two real data sets as well as simulated data.

Alpha Power Exponentiated Inverse Rayleigh (APEIR) distribution
By applying the cumulative distribution function of the Exponentiated inverse Rayleigh distribution to the ALPF, we obtained the following Cdf and Pdf for the APEIR and is given by  The hazard and survival function corresponding to the probability density function are as follows Fig 2 shows the hazard rate function and survival function of APEIR distribution with various values of parameters. Clearly, the hazard rate function of APEIR distribution is unimodal and positively skewed for α > 1.
Lemma 1: If α < 1 then f(x) is a decreasing function, this implies that f APEIR (x) is decreasing function.
Proof: If f(x) is a differentiable function and if its first order derivative or d dx logf x ð Þ < 0 for x in (α, β, θ) then f(x) is a decreasing function and vice versa.
Taking the first derivative of logf APEIR (x) i.e.
Taking second order derivative of Eq (11), we get

Quantile function
Let X *APEIR (α β, θ) then its Quantile function is given below; Median of APEIR distribution is obtained by substituting u = 1/2 in Eq (13), we get

Mode
Mode of APEIR distribution is that point by which the distribution reaches its maximum point and it is obtained by solving the following equation We finally, obtained the result

Rth moment of APEIR distribution
Let X *APEIR (α β, θ), then the expression of its r th moment is as follows; Using y = x −2 and series notation a À z ¼ P 1 k¼0 ðÀ logaÞ k k! ðÀ zÞ k and logz ¼ À ðÀ 1Þ j ðÀ 1þzÞ j j , we get the final result as where l ¼ À r 2 . which is the required result.

Moment Generating Function (MGF)
Let X~APEIR(α, β, θ) then the expression for its MGF is as follows; Using t r x r r! and the series representation a À z ¼ The MGF of APEIR distribution has the following form where l ¼ À r 2 .

Mean residual life function
Let X be the survival time of an object having pdf "f(x)" and survival function specified in Eq (10), the mean residual life function is the average remaining lifespan, which is a component survived up to time t. The mean residual life function, say, μ(t) has the following expression.

Order statistics
Let X 1 , X 2 , X 3 , . . ., X n be a random sample of size n from APEIR distribution and let X (1) � X (2) � . . . � X (n) denote the order statistics. Let X i:n denotes the i th order statistics, then the Probability Density function of X i:n is given by putting Eqs (7) and (8) of APEIR distribution in (24), we obtain the pdf of i th order statistic for x > 0, as is mentioned below by inserting i = 1 in Eq (25), we obtain the smallest order statistic as follows: For largest order statistic insert i = n in Eq (25), we get Put i = n /2 in Eq (25), to obtain the distribution of median, we have x 2 a e À by x 2 a e À by
using Eqs (7) and (8) of APEIR distribution then Stress Strength Parameter is given as; after simplification, we finally obtained the equation for Stress-Strength Parameter.
Lemma 3: Let X~APEIR(α, θ, β), then its Renyi entropy is defined by where l ¼ 3ðvÀ 1Þ 2 : Proof: For APEIR distribution, Renyi entropy has the following expression; The result can be obtained easily by substituting a À z ¼ P 1 k¼0 ðÀ logaÞ k k! ðÀ zÞ k : Lemma4: The Mean Waiting Time say � mðtÞ of APEIR distribution is as follows; Proof: For APEIR, the mean waiting time is given by

Simulations study
The parameter estimates of APEIR distribution, their Mean Square Error (MSE) as well as the bias measure are computed using a simulation study with 1000 replications each with a sample of size n = 30, 70, 130 and 170. A simulated data is generated from APEIR distribution using the following expression ; where U follows a standard uniform distribution. The average bias and MSE are computed by using the mathematical formulae as under β, θ). The average bias and MSEs are given in Table 1. It has been observed that the MSEs and bias of the estimates are decreasing for all parameter combinations with the increase in the sample of size n.

Applications
In this section, we provide two applications of the proposed distribution to the lifetime data. The performance of the suggested model is checked by the goodness of fit criteria including they are the AIC, CAIC, BIC, HQIC, and the P-value. For more details of the goodness of fit criteria, we refer to see [14][15][16][17][18][19]. In general, with fever values of these statistics, a probability model would perform better than others. The proposed model is compared with Exponentiated Inverse Rayleigh distribution by Rehman and Sajjad [12], Weibull Rayleigh distribution by Merovci and Elbatal [20], Generalized Rayleigh distribution by Raqab and Madi [21], two parameter Rayleigh distribution by Dey et.al [22], Transmuted inverse Rayleigh distribution by Afaq et al [23] and modified inverse Rayleigh distribution by Muhammad [24]. The probability functions of these distributions are given by x 2 y; b; X > 0: a; y; b; X > 0: • Generalized Rayleigh (GR) Distribution f ðxÞ ¼ 2ag 2 xe À ðgxÞ 2 ð1 À e À ðgxÞ 2 Þ aÀ 1 a; g; X > 0: • Two Parameter Rayleigh (TPR) Distribution f ðxÞ ¼ 2aðx À mÞe À aðxÀ mÞ 2 x > m; a > 0: a; y; X > 0: • Transmuted Inverse Rayleigh Distribution.
f ðxÞ ¼ 2y  Table 2 describes the parameter values of the probability models and also describes the goodness of fit measures. It is evident that the goodness of fit measures has fever values for the proposed model and hence it is concluded that the proposed model increased the flexibility of the model.
In Fig 4, the histogram represents the theoretical densities of the Alpha Power Exponentiated Inverse Rayleigh (APEIR), Two Parameter Rayleigh (TPR) and Exponentiated Inverse Rayleigh (EIR) by continuous red color line, dotted blue line and dotted green line respectively. It is evident from the above figure that the Alpha Power Exponentiated Inverse Rayleigh (APEIR) is leptokurtic and positively skewed as compared to other densities. Furthermore, the graph suggests that the Alpha Power Exponentiated Inverse Rayleigh (APEIR) distribution is less thick as compared to Two Parameter Rayleigh (TPR) distribution and thicker than Exponentiated Inverse Rayleigh (EIR) in the tail.  If the plot of empirical against the theoretical CDFs is observed, then Alpha Power Exponentiated Inverse Rayleigh (APEIR) provides a better fit as compared to Exponentiated Inverse Rayleigh (EIR) and Two Parameter Rayleigh (TPR).  Table 3 describes the MLE of the probability models and describe the goodness of fit measures. Again, it is concluded that by increasing another parameter, we get a more significant result as compared to others.

Conclusion
The paper presents a new probability distribution called Alpha Power Exponentiated Inverse Rayleigh (APEIR) distribution. The objective of the proposed distribution is to model the data with both monotonic and non-monotonic hazard rate shapes. The proposed distribution is of keen interest due its desirable properties. To estimate the parameters of the new distribution, Maximum likelihood estimation procedure is used. Furthermore, to evaluate the performance of the proposed distribution, it was fitted to two real data sets. The results showed that the new distribution provides a better fit to these data sets as compared to other versions of the Rayleigh distributions. Future researchers may derive new flexible distributions by using transmutation technique, or by increasing the scale or shape parameter to the proposed distributions in this paper. Further one can study the Bayesian analysis by choosing informative and noninformative priors.