Next Article in Journal
An Interactive Consensus Model in Group Decision Making with Heterogeneous Hesitant Preference Relations
Next Article in Special Issue
Some Probabilistic Generalizations of the Cheney–Sharma and Bernstein Approximation Operators
Previous Article in Journal
Research Agenda on Multiple-Criteria Decision-Making: New Academic Debates in Business and Management
Previous Article in Special Issue
An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mixture of Akash Distributions: Estimation, Simulation and Application

by
Anum Shafiq
1,
Tabassum Naz Sindhu
2,
Showkat Ahmad Lone
3,
Marwa K. H. Hassan
4 and
Kamsing Nonlaopon
5,*
1
School of Mathematics and Statistics, Nanjing University of Information Science and Technology, Nanjing 210044, China
2
Department of Statistics, Quaid-i-Azam University, Islamabad 44000, Pakistan
3
Department of Basic Sciences College of Science and Theoretical Studies, Saudi Electronic University, Riyadh 11673, Saudi Arabia
4
Department of Mathematics, Faculty of Education, Ain Shams University, Cairo 11566, Egypt
5
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(10), 516; https://doi.org/10.3390/axioms11100516
Submission received: 4 September 2022 / Revised: 20 September 2022 / Accepted: 21 September 2022 / Published: 29 September 2022

Abstract

:
In this paper, we propose a two-component mixture of Akash model (TC-MAM). The behavior of TC-MAM distribution has been presented graphically. Moment-based measures, including skewness, index of dispersion, kurtosis, and coefficient of variation, have been determined and hazard rate functions are presented graphically. The probability generating function, Mills ratio, characteristic function, cumulants, mean time to failure, and factorial moment generating function are all statistical aspects of the mixed model that we explore. Furthermore, we figure out the relevant parameters of the mixture model using the most suitable methods, such as least square, weighted least square, and maximum likelihood mechanisms. Findings of simulation experiments to examine behavior of these estimates are graphically presented. Finally, a set of data taken from the real world is examined in order to demonstrate the new model’s practical perspectives. All of the metrics evaluated favor the new model and the superiority of proposed distribution over mixture of Lindley, Shanker, and exponential distributions.

1. Introduction

In most reliability scenarios, data are modelled using a single parametric model. However, in certain circumstances, a population can be split into many subgroups, each showing a particular category of collapse. Finite mixture models serve a significant role in modelling such diverse data. Biology, business, engineering, healthcare, genetics, marketing, real-world applications, and social sciences all benefit from finite mixture models. Mixture models are created by varying the proportions of two or more models to generate a new distribution with novel properties. Consequently, it is essential to examine the statistical characteristics of the suggested mixture model and employ the suitable methods for estimating the unexplained parameters. Mixture models are used in a diversity of applications, such as clustering and classification [1,2,3,4]. Sultan et al. [5] proposed a mix of inverse Weibull models and utilized density and hazard function graphs to study some of its features. The conventional characteristics of the concoction of Burr XII and Weibull distributions were examined by [6]. Recently, the authors [7,8], Ateya [9], Mohammadi et al. [10], and Al-Moisheer et al. [11] are among the scientists who study mixture modelling in a variety of contexts.
Many applied sciences, including medical, engineering, insurance, and finance, rely on lifetime data modelling and analysis. Some of the continuous distributions used to explain lifetime data are Weibull, gamma, exponential, lognormal, and Lindley, as well as their generalizations. Since many investigators have employed Lindley distribution to predict lifetime data, and Hussain [12] has demonstrated that Lindley model is effective for stress-strength dependability modelling, Lindley model may not be suitable for describing real-world data in many cases. Shanker [13] developed a novel model by using a two-component concoction of an exponential model ϑ and a gamma model 3 , ϑ to have a unique distribution that is more flexible than Lindley and exponential distribution for modelling lifespan data in terms of dependability and hazard rate shapes. Shanker et al. [14] have devised and addressed the concept of modelling lifetime data using one parameter families of distributions, such as Akash, exponential, and Lindley distributions. Many lifetime datasets are employed to exhibit its adaptability over the exponential distribution. Shanker and Shukla [15] examined the two-parameter Akash model and determined its statistical characteristic, estimation problem, and application to it. As a reason, the Akash distribution can be used as an alternate lifetime model in reliability analysis.
The maximal likelihood estimation (MLE) is well-known estimation approach. Despite the fact that MLE is efficient and has strong conceptual features, there is confirmation that it does not work well, especially with small samples. As a result, different estimation approaches have been offered in the studies as options to conventional method. The weighted least-squares estimation (WLSE), L-moments estimator (LME), percentile estimator (PCE) and least squares estimator (LSE) are among the most frequently recommended. These approaches, in general, do not possess desirable theoretical features, but they can offer better estimates of unknown parameters in specific instances than the MLE. Various estimating approaches for many models have been investigated in the studies, as illustrations [16,17,18,19,20,21]. The goal of this research is to give a mechanism for expert statisticians to choose the best evaluation method for the Two-Component Mixture of Akash Model (TC-MAM). In this investigation, we estimate the TC-MAM using LSE and WLSE, in conjunction to MLE.
Our goal in this investigation is to develop a novel mixture model for modeling real lifespan datasets from various disciplines of knowledge that is better fitting than mixture of Shanker, exponential and Lindley distributions. The TC-MAM model has an advantage over the Shanker and exponential models because the exponential distribution has a constant hazard rate function and the Shanker model has an increasing hazard rate function, whereas the failure rate function for a TC-MAM model exhibits monotonically increasing, modified declining, decreasing–increasing–decreasing (DID), declining–increasing (DI), and upside-down bathtub behavior. The novel TC-MAM is being developed in particular to offer a novel flexible parametric model for modeling complex data that emerges in dependability research, investigation of lifespan, quality control, statistical mechanics, economics, biological investigations, and other fields. The purpose is to provide a novel model for lifespan analysis that can handle various types of failure rates, as well as various close form features of novel model with simple physical interpretations.
The originality of this research is due to the fact that we present a thorough explanation of the statistical aspects of TC-MAM in the hopes of attracting more applications in lifespan analysis. Additionally, as far as we know, no investigation has been performed to evaluate all of these estimators of the TC-MAM, as well as their mathematical and statistical features and assessment methods to estimate of unexplained parameters of TC-MAM. For various sample sizes and parametric values, we demonstrate how alternative frequentist estimators of the suggested distribution work.

2. The Two-Component Mixture of Akash Model

A r . v . T is stated to have a TC-MAM if its PDF and CDF can be integrated as:
f ( t ϖ ) = δ f 1 ( t λ 1 ) + δ ˜ f 2 ( t λ 2 ) , δ ˜ = 1 δ , t , λ i > 0 , i = 1 , 2 .
f ( t ϖ ) = δ λ 1 3 λ 1 2 + 2 ( 1 + t 2 ) exp λ 1 t + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t 2 ) exp λ 2 t ,
and
F ( t ϖ ) = δ F 1 ( t λ 1 ) + δ ˜ F 2 ( t λ 2 ) ,
F ( t ϖ ) = δ 1 1 + λ 1 t λ 1 t + 2 λ 1 2 + 2 exp λ 1 t + δ ˜ 1 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t ,
where ϖ = λ 1 , λ 2 , δ and δ is a positive mixing parameter, whereas λ i are positive scale parameters.

2.1. Mode

By tackling the given non-linear equation with respect to t, the mode of the TC-MAM( ϖ ) is derived
δ λ 1 3 λ 1 2 + 2 exp λ 1 t t ( 2 λ 1 t ) λ 1 + δ ˜ λ 2 3 λ 2 2 + 2 exp λ 2 t t ( 2 λ 2 t ) λ 2 = 0 .

2.2. Median

The median of TC-MAM ϖ is given here. Let F ( t ϖ ) be CDF of TC-MAM ϖ the median is at 50th quantiles that is Q 0.5 . The median ( t * ) can, therefore, be determined by resolving given equation for t.
δ 1 1 + λ 1 t λ 1 t + 2 β 1 2 + 2 exp λ 1 t + δ ˜ 1 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t = 0.5 ,
δ 1 + λ 1 t λ 1 t + 2 β 1 2 + 2 exp λ 1 t + δ ˜ 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t = 0.5 ,
Numerical strategies like Newton–Raphson approach can be utilised to find t * from Equation (7).
Several graphs of PDF and CDF of TC-MAM, as well as both component densities, for various parametric values are shown in Figure 1 and Figure 2. It should be indicated that input parameters were selected at random until a wide range of patterns could be examined. The PDF exemplifies its adaptability. The PDF curves of TC-MAM( ϖ ) indicate that it can be monotonically decreasing, positively skewed, inverted U, and declining–increasing–decreasing (DID), as well as modified monotonically decreasing with platykurtic, mesokurtic, and leptokurtic curves. As a result, it can be used to model a diverse set of data.

2.3. mth Moments about Origin

For a r.v. T, the mth moments of TC-MAM ϖ are as:
μ ˘ m = E ( T m ) = 0 t m f ( t Δ ˘ ) d t = 0 t m δ ( 1 + t 2 ) λ 1 3 λ 1 2 + 2 exp λ 1 t + δ ˜ ( 1 + t 2 ) λ 2 3 λ 2 2 + 2 exp λ 2 t d t ,
E ( T m ) = δ m ! λ 1 2 + m + 1 m + 2 λ 1 m λ 1 2 + 2 + δ ˜ m ! λ 2 2 + m + 1 m + 2 λ 2 m λ 2 2 + 2 , m = 1 , 2 .
The mean of the TC-MAM ϖ is:
μ ˘ 1 = δ λ 1 2 + 6 λ 1 λ 1 2 + 2 + δ ˜ λ 1 2 + 6 λ 1 λ 1 2 + 2 = μ ,
while the variance is given by
σ ˘ 2 = δ λ 1 4 + 16 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 2 + δ ˜ λ 2 4 + 16 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 2 .
Graphs of the mean and variance of TC-MAM ( ϖ ) for a variety of parameter values that can be identified in Figure 3 and Figure 4. The mean of TC-MAM ( ϖ ), shows a monotonically decreasing behavior for fixed value of λ 1 and δ and varying values of λ 2 (see Figure 3a). The escalating changes of mixing parameter δ enhanced the mean concentration, according to this analysis. We draw mean graphs (see Figure 3b) to demonstrate the behavior of the mean for fixed values of λ 2 and δ and varying values of λ 1 . It reveals the characteristics of component parameter λ 1 in relation to a mean profile. The boosting attitude of δ reduce the concentration of mean for the all varying values of λ 1 . The significances of the mixing parameter versus the mean profile are defined in Figure 3c for various levels of λ 1 . From these drawn lines, it can be deduced that the concentration of the mean profile is a deteriorating function for parameter λ 1 . The variance exhibits the same behavior as the mean in all scenarios (see Figure 4a–c).
In particular moments about origin
μ ˘ 2 = δ 2 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 + δ ˜ 2 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 ,
μ ˘ 3 = δ 6 λ 1 2 + 20 λ 1 3 λ 1 2 + 2 + δ ˜ 6 λ 2 2 + 20 λ 2 3 λ 2 2 + 2 ,
μ ˘ 4 = δ 24 λ 1 2 + 30 λ 1 4 λ 1 2 + 2 + δ ˜ 24 λ 2 2 + 30 λ 2 4 λ 2 2 + 2 ,
and the moments about mean of the TC-MAM ϖ are:
μ 2 = δ λ 1 4 + 16 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 2 + δ ˜ λ 2 4 + 16 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 2 ,
μ 3 = δ 2 λ 1 6 + 30 λ 1 4 + 36 λ 1 2 + 24 λ 1 3 λ 1 2 + 2 3 + δ ˜ 2 λ 2 6 + 30 λ 2 4 + 36 λ 2 2 + 24 λ 2 3 λ 2 2 + 2 3 ,
μ 4 = δ 3 3 λ 1 8 + 128 λ 1 6 + 408 λ 1 4 + 576 λ 1 2 + 240 λ 1 4 λ 1 2 + 2 4 + δ ˜ 3 3 λ 2 8 + 128 λ 2 6 + 408 λ 2 4 + 576 λ 2 2 + 240 λ 2 4 λ 2 2 + 2 4 .
The φ ˘ C V (Coefficient of Variation), Ψ ˘ S k (Skewness) and ψ ˘ K (Kurtosis) of TC-MAM ϖ are:
φ ˘ C V = δ λ 1 4 + 16 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 2 + δ ˜ λ 2 4 + 16 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 2 δ λ 1 2 + 6 λ 1 λ 1 2 + 2 + δ ˜ λ 1 2 + 6 λ 1 λ 1 2 + 2 ,
Ψ S k = δ 2 λ 1 6 + 30 λ 1 4 + 36 λ 1 2 + 24 λ 1 3 λ 1 2 + 2 3 + δ ˜ 2 λ 2 6 + 30 λ 2 4 + 36 λ 2 2 + 24 λ 2 3 λ 2 2 + 2 3 δ λ 1 4 + 16 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 2 + δ ˜ λ 2 4 + 16 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 2 3 / 2 ,
ψ ˘ K = δ 3 3 λ 1 8 + 128 λ 1 6 + 408 λ 1 4 + 576 λ 1 2 + 240 λ 1 4 λ 1 2 + 2 4 + δ ˜ 3 3 λ 2 8 + 128 λ 2 6 + 408 λ 2 4 + 576 λ 2 2 + 240 λ 2 4 λ 2 2 + 2 4 δ λ 1 4 + 16 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 2 + δ ˜ λ 2 4 + 16 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 2 2 ,
and Index of Dispersion I D is
I D = δ λ 1 4 + 16 λ 1 2 + 12 λ 1 2 λ 1 2 + 2 2 + δ ˜ λ 2 4 + 16 λ 2 2 + 12 λ 2 2 λ 2 2 + 2 2 δ λ 1 2 + 6 λ 1 λ 1 2 + 2 + δ ˜ λ 1 2 + 6 λ 1 λ 1 2 + 2 .
The TC-MAM ϖ is readily explained to be over-distributed when μ 2 > μ , equi-dispersed μ 2 = μ , as well as under-dispersed μ 2 < μ .
Graphs of the ID of TC-MAM ( ϖ ) for various parameter settings are illustrated in Figure 5. The boosting attitude of δ rise the concentration of mean for the all varying values of λ 2 (see Figure 5a). However, boosting attitude of δ reduces the concentration of mean for the all varying values of λ 1 (see Figure 5b). The effects of δ against the concentration of the mean profile are shown in Figure 5c. The concentration of the mean profile is a decreasing function for parameter λ 1 according to these depicted lines. Figure 6 and Figure 7 explain the nature of Ψ S k and ψ ˘ K in relation to λ 1 , λ 2 and δ . The coefficient of skewness and kurtosis of TC-MAM ( ϖ ), shows a decreasing behavior for fixed value of λ 1 and δ and varying values of λ 2 (see Figure 6a and Figure 7a).
To expose the behavior of Ψ S k and ψ ˘ K for fixed value of λ 2 and δ and varying values of λ 1 (see Figure 6b and Figure 7b). The escalating changes of mixing parameter δ enhanced Ψ S k and ψ ˘ K concentration, according to this analysis. The significances of the mixing parameter versus the coefficient of skewness and kurtosis profile are defined in Figure 6c and Figure 7c for various levels of λ 1 . From these drawn lines, it can be deduced that the concentration of the skewness and kurtosis profile is a deteriorating function for parameter λ 1 .

2.4. Moment Generating Function

The MGF of TC-MAM ϖ is specified as:
M ˜ t υ = E e t υ = 0 e t υ δ λ 1 3 λ 1 2 + 2 ( 1 + t 2 ) exp λ 1 t + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t 2 ) exp λ 2 t d t ,
M ˜ t υ = δ λ 1 3 λ 1 2 + 2 1 λ 1 υ + 2 λ 1 υ 3 + δ ˜ λ 2 3 λ 2 2 + 2 1 λ 2 υ + 2 λ 2 υ 3 .
M ˜ t υ = δ k = 0 λ 1 2 + k + 1 k + 2 λ 1 2 + 2 υ λ 1 k + δ ˜ k 1 = 0 λ 2 2 + k 1 + 1 k 1 + 2 λ 2 2 + 2 υ λ 2 k 1 .

2.5. Cumulants

The cumulants (CF), ξ ˘ υ = E exp i υ t of TC-MAM ϖ is derived by plugging υ with ‘ i υ ’ in Equation (22), the following formula can be used to obtain the CF:
ξ ˘ υ = δ k = 0 λ 1 2 + k + 1 k + 2 λ 1 2 + 2 i υ λ 1 k + δ ˜ k 1 = 0 λ 2 2 + k 1 + 1 k 1 + 2 λ 2 2 + 2 i υ λ 2 k 1 ,
where the complex unit i = 1 .

2.6. Probability Generating Function (PGF)

In Equation (22), the PGF by plugging υ with “ln( ω )” is:
P t ω = E ω t = δ k = 0 λ 1 2 + k + 1 k + 2 λ 1 2 + 2 l n ( ω ) λ 1 k + δ ˜ k 1 = 0 λ 2 2 + k 1 + 1 k 1 + 2 λ 2 2 + 2 l n ( ω ) λ 2 k 1 .

2.7. Factorial Moment Generating Function

By plugging υ with ‘ln 1 + ϕ ’ in Equation (22), the FMGF can be shown as
F ˘ t ω   = E e t ln 1 + ϕ = δ k = 0 λ 1 2 + k + 1 k + 2 λ 1 2 + 2 l n 1 + ϕ λ 1 k + δ ˜ k 1 = 0 λ 2 2 + k 1 + 1 k 1 + 2 λ 2 2 + 2 l n 1 + ϕ λ 2 k 1 .

3. Reliability Measures

In reliability framework, lifetime models are classified using the reliability/survival function and the failure/hazard rate function. TC-MAM ϖ is currently being studied for its reliability properties.

3.1. Reliability Function

The reliability function R t ϖ of TC-MAM ϖ is.
R t ϖ = δ 1 + λ 1 t λ 1 t + 2 λ 1 2 + 2 exp λ 1 t + δ ˜ 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t

3.2. Hazard Function

