A New Truncated Muth Generated Family of Distributions with Applications

In recent years, theMuth distribution has been used for the construction of accurate statistical models, with applications in various applied fields. In this paper, we use a truncated-composed scheme to create a new unit Muth distribution, fromwhich we motivate a more general family of continuous distributions called the truncated Muth generated family. )e key benefits of this family are its analytical simplicity, connections with the exponential generated family, and flexibility conferred on any parental distribution. In particular, it improves the capability of the functions of the parental distribution, enhancing their peak, asymmetry, tail, and flatness levels, among others.)e characteristics of quantile andmoment measures and functions of the truncatedMuth generated family are described in detail. As a concrete example, a particular distribution that extends theWeibull distribution is highlighted. In an applied part, the parameters are calculated using the maximum likelihood procedure. We use a comprehensive simulation analysis to demonstrate the accuracy of the derived estimates. )e revised Weibull model is then used to fit two real-world datasets. )e new model is shown to be more suited to these datasets than other competing models.


Introduction
e Muth (M) distribution, in its most basic form, is a lifetime distribution introduced by [1]. It is mathematically defined by the following one-parameter cumulative distribution function (cdf ): with α ∈ (0, 1], and F M (x) � 0 for x < 0. e following properties are satisfied. As α approaches 0, the M distribution transforms into a classical exponential distribution with the parameter equal to one. Its right tail is less weighted than that of some other lifetime distributions, and it has enough versatility to suit a large panel of lifetime data sets properly, especially those resulting from reliability experiments. Many of these properties are developed in [1][2][3]. e M distribution was later generalized by Jodra et al. [4] using the power transform. It is also possible to refer to the exponentiated Tessier distribution by Sharma et al. [5], where it provides an in-depth analytical connection with the exponentiated version of the M distribution.
Also, an innovative approach has sought to develop generic families of continuous distributions that generalize or expand a parental distribution using specific mathematical schemes. By combining the type II of the T-X scheme by Alzaatreh et al. [6] and the M distribution, in [7], the authors proposed the M generated (M-G) class of distributions with the perspective of creating attractive distributions for statisticians. Concretely, the M-G class is defined by the following cdf: (2) where G(x) refers to the cdf of a continuous distribution, that we call the parental distribution. en, to demonstrate the importance of the M-G class, in [7], the authors looked at the M uniform, M Rayleigh, M Lomax, M exponential, and M Weibull distributions. e graphs show that the main functions, and the probability density functions (pdfs) and hazard rate functions (hrfs) in particular, have a variety of curvatures, demonstrating their ability to model heterogeneous phenomena. is is shown in [7] using the M Weibull distribution as a parangon of the M-G class and some data by Abouelmagd et al. [8] on aircraft windshield failure times. In fact, it is shown that the M Weibull model adjusts these data better than other extensions of the Weibull models. An extension of the M-G family was proposed by [9] through the transmuted scheme. As an important remark made in [9], one can take α such that α ∈ (∞, 1]/ 0 { }, a condition that significantly improves the modeling ability of the M-G family, among others.
is condition on α will be used throughout the rest of the study. To conclude this paragraph, we mention that another type of M-G family has been established in [10], also based on the T-X transformation.
In this paper, we contribute to the capability of the M distribution to generate flexible distributions by considering the truncated-composed scheme. A list of recent families of distributions based on this scheme is the following: truncated Fréchet generated (TFG) family by Abid and Abdulrazak [11], truncated Weibull generated (TWG) family by Najarzadegan et al. [12], truncated inverted Kumaraswamy generated (TIKG) family by Bantan et al. [13], type II truncated Fréchet generated (TIITFG) family by Aldahlan [14], truncated Cauchy power generated (TCPG) family by Aldahlan et al. [15], exponentiated truncated inverse Weibull generated (ETIWG) family by Almarashi et al. [16], truncated Burr generated (TBG) family by Jamal et al. [17], truncated generalized Fréchet generated (TGFG) family by ZeinEldin et al. [18], truncated inverse Lomax generated (TILG) family by Algarni et al. [19], and truncated Burr X generated (TBX) family by Bantan et al. [20]. Despite this extensive literature, there is no work on what can be called the new truncated M generated (NTM-G) distribution. e main interest of this family is to provide a simplified alternative to the M-G family, while using the functionalities of the M distribution to extend the modeling features of any parental distribution. e following plan is used to illustrate this claim. Section 2 is devoted to a special truncated M distribution belonging to the family of unit distributions. It is new and also is the key ingredient in defining our main general family. We will look at some of its most interesting properties, with a focus on quantitative analysis and moment analysis. e considered TM-G family is defined in Section 3. After the investigation of its related functions, we perform a quantile analysis, followed by a moment analysis. In addition, a numerical work is given on some moment measures. Section 4 is the applied part; we show concretely how the TM-G family can be used to estimate a distribution from data. e maximum likelihood (ML) procedure is employed, and reinforced by a simulation study to guarantee the effectiveness of the method. Real data are then analysed through comparable models, showing that our model is competitive and more accurate than some models derived from the M-G family. A conclusion is given in Section 5.

Motivation.
e creation of unit interval distributions is growing rapidly in the literature. As the main explanation, these distributions are useful for modeling proportions, percentages, and rates that are defined at a unit interval. Modern applications are numerous in the fields of psychology, economics, biology, and engineering. For more information on this topic, see [21][22][23][24] as well as the references therein. is section is dedicated to the unit truncated M distribution, which will be the main ingredient to the NTM-G family. It is worth noting that this unit distribution is not documented in the literature to our knowledge.

Presentation.
Based on the cdf of the M distribution as described in equation (1), we propose a new unit distribution, called the new truncated M (NTM) distribution. It is defined by the following cdf: and the following adjustments at the boundaries: F NTM (x) � 0 for x < 0 and F NTM (x) � 1 for x > 1, that is, for x ∈ (0, 1), we have To our knowledge, there is no reference to the NTM distribution in the literature. e NTM distribution, like the other unit distributions in the abovementioned references, can be used to interpret any proportional-like data. As a result, it offers an alternative to some useful one-parameter unit distributions, such as the power, beta, and the Topp-Leone distributions by Topp and Leone [25].
As a basic property, when α decreases to 0, the NTM distribution has a unit truncated exponential distribution as the limit distribution.
In complement to the cdf, note that the pdf of the NTM distribution is specified by with f NTM (x) � 0 for x < 0 and f NTM (x) � 1 for x > 1, and the hrf is obtained as 2 Complexity with h NTM (x) � 0 for x < 0 and h NTM (x) � 1 for x > 1. e shape behavior of these two functions reveals a flexible model that may be useful for a variety of statistical applications. We, however, omit the details of this aspect and prefer to focus on the quantile analysis and moment analysis of the NTM distribution, which will have important uses for the coming NMTG family.

Quantile
Analysis. e quantile function (qf ) of the NTM distribution is described in the following result. As for the classical M distribution, the Lambert function plays a central role.

Proposition 1.
e qf of the NTM distribution is Proof.
e following equation, F NTM (x) � u with respect to x, is used to determine the qf. is equation is equivalent to and by applying (Corollary 2 in [2]) with "ud α " instead of "u," we get the stated result.
□ e concept of qf has been used in a variety of publications, both theoretical and practical. Essentially, we can use it to define the three main quartiles of a distribution; for the NTM distribution, they are given by Q 1 � Q NTM (1/4), Q 2 � Q NTM (1/2), and Q 3 � Q NTM (3/4). e median of the NTM distribution is simply Q 2 , that is, e qf is also the basis of some distributional measures, such as the skewness measure of Bowley expressed as and the kurtosis measure of Moors is defined by where O i denotes the i th octile of the NTM distribution, i.e., O i � Q NTM (i/8), i � 1, 3, 5 and 7. As a last remark, we can generate values from the NTM distribution by exploiting the fact that, for any random variable U with the unit uniform distribution, Q NTM (U) has the NTM distribution. In other words, using the quantile transformation, values generated by the unit uniform distribution are converted to values created by the NTM distribution. More information on the topic of quantitative analysis can be found in the book of [26].

Moment Analysis.
In this section, we perform a moment analysis. With this aim, we introduce the following special integral function: for a, b, c, and d in such a way that this special integral exists. Connections with well-established integral functions, including gamma and general exponential integral functions, exist. We can also view this integral as a truncated version of the exponential transform of the power-log function u(x) � ( ± x) c (log( ± x)) d . e moment generating function (mgf ) of the NTM distribution is described in the following result. Proposition 2. Let X be a random variable with the NTM distribution. e mgf of X is defined by M(t) � E(e tX ) with t ∈ R; it always exists and can be expressed as where I(a, b, c, d, ± ) is the integral function defined in equation (12).
Proof. By changing the variable y � e αx /α, we have Complexity If α < 0, we can write On the contrary, if α ∈ (0, 1], we can write is ends the proof of the proposition. □ Remark 1. For the case α ∈ (0, 1], by considering the incomplete upper gamma function specified by Γ(a, b) � +∞ b y a− 1 e − y dy for a > 0 and b ≥ 0 and the fol- Proposition 3. Let s be a positive integer and X be a random variable with the NTM distribution. en, the sth ordinary moment of X defined by m(s) � E(X s ) always exists and can be expressed as Proof. We proceed in the same way as the proof of Proposition 2. By making the change of variable y � e αx /α, we have If α < 0, as a result of the binomial formula, we have On the contrary, if α ∈ (0, 1], by invoking the same arguments, we get e proof is now completed. □ Based on Proposition 3, the mean m(1) and standard deviation σ of X can be derived. Also, by using standard relations, we can determine the sth moment of X about the mean given as e moments skewness and kurtosis coefficients of X are given by S � m † (3)/σ 3 and K � m † (4)/σ 4 .
While adjusting the settings of the parameter α, Table 1 displays some numerical values for the measures above.
We can note from Table 1 that, for the considered values only, when α increases, m(1) and S are increased. But K increases when α < 0 and decreases when α > 0. However, the rigorous monotonicity of these measures needs further investigation that we omit here.

NTM-G Family
e NTM-G family is now the focus of attention.

Presentation.
e truncated-composition scheme is applied to define the NTM-G family. By considering a parental cdf denoted by G(x) and the NTM distribution, we define the NTM generated (NTM-G) family by the cdf obtained by the following composition: e main motivations for defining this family are (i) to transpose the functionalities of the NTM distribution to increase those of the parental distribution defined by G(x) and (ii) offer a more simple alternative to the former M-G family introduced in [7]. In some senses, it is more connected with the true nature of the M distribution. e pdf is obtained as where g(x) refers to the pdf derived to G(x). e hrf is given as ese functions are modulable with respect to α, and the definition of the parental distribution is governed by G(x) and g(x). We thus extend the scope of this parental distribution through the use of the NTM strategy.

An Extended Weibull Distribution.
Here, we aim to extend the Weibull distribution through the NTM scheme. To begin, the Weibull distribution with parameters c > 0 and k > 0 is specified by the following cdf and pdf: with G W (x) � 0 for x ≤ 0, and with g W (x) � 0 for x ≤ 0, respectively. By substituting these functions for those defining the NTM-G family, we introduce the NTMW distribution with the following cdf and pdf: and F NTMW (x) � 0 for x < 0, and and f NTMW (x) � 0 for x ≤ 0, respectively. Also, the hrf of the NTMW distribution is given by and h NTMW (x) � 0 for x ≤ 0. us, the NTMW distribution offers a new three-parameter distribution that modifies the analytical structure of the Weibull distribution by adding more tuning parameters. As a visual approach, Figure 1 depicts some plots showing different shapes of the pdf and hrf of the NTMW distribution. From Figure 1, we observe that the pdf of the NTMW distribution has monotonic and nonmonotonic shapes, including diverse decreasing and unimodal shapes, with left and right skewed characteristics. e shape panel of the hrf is quite extensive; we see increasing (concave or convex), decreasing, almost constant, and upside-down shapes. All of these observations demonstrate the flexibility of the NTMW model, which justifies its use in applied statistical scenarios. e pdf and hrf of the MW distribution display in [7] do not have as many shapes. e possible applications of this model will be found in the concrete examples in the applied section of the paper.

Quantile Analysis.
e quantile features of the NTM-G family can be easily obtained using the quantile analysis of the NTM distribution performed in Section 2.3. e qf of the NTM-G family is discussed in detail in the following proposition.

Proposition 4. e qf of the NTM-G family is
where Q G (x) denotes the qf associated to G(x).
Proof. e proof is a consequence of the facts that F NTM− G (x) � F NTM (G(x)), x ∈ R, and the application of Proposition 1.

□
Based on this qf, the comments formulated on the qf for the NTM distribution in Section 2.3 can be transposed.
For the Weibull distribution, we have is qf can be used for further quantile analysis of the NTMW distribution, beyond the scope of this paper.

Moment Analysis.
e moment analysis of the NTM-G family is now being considered. Two approaches are proposed: an integral approach and a series decomposition approach.

Integral Approach.
e following result expresses an integral representation of the ordinary moments of the NTM-G family.
Proposition 5. Let s be a positive integer and X be a random variable with the pdf of the NTM-G family. en, the sth ordinary moment of X can be expressed as provided that it exists.
Proof. We proceed as the proof of Proposition 2. By making the incremental changes of variable y � G(x), then z � e αy /α, and we have e demonstration of the proposition is now complete.

□ 6 Complexity
According to the complexity of Q G (x), the integral in Proposition 5 may have further development. Based on ordinary moments, the well-known moment measures as presented in Section 2.4 can be defined in a similar way; mean m G (1), standard deviation σ, moments skewness coefficient S, and moment kurtosis coefficient K, among others.
In the setting of the NTMW distribution, we have m G (s) � C α αe 1/α c s is integral seems to be not reported in the literature. It can, however, be computed by using mathematical software.

Series Approach.
e following result proposes another approach for the calculation of ordinary moments in the form of a series expansion. is may be of interest for indepth mathematical or computational manipulations.

Proposition 6. Let s be a positive integer and X be a random variable with the pdf of the NTM-G family. Under the condition that sums and integral can be inverted, the sth ordinary moment of X can be expanded as
where A k,ℓ � C α e 1/α (− 1) k α ℓ− k ((k + 2) ℓ − α(k + 1) ℓ )/(k!ℓ!) and Proof. First, we recall that m G (s) � +∞ − ∞ x s f NTM− G (x)dx. We can apply the Taylor series to the exponential function several times after rewriting f NTM− G (x) in a suitable manner, yielding the following expansions in turn: Complexity e sth ordinary moment of X is expanded under the condition that sums and integrals can be inverted, which is primarily determined by the definition of G(x). In this case, we have (38) e stated result is obtained. □ e integral term T ℓ;G (s) can be expressed in the context of the NTMW distribution, and the series expansion becomes explicit. In addition, by taking an integer K large enough, we can approximate m G (s) as m G (s) ≈ K k,ℓ�0 A k,ℓ T ℓ;G (s).

Numerical Study.
We now provide a numerical illustration of the moment analysis. In this regard, values of moment measures of the NTMW distribution are calculated for selected values of the parameters. ese values are presented in Tables 2-4 .
We can see from Tables 2-4 that, for the considered values, when α increases, m(1) tends to increase, but S and K tend to decrease. Additionally, when k increases, S and K tend to decrease.

Applied Study
In this section, we use the ML procedure to suggest some estimates for the unknown parameters of the NTMW model. Simulation experiments are used to evaluate the efficiency of these estimates. en, two real-life datasets are used to suit the proposed distribution, with comparisons to valuable competitors.

ML Procedure.
Let us now describe the ML procedure employed to estimate the NTMW model parameters. We denote by x 1 , . . . , x n some generic data that are supposed to be observed independently from a random variable with a distribution belonging to the NTM-G family. With this data, the total likelihood function is defined by where ζ denotes the possible parameters involved in G(x) and thus in g(x). In an expanded form, the total likelihood function becomes e total log-likelihood function is derived as ℓ(α, ζ) � log[L(α, ζ)], that is, in an expanded form, e ML estimates (MLEs) of (α, ζ) are defined as (α, ζ) � argmax (α,ζ) L(α, ζ) or equivalently (α, ζ) � argmax (α,ζ) ℓ (α, ζ).
us, when ℓ(α, ζ) is differentiable, (α, ζ) can be obtained by solving the following equations simultaneously according to the parameters: zℓ(α, ζ)/zα � 0 and zℓ(α, ζ)/zζ � 0, where with Due to the complex nature of these equations, there are few chances to obtain closed-form expressions for (α, ζ). Numerical techniques such as the (quasi) Newton-Raphson method can be used to find precise numerical solutions. An advantage of using the ML procedure is that, under wellidentified conditions, the underlying asymptotic distribution of (α, ζ) can be approximated by a multivariate normal distribution with a mean vector as (α, ζ) and the following matrix of covariance: [− z 2 ℓ(α, ζ)/zξ zξ t ] − 1 | (α,ζ)�(α,ζ) with ξ � (α, ζ). is matrix can be computed with the help of mathematical software. anks to this standard multivariate distribution, we can construct approximate confidence intervals (CIs) for the model parameters at a certain level, say 8 Complexity 100p% with p ∈ (0, 1); the formula for the associated lower bounds (LBs) and upper bounds (UBs) remains quite standard. All technical details and results of the ML procedure can be found in [27].

Simulation
Work. e methodology of the ML procedure presented above can be applied to the NTMW distribution by taking G(x) � G W (x) and g(x) � g W (x). In this case, ξ � (α, c, k) becomes the unknown vector of parameters.
us, it is natural to think that the MLEs are adequate at estimating the parameters α, c, and k. We illustrate this claim by providing a simulation study to assess the behavior of these MLEs by considering the mean square errors (MSEs), lower bounds (LBs) of CIs, upper bounds (UBs) of CIs, and the corresponding average lengths (ALs) with the levels chosen as 90% or 95%. e results are obtained using the R software (see [28]). Our simulation process is designed as follows :  Tables 5-8 e values in Tables 5-8 for the four considered sets of parameters indicate that when n increases, the MLEs are near equal to the expected values. Also, the MSEs and ALs decrease and tend to 0 as n becomes consequent. Overall, the simulation work found that estimating the model parameters using the ML procedure is appropriate.

Application to Real Datasets.
In this section, the fits of the NTMW model are compared to those of modified Weibull models. Precisely, we consider the beta Weibull (BW) model by Famoye et al. [29], M Weibull (MW) model by Almarashi and Elgarhy [7], exponentiated Weibull (EW) model by Mudholkar and Srivastava [30], and classical M model. e cdf related to these models is described below.
(i) e underlying cdf of the BW model is given by     Complexity 9 and F BW (x) � 0 for x ≤ 0, with α, β, λ, and θ > 0, where B(α, β) denotes the standard beta function, i.e., B(α, β) � 1 0 t α− 1 (1 − t) β− 1 dt (ii) e underlying cdf of the MW model is given by and F MW (x) � 0 for x ≤ 0, with α, β, and c > 0 (iii) e underlying cdf of the EW model is given by and F EW (x) � 0 for x ≤ 0, with α, β, and c > 0 (iv) e underlying cdf of the M model is given by equation (1), with α ∈ (− ∞, 1)/ 0 { } e first real dataset, called dataset I, corresponds to the breaking stress of carbon fibres (in Gba). It was researched in [31]. It contains n � 66 data values. e second dataset, called dataset II, was derived from studies in [8] and represents n � 84 aircraft windshield failure times. e R software is again employed. Tables 9 and 10 show the MLEs for the parameters in each of the five models, along with their standard errors (SEs), for datasets I and II, respectively.
Classically, the estimated negative maximum log-likelihood function (ρ), Akaike information criterion (AIC), Bayesian information criterion (BIC), Cramér von-Mises (W) statistic, Anderson-Darling (A) statistic, and Kolmogorov-Smirnov (KS) statistic with p value (p value) are calculated in Tables 11  and 12 for datasets I and II, respectively. e model with the smallest AIC, BIC, W, A, and KS, as well as the greatest p value for the KS test, is the overall best model using these criteria.
Based on Tables 11 and 12, all of the models appear to be competitive in fitting the datasets, except the M model, which has very small p values for the KS test. e proposed NTMW model, however, outperforms the competitors on all criteria and is the best model for both datasets. Furthermore, the obtained p values are all greater than 0.9, which is quite satisfactory in terms of adequacy.
In order to visualize the efficiency of the modeling, the fitted pdfs and cdfs of each of the five models are plotted in Figures 2 and 3 , for datasets I and II, respectively. e estimated pdfs and cdfs of the NTMW model have particularly well captured the form of the corresponding empirical objects, as shown in Figures 2 and 3. In particular, for the histograms of the datasets, the maximum and tails have received better treatment for the NTMW model than for the competitive models.

Conclusion
e Muth distribution is a versatile distribution that can serve various statistical objectives. In this study, we have briefly discussed the possibility of a unit truncated version of this distribution.
en, we exploited it to construct the new truncated Muth generated (NTM-G) family of distributions. e mathematical properties of this family have been thoroughly examined. We have outlined a special distribution of the NTM-G family, providing a motivated three-parameter modification of the Weibull distribution. It is called the new truncated Muth Weibull (NTMW) distribution. A rule-ofthumb simulation research ensures the effectiveness of the maximum likelihood procedure, which is employed to estimate the model parameters. e flexibility of the model for data fitting has been demonstrated using two realistic datasets. e NTMW model was shown to be the best for the given datasets, when compared to various competing models.
Future works of the proposed family include the applications of the distributions in important applied scenarios, such as those considered in [32,33], a bivariate extended family based on the method of [34], and a discrete version following the approach of [35].

Data Availability
e data used to support the findings of the study are available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.