Next Article in Journal
Partial Sums of the Normalized Le Roy-Type Mittag-Leffler Function
Next Article in Special Issue
A New Probability Distribution: Model, Theory and Analyzing the Recovery Time Data
Previous Article in Journal
Approximation Relations on the Posets of Pseudoultrametrics
Previous Article in Special Issue
Risk Analysis: Changing the Story with the Statistical Stochastic Process and VaR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Markov Bernoulli Lomax with Applications Censored and COVID-19 Drought Mortality Rate Data

by
Bahady I. Mohammed
1,
Yusra A. Tashkandy
2,
Mohmoud M. Abd El-Raouf
3,
Md. Moyazzem Hossain
4 and
Mahmoud E. Bakr
5,*
1
Department of Mathematics, Faculty of Science, Al-Azhar University, Nasr City 11884, Cairo, Egypt
2
Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
3
Basic and Applied Science Institute, Arab Academy for Science, Technology and Maritime Transport (AASTMT), Alexandria P.O. Box 1029, Egypt
4
School of Mathematics, Statistics & Physics, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
5
Giza Engineering Institute GEI, Basic Sciences Department, Giza 12519, Cairo, Egypt
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(5), 439; https://doi.org/10.3390/axioms12050439
Submission received: 10 February 2023 / Revised: 15 April 2023 / Accepted: 24 April 2023 / Published: 28 April 2023

Abstract

:
In this article, we present a Markov Bernoulli Lomax (MB-L) model, which is obtained by a countable mixture of Markov Bernoulli and Lomax distributions, with decreasing and unimodal hazard rate function (HRF). The new model contains Marshall- Olkin Lomax and Lomax distributions as a special case. The mathematical properties, as behavior of probability density function (PDF), HRF, rth moments, moment generating function (MGF) and minimum (maximum) Markov-Bernoulli Geometric (MBG) stable are studied. Moreover, the estimates of the model parameters by maximum likelihood are obtained. The maximum likelihood estimation (MLE), bias and mean squared error (MSE) of MB-L parameters are inspected by simulation study. Finally, a MB-L distribution was fitted to the randomly censored and COVID-19 (complete) data.

1. Introduction

The emergence of modern applications and many developments in various fields, in addition to the limitations of some well-known distributions, lead us to create other distributions that are more suitable for modern applications and free from restrictions.
Gharib et al. [1] used a countable mixture with a Markov Bernoulli geometric model to introduced a new family, which its survival function (SF) is given by:
G ¯ x , α , ρ = α F ¯ x 1 ρ F ¯ x 1 1 1 ρ α F ¯ x , x ,   α > 0 ,   0 ρ < 1
When ρ = 0 , the SF G ¯ x , α , ρ in (1) reduces to the SF G ¯ x , α which is introduced by Marshall and Olkin [2]. Moreover, if, ρ = 1 1 α then G ¯ x , α , ρ in (1) reduces to:
G ¯ x , α = F ¯ x α F x + F ¯ x . x ,   α > 1
For ρ = 0 and α = 1 the SF G ¯ x , α , ρ reduces to F ¯ x .
The SF of Lomax or (Pareto Type-II) distribution is:
F ¯ x = ( 1 + β x ) θ , x > 0 , β , θ > 0 .
The probability density function (PDF) and hazard rate function (HRF) are:
f x = β θ 1 + β x θ 1 ,
h F r = β θ 1 + β x 1 .
Johnson et al. [3] used the Lomax model in practical and theoretical fields as, economics and biological. Moreover, Harris [4], Bryson [5], Cordeiro et al. [6] and Bhagwati Devi [7] are, respectively, used this model in reliability & life testing, income and wealth data, firm size data and Entropy.
In the past few years, several authors have expanded the Lomax distribution due to its importance in life time distributions as: generalized Lomax (Raj Kamal Maurya et al. [8]), Poisson-Lomax (Mohammed et al. [9]), Marshall–Olkin Power Lomax (Muhammad Ahsanul Haq et al. [10]), the type II Topp Leone-Power Lomax (Sirinapa Aryuyuen and Winai Bodhisuwan [11]), new weighted Lomax (Huda M. Alshanbari et al. [12]) and reflected-shifted-truncated Lomax Distribution (Sanku Dey et al. [13]), there are other generalizations of the Lomax distribution, and they different in terms of the form of the PDF and the behavior of the HRF; see, Ghitany et al. [14], Lemonte and Cordeiro [15], Cordeiro et al. [16], Al-Zahrani and Sagor [17,18], Tahir et al. [19], El-Bassiouny et al. [20], Rady et al. [21] and Cooray et al. [22], Wael S. Abu El Azm et al. [23], Hassan Alsuhabi et al. [24] and Adebisi A. Ogunde [25].
If we put F ¯ x = ( 1 + β x ) θ , x > 0 , β , θ > 0 , which is the SF of the Lomax distribution, in (1), we have the MB-L (α, β, θ, ρ ) model is:
G ¯ x , α , β , θ , ρ = α ( 1 + β x ) θ 1 ρ ( 1 + β x ) θ 1 1 1 ρ α ( 1 + β x ) θ , x > 0 ,   α , λ > 0 ,   0 ρ < 1

2. The PDF of the MB-L Model