The failure rate function h ( t ϖ ) of the TC-MAM( ϖ ) is described as follows.
h ( t ϖ ) = δ λ 1 3 λ 1 2 + 2 ( 1 + t 2 ) exp λ 1 t + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t 2 ) exp λ 2 t δ 1 + λ 1 t λ 1 t + 2 λ 1 2 + 2 exp λ 1 t + δ ˜ 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t .
In Figure 8, the HRF of TC-MAM ϖ shows monotonically increasing, modified decreasing, decreasing–increasing–decreasing (DID), decreasing–increasing (DI), and upside down bathtub behavior. Figure 8a,b,d signifies that the reduction in failure rate function profile and is noted by enlarging the value of mixing parameter δ , and for λ 1 < λ 2 . Figure 8c,e,f exhibits the diversion in the failure rate function for various values of δ . It is found in Figure 8 that the failure rate distribution is expanding due to higher the value of mixing parameter δ , and for λ 1 > λ 2 .

3.3. Mills Ratio

Mills ratio is an another method of quantifying reliability due to its relation to failure rate. Mills ratio Υ ˇ t ϖ of TC-MAM( ϖ ) is
Υ ˇ t ϖ = R t ϖ f t ϖ = δ 1 + λ 1 t λ 1 t + 2 λ 1 2 + 2 exp λ 1 t + δ ˜ 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t δ λ 1 3 λ 1 2 + 2 ( 1 + t 2 ) exp λ 1 t + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t 2 ) exp λ 2 t .

3.4. Cumulative Hazard Rate Function

The CHRF of TC-MAM ϖ is
H t ϖ = 0 t h ( y ϖ ) d y = log R ( t ϖ ) .
It is a risk indicator: the stronger the H t ϖ estimate, the greater the chance of failure by t-time. It must be stated that
R ( t ϖ ) = e H t ϖ and f t ϖ = h t ϖ e H t ϖ .
So,
H ( t ϖ ) = log δ 1 + λ 1 t λ 1 t + 2 λ 1 2 + 2 exp λ 1 t + δ ˜ 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t .

3.5. Reversed Hazard Rate Function

The RHRF of a random life of TC-MAM ϖ is defined as
h ˘ ( t ϖ ) = δ λ 1 3 λ 1 2 + 2 ( 1 + t 2 ) exp λ 1 t + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t 2 ) exp λ 2 t 1 δ 1 + λ 1 t λ 1 t + 2 λ 1 2 + 2 exp λ 1 t δ ˜ 1 + λ 2 t λ 2 t + 2 λ 2 2 + 2 exp λ 2 t .

3.6. Mean Time to Failure (MTTF)

The expected time for which the device performs efficiently is given by the mean time to failure (MTTF). If TC-MAM ϖ then reliability function is used to express MTTF, which is as follows:
M ˘ t ϖ = 0 + R t ϖ d t ,
and R t is provided in Equation (28). Thus
M ˘ t ϖ = δ λ 1 2 + 6 λ 1 λ 1 2 + 2 + δ ˜ λ 1 2 + 6 λ 1 λ 1 2 + 2 .

4. Estimation Inference via Simulation

Given that the parametric vector ϖ is undetermined, certain statistical properties of the TC-MAM( ϖ ) are presented to this section. The evaluation of parametric vector ϖ is accomplished by three widely known estimation mechanisms, such as MLE, LSE, and WLSE. From now, t 1 , t 2 , , t n signify n determined values from T and their ascending sorting values t 1 t 2 t n .

4.1. Maximum Likelihood Estimation (MLE)

The MLE method is the best methodology for parameter assessment. The popularity of the approach stems from its many advantageous characteristics, such as consistency, normality, and asymptotic efficiency. Let t 1 , t 2 , , t n be n determined values from the Equation (2) and ϖ be the vector of undetermined parameters. The evaluations of MLEs of ϖ can be given by optimizing the likelihood function with respect to λ 1 , λ 2 , and δ given by L t ϖ = i = 1 n f ( t i ; ϖ ) or likewise the log-likelihood function for ϖ is
l t ϖ = ln i = 1 n f ( t i ; ϖ )
l t ϖ = i = 1 n ln δ λ 1 3 λ 1 2 + 2 ( 1 + t i 2 ) exp λ 1 t i + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t i 2 ) exp λ 2 t i .
So, by partially differentiating l t ϖ in terms of each parameter ( λ 1 , λ 2 , δ ) and placing the results to zero, the MLEs of the relevant parameters are determined as
l t ϖ λ 1 = i = 1 n δ ( 1 + t i 2 ) λ 1 2 exp λ 1 t i 3 λ 1 2 + 2 λ 1 t i λ 1 2 + 2 2 λ 1 2 λ 1 2 + 2 2 δ λ 1 3 λ 1 2 + 2 ( 1 + t i 2 ) exp λ 1 t i + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t i 2 ) exp λ 2 t i ,
l t ϖ λ 2 = i = 1 n δ ˜ ( 1 + t i 2 ) λ 2 2 exp λ 2 t i 3 λ 2 2 + 2 λ 2 t i λ 2 2 + 2 2 λ 2 2 λ 2 2 + 2 2 δ λ 1 3 λ 1 2 + 2 ( 1 + t i 2 ) exp λ 1 t i + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t i 2 ) exp λ 2 t i ,
l t ϖ δ = i = 1 n λ 1 3 λ 1 2 + 2 ( 1 + t i 2 ) exp λ 1 t i λ 2 3 λ 2 2 + 2 ( 1 + t i 2 ) exp λ 2 t i δ λ 1 3 λ 1 2 + 2 ( 1 + t i 2 ) exp λ 1 t i + δ ˜ λ 2 3 λ 2 2 + 2 ( 1 + t i 2 ) exp λ 2 t i .
As a consequence, the MLE is found by evaluating this non-linear set of equations. However such equations cannot be handled analytically, we can use statistical software to solve them using an iterative methodology namely the Newton method or fixed point iteration methods.

4.2. Least Square Estimators (LSE)

The ordinary least square approach [22] is widely used for assessing undetermined parameters. The LSEs of λ 1 , λ 2 , and δ , indicated by λ ˜ 1 L S E , λ ˜ 2 L S E , and δ ˜ L S E , can be determined by minimizing Equation (42)
L S ( ϖ ) = i = 1 n F ( t ( i ) ϖ ) i n + 1 2 ,
with respect to λ 1 , λ 2 , and δ , where F ( ·)is given by Equation(4). They may be determined in the similar way by solving the non-linear equations below:
L S ( ϖ ) λ 1 = i = 1 n F ( t ( i ) ϖ ) i n + 1 Ψ ˘ 1 ( t ( i ) λ 1 ) = 0 ,
L S ( ϖ ) λ 2 = i = 1 n F ( t ( i ) ϖ ) i n + 1 Ψ ˘ 2 ( t ( i ) λ 2 ) = 0 ,
and
L S ( ϖ ) δ = i = 1 n F ( t ( i ) ϖ ) i n + 1 Ψ ˘ 3 ( t ( i ) δ ) = 0 ,
where
Ψ ˜ 1 ( t ( i ) λ 1 ) = δ t ( i ) λ 1 2 exp λ 1 t ( i ) 6 + λ 1 2 + 2 λ 1 t ( i ) + 2 + λ 1 2 t ( i ) λ 1 2 + 2 2 ,
Ψ ˜ 2 ( t ( i ) λ 2 ) = δ ˜ t ( i ) λ 2 2 exp λ 2 t ( i ) 6 + λ 2 2 + 2 λ 2 t ( i ) + 2 + λ 2 2 t ( i ) λ 2 2 + 2 2 ,
Ψ ˜ 3 ( t ( i ) δ ) = 1 + λ 2 t ( i ) 2 + λ 2 t ( i ) λ 2 2 + 2 exp λ 2 t ( i ) 1 + λ 1 t ( i ) 2 + λ 1 t ( i ) λ 1 2 + 2 exp λ 1 t ( i ) .

4.3. Weighted Least Squares Estimators (WLSE)

Take a look at the following weighted function (see [23])
κ ˘ i = n + 1 2 n + 2 i n i + 1 .
The WLSEs λ ˜ 1 W L S E , λ ˜ 2 W L S E , and δ ˜ W L S E , can be obtained by minimizing Equation (50)
W L S ( ϖ ) = i = 1 n n + 1 2 n + 2 i n i + 1 F ( t ( i ) ϖ ) i n + 1 2 ,
One can also obtain these estimators by solving:
W L S ( ϖ ) λ 1 = i = 1 n n + 1 2 n + 2 i n i + 1 F ( t ( i ) ϖ ) i n + 1 Ψ ˘ 1 ( t ( i ) λ 1 ) = 0 ,
W L S ( ϖ ) λ 2 = i = 1 n n + 1 2 n + 2 i n i + 1 F ( t ( i ) ϖ ) i n + 1 Ψ ˘ 2 ( t ( i ) λ 2 ) = 0 ,
and
W L S ( ϖ ) δ = i = 1 n n + 1 2 n + 2 i n i + 1 F ( t ( i ) ϖ ) i n + 1 Ψ ˘ 3 ( t ( i ) δ ) = 0 ,
where Ψ ˜ 1 ( t ( i ) λ 1 ) , Ψ ˜ 2 ( t ( i ) λ 2 ) and Ψ ˜ 3 ( t ( i ) δ ) are given in Equations (46)–(48).

4.4. Simulation Study

The simulation study is used to evaluate the various estimating methodologies outlined in the preceding subsection. Monte Carlo simulations are performed with a variety of mixing proportion δ and distribution parameters. The performance of MLE, LSEs, and WLSEs of the TC-MAM( ϖ ) parameters is evaluated using four simulation experiments. The proficiency of the MLEs, LSEs, and WLSEs is discussed using the bias and MSE indicators. In terms of n, the efficiency of each parameter estimation strategy for the TC-MAM( ϖ ) model is examined. The simulation algorithm is subdivided into six steps:
  • By adjusting the mixing proportion and model parameters λ 1 , λ 2 , δ = Set-I 0.15 , 0.30 , 0.40 , Set-II 0.25 , 0.5 , 0.6 and Set-III 0.5 , 0.2 , 0.6 , generate random samples of sizes 25 , 30 , , 500 from TC-MAM( ϖ ). The random samples for the simulation are obtained as specified in the upcoming stage.
  • Generate a random variable u from uniform distribution U ( 0 , 1 ) , employing R uniform generator (runif).
  • If u δ , then generate a random variable from the first component, which is a Akash distribution λ 1 . If u > δ , the second component, Akash distribution λ 2 , is utilized to produce a random variate.
  • Follow (2 and 3) till you have the prescribed sample size n.
  • Employing 1000 iterations, continue steps 1–4 each time. Evaluate MLEs, LSEs, and WLSEs for the 1000 samples, say κ ˜ j for j = 1 , 2 , , 1000 , having optima function and the Nelder–Mead technique in R to compute estimates.
  • Determine biases and MSEs. The two metrics are utilized to meet these targets:
    B i a s κ n = 1 1000 j = 1 1000 κ ˜ j κ ,
    M S E κ n = 1 1000 j = 1 1000 κ ˜ κ 2 ,
    where κ = λ 1 , λ 2 , δ .
The empirical findings are depicted in Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14. These results suggest that the proposed estimation methods are effective at estimating the TC-MAM parameters. We can deduce that the estimators display asymptotic unbiasedness because the bias goes to zero as n rises. On the other hand, MSE behavior implies consistency because the errors trend to zero as n increases. The following conclusions can be drawn from Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14.
  • Under all three estimation procedures, the estimated bias of parameters λ 1 , λ 2 , δ reduces as n grows.
  • For parametric set-I, the estimated bias of parameter λ 1 under LSE, WLSE is negative, for Set II λ 1 and δ under all three estimation methods and for set III, the estimated bias of parameter δ is negative (see Figure 9, Figure 11 and Figure 13).
  • Figure 9 and Figure 10 show the bias and MSE of λ ˜ 1 , λ ˜ 2 , and δ ˜ , for parametric Set-I and the WLSE always has the smallest value of bias and MSE of all estimators.
  • In the second scenario, the MLE estimators of λ 2 is over-estimated, while MSE of λ 2 is highest among the three considered estimators (see Figure 11 and Figure 12).
  • The estimators of λ 1 are over-estimated in the third scenario, however the over and under estimation of λ 2 and δ are seen among the three investigated estimators, and the WLSE always has the minimum value of bias of all estimators (see Figure 13).
  • Among the three estimators evaluated, the MSE of λ 1 is the greatest (see Figure 14).
  • The MSE of λ ˜ 2 is strongly stimulated and higher under MLE and LSE estimation methods when n < 50 (see Figure 10).
  • Figure 13 and Figure 14 demonstrate the influence of choice of parameters on the estimation approaches, here bias and MSEs are comparatively low among the selected set of parameters.
  • Some big shifts in MSEs of considered estimators under MLE, LSE, and WLSE are observed when n < 50 .
  • In terms of bias, the WLSE’s performance is relatively favorable.
  • Furthermore, when n increases, the MSE for all three estimating strategies decreases, satisfying the consistency criteria (Figure 10, Figure 12 and Figure 14).
  • In all estimating methodologies, the difference between estimates and stated parameters reduces as n rises.
  • As n approaches infinity, WLSE estimation is frequently better in terms of bias and MSE when likened to other estimation methods for all given parameter values.
  • The estimated MSEs of parameters λ ˜ 1 , λ ˜ 2 , and δ ˜ under the MLE estimation technique decrease quickly as n increases, demonstrating the effectiveness of the MLE procedure.
The final conclusion drawn from the foregoing figures is that, as n rises, estimated bias and MSE graphs for estimators λ ˜ 1 , λ ˜ 2 , and δ ˜ finally approach zero for all estimating methods. This demonstrates the accuracy of both the estimating methods and the numerical computations for the TC-MAM parameters.

5. Applications

We demonstrate the flexibility of the TC-MAM in this section by examining a real dataset. The TC-MAM distribution is compared to competing models, such as the two component mixture of Shanker distribution (2C-MSM), the two component mixture of exponential model (2-CMEM), and the two component mixture of Lindley distribution (2C-MLM) using the R function maxLik(). The -Log-likelihood (-LL), the AIC, BIC, and AICC have all been used to compare these models. The model having the least quantities of above-mentioned goodness-of-fit (GoF) measures may be the best fit for the real dataset.
Dataset: There are 56 observations in this dataset pertaining to the burning velocity of various chemical substances. The laminar flame speed at the specified composition, temperature, and pressure circumstances is the burning speed/velocity. It lowers as the inhibitor concentration rises, and it may be observed directly by analysing the pressure distribution in the spherical vessel and monitoring the flame propagation. We consider a real-life dataset which represents the burning velocity (cm/s) of several chemical compounds to show the TC-MAM distribution’s suitability. This dataset is extracted from https://www.cheresources.com/mists.pdf (accessed on 4 September 2022) and the data are as follows: 68, 61, 64, 55, 51, 68, 44, 50, 82, 60, 89, 61, 54, 166, 66, 50, 87, 48, 42, 58, 46, 67, 46, 46, 44, 48, 56, 47, 54, 47, 89, 38, 108, 46, 40, 44, 312, 41, 31, 40, 41, 40, 56, 45, 43, 46, 46, 46, 46, 52, 58, 82, 71, 48, 39, and 41 [24,25,26,27] contains further data applications. The MLEs for the TC-MAM and GoF measures are shown in Table 1. The TC-MAM clearly outperforms the 2-CMSM, 2-CMEM, and 2-CMLM, as shown in Table 1. The profiles of the log-likelihood function (PLLF) based on the dataset that confirm the conclusions of Table 1 are shown in Figure 15. Figure 15 and Figure 16 show a graphical illustration of MLE existence and uniqueness, respectively. To summarize, the TC-MAM emerges as the better model for the dataset, indicating its usefulness in a real-world setting. We can deduce from this graphical representation and results obtain from Table 1 that the TC-MAM is a better fit for the dataset in consideration.

6. Conclusions

In this investigation, we used three estimated techniques: MLE, LSE, and WLSE to work on two component mixtures of Akash models. In particular, the Akash mixing model’s statistical and reliability features were achieved, such as central moments, Cumulants, Cumulant Generating Function, Probability Generating Function, Mean Time to Failure, Factorial Moment Generating Function, Coefficient of variation, Mills ratio, skewness and kurtosis, Reversed Hazard Rate Function, and Mean Residual Life. To investigate and assess the estimating approaches’ performance, a simulation study with 1000 iterations was performed and it was noted that when n increases, the estimated MSEs of parameters λ ˜ 1 , λ ˜ 2 , and δ ˜ under the MLE estimation technique rapidly decrease, illustrating the efficiency of the MLE procedure. As a result, we found that estimating model unknown parameters with regards of accuracy and consistency, the MLE approach surpassed the rest. Furthermore, we used real datasets to explain the utility of the underlying mixture model.

Author Contributions

Author Contributions: Conceptualization, A.S., T.N.S. and S.A.L.; methodology, A.S., M.K.H.H. and T.N.S.; validation, S.A.L. and T.N.S.; formal analysis, A.S., T.N.S. and K.N.; investigation, M.K.H.H., T.N.S. and S.A.L.; writing—original draft preparation, A.S., T.N.S. and S.A.L.; writing—review and editing, A.S., M.K.H.H. and K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All the data are available in this manuscript.

Conflicts of Interest

Authors declare that they do not have conflicts of interest.

Nomenclature

Symbols
Υ ˇ t ϖ Mills Ratio
P t ϖ PGF
H ( t ϖ ) CHRF
f ( t ϖ ) PDF
M ˜ t ϖ MGF
R ( t ϖ ) RF
K ˘ ϖ CGF
M ˘ t ϖ MTTF
F ( t ϖ ) CDF
h ( t ϖ ) HRF
Q ( q ; ϖ ) QF
ξ ˘ t υ CF
F ˘ t ω FMGF
h t ϖ RHRF
M ˘ R t ϖ MRL
Abbreviations
CHRFCumulative Hazard Rate Function
MGFMoment Generating Function
PDFProbability Density Function
CGFCumulant Generating Function
CDFCumulative Distribution Function
FMGFFactorial Moment Generating Function
RHRFReversed Hazard Rate Function
PGFProbability Generating Function
WLSEWeighted Least Square Estimator
MLEMaximum likelihood Estimator
AICCAkaike Information Criterion Corrected
AICAkaike Information Criterion
TTFTime-To-Failure
QFQuantile Function
RFReliability Function
CFCharacteristic Function
MSEMean square error
MRLMean Residual Life
LSELeast Square Estimator
r.v.Random Variable
HRFHazard Rate Function
MTTFMean Time to Failure
GoFGoodness-of-Fit
BICBayesian Information Criterion

References

  1. Everitt, B.S. A finite mixture model for the clustering of mixed-mode data. Stat. Probab. Lett. 1988, 6, 305–309. [Google Scholar] [CrossRef]
  2. Lindsay, B.G. Mixture models: Theory, geometry and applications. In NSF-CBMS Regional Conference Series in Probability and Statistics; Institute of Mathematical Statistics: New York, NY, USA, 1995; p. i-163. Available online: https://www.jstor.org/stable/4153184 (accessed on 4 September 2022).
  3. McLachlan, G.J.; Basford, K.E. Mixture Models: Inference and Applications to Clustering; M. Dekker: New York, NY, USA, 1988; Volume 38. [Google Scholar]
  4. McLachlan, G.; Peel, D. Finite Mixture Models; John Wiley & Sons: New York, NY, USA, 2000. [Google Scholar]
  5. Sultan, K.S.; Ismail, M.A.; Al-Moisheer, A.S. Mixture of two inverse Weibull distributions: Properties and estimation. Comput. Stat. Data Anal. 2007, 51, 5377–5387. [Google Scholar] [CrossRef]
  6. Mohammad, D.; Muhammad, A. On the Mixture of BurrXII and Weibull Distribution. J. Stat. Appl. Probab. 2014, 3, 251–267. [Google Scholar]
  7. Sindhu, T.N.; Feroze, N.; Aslam, M.; Shafiq, A. Bayesian inference of mixture of two Rayleigh distributions: A new look. Punjab Univ. J. Math. 2020, 48, 49–64. [Google Scholar]
  8. Sindhu, T.N.; Hussain, Z.; Aslam, M. Parameter and reliability estimation of inverted Maxwell mixture model. J. Stat. Manag. Syst. 2019, 22, 459–493. [Google Scholar] [CrossRef]
  9. Ateya, S.F. Maximum likelihood estimation under a finite mixture of generalized exponential distributions based on censored data. Stat. Pap. 2014, 55, 311–325. [Google Scholar] [CrossRef]
  10. Mohammadi, A.; Salehi-Rad, M.R.; Wit, E.C. Using mixture of Gamma distributions for Bayesian analysis in an M/G/1 queue with optional second service. Comput. Stat. 2013, 28, 683–700. [Google Scholar] [CrossRef]
  11. Al-Moisheer, A.S.; Daghestani, A.F.; Sultan, K.S. Mixture of Two One-Parameter Lindley Distributions: Properties and Estimation. J. Stat. Theory Pract. 2021, 15, 11. [Google Scholar] [CrossRef]
  12. Hussain, E. The Non-Linear Functions of Order Statistics and Their Properties in Selected Probability Models. Ph.D. Thesis, Department of Statistics, University of Karachi, Karachi, Pakistan, 2006. [Google Scholar]
  13. Shanker, R. Akash distribution and its applications. Int. J. Probab. Stat. 2015, 4, 65–75. [Google Scholar]
  14. Shanker, R.; Hagos, F.; Sujatha, S. On modeling of lifetime data using one parameter Akash, Lindley and exponential distributions. Biom. Biostat. Int. J. 2016, 3, 1–10. [Google Scholar] [CrossRef]
  15. Shanker, R.; Shukla, K.K. On two-parameter Akash distribution. Biom. Biostat. Int. J. 2017, 6, 416–425. [Google Scholar] [CrossRef]
  16. Dey, S.; Kumar, D.; Ramos, P.L.; Louzada, F. Exponentiated Chen distribution: Properties and estimation. Commun. Stat.-Simul. Comput. 2017, 46, 8118–8139. [Google Scholar] [CrossRef]
  17. Dey, S.; Alzaatreh, A.; Zhang, C.; Kumar, D. A new extension of generalized exponential distribution with application to ozone data. Ozone Sci. Eng. 2017, 39, 273–285. [Google Scholar] [CrossRef]
  18. Rodrigues, G.C.; Louzada, F.; Ramos, P.L. Poisson exponential distribution: Different methods of estimation. J. Appl. Stat. 2018, 45, 128–144. [Google Scholar] [CrossRef]
  19. Dey, S.; Moala, F.A.; Kumar, D. Statistical properties and different methods of estimation of Gompertz distribution with application. J. Stat. Manag. Syst. 2018, 21, 839–876. [Google Scholar] [CrossRef]
  20. Dey, S.; Josmar, M.J.; Nadarajah, S. Kumaraswamy distribution: Different methods of estimation. Comput. Appl. Math. 2018, 37, 2094–2111. [Google Scholar] [CrossRef]
  21. Shafiq, A.; Batur Çolak, A.; Naz Sindhu, T.; Ahmad Lone, S.; Alsubie, A.; Jarad, F. Comparative study of artificial neural network versus parametric method in COVID-19 data analysis. Results Phys. 2022, 38, 105613. [Google Scholar] [CrossRef]
  22. Swain, J.J.; Venkatraman, S.; Wilson, J.R. Least-squares estimation of distribution functions in Johnson’s translation system. J. Stat. Comput. Simul. 1988, 29, 271–297. [Google Scholar] [CrossRef]
  23. Gupta, R.D.; Kundu, D. Generalized exponential distribution: Different method of estimations. J. Stat. Simul. 2001, 69, 315–337. [Google Scholar] [CrossRef]
  24. Sindhu, T.N.; Shafiq, A.; Al-Mdallal, Q.M. On the analysis of number of deaths due to Covid- 19 outbreak data using a new class of distributions. Results Phys. 2021, 21, 103747. [Google Scholar] [CrossRef]
  25. Lone, S.A.; Sindhu, T.N.; Shafiq, A.; Jarad, F. A novel extended Gumbel Type II model with statistical inference and COVID-19 applications. Results Phys. 2022, 35, 105377. [Google Scholar] [CrossRef]
  26. Sindhu, T.N.; Hussain, Z.; Alotaibi, N.; Muhammad, T. Estimation method of mixture distribution and modeling of COVID-19 pandemic. Aims Math. 2022, 7, 9926–9956. [Google Scholar] [CrossRef]
  27. Lone, S.A.; Sindhu, T.N.; Jarad, F. Additive Trinomial Fréchet distribution with practical application. Results Phys. 2022, 33, 105087. [Google Scholar] [CrossRef]
Figure 1. Behavior of f 1 ( t λ 1 ) (first component density), f 2 ( t λ 2 ) (second component density) and density of TC-MAM f m ( t ϖ ) with δ against t.
Figure 1. Behavior of f 1 ( t λ 1 ) (first component density), f 2 ( t λ 2 ) (second component density) and density of TC-MAM f m ( t ϖ ) with δ against t.
Axioms 11 00516 g001
Figure 2. Behavior of F 1 ( t λ 1 ) (first component CDF), F 2 ( t λ 2 ) (second component CDF) and CDF of TC-MAM ( F m ( t ϖ ) ) with δ against t.
Figure 2. Behavior of F 1 ( t λ 1 ) (first component CDF), F 2 ( t λ 2 ) (second component CDF) and CDF of TC-MAM ( F m ( t ϖ ) ) with δ against t.
Axioms 11 00516 g002
Figure 3. Variations of Mean of TC-MAM.
Figure 3. Variations of Mean of TC-MAM.
Axioms 11 00516 g003
Figure 4. Variations of Variance of TC-MAM.
Figure 4. Variations of Variance of TC-MAM.
Axioms 11 00516 g004
Figure 5. Variations of index of dispersion of TC-MAM ϖ .
Figure 5. Variations of index of dispersion of TC-MAM ϖ .
Axioms 11 00516 g005
Figure 6. Variations of coefficient of skewness S K , of TC-MAM ϖ .
Figure 6. Variations of coefficient of skewness S K , of TC-MAM ϖ .
Axioms 11 00516 g006
Figure 7. Behavior of coefficient of kurtosis K of TC-MAM ϖ .
Figure 7. Behavior of coefficient of kurtosis K of TC-MAM ϖ .
Axioms 11 00516 g007
Figure 8. Variations in h ( t ϖ ) for λ 1 , λ 1 and δ .
Figure 8. Variations in h ( t ϖ ) for λ 1 , λ 1 and δ .
Axioms 11 00516 g008
Figure 9. Behavior of bias of estimators with different methods under parametric set I against n.
Figure 9. Behavior of bias of estimators with different methods under parametric set I against n.
Axioms 11 00516 g009
Figure 10. Behavior of MSE of estimators with different methods under parametric set I against n.
Figure 10. Behavior of MSE of estimators with different methods under parametric set I against n.
Axioms 11 00516 g010
Figure 11. Behavior of bias of estimators with different methods under parametric set II against n.
Figure 11. Behavior of bias of estimators with different methods under parametric set II against n.
Axioms 11 00516 g011
Figure 12. Behavior of MSE of estimators with different methods under parametric set II against n.
Figure 12. Behavior of MSE of estimators with different methods under parametric set II against n.
Axioms 11 00516 g012
Figure 13. Behavior of bias of estimators with different methods under parametric set III against n.
Figure 13. Behavior of bias of estimators with different methods under parametric set III against n.
Axioms 11 00516 g013
Figure 14. Behavior of MSE of estimators with different methods under parametric set III against n.
Figure 14. Behavior of MSE of estimators with different methods under parametric set III against n.
Axioms 11 00516 g014
Figure 15. The plots of PLLF for Dataset.
Figure 15. The plots of PLLF for Dataset.
Axioms 11 00516 g015
Figure 16. The graphs of score functions cross the horizontal axis at λ ˜ 1 , λ ˜ 2 , and δ ˜ of Dataset.
Figure 16. The graphs of score functions cross the horizontal axis at λ ˜ 1 , λ ˜ 2 , and δ ˜ of Dataset.
Axioms 11 00516 g016
Table 1. MLEs, and GoF statistics for the Dataset I.
Table 1. MLEs, and GoF statistics for the Dataset I.
DistributionsMLEs L L AICBICAICC
TC-MAM λ ˜ 1 0.011647260.0557526.1114532.1875526.5729
λ ˜ 2 0.053562
δ ˜ 0.024836
2C-MSM λ ˜ 1 0.008656269.2051544.4102550.4863544.8717
λ ˜ 2 0.035161
δ ˜ 0.023762
2C-MLM λ ˜ 1 0.008684270.0022546.0045552.0805546.4659
λ ˜ 2 0.034539
δ ˜ 0.023282
2C-MEM λ ˜ 1 0.016401286.1761578.3523584.4283578.8137
λ ˜ 2 0.016396
δ ˜ 0.746712
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shafiq, A.; Sindhu, T.N.; Lone, S.A.; Hassan, M.K.H.; Nonlaopon, K. Mixture of Akash Distributions: Estimation, Simulation and Application. Axioms 2022, 11, 516. https://doi.org/10.3390/axioms11100516

AMA Style

Shafiq A, Sindhu TN, Lone SA, Hassan MKH, Nonlaopon K. Mixture of Akash Distributions: Estimation, Simulation and Application. Axioms. 2022; 11(10):516. https://doi.org/10.3390/axioms11100516

Chicago/Turabian Style

Shafiq, Anum, Tabassum Naz Sindhu, Showkat Ahmad Lone, Marwa K. H. Hassan, and Kamsing Nonlaopon. 2022. "Mixture of Akash Distributions: Estimation, Simulation and Application" Axioms 11, no. 10: 516. https://doi.org/10.3390/axioms11100516

APA Style

Shafiq, A., Sindhu, T. N., Lone, S. A., Hassan, M. K. H., & Nonlaopon, K. (2022). Mixture of Akash Distributions: Estimation, Simulation and Application. Axioms, 11(10), 516. https://doi.org/10.3390/axioms11100516

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