From Equation (3), we have:
g x = α β 1 + x β 1 θ θ 1 + x β 2 θ 1 + α + 2 1 + x β θ ρ + α ρ 2 1 + α + 1 + x β θ α ρ 2 , x > 0 , α > 0 , λ > 0 , 0 ρ < 1
When α = 0, G ¯ x , α , β , θ , ρ and g x , reduce to corresponding the MOEL distribution Ghitany et al. [14].
Also, when ρ = 0, α = 0, the G ¯ x , α , β , θ , ρ and g x , reduce to the Lomax distribution Lomax [26]. Figure 1 and Figure 2 give the graph MB-L PDF for values of, α ,   β ,   θ and ρ .
Figure 1 and Figure 2 show different shapes of the PDF while it gives, a monotonic increasing, decreasing, constant and unimodal shapes, so we can conclude that the MB-L model is a very flexible distribution in modeling various type of data.
The next theorem gives the behavior of the MB-L PDF.
Theorem 1. 
For the MB-L ( α , β , θ , ρ ) model, The PDF given by (4) is decreasing if 0 < θ < α ( 1 ρ ) ( 1 + α ρ ) 2 + α 2 ρ 1 ρ + α ( 1 + ρ ) , independent of β and is unimodal if θ > α 2 ( ρ 1 ) 2 + α 2 ρ 1 α ( 1 + ρ ) .
Proof. 
We can rewrite Equation (4) as:
g x , α , β , θ , ρ = α   β   θ ( 1 + x   β 2 θ 1 + α + 2 1 + x   β θ ρ + α   ρ 2 ) 1 + x   β θ + 1 ( 1 + α + 1 + x   β θ α   ρ ) 2
then,
g ` x = α β 2 θ ( 1 + x β ) 2 ( 1 + α ( 1 ρ ) + ( 1 + x β ) θ ) 3 Φ x , x > 0
where,
Φ x = α β 2 θ ( 2 θ ( ( 1 + x β ) θ ρ ) ( 1 + α + ( 1 + x β ) θ α ρ ) 2 θ ( ( 1 + x β ) 2 θ ( 1 + α + 2 ( 1 + x β ) θ ) ρ + α ρ 2 ) ( 1 + x β ) θ ( 1 + θ ) ( 1 + α + ( 1 + x β ) θ α ρ ) ( ( 1 + x β ) 2 θ ( 1 + α + 2 ( 1 + x β ) θ ) ρ + α ρ 2 ) )
if Φ 0 = α β 2 θ ( ρ 1 ) ( 2 θ + α 2 1 + θ 1 ρ ρ + α ( 1 + θ ( 1 + ρ ) + ρ ) ) 0 , then θ α ( 1 ρ ) ( 1 + α ρ ) 2 + α 2 ρ 1 ρ + α ( 1 + ρ ) , g ` x < 0 , then g ( x ) is decreasing.
if Φ 0 > 0 , then θ < α ( 1 ρ ) ( 1 + α ρ ) 2 + α 2 ρ 1 ρ + α ( 1 + ρ ) ,   lim x g ( x ) = 0 , and lim x 0 g x = β θ ( 1 ( 1 + α ) ρ + α ρ 2 ) α ( 1 ρ ) 2 then g(x), first increases than decrease to zero and hence has a mode x m o d given by: Φ x m o d = 0 . Moreover, this mode is unique Dharmadhikari et al. [27]. □
Remark 1.
  • For  ρ = 0 , α = 1 ,  g(x) is decreasing (unimodal) if  θ 1 ( θ > 1 which is the well-known result for the Lomax distribution.
  • For ρ = 0 , g(x) is decreasing if  θ < α 2 α (i.e., α + ( 2 α ) θ 0 ) and is unimodal if  θ > α 2 α (i.e., α + ( 2 α ) θ < 0 ) which is the well-known result for the MOEL distribution (Ghitany et al. [14]).
The  r t h  moment of MB-L model
For the MB-L ( α , β , θ , ρ ) model, the r t h moment E X r , r 1 , is:
E X r = r 0 x r 1 G ¯ x d x = r 0 x r 1 α 1 + β x θ 1 ρ 1 + β x θ 1 1 1 ρ α 1 + β x θ d x = r β r + 2 θ u = 0 1 1 ρ α u 0 1 y u 1 θ y 1 θ 1 r 1 ρ y d y = r β r + 2 θ u = 0 1 1 ρ α u ( 1 ) r θ   Gamma 1 + r
The MGF of MB-L model
We present the MGF M (t, α , β , θ , ρ ) of the MB-L model. Using Equation (4), the substitution u = ( 1 + β x ) and Maclaurin expansion of e x for all x, we get the following:
M t , α , β , θ , ρ = E e t X = 0 e t x α   β   θ ( 1   +   x   β 2 θ       1   +   α   +   2 1   +   x   β θ ρ   +   α   ρ 2 ) 1   +   x β θ   +   1 (   1   +   α   +   1   +   x   β θ     α   ρ ) 2 d x = 1 e t β ( u     1 ) α   θ ( u 2 θ       1   +   α   +   2 u θ ρ   +   α   ρ 2 ) u θ   +   1 (   1   +   α   +   u θ     α   ρ ) 2 d u = α   θ   e t β 1 e t u β α   θ ( u 2 θ       1   +   α   +   2 u θ ρ   +   α   ρ 2 ) u θ   +   1 (   1   +   α   +   u θ     α   ρ ) 2 = α   θ   e t β i   =   0 j   =   0 ( α   ρ ) j t β α   ρ θ j   +   2     i   +   1 θ j   +   1     i     1 i ! j ! .
where 0 < α   ρ < 1 . For more details of MGF see BS [28] and YF and SY [29].
Now we have the following results:
Theorem 2. 
For the SF (3) of MB-L model, if  X i , i = 1 , 2 , , n  are independent identically distributed (i.i.d), then  U N = min X 1 , X 2 , , X n ,  has the SF:
G ¯ U N x , α , β , θ , ρ = α 1 ρ   ( 1 + β x ) θ ( 1 + β x ) θ 1 1 ρ   α n .
Proof. 
We know that,
G ¯ x , α , β , θ , ρ = α   ( 1 + β x ) θ 1 ρ   ( 1 + β x ) θ 1 1 1 ρ   α   ( 1 + β x ) θ
G ¯ U N x , α , β , θ , ρ = P min X 1 , X 2 , , X n > x = i = 1 n G ¯ x , α , β , θ , ρ = α 1     ρ ( 1   +   β x ) θ ( 1   +   β x ) θ     1     1     ρ α n .
Theorem 3. 
For the SF (3) of MB-L model, if  X i , i = 1 , 2 , , N  are i.i.d,  N  be a Markov Bernoulli geometric distribution with parameters  p , ρ  such that:
P N = n = p , n = 1 1 p a 1 1 a 1 n 2 , n 2 ,
which is independent of X i for all I = 1,2, …, N. Then, U N = min 1 i N X i , is distributed MB-L if X i is distributed as Lomax distribution.
Proof. 
Suppose that,
G ¯ x = Pr U N > x = P r X i > x i , i = 1 , 2 , , N = n = 1 F ¯ n x Pr N = n = p F ¯ x + 1 p a 1 F ¯ 2 x n = 2 F ¯ n 2 x ( 1 a 1 ) n 2 = p F ¯ x + 1 p a 1 F ¯ 2 x 1 [ 1     ( 1     a 1 ) F ¯ x ] = F ¯ x [ p     p F ¯ x   +   a 1 F ¯ x ] [ 1     ( 1     a 1 ) F ¯ x ] ,   w h e r e   a 1 = 1 ρ p
hence,
G ¯ x = p [ 1 ρ 1 + x β θ ] 1 + x β θ ( 1 a 1 )
Which is MB-L mdel with p = α . □

3. The HRF of the MB-L Model

From Equation (3) we have:
h x = β θ ( ( 1 + x β ) 2 θ ( 1 + α + 2 ( 1 + x β ) θ ) ρ + α ρ 2 ) ( 1 + x β ) ( ( 1 + x β ) θ ρ ) ( 1 + α + ( 1 + x β ) θ α ρ ) .
For all θ , α > 0 , 0 < ρ < 1 , then,
lim x 0 h x = β   θ ( 1 + α   ρ ) α ( 1 + ρ ) ,   lim x h x = 0 .
Figure 3 and Figure 4 show different shapes of the HRF while it gives, a monotonic increasing, decreasing, constant and unimodal shapes, so we can conclude that the HRF is a very flexible distribution in modeling various type of data.
Now we will study the behavior of HRF according to the following theorem:
Theorem 4.
For the MB-L ( α , β , θ , ρ ) model, The HRF (3) is decreasing (unimodal) if  θ α 1 ρ α ρ 1 1 α 2 ρ α 1 ρ θ < α 1 ρ ( α ρ 1 ) ( 1 α 2 ρ α ( 1 ρ ) )  independent of  β .
Proof. 
From Equation (3) we have:
h ` x = β 2   θ   ( 1 + x β ) θ ( 1 + x β 2 1 + x β θ ρ 2 1 + α + 1 + x β θ α   ρ 2 ) ) ψ ( x )
where, ψ ( x ) = θ ( ( 1 + x β ) 2 θ ρ ) ( 1 + ρ ) α 2 ( 1 + θ ρ ( 1 + x β ) θ ) ( 1 + ρ ) 2 ρ +   ( ( 1 + x β ) θ + 1 ) ( ( 1 + x β ) 3 θ + ( ( 1 + x β ) θ 3 ( 1 + x β ) 2 θ ) ρ + ( 1 + 2 ( 1 + x β ) θ ) ρ 2 ) +   α 1 ( 1 + ρ ) ( ( 1 + x β ) 2 θ + 2 ( 1 + 2 ( 1 + x β ) θ ) ρ + ( 2 ( 1 + x β ) θ 3 ) ρ 2 + θ ( ( 1 + x β ) 2 θ   2 ρ + ρ 2 ) ) ]
The proof is as in the Theorem 1. □
Remark 2.
  • For  ρ = 0 , α = 1 ,   θ 1 ( θ < 1 ) h(x) is decreasing (unimodal). It is the same result for the Lomax distribution.
  • For  ρ = 0 ,  h(x) is decreasing if  θ α 1 α (i.e.,  α + ( 1 α ) θ 0 ) and is unimodal if  θ < α 1 α (i.e., α + ( 1 α ) θ < 0 ) which is the well-known result for the MOEL distribution.

4. Estimation of MB-L Parameters and Asymptotic Confidence Intervals (CI)

Here, the MLE for the MB-L Parameters are developed. Asymptotic confidence intervals of ψ ^ = ( α ^ , β ^ , θ ^ , ρ ^ ) are obtained using the inverse Fisher’s information matrix elements. Simulation studies are carried out to investigate the accuracy of the estimates of the model’s parameter.
Suppose that ( t 1 , δ 1 ), ( t 2 , δ 2 ), …, ( t n , δ n ) be a random sample from the MB-L ( α , β , θ , ρ ) model, where δ i = 0 or δ i = 1 if t i are censored or complete observations, respectively.
The log-likelihood function for the MB-L ( α , β , θ , ρ ) model is:
l n t i , α , β , θ , ρ = i = 1 n { δ i log α   β   θ ( 1 + x   β 2 θ 1 + α ( 1 + ρ 2 ) + 2 1 + x   β θ ρ ) 1 + x β θ + 1 ( 1 + α ( 1 ρ ) + 1 + x   β θ ) 2 + 1 δ i log α ( 1 + β x ) θ 1 ρ ( 1 + β x ) θ 1 1 1 ρ α ( 1 + β x ) θ }
The first derivative of l n t i , α , β , θ , ρ with respect to α ,   β ,   θ and ρ , respectively, are given by
l n α = 0 , l n β = 0 , l n θ = 0 , l n ρ = 0 .
The MLE ψ ^ = ( α ^ , β ^ , θ ^ , ρ ^ ) can be obtained numerically using these equations. For the asymptotic CI, the normal approximation of the MLE can be used to construct asymptotic CIs for the parameters ψ when the sample size is large enough. A two-sided (1−α) 100% CIs for ψ are ( ψ ^ ± Z α / 2 V a r ( ψ ^ ) ) , where V a r ( ψ ^ ) are the asymptotic variances of ψ ^ .
To compare the MB-L model with MOEL and Lomax distributions, we used the likelihood ratio test (LRT) as:
First: the null hypothesis H 10 : α = 1 , ρ = 0 (Lomax distribution). Under H 10 the likelihood ratio statistic: A 1 = 2 ln ( 1 , β ^ , θ ^ , 0 ) ln ( α ^ , β ^ , θ ^ , ρ ^ ) , it has a χ 2 distribution with 2 degrees of freedom
Second: the null hypothesis H 20 : ρ = 0 (MOEL distribution). Under H 20 the likelihood ratio statistic: A 2 = 2 ln ( α ^ , β ^ , θ ^ , 0 ) ln ( α ^ , β ^ , θ ^ , ρ ^ ) , which has an asymptotic χ 2 distribution with 1 degree of freedom.
Also, the model selection: Akiake information criterion (AIC) (Akiake [30]), Bayesian information criterion (BIC) and Consistent Akaike Information Criteria C A I C defined as:
A I C = log l i k e l i h o o d 2 k ,
B I C = l o g l i k e l i h o o d k 2 log n ,
C A I C   =   2   l o g   l i k e l i h o o d 2 k n n k 1 .
where, k is the number of model parameters and n is the sample size. The model with higher AIC, CAIC and BIC is the one that better fits the data.

5. Simulation

The calculation of the estimation is based on N = 10,000 simulated samples from the MB-L model. The sample sizes are 50, 100, 200 and 300 and the parameter values are ψ = α , β , θ , ρ = ( 0.7 , 1.2 , 0.05 , 0.03 ) and ( 0.3 , 0.3 , 0.1 , 0.2 ) . The validity of the estimate of ψ is studied by the following measures:
  • Bias of ψ ( α , β , θ , ρ ) is 1 N i = 1 N ( ψ ^ ψ )
  • Mean square error (MSE) of ψ ( α , β , θ , ρ ) is 1 N i = 1 N ( ψ ^ ψ ) 2 .
  • Coverage probability (CP) of the N simulated confidence intervals.
When the n is large, the values of ψ ^ are close to the initial values of ψ see Table 1.

6. Applications

6.1. Censored Data

Lee and Wang [31] P. 231 obtained the data which represent the 137-bladder cancer patient. This data has completed at 0.08 to 79.05 months and censored at 0.87, 3.02, 4.33, 4.65, 4.70, 8.60, 10.86, 19.36, 24.80 months. Table 2 shows values of MLE, Log-likelihood, AIC, BIC and CAIC of MB-L model with other models.
We note that, AIC, BIC and CAIC of MB-L model more than the corresponding of the MOEL and Lomax distributions which means that MB-L model is better to fit for the given data. Moreover, the approximate 95% two-sided CI of the parameters α , β , θ and ρ are given respectively as [ 2.075 , 3.88 ], [ 1.275 , 4.838 ], [ 0.147 , 0.342 ] and [ 0.006 , 0.179 ].
For the given censored data, under H 10 thus X L 1 = 2 [ 422.348 + 417.228 ] = 10.24 , then X L 1 > χ 2 , 0.05 2 = 5.991 . Also, under H 20 thus X 6 L 2 = 2 [ 419.842 + 417.228 ) ] = 5.228 , then X L 2 > χ 1 , 0.05 2 = 3.84 . So, the LRT rejects the null hypothesis that the Lomax and MOEL models is proper for the specific data.
Now, suppose that t i and δ i ,   i = 1 , 2 , , n are, respectively, the ordered survival times and corresponding censoring indicators. The product-limit estimator or Kaplan-Meier estimator (KME) (Kaplan-Meier [32]) of a SF is:
G ¯ n t = t : t i t i = 1 , , n 1 δ i n i + 1 , t > 0 .
Figure 5, Figure 6 and Figure 7 show the probability-probability (P-P) plot of the KME versus the fitted Lomax, MOEL and MB-L SFs for 137 censored data.
From the previous figures, we notice that the drawn points for the fitted MB-L SF are close to the 45° line, indicating good fit as comparing with the fitted MB-L, MOEL and Lomax SFs.
Since α = 2.9782 , θ = 3.0578 , ρ = 0.2449 and β = 0.0927 , then the estimated hazard rate function ((a) MB-L model, (b) Lomax distribution) is as shown in the Figure 8.

6.2. COVID-19 Data

Applying the data from Huda M. Alshanbari et al. [33], which shows the two complete data of COVID-19 which represent a mortality rate.
First: COVID-19 data obtained through 37 days, from 27 June to 2 August 2021 (Saudi Arabia). The data and its measures show in Table 3 and Table 4.
From this table, AIC, BIC and CAIC of MB-L model more than the corresponding of the MOEL and Lomax distributions which means that MB-L model is better to fit for the given data. Moreover, the approximate 95% two-sided CI of the parameters α , β , θ and ρ are given respectively as [ 0.0049 , 0.0038 ], [ 0.068 , 0.563 ], [ 0.868 , 7.358 ] and [ 0.068 , 0.307 ].
For the given 37 COVID-19 data, under H 10 thus X L 1 = 2 [ 75.905 81.646 ] = 11.482 , then X L 1 > χ 2 , 0.05 2 = 5.991 . Also, under H 20 thus X L 2 = 2 [ 79.529 81.646 ] = 4.234 , then X L 2 > χ 1 , 0.05 2 = 3.84 . So, the LRT rejects the null hypothesis that the Lomax and MOEL models is proper for the specific data. The estimated hazard rate function ((a) MB-L model, (b) Lomax distribution) is as shown in the Figure 9.
Second: COVID-19 data obtained during 172 days from the first of 1 March to 20 August 2020 (Italy). The data and its measures show in Table 5 and Table 6.
From Table 6, AIC, BIC and CAIC of MB-L model more than the corresponding of the MOEL and Lomax models which means that MB-L Distribution is better to fit for the given data. In addition to, the approximate 95% two-sided CI of the parameters α , β , θ and ρ are given respectively as [ 0.219 , 0.993 ], [ 0.188 , 0.776 ], [ 3.203 , 11.916 ] and [ 0.111 , 0.316 ].
For the given 172 COVID-19 data, under H 10 thus X L 1 = 2 73.756 154.215 = 160.918 , then X L 1 > χ 2 , 0.05 2 = 5.991 . Also, under H 20 thus X L 2 = 2 [ 134.216 154.215 ] = 39.998 , then X L 1 > χ 1 , 0.05 2 = 3.84 . So, the LRT rejects the null hypothesis that the Lomax and MOEL models is proper for the specific data. The estimated hazard rate function ((a) MB-L model, (b) Lomax distribution) is as shown in the Figure 10.

7. Conclusions

In this paper, we introduced a four-parameter continuous distribution that generalizes MOEL and Lomax distributions. The new model is referred to as MB-L distribution, the derived properties including PDF, HRF, moments, MGF and minimum (maximum) MBG stable. The MLE procedure is straightforward. The active fitting of MB-L distribution is shown on bladder cancer and COVID-19 applications. We note that the value of the selected model choices is higher for MB-L distribution than for the MOEL and Lomax distributions. Furthermore, the relationship between the empirical and fitted SFs for MB-L distribution is higher compared to MOEL and Lomax distributions. All previous results indicate the advantage of MB-L distribution for bladder cancer and COVID-19 data.

Author Contributions

Conceptualization, B.I.M.; methodology, B.I.M.; software, B.I.M.; validation, Y.A.T., M.E.B. and M.M.H.; investigation, Y.A.T., M.E.B. and M.M.A.E.-R.; resources, B.I.M., Y.A.T., M.E.B., M.M.A.E.-R. and M.M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research project was supported by the Researchers Supporting Project Number (RSP2023R488), King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

The data was mentioned along the paper.

Conflicts of Interest

The authors declare there is no conflict of interest.

References

  1. Gharib, M.; Mohammed, B.I.; Al-Ajmi, K.A.H. A New Method for Adding Two Parameters to a Family of Distributions with Application. J. Stat. Appl. Pro. 2017, 6, 487–497. [Google Scholar] [CrossRef]
  2. Marshall, A.W.; Olkin, I. A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 1997, 84, 641–652. [Google Scholar] [CrossRef]
  3. Johnson, N.L.; Kotz, S.; Balakrishnan, N. Contineous Univariate Distributions, 2nd ed.; Wiley: New York, NY, USA, 1994; Volume 1. [Google Scholar]
  4. Harris, C.M. The Pareto Distribution as a Queue Service Discipline. Oper. Res. 1968, 16, 307–313. [Google Scholar] [CrossRef]
  5. Bryson, M. Heavy tailed distributions: Properties and tests. Technometrics 1974, 16, 61–68. [Google Scholar] [CrossRef]
  6. Cordeiro, G.M.; Ortega, E.M.; Popović, B.V. The gamma-Lomax distribution. J. Stat. Comput. Simul. 2013, 85, 305–319. [Google Scholar] [CrossRef]
  7. Devi, B.; Kumar, P.; Kour, K. Entropy of Lomax Probability Distribution and its Order Statistic. Int. J. Stat. Syst. 2017, 12, 175–181. [Google Scholar]
  8. Maurya, R.K.; Tripathi, Y.M.; Lodhi, C.; Rastogi, M.K. On a generalized Lomax distribution. Int. J. Syst. Assur. Eng. Manag. 2019, 10, 1091–1104. [Google Scholar] [CrossRef]
  9. Mohammed, B.I.; Abu-Youssef, S.E.; Sief, M.G. A New Class with Decreasing Failure Rate Based on Countable Mixture and Its Application to Censored Data. J. Test. Eval. 2019, 48, 273–288. [Google Scholar] [CrossRef]
  10. Haq, M.A.U.; Rao, G.S.; Albassam, M.; Aslam, M. Marshall–Olkin Power Lomax distribution for modeling of wind speed data. Energy Rep. 2020, 6, 1118–1123. [Google Scholar] [CrossRef]
  11. Aryuyuen, S.; Bodhisuwan, W. The Type II Topp Leone-Power Lomax Distribution with Analysis in Lifetime Data. J. Stat. Theory Pract. 2020, 14, 31. [Google Scholar] [CrossRef]
  12. Alshanbari, H.M.; Ijaz, M.; Asim, S.M.; Hosni El-Bagoury, A.A.-A.; Dar, J.G. New Weighted Lomax (NWL) Distribution with Applications to Real and Simulated Data. Math. Probl. Eng. 2021, 2021, 8558118. [Google Scholar] [CrossRef]
  13. Dey, S.; Altun, E.; Kumar, D.; Ghosh, I. The Reflected-Shifted-Truncated Lomax Distribution: Associated Inference with Applications. Ann. Data Sci. 2021, 2021, 1. [Google Scholar] [CrossRef]
  14. Ghitany, M.E.; Al-Awadhi, F.A.; Alkhalfan, L.A. Marshall–Olkin Extended Lomax Distribution and Its Application to Censored Data. Commun. Stat. Theory Methods 2007, 36, 1855–1866. [Google Scholar] [CrossRef]
  15. Lemonte, A.J.; Cordeiro, G.M. An extended Lomax distribution. Statistics 2013, 47, 800–816. [Google Scholar] [CrossRef]
  16. Cordeiro, G.; Alizadeh, M.; Marinho, P. The type I half-logistic family of distributions. J. Stat. Comput. Simul. 2016, 86, 707–728. [Google Scholar] [CrossRef]
  17. Al-Zahrani, B.; Sagor, H. Statistical analysis of the Lomax–Logarithmic distribution. J. Stat. Comput. Simul. 2014, 85, 1883–1901. [Google Scholar] [CrossRef]
  18. Al-Zahrani, B.; Sagor, H. The Poisson-Lomax Distribution. Rev. Colomb. Estadística 2014, 37, 225–245. [Google Scholar] [CrossRef]
  19. Tahir, M.H.; Cordeiro, G.M.; Mansoor, M.; Zubair, M. Weibull-Lomax distribution: Properties and applications. Hacet. J. Math. Stat. 2015, 44, 455–474. [Google Scholar] [CrossRef]
  20. El-Bassiouny, A.; Abdo, N.F.; Shahen, H.S. Exponential lomax distribution. Int. J. Comput. Appl. 2015, 121, 24–29. [Google Scholar]
  21. Rady, E.-H.A.; Hassanein, W.A.; Elhaddad, T.A. The power Lomax distribution with an application to bladder cancer data. Springerplus 2016, 5, 1838. [Google Scholar] [CrossRef]
  22. Cooray, K. Analyzing lifetime data with long-tailed skewed distribution: The logistic-sinh family. Stat. Model. 2005, 5, 343–358. [Google Scholar] [CrossRef]
  23. Abu El Azm, W.S.; Almetwally, E.M.; Al-Aziz, S.N.; El-Bagoury, A.A.-A.H.; Alharbi, R.; Abo-Kasem, O.E. A New Transmuted Generalized Lomax Distribution: Properties and Applications to COVID-19 Data. Comput. Intell. Neurosci. 2021, 2021, 5918511. [Google Scholar] [CrossRef]
  24. Alsuhabi, H.; Alkhairy, I.; Almetwally, E.M.; Almongy, H.M.; Gemeay, A.M.; Hafez, E.; Aldallal, R.; Sabry, M. A superior extension for the Lomax distribution with application to Covid-19 infections real data. Alex. Eng. J. 2022, 61, 11077–11090. [Google Scholar] [CrossRef]
  25. Ogunde, A.A.; Chukwu, A.U.; Oseghale, I.O. The Kumaraswamy Generalized Inverse Lomax distribution and applications to reliability and survival data. Sci. Afr. 2023, 19. [Google Scholar] [CrossRef]
  26. Lomax, K.S. Business failures: Another example of the analysis of failure data. J. Am. Stat. Assoc. 1954, 45, 21–29. [Google Scholar] [CrossRef]
  27. Dharmadhikari, S.; Joag-dev, K. Unimodality, Convexity, and Applications; Academic Press: Cambridge, MA, USA, 1998. [Google Scholar]
  28. Simsek, B. Formulas Derived from Moment Generating Functions And Bernstein Polynomials, Applicable Analysis and Discrete Mathematics. Appl. Anal. Discret. Math. 2019, 13, 839–848. [Google Scholar] [CrossRef]
  29. Yalcin, F.; Simsek, Y. Formulas for characteristic function and moment generating functions of beta type distribution. Rev. Real Acad. Cienc. Exactas Físicas Y Naturales. Ser. A Matemáticas 2022, 116, 86. [Google Scholar] [CrossRef]
  30. Akaike, H. Fitting Autoregressive Models for Prediction. Ann. Inst. Stat. Math. 1969, 21, 243–247. [Google Scholar] [CrossRef]
  31. Lee, E.T.; Wang, J.W. Statistical Methods for Survival Data Analysis, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2003; Volume 476. [Google Scholar]
  32. Kaplan, E.L.; Meier, P. Nonparametric Estimation from Incomplete Observations. J. Am. Stat. Assoc. 1958, 53, 457–481. [Google Scholar] [CrossRef]
  33. Alshanbari, H.M.; Odhah, O.H.; Almetwally, E.M.; Hussam, E.; Kilai, M.; El-Bagoury, A.-A.H. Novel Type I Half Logistic Burr-Weibull Distribution: Application to COVID-19 Data. Comput. Math. Methods Med. 2022, 2022, 1444859. [Google Scholar] [CrossRef]
Figure 1. Draw of decreasing MB-L PDF for values of α ,   β ,   θ and ρ .
Figure 1. Draw of decreasing MB-L PDF for values of α ,   β ,   θ and ρ .
Axioms 12 00439 g001
Figure 2. Draw of increasing-decreasing MB-L PDF for values of α ,   β ,   θ and ρ .
Figure 2. Draw of increasing-decreasing MB-L PDF for values of α ,   β ,   θ and ρ .
Axioms 12 00439 g002
Figure 3. Draw of decreasing HRF for values of α ,   β ,   θ and ρ .
Figure 3. Draw of decreasing HRF for values of α ,   β ,   θ and ρ .
Axioms 12 00439 g003
Figure 4. Draw of increasing-decreasing MB-L HRF for values of α ,   β ,   θ and ρ .
Figure 4. Draw of increasing-decreasing MB-L HRF for values of α ,   β ,   θ and ρ .
Axioms 12 00439 g004
Figure 5. P-P plot of KME and fitted Lomax SF.
Figure 5. P-P plot of KME and fitted Lomax SF.
Axioms 12 00439 g005
Figure 6. P-P plot of KME and fitted MB-L SF.
Figure 6. P-P plot of KME and fitted MB-L SF.
Axioms 12 00439 g006
Figure 7. P-P plot of KME and fitted MOEL SF.
Figure 7. P-P plot of KME and fitted MOEL SF.
Axioms 12 00439 g007
Figure 8. The estimated hazard rate function of MB-L model from the given censored data.
Figure 8. The estimated hazard rate function of MB-L model from the given censored data.
Axioms 12 00439 g008
Figure 9. The estimated hazard rate function of MB-L model based on COVID-19 data (Saudi Arabia).
Figure 9. The estimated hazard rate function of MB-L model based on COVID-19 data (Saudi Arabia).
Axioms 12 00439 g009
Figure 10. The estimated hazard rate function of MB-L model based on COVID-19 data (Italy).
Figure 10. The estimated hazard rate function of MB-L model based on COVID-19 data (Italy).
Axioms 12 00439 g010
Table 1. The MLE, bias, MSE and CP values for the MB-L ( α , β , θ , ρ ) model.
Table 1. The MLE, bias, MSE and CP values for the MB-L ( α , β , θ , ρ ) model.
nParameterInitialMLEBiasMSECPInitialMLEBiasMSECP
50α0.70.70700.00700.00310.98910.30.30410.004110.00080.9782
β1.21.1934−0.00660.00090.95350.30.2981−0.00190.00020.9643
θ0.050.63900.58900.34830.98520.10.43270.33270.11140.9665
ρ0.030.03370.00390.00070.95520.20.23310.03310.02070.9663
100α0.70.70080.00080.00090.98850.30.30470.00470.00040.9757
β1.21.2030.00310.00040.95350.30.2975−0.00250.00010.9527
θ0.050.64150.59140.35060.97640.10.43470.33470.11240.9791
ρ0.030.0296−0.00040.00030.95740.20.22680.02680.00970.9546
200α0.70.70450.00450.00090.99990.30.30080.00080.00010.9887
β1.21.1991−0.00090.00020.95730.30.2997−0.00030.000050.9592
θ0.050.63680.58680.34470.98630.10.43760.33760.11420.9854
ρ0.030.03220.00220.00030.95830.20.20520.00520.00380.9575
300α0.70.70480.00480.00060.99920.30.2997−0.00030.00010.9985
β1.21.1976−0.00240.00010.95920.30.30020.00020.000030.9587
θ0.050.63630.58630.34400.98450.10.43680.33680.113540.9863
ρ0.030.03260.00260.00020.95470.20.20290.00290.00300.9638
Table 2. MLE, standrd error (S.E), log-likelihood, AIC, BIC for MB-L, MOEL and Lomax models from the 137- censored data.
Table 2. MLE, standrd error (S.E), log-likelihood, AIC, BIC for MB-L, MOEL and Lomax models from the 137- censored data.
ModelParameterMLES.ELog-LikelihoodAICBICCAIC
MB-Lα2.97820.433 417.228 426.229 428.069 424.467
θ3.05780.231
ρ0.24490.031
β0.09270.005
MOELα2.9590.087−419.842 426.287 428.127 426.590
θ4.1150.012
β0.0580.004
Lomaxθ3.7330.036 422.348 430.347 432.188 430.651
β0.03250.006
Table 3. 37- COVID-19 data.
Table 3. 37- COVID-19 data.
37 COVID-19 Data
0.01950.02130.02140.02170.02310.02330.02350.0235
0.02380.02390.02450.02600.02640.02680.02700.0271
0.02750.02780.02780.02820.02820.02850.02870.0294
0.02960.03000.03010.03090.03100.03130.03140.0315
0.03240.03250.03280.03320.0358
Table 4. MLE, S.E, log-likelihood, AIC, CAIC and BIC for MB-L, MOEL and Lomax models from COVID-19 data (Saudi Arabia).
Table 4. MLE, S.E, log-likelihood, AIC, CAIC and BIC for MB-L, MOEL and Lomax models from COVID-19 data (Saudi Arabia).
ModelParameterMLES.ELog-LikelihoodAICBICCAIC
MBELα0.01030.033 81.646 73.646 74.425 72.396
θ4.11200.131
ρ0.40390.012
β0.05550.006
MOELα0.10430.092 79.529 73.529 74.113 72.802
θ0.49010.211
β9.2670.005
Lomaxθ0.82130.321 75.905 71.905 72.294 71.552
β86.0390.432
Table 5. 172- COVID-19 data.
Table 5. 172- COVID-19 data.
172 COVID-19 Data
0.01070.04900.06010.04600.0533 0.06300.02970.0885
0.05400.17200.08470.07130.09890.04950.10250.1079
0.09840.11240.08070.10440.12120.11670.12550.1416
0.13150.10730.16290.14850.14530.20000.20700.1520
0.16280.16660.14170.12210.17670.19870.14080.1456
0.14430.13190.10530.17890.20320.21670.13870.1646
0.13750.14210.20120.19570.12970.17540.13900.1761
0.11190.19150.18270.15480.15220.13690.24950.1253
0.15970.21950.25550.19560.18310.17910.20570.2406
0.12270.21960.26410.30670.17490.21480.21950.1993
0.24210.24300.19940.17790.09420.30670.19650.2003
0.11800.16860.26680.21130.33710.17300.22120.4972
0.16410.26670.26900.23210.27920.35150.13980.3436
0.22540.13020.08640.16190.13110.19940.31760.1856
0.10710.10410.15930.05370.11490.11760.04570.1264
0.04760.16200.11540.14930.06730.08940.03650.0385
0.21900.07770.05610.04350.03720.03850.07690.1491
0.08020.08700.04760.05620.01380.06840.11720.0321
0.03270.01980.01820.01970.02980.05450.02080.0079
0.02370.01690.03360.07550.02630.02600.01500.0054
0.03750.00430.01540.01460.02100.01150.00520.2512
0.00840.01250.01250.01090.0071
Table 6. MLE, log-likelihood, AIC, CAIC and BIC for MB-L, MOEL and Lomax models from the COVID-19 data (Italy).
Table 6. MLE, log-likelihood, AIC, CAIC and BIC for MB-L, MOEL and Lomax models from the COVID-19 data (Italy).
ModelParameterMLES.ELog-LikelihoodAICBICCAIC
MB-Lα0.38770.221 154.215 146.215 143.920 145.976
θ7.56060.431
ρ0.10280.023
β0.48240.321
MOELα0.00130.006 134.216 128.216 126.495 128.073
θ1.00210.453
β0.01360.051
Lomaxθ0.47510.254 73.756 69.756 68.608 69.6845
β71.7020.534
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mohammed, B.I.; Tashkandy, Y.A.; El-Raouf, M.M.A.; Hossain, M.M.; Bakr, M.E. The Markov Bernoulli Lomax with Applications Censored and COVID-19 Drought Mortality Rate Data. Axioms 2023, 12, 439. https://doi.org/10.3390/axioms12050439

AMA Style

Mohammed BI, Tashkandy YA, El-Raouf MMA, Hossain MM, Bakr ME. The Markov Bernoulli Lomax with Applications Censored and COVID-19 Drought Mortality Rate Data. Axioms. 2023; 12(5):439. https://doi.org/10.3390/axioms12050439

Chicago/Turabian Style

Mohammed, Bahady I., Yusra A. Tashkandy, Mohmoud M. Abd El-Raouf, Md. Moyazzem Hossain, and Mahmoud E. Bakr. 2023. "The Markov Bernoulli Lomax with Applications Censored and COVID-19 Drought Mortality Rate Data" Axioms 12, no. 5: 439. https://doi.org/10.3390/axioms12050439

APA Style

Mohammed, B. I., Tashkandy, Y. A., El-Raouf, M. M. A., Hossain, M. M., & Bakr, M. E. (2023). The Markov Bernoulli Lomax with Applications Censored and COVID-19 Drought Mortality Rate Data. Axioms, 12(5), 439. https://doi.org/10.3390/axioms12050439

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop