Next Article in Journal
Klein-Gordon Theory in Noncommutative Phase Space
Next Article in Special Issue
On Statistical Modeling Using a New Version of the Flexible Weibull Model: Bayesian, Maximum Likelihood Estimates, and Distributional Properties with Applications in the Actuarial and Engineering Fields
Previous Article in Journal
Organ Patterning at the Shoot Apical Meristem (SAM): The Potential Role of the Vascular System
Previous Article in Special Issue
The Arcsine Kumaraswamy-Generalized Family: Bayesian and Classical Estimates and Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Amount of Carbon Dioxide Emissions Application: New Modified Alpha Power Weibull-X Family of Distributions

Department of Statistics and Operations Research, Faculty of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(2), 366; https://doi.org/10.3390/sym15020366
Submission received: 5 November 2022 / Revised: 20 December 2022 / Accepted: 14 January 2023 / Published: 30 January 2023

Abstract

:
The use of statistical distributions to model life phenomena has received considerable attention in the literature. Recent studies have shown the potential of statistical distributions in modeling data in applied sciences, especially in environmental sciences. Among them, the Weibull distribution is one of the most well-known models that can be used very effectively for modeling data in the fields of pollution and gas emissions, to name a few. In this paper, we introduce a family of distributions, which we call the modified Alpha-Power Weibull-X family of distributions. Based on the proposed family, we introduce a new model with five parameters, the modified Alpha-Power Weibull–Weibull distribution. Some mathematical properties were determined. Bayesian and maximum likelihood estimates for the model parameters were derived. The MLEs, bootstrap and Bayesian HPD credibility intervals for the unknown parameters were performed. A Monte Carlo simulation study was performed to evaluate the performance of the estimates. A simulation study was performed based on the parameters of the proposed model. An application to the carbon dioxide emissions dataset was performed to predict unique symmetric and asymmetric patterns and illustrate the applicability and potential of the model. For this data set, the proposed model is compared with the modified alpha power Weibull exponential distribution and the two-parameter Weibull distribution. To show which of the competing distributions is the best, we draw on certain analytical tools such as the Kolmogorov–Smirnov test. Based on these analytical measures, we found that the new model outperforms the competing models.

1. Introduction

Increased burning of fossil fuels, deforestation, soil degradation, and various industrial practices have increased carbon dioxide levels in the atmosphere, raising global concerns about climate change and its impact on the environment. In the field of Big Data science and other related fields, the best possible description of real-world phenomena is an important research topic. We refer the reader to [1,2], for more information. Global carbon dioxide emissions from fossil fuels have increased substantially since 1900. Since 1970, carbon dioxide emissions have increased by about 90%, with increasing emissions from fossil fuel combustion and industry accounting for about 78% of total greenhouse gas emissions from 1970 to 2011. Agriculture, deforestation, and other land use changes accounted for the second largest share. We refer the reader to [3,4,5] for more background on recent reports on carbon dioxide emissions.
The quality of the statistical analysis depends on the statistical distribution chosen to model the data. Since various distributions have been used to represent the data and the well-known classical distributions are not sufficient to explain the actual behavior of the data, many transformed, augmented, composite, and mixed distributions have been developed and applied in various fields. However, there are still many important problems that cannot be explained by the current distributions, so we need more flexible and consistent distributions for these problems; see [6,7,8]. One of the most important and recent problems that has piqued our interest is carbon dioxide emissions. For more information on statistical modeling with different models, see [9,10,11].
The problem of complex behavior of data related to carbon dioxide emissions, which do not have a fixed shape or behavior, but rather fluctuate and are unstable. Therefore, we have taken it upon ourselves to find a suitable probability distribution to explain the behavior of carbon dioxide emissions as they are an environmental and health problem.
The objectives of this research are:
  • To introduce a new probability distribution that is more flexible and suitable for modeling real data by adding three additional parameters to the Weibull distribution function.
  • Derive general mathematical properties of the new distribution.
  • Estimate the parameters of the probability distribution of the complete data using the maximum likelihood and Bayesian estimation methods, and compare them using Monte Carlo simulation estimates, biases, and expected errors.
  • Apply experimentally obtained results to the study of carbon dioxide emissions. To show the adequacy of the distribution, it is compared with some other special and standard distributions.
Now, we introduce the suggested family of distributions called the new modified alpha power Weibull-X family. Let ν ( τ ) be the probability density function (PDF) of a random variable (RV) T, where T ( τ 1 , τ 2 ) for < τ 1 < τ 2 < , and let Υ ( G ( x ) ) be a function of G ( x ) , and G ( x ) is the cumulative distribution function (CDF) of an RV X, satisfying the conditions
(1)
Υ ( G ( x ) ) ( τ 1 , τ 2 ) ,
(2)
Υ ( G ( x ) ) is differentiable and monotonically increasing, and
(3)
Υ ( G ( x ) ) τ 1 as x and Υ ( G ( x ) ) τ 2 as x .
Alzaatreh et al. [12] define the cdf of the T-X family of distributions as
F ( x ) = τ 1 Υ ( G ( x ) ) ν ( τ ) d τ , x R
The corresponding PDF is
f ( x ) = ν ( Υ ( G ( x ) ) ) x Υ ( G ( x ) ) , x R .
Based on the alpha power transformation method, the modified alpha power Weibull distribution was proposed by Chettri et al. [13]. The modified alpha power Weibull distribution CDF and PDF are respectively given by
F ( x , α , β , γ ) = α 1 e ( β x ) γ 1 α 1 if α 1 1 e ( β x ) γ if α = 1 ,
and
f ( x , α , β , γ ) = L o g [ α ] α 1 γ β γ x γ 1 e ( β x ) γ α 1 e ( β x ) γ if α 1 γ β γ x γ 1 e ( β x ) γ if α = 1 ,
where x > 0 ; α , β , γ > 0 .
The remainder of this article is organised as follows: Section 2 introduces the new modified Alpha-Power Weibull-X family. Section 3 generates the modified alpha-power Weibull–Weibull distribution. Section 4 discusses the estimation of unknown parameters under the quadratic error, LINEX and the general entropy loss function. A Monte Carlo simulation study is discussed in Section 5. In Section 6, the application of carbon dioxide emissions is carried out to evaluate the efficiency of the proposed model. In Section 7, the discussion and some future frameworks are presented. Finally, brief conclusions are drawn in Section 7.

2. The New Modified Alpha Power Weibull-X Family

Let T e x p ( 1 ) ; then, let its CDF and PDF be given by V ( τ ) = 1 e τ , τ 0 , and ν ( τ ) = e τ , τ 0 , respectively. Considering the modified alpha power Weibull distribution as a generator, we obtain the modified alpha power Weibull X family of distributions by replacing x with G ( x , ζ ) in the alpha power Weibull distribution; this is due to its unique symmetric and asymmetric patterns. Now, let
Υ ( G ( x , ζ ) ) = log 1 α 1 e ( β G ( x , ζ ) ) γ 1 α 1 e β γ 1 if α 1 ( β G ( x , ζ ) ) γ if α = 1 .
where ζ is a vector of parameters. By using (1), we define the cdf of the modified alpha power Weibull (APMW)-X family by
F ( x , α , β , γ , ζ ) = α 1 e ( β G ( x , ζ ) ) γ 1 α 1 e β γ 1 if α 1 1 e ( β G ( x , ζ ) ) γ if α = 1 .
The APMW-X density function is
f ( x , α , β , γ , ζ ) = γ β γ L o g [ α ] e ( β G ( x , ζ ) ) γ g ( x , ζ ) α e β γ e ( β G ( x , ζ ) ) γ G ( x , ζ ) γ 1 if α 1 γ β γ G ( x , ζ ) γ 1 e ( β G ( x , ζ ) ) γ g ( x , ζ ) if α = 1 ,
The survival and hazard functions of APMW-X are given, respectively, by
S ( x , α , β , γ , ζ ) = 1 α 1 e ( β G ( x , ζ ) ) γ 1 α 1 e β γ 1 if α 1 e ( β G ( x , ζ ) ) γ if α = 1 ,
and
H ( x , α , β , γ , ζ ) = γ β γ L o g [ α ] e ( β G ( x , ζ ) ) γ g ( x , ζ ) α 1 e ( β G ( x , ζ ) ) γ ( G ( x , ζ ) ) γ 1 α 1 e β γ 1 1 α 1 e ( β G ( x , ζ ) ) γ 1 α 1 e β γ 1 if α 1 γ β γ G ( x , ζ ) γ 1 g ( x , ζ ) if α = 1 ,
The main motivations for using the APMW-X family in practice are the following:
(i)
To develop the flexibility and properties of the basic models;
(ii)
To provide a suitable procedure for adding additional parameters in extended models with strong outliers, which are very useful in gas emission modeling;
(iii)
Introduce the extended version of a basic model with closed forms for the cdf and hazard rate function, where the special submodels of this family can be used in the analysis of censored data sets;
(iv)
Compared to existing competing models, the special cases of the APMW-X approach are able to model data sets with high tail content.
Figure 1 plots different SF for the APMW-X ( α , β , γ ) family.

2.1. Quantile Function

Supposing p U n i f o r m ( 0 , 1 ) , we have solved the following equation for the quantile function Q ( p ) :
p = α 1 e ( F ( Q ( p ) ) β ) γ 1 α 1 e β γ 1 if α 1 1 e ( β F ( Q ( p ) ) ) γ if α = 1 .
Letting y = F ( Ω ( p ) ) , we have
p = α 1 e ( y β ) γ 1 α 1 e β γ 1 if α 1 1 e ( β y ) γ if α = 1 .
By solving Equation (11) for y, we have
F ( Q ( p ) ) = 1 β L o g 1 L o g 1 + p α 1 e β γ 1 L o g [ α ] 1 γ if α 1 1 β L o g 1 p 1 γ if α = 1 .
Thus,
Q ( p ) = F 1 1 β L o g 1 L o g 1 + p α 1 e β γ 1 L o g [ α ] 1 γ if α 1 F 1 1 β L o g 1 p 1 γ if α = 1 .
where F 1 is the inverse cumulative of the baseline distribution.

2.2. The Likelihood Function of the APMW-X Family

Let { X i , i = 1 , , n } be the observed sample. The likelihood function based on the APMW-X ( α , β , γ , ζ ) for α 1 is given by
l = β γ γ L o g [ α ] α 1 e β γ 1 n i = 1 n e β G x i , ζ γ g x i , ζ α 1 e β G x i , ζ γ G x i , ζ γ 1 .
The corresponding log-likelihood function is given by
L = n L o g [ β γ L o g [ α ] ] n L o g α 1 e β γ 1 i = 1 n β G x i , ζ γ + i = 1 n L o g g x i , ζ + i = 1 n 1 e β G x i , ζ γ L o g [ α ] + ( γ 1 ) i = 1 n L o g β G x i , ζ .
Let G ( x i , ζ ) = G i and g ( x i , ζ ) = g i be the baseline CDF and PDF, respectively. The first partial derivatives of (15) with respect to α , β , γ are given by
L α = 1 e β γ n α e β γ α 1 e β γ 1 + n α L o g [ α ] + i = 1 n 1 e β G i γ α ,
L β = n β + n ( γ 1 ) β e β γ n α 1 e β γ β γ 1 γ L o g [ α ] α 1 e β γ 1 i = 1 n γ G i β G i γ 1 + i = 1 n e β G i γ γ L o g [ α ] G i β G i γ 1 ,
and
L γ = n γ e β γ n α 1 e β γ β γ L o g [ α ] L o g [ β ] α 1 e β γ 1 + i = 1 n L o g β G i i = 1 n L o g β G i β G i γ + i = 1 n e β G i γ L o g [ α ] L o g β G i β G i γ .

3. The Modified Alpha Power Weibull–Weibull (APMW-W) Distribution

The aim of this section is to propose a new composite probability distribution model (the modified alpha power Weibull–Weibull distribution) using the alpha power Weibull-X family to obtain a more convenient and flexible distribution for modeling observations. The main motive for studying and applying the APMW-X method to the Weibull distribution is the following:
  • The APMW method is an effective way to add more than three parameters to the distribution family;
  • The APMW method makes the distribution richer and more flexible;
  • The APMW method provides models that can model both monotonic and non-monotonic hazard rate function (HRF);
  • The APMW method gives us a better fit than other modified models with the same or fewer parameters.
Let X be a random variable (R.V.) that follows the two-parameters Weibull distribution ( λ , μ ) , then its CDF, denoted by G ( x ; λ , μ ) , is given by
G ( x ; λ , μ ) = 1 e x μ λ , x 0 ; λ , μ > 0 .
Here, x 0 , λ > 0 and μ > 0 are the shape and the scale parameter, respectively. The corresponding PDF, denoted by g ( x ; λ , μ ) , is given by
g ( x ; λ , μ ) = λ μ x μ λ 1 e x μ λ , x 0 ; λ , μ > 0 .
by taking G ( x ) and g ( x ) to be G ( x ; λ , μ ) and g ( x ; λ , μ ) , respectively. The CDF and PDF of the APMW-W distribution are given, respectively, by
F ( x ; α , β , γ , λ , μ ) = α 1 e x p β 1 e x μ λ γ 1 α 1 e β γ 1 if α 1 1 e x μ λ if α = 1 ,
f ( x ; α , β , γ , λ , μ ) = γ λ β γ L o g [ α ] x α 1 e β γ 1 e β 1 e x μ λ γ x μ λ α 1 exp β 1 e x μ λ γ × 1 e x μ λ γ 1 x μ λ if α 1 λ μ e x μ λ x μ λ 1 if α = 1 ,
The survival and hazard functions of APMW-W are given, respectively, by
S ( x ; α , β , γ , λ , μ ) = 1 α 1 e x p β 1 e x μ λ γ 1 α 1 e β γ 1 if α 1 e x μ λ if α = 1 ,
H ( x ; α , β , γ , λ , μ ) = γ λ L o g [ α ] x μ λ exp β 1 e x μ λ γ α e β γ β 1 e x μ λ γ x e x μ λ 1 α e β γ α exp β 1 e x μ λ γ if α 1 λ μ x μ λ 1 if α = 1 ,
Figure 2, Figure 3 and Figure 4 plot different SF, PDF and HRF, respectively, for the APMW-W ( α , β , γ , λ , μ ) .
By substituting the CDF of the APMW-W distribution (21) in (13), the quantile function Ω ( p ) is
Q ( p ) = F 1 1 β L o g 1 L o g 1 + p α 1 e β γ 1 L o g [ α ] 1 γ if α 1 F 1 1 β L o g 1 p 1 γ if α = 1 .
if F 1 is the quantile of the baseline distribution APMW-W ( α , β , γ , λ , μ ) . We can write
Ω ( p ) = μ L o g β β L o g 1 L o g 1 + p α 1 e β γ 1 L o g [ α ] 1 γ 1 1 λ .
In addition, the effects of shape parameters on skewness and kurtosis can be determined using quantile measures. We obtain skewness and kurtosis measures of APMW-W ( α , β , γ , λ , μ ) . The skewness measure (SK) of APMW-W ( α , β , γ , λ , μ ) (see Bowley [14]) of X is given by
S K ( α , β , γ , λ , μ ) = 2 Q 1 / 2 Q 3 / 4 Q 1 / 4 Q 1 / 4 Q 3 / 4 ,
and the kurtosis (K) (see Moor [15]) is given by
K ( α , β , γ , λ , μ ) = Q 1 / 8 Q 3 / 8 + Q 5 / 8 Q 7 / 8 Q 2 / 8 Q 6 / 8 .
Some quantile values for α = 3.1 , β = 2.2 , γ = 1 , λ = 0.5 , and μ = 1 are shown in Figure 5. The skewness obtained is 0.592713, and the kurtosis is 2.19598 for the same case shown in Figure 5. Table 1 shows some quantile values for the same case. Figure 6, Figure 7 and Figure 8 show the SK and K for some cases of APMW-W ( α , β , γ , λ , μ ) .

4. Estimation of the Parameters

We obtain estimators of the model parameters of the APMW-W ( α , β , γ , λ , μ ) distribution in this section.

4.1. The Maximum Likelihood Estimation

Here, we discuss the maximum likelihood estimators (MLEs) of the model parameters of the APMW-W ( α , β , γ , λ , μ ) distribution. The first partial derivatives of (15) with respect to λ , μ are given by
L λ = i = 1 n g i ( λ ) g i + ( γ 1 ) i = 1 n G i ( λ ) G i i = 1 n γ β γ G i γ 1 G i ( λ ) + γ L o g [ α ] β γ i = 1 n e β G i γ G i γ 1 G i ( λ ) ,
and
L μ = i = 1 n g i ( μ ) g i + ( γ 1 ) i = 1 n G i ( μ ) G i i = 1 n γ β γ G i γ 1 G i ( μ ) + γ L o g [ α ] β γ i = 1 n e β G i γ G i γ 1 G i ( μ ) ,
where
G i λ = e x μ λ x μ λ L o g x μ ,
G i μ = λ μ 2 x e x μ λ x μ λ 1 ,
g i λ = 1 x e x μ λ x μ λ 1 λ 1 + x μ λ L o g x μ ,
and
g i μ = λ 2 μ 2 e x μ λ 1 + x μ λ x μ λ 1 ,
The MLEs of the parameters α , β , γ , λ , and μ are obtained by equating Equations (16)–(18), (26) and (27) to zero and solving the above equation simultaneously. However, it is difficult to solve these equations to obtain the estimates of the unknown parameters in explicit form. Therefore, a numerical technique can be used to solve these nonlinear equations.

4.2. Bayesian Estimation

Bayesian inference is a suitable method to work with the full samples of APMW-W ( α , β , γ , λ , μ ) . Prior predictive distributions can be used to check the reasonableness of a prior for a given situation before observing sample data. Gamma distribution is one of the most commonly used distributions as a pre-distribution and gives good experimental results. We assume that α , β , γ , λ , and μ are R.V.s that follow the prior PDFs Gamma ( α ; a 1 , b 1 ) , Gamma ( β ; a 2 , b 2 ) , Gamma ( γ ; a 3 , b 3 ) , Gamma ( λ ; a 4 , b 4 ) , and Gamma ( μ ; a 5 , b 5 ) , respectively. Then, the posterior density of α , β , λ , μ , ν and the data are given by
l = β γ γ L o g [ α ] α 1 e β γ 1 n i = 1 n e β G x i , ζ γ g x i , ζ α 1 e β G x i , ζ γ G x i , ζ γ 1 .
π * ( α , β , γ , λ , μ | x ) = J 1 L o g [ α ] α 1 e β γ 1 n α n + a 1 1 i = 1 n e b 1 α + b 2 β + b 3 γ + b 4 λ + b 5 μ β G x i , ζ γ × β γ n + a 2 1 γ n + a 3 1 λ a 4 1 μ a 5 1 e i = 1 n β G x i , ζ γ i = 1 n g x i , ζ G x i , ζ γ 1 ,
where J is the normalizing constant.

5. Monte Carlo Simulation Study

This section is concerned with evaluating the performance of the maximum likelihood and Bayesian estimators of the APMW-W ( α , β , γ , λ , μ ) distribution through a Monte Carlo simulation study. The simulation is performed for the parameters α = 1.025 , β = 2 , γ = 3.27 , λ = 0.52 , and μ = 1.16 of the APMW-W distribution model.

5.1. MLE Monte Carlo Simulation

The simulation study is conducted as follows.
  • Random samples of size n = 25 , 50 , , 400 are generated from the APMW-W ( 1.025 , 2 , 3.27 , 0.52 , 1.16 ) distribution.
  • Model parameters were estimated using the maximum likelihood method;
  • One-thousand replicates were performed to calculate the biases and expected errors (ERs) of these estimators;
  • The formulas used to calculate the estimate, biases, and ERs are as follows:
    α ^ = 1 1000 i = 1 1000 α ^ i ,
    B i a s α ^ = 1 1000 i = 1 1000 α ^ i α ,
    and
    E R α ^ = 1 1000 i = 1 1000 α ^ i α 2 ;
  • Step (4) is also repeated for the parameters β , γ , λ , and μ .

5.2. The Bootstrap Confidence Intervals: Boot-p Algorithm

Next, obtain the bootstrap confidence intervals for boot-p for the unknown parameters δ = ( α , β , γ , λ , μ ), we apply the following algorithms:
  • Generate sample { x i } of size n from the APMW-W ( α , β , γ , λ , μ ) and estimate a δ ^ ;
  • Generate another sample { x i * } of size n using δ ^ . Then, estimate δ ^ * ;
  • Repeat step 2 B times;
  • Via F ^ ( x ) = P ( δ ^ * x ) , that is, the CDF of δ ^ * , the 100 ( 1 ε ) % C.I. of δ is given by
    δ ^ B o o t p ( ϵ 2 ) , δ ^ B o o t p ( 1 ϵ 2 ) ,
    where δ ^ B o o t p ( κ ) = F ^ 1 ( κ ) and x is prefixed.
For more details about the bootstrap confidence intervals, one may refer to Kundu and Joarder [16].

5.3. Bayesian Monte Carlo Simulation Study

We assume that α , β , γ , λ , and μ has the prior PDFs Gamma ( α ; 0.3 , 0.5 ) , Gamma ( β ; 0.09 , 0.01 ) , Gamma ( γ ; 0.3 , 0.3 ) ,Gamma ( λ ; 0.2 , 0.8 ) and Gamma ( μ ; 0.9 , 0.3 ) , respectively. We use the Metropolis–Hastings procedure as:
  • Set start values α ( 0 ) = 1.025 , β ( 0 ) = 2 , γ ( 0 ) = 3.27 , λ ( 0 ) = 0.52 and μ ( 0 ) = 1.16 . Then, simulate sample of size n from A P M W W ( α ( 0 ) , β ( 0 ) , γ ( 0 ) , λ ( 0 ) , μ ( 0 ) ) , next set l = 1 ;
  • Simulate α ( * ) , β ( * ) , γ ( * ) , λ ( * ) and μ ( * ) . using the proposal distributions N ( α ( l 1 ) , V ( α ^ ) ) , N ( β ( l 1 ) , V ( β ^ ) ) , N ( γ ( l 1 ) , V ( λ ^ ) ) , N ( λ ( l 1 ) , V ( λ ^ ) ) and N ( μ ( l 1 ) , V ( μ ^ ) ) ;
  • Calculate r = m i n π * ( α ( * ) , β ( * ) , γ ( * ) , λ ( * ) , μ ( * ) ) π * ( α ( l 1 ) , β ( l 1 ) , γ ( l 1 ) , λ ( l 1 ) , μ ( l 1 ) ) , 1 ;
  • Simulate U from Uniform(0, 1);
  • If U < r , then α ( l ) , β ( l ) , γ ( l ) , λ ( l ) , μ ( l ) = α ( * ) , β ( * ) , γ ( * ) , λ ( * ) , μ ( * ) ;
    If U r , then α ( l 1 ) , β ( l 1 ) , γ ( l 1 ) , λ ( l 1 ) , μ ( l 1 ) = α ( * ) , β ( * ) , γ ( * ) , λ ( * ) , μ ( * ) ;
  • Set l = l + 1 ;
  • Iterate Steps 2–6, M repetitions, and obtain α ( l ) , β ( l ) , γ ( l ) , λ ( l ) and μ ( l ) for l = 1 , , M .
Suppose the squared error loss function, given by L S E ( δ , δ ^ ) = ( δ δ ^ ) 2 , by using the generated random samples from the M-H technique, and N is the nburn. Then, the Bayes estimator of δ against the squared error S E loss function is given by
δ ^ S E = E δ δ | x = 1 M N l = N + 1 M δ ( l ) .
Next, suppose the LINEX ( L E ) loss function, given by
L L E δ , δ ^ = exp ρ δ δ ^ ρ δ δ ^ 1 , ρ 0 .
The approximate Bayes estimate of δ under L E loss function is given by
δ ^ L E = 1 ρ log E δ exp ρ δ | x = 1 ρ log l = N + 1 M e x p ρ δ ( l ) M N .
The parameter ρ in LINEX is chosen as 0.2 ( L E 1 ), and 0.8 ( L E 2 ). Finally, suppose the general entropy (GE) loss function, given by
L G E δ , δ ^ = δ ^ δ ε ε log δ ^ δ 1 .
The parameter ε in GE is chosen as 0.6. The approximate Bayes estimate of the parameters, given by
δ ^ G E = E δ δ ε x 1 ε = 1 M N l = N + 1 M δ ( l ) ε 1 ε ,

5.4. MCMC HPD Credible Interval Algorithm

  • Arrange α ( * ) , β ( * ) , γ ( * ) , λ ( * ) and μ ( * ) in rising values;
  • The lower bounds of α , β , γ , λ , and μ are in the rank ( M N ) ε / 2 ;
  • The upper bounds of α , β , γ , λ , and μ is in the rank ( M N ) ( 1 ε / 2 ) ;
  • Iterate the previous steps M times. Obtain the average value of the lower and upper bounds of α , β , γ , λ , and μ .
The point and interval simulation results of the APMW-W distribution for α = 1.025 , β = 2 , γ = 3.27 , λ = 0.52 , and μ = 1.16 are, respectively, presented in Table 2 and Table 3.

6. The Carbon Dioxide Emissions Application

The data set of 50 carbon dioxide emissions for the period (1970–2019) given by Albank Aldawli for Saudi Arabia is considered as an application of the APMW-W distribution. Table 4 shows the descriptive statistics of the data for carbon dioxide emissions ( × 10 5 ) . The boxplot and Q-Q plot are shown in Figure 9. Figure 10 shows the fitted PDF of APMW-W and CDF. Figure 11 shows the PP plot and the Kaplan–Meier survival function of APMW-W.
The goodness-of-fit of APMW-W is compared with some other models, including the modified alpha-power Weibull exponential distribution (APMW-E = A P M W X ( α , β , γ , λ , 1 ) ) and the two-parameter Weibull distribution (TW-D ( λ , μ ) ) (Equation (19)). The estimated values of the APMW-W and the competing models for the given data set of 50 carbon dioxide emissions are shown in Table 4. The distribution functions of these competitive distributions are given by:
  • APMW-W distribution:
    f ( x ; α , β , γ , λ , μ ) = γ λ β γ L o g [ α ] x α 1 e β γ 1 e β 1 e x μ λ γ x μ λ α 1 exp β 1 e x μ λ γ × 1 e x μ λ γ 1 x μ λ ,
  • APMW-E distribution:
    f ( x , α , β , γ , μ ) = γ β γ μ L o g [ α ] e β 1 e x μ γ e x μ α e β γ e β 1 e x μ γ 1 e x μ γ 1 ,
  • TW-d distribution:
    f ( x ; λ , μ ) = λ μ x μ λ 1 e x μ λ , x 0 ; λ , μ > 0 .
Table 5 shows the estimated parameter values of APMW-W and the competing models. In Table 6, the Kolmogorov–Smirnov test is performed.

7. Discussion and Future Framework

The addition of the three parameters to the family shows a significant effect on the diversity of SF as in Figure 1. The Weibull model became very flexible after the APMW-X family parameters were added. Sometimes, it resembles a bell curve with some torsion and at other times it seems to have strong swings as seen in Figure 3, which depends on the specific values of the parameters. The proposed model is a good candidate for data modeling in various financial, industrial, medical and other applications. However, it can be seen from Figure 3 that APMW-W has simple features and an elastic failure rate. The simple hazard rate and flexible features are another superiority of the proposed model along with its heavy-tail behavior.
The Kolmogorov–Smirnov (KS) statistics for one sample with p-values are given in Table 6. From the results in Table 6, it can be seen that the APMW-W model could be selected as the best model among the fitted models. In the future, it is possible to expand the scientific aspects associated with the application, such as the medical, technical, and industrial aspects.
We empirically show that the new five-parameter expansion of the Weibull distribution provides the best fit to the carbon dioxide emission data than the competing distributions. The practical example shows that the proposed model is a suitable alternative distribution for modeling carbon dioxide emission data.

8. Conclusions

The main objective of this study is to instruct a new flexible modification of the Weibull model by introducing three additional parameters. The introduction of the additional parameters leads to greater flexibility to improve the goodness of fit to the reliability data. We determined the maximum likelihood estimators for the intended model parameters and performed a Bayesian Monte Carlo simulation study. The performance of the Bayesian estimators is better than that of the corresponding ML estimators. The expected errors support the Bayesian estimator in most cases. The width of the intervals of the Bayesian estimator is shorter than that of the maximum likelihood estimator at the same confidence level. The loss functions of LINEX and general entropy behave better and are close in terms of variances. The biases and mean squared errors decrease with increasing sample size. It is clear that the proposed model fits well with the estimated PDF and CDF plots. The proposed model fits the Kaplan–Meier survival plot very well. Based on the Kolmogorov–Smirnov one-sample test, the new model provides a better fit than other competing models.

Author Contributions

Formal analysis, W.E.; Funding acquisition, Y.T.; Methodology, W.E.; Software, W.E.; Supervision, W.E.; Writing—original draft, W.E.; Writing—review and editing, Y.T. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Researchers Supporting Project number (RSP2023R488), King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

All the datasets used in this paper are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gomez-Deniz, E.; Calder in-Ojeda, E. On the usefulness of the logarithmic skew normal distribution for describing claims size data. Math-Ematical Probl. Eng. 2020, 2020, 1420618. [Google Scholar] [CrossRef]
  2. Ahmad, Z.; Mahmoudi, E.; Hamedani, G. A class of claim distributions: Properties, characterizations and applications to insurance claim data. Commun. Stat.-Deory Methods 2020, 49, 2183–2208. [Google Scholar] [CrossRef]
  3. Emissions due to Agriculture Global, Regional and Country Trends 2000–2018. Available online: https://www.fao.org (accessed on 15 September 2021).
  4. 2019 UK Greenhouse Gas Emissions, Final Figures. Available online: https://assets.publishing.service.gov.uk (accessed on 15 September 2021).
  5. Shafiq, A.; Lone, S.A.; Sindhu, T.N.; Khatib, Y.E.; Al-Mdallal, Q.M.; Muhammad, T. A new modified Kies Fréchet distribution: Applications of mortality rate of COVID-19. Results Phys. 2021, 28, 104638. [Google Scholar] [CrossRef] [PubMed]
  6. Ahmad, Z.; Mahmoudi, E.; Dey, S. A new family of heavy tailed distributions with an application to the heavy tailed insurance loss data. Commun. Stat.-Simul. Comput. 2020, 49, 4372–4395. [Google Scholar] [CrossRef]
  7. Klakattawi, H.S.; Aljuhani, W.H. A new technique for generating distributions based on a combination of two techniques: Alpha power transformation and exponentiated TX distributions family. Symmetry 2021, 13, 412. [Google Scholar] [CrossRef]
  8. Mansoor, M.; Tahir, M.H.; Cordeiro, G.M.; Alzaatreh, S.B. The Marshall-Olkin logistic-exponential distribution. Commun. Stat.-Theory Methods 2019, 48, 220–234. [Google Scholar] [CrossRef]
  9. El-Khatib, Y.; Hatemi-J, A. Computations of price sensitivities after a financial market crash. In Electrical Engineering and Intelligent Systems; Springer: New York, NY, USA, 2013; pp. 239–248. [Google Scholar]
  10. El-Khatib, Y.; Al-Mdallal, Q.M. Numerical simulations for the pricing of options in jump diffusion markets. Arab. J. Math. Sci. 2012, 18, 199–208. [Google Scholar] [CrossRef] [Green Version]
  11. Foss, S.; Korshunov, D.; Zachary, S. An Introduction to Heavy-Tailed and Subexponential Distributions; Springer Science & Business Media: Berlin, Germany, 2013. [Google Scholar]
  12. Alzaatreh, A.; Lee, C.; Famoye, F. A new method for generating families of continuous distributions. Metron 2013, 71, 63–79. [Google Scholar] [CrossRef] [Green Version]
  13. Chettri, S.; Das, B.; Chakraborty, S. A New Modified Alpha Power Weibull Distribution: Properties, Parameter Estimation and Application. J. Indian Soc. Probab. Stat. 2021, 22, 417–449. [Google Scholar] [CrossRef]
  14. Bowley, A.L. Elements of Statistics, 4th ed.; Charles Scribner’s Sons: New York, NY, USA, 1920. [Google Scholar]
  15. Moors, J. The meaning of kurtosis: Darlington re-examined. Am. Stat. 1986, 40, 283–284. [Google Scholar]
  16. Kundu, D.; Joarder, A. Analysis of Type-II progressively hybrid censored data. Comput. Stat. Data Anal. 2006, 50, 2509–2528. [Google Scholar] [CrossRef]
Figure 1. Different SF for the APMW-X ( α , β , γ ) .
Figure 1. Different SF for the APMW-X ( α , β , γ ) .
Symmetry 15 00366 g001
Figure 2. Different SF for the APMW-X ( α , β , γ , λ , μ ) .
Figure 2. Different SF for the APMW-X ( α , β , γ , λ , μ ) .
Symmetry 15 00366 g002
Figure 3. Different PDF for the APMW-X ( α , β , γ , λ , μ ) .
Figure 3. Different PDF for the APMW-X ( α , β , γ , λ , μ ) .
Symmetry 15 00366 g003
Figure 4. Different HRF for the APMW-X ( α , β , γ , λ , μ ) .
Figure 4. Different HRF for the APMW-X ( α , β , γ , λ , μ ) .
Symmetry 15 00366 g004
Figure 5. Different quantile functions for the APMW-X ( α , β , γ , λ , μ ) .
Figure 5. Different quantile functions for the APMW-X ( α , β , γ , λ , μ ) .
Symmetry 15 00366 g005
Figure 6. Plots for the S K ( α , β , 1 , 0.5 , 1.3 , 0.1 ) and K ( α , β , 1 , 0.5 , 1.3 , 0.1 ) .
Figure 6. Plots for the S K ( α , β , 1 , 0.5 , 1.3 , 0.1 ) and K ( α , β , 1 , 0.5 , 1.3 , 0.1 ) .
Symmetry 15 00366 g006
Figure 7. Plots for the S K ( 3.1 , 2.2 , λ , μ , 1.3 , 0.1 ) and K ( 3.1 , 2.2 , λ , μ , 1.3 , 0.1 ) .
Figure 7. Plots for the S K ( 3.1 , 2.2 , λ , μ , 1.3 , 0.1 ) and K ( 3.1 , 2.2 , λ , μ , 1.3 , 0.1 ) .
Symmetry 15 00366 g007
Figure 8. Plots for the S K ( 3.1 , 2.2 , 1 , 0.5 , μ ) and K ( 3.1 , 2.2 , 1 , 0.5 , μ ) .
Figure 8. Plots for the S K ( 3.1 , 2.2 , 1 , 0.5 , μ ) and K ( 3.1 , 2.2 , 1 , 0.5 , μ ) .
Symmetry 15 00366 g008
Figure 9. The boxplot and Q-Q plot of the carbon dioxide emissions data.
Figure 9. The boxplot and Q-Q plot of the carbon dioxide emissions data.
Symmetry 15 00366 g009
Figure 10. The fitted PDF and CDF of the APMW-W distribution.
Figure 10. The fitted PDF and CDF of the APMW-W distribution.
Symmetry 15 00366 g010
Figure 11. The PP plot and the Kaplan–Meier survival function of the APMW-W distribution.
Figure 11. The PP plot and the Kaplan–Meier survival function of the APMW-W distribution.
Symmetry 15 00366 g011
Table 1. Some quantile values for α = 3.1 , β = 2.2 , γ = 1 , λ = 0.5 , and μ = 1 .
Table 1. Some quantile values for α = 3.1 , β = 2.2 , γ = 1 , λ = 0.5 , and μ = 1 .
xQ(x)
0.10.0051
0.20.0223
0.30.0557
0.40.1125
0.50.2062
0.60.3630
0.70.6420
0.81.2066
0.92.7380
Table 2. Point estimation of the APMW-W parameters.
Table 2. Point estimation of the APMW-W parameters.
Point
nPar. ML SE LE 1 LE 2 GE
25 α 0.96920.75060.76860.68260.5059
−0.0562−0.2748−0.2568−0.3428−0.5196
0.28380.37680.37550.3920.5872
β 1.53031.76431.77641.71691.6556
−0.4718−0.2378−0.2257−0.2852−0.3465
0.07980.30550.29130.36320.524
γ 1.51553.25143.25283.24593.248
−1.7589−0.023−0.0216−0.0284−0.0264
0.01960.04540.04530.0460.046
λ 1.72760.58050.58130.57750.5704
1.20590.05880.05960.05580.0487
0.01680.0080.00820.00760.0068
μ 1.1181.15961.16821.12571.0898
−0.04040.00120.0098−0.0327−0.0686
0.06130.31220.31450.30410.3352
50 α 1.02830.74310.77160.64620.5068
0.0029−0.2823−0.2538−0.3792−0.5186
0.13590.53820.56510.48160.6411
β 1.60512.02792.03252.00952.0089
−0.39690.02580.03050.00750.0069
0.13040.11040.11090.10930.1118
γ 1.41533.15053.15723.12353.1312
−1.8591−0.1238−0.1171−0.1508−0.1431
0.00690.29890.29470.3170.316
λ 1.54910.61180.61240.60910.6038
1.02740.090.09070.08740.0821
0.02760.02250.02270.02150.0202
μ 1.10251.19481.20221.1651.1439
−0.05590.03650.04380.0066−0.0145
0.04380.26130.26580.24470.253
100 α 0.89420.6460.65870.59860.4777
−0.1312−0.3794−0.3667−0.4268−0.5477
0.08810.31930.31510.33850.4707
β 1.65461.93591.93921.92281.9196
−0.3475−0.0661−0.0628−0.0792−0.0824
0.01490.09430.0930.09920.1072
γ 1.48173.34433.34743.33223.337
−1.79260.070.0730.05780.0627
0.01360.11980.12050.11720.1188
λ 1.34930.56990.57020.5690.5667
0.82760.04820.04850.04730.045
0.0060.00650.00650.00630.0061
μ 1.11811.29641.30431.26571.2507
−0.04020.1380.14590.10730.0924
0.02680.21390.22070.18980.1938
200 α 1.04921.08781.09961.04030.9805
0.02380.06240.07420.0149−0.0449
0.06270.35920.36850.32490.3654
β 1.76631.91871.92111.90931.9087
−0.2357−0.0833−0.0809−0.0928−0.0933
0.03090.04560.04540.04660.0474
γ 1.46983.19433.19713.18343.1869
−1.8046−0.0801−0.0773−0.091−0.0875
0.00910.08380.08230.08960.0883
λ 1.25580.57410.57430.57340.5717
0.73410.05240.05260.05170.05
0.00260.00580.00590.00570.0054
μ 1.13571.00351.00710.98980.9744
−0.0227−0.1548−0.1513−0.1686−0.184
0.02220.1160.11540.11890.1284
400 α 0.94481.08081.08651.05831.0335
−0.08060.05540.06110.03290.0081
0.02990.17580.17760.16930.1815
β 1.75592.02582.02682.02172.0218
−0.24620.02370.02470.01960.0197
0.01630.02760.02770.02710.0272
γ 1.45013.24643.24793.24063.2427
−1.8243−0.0279−0.0265−0.0338−0.0317
0.00590.04130.0410.04240.0421
λ 1.21780.57070.57080.57040.5697
0.69610.0490.04910.04870.048
0.00150.00370.00370.00360.0035
μ 1.1811.14921.15151.14031.1341
0.0227−0.0092−0.0069−0.0181−0.0243
0.01360.06020.06060.05870.0593
The first line represents estimate, the second line represents bias, and the third line represents ER.
Table 3. Interval estimation of the APMW-W parameters.
Table 3. Interval estimation of the APMW-W parameters.
nPar. ML Boot HPD S HPD LE 1 HPD LE 2 HPD GE
25 α −0.075 2.01340.001 4.9810.112 2.010.1122 2.0290.1116 1.94090.0054 1.9402
2.08834.981.8981.91681.82931.9348
0.0902 1.84820.003 3.9440.144 1.8290.1476 1.85490.1269 1.81640.0233 1.8152
1.7583.9411.6851.70741.68951.7919
β 0.9767 2.08390.0001 6.2830.549 2.4760.5829 2.48870.4167 2.42140.0444 2.4298
1.10726.2831.9271.90582.00472.3854
1.0642 1.99630.001 5.1320.738 2.3980.7632 2.40720.5393 2.35570.2286 2.3607
0.93215.1311.661.6441.81652.1321
γ 1.241 1.790.177 5.7822.762 3.642.7673 3.64482.7296 3.62032.7399 3.6292
0.54915.6050.8780.87750.89070.8893
1.2844 1.74660.227 5.0462.846 3.5792.8469 3.57932.8438 3.57782.8445 3.5783
0.46224.8190.7330.73250.7340.7338
λ 1.4735 1.98180.491 5.1990.46 0.7270.4607 0.7280.4591 1.5895099780.4549 0.7174
0.50834.7080.2670.26730.26460.2625
1.5137 1.94160.529 4.3340.475 0.7060.4751 1.5286108040.4744 0.69630.4721 0.6919
0.42793.8050.2310.23360.22190.2199
μ 0.6329 1.60310.101 3.7980.237 2.2970.2379 2.31190.2321 2.23520.2029 2.2408
0.97023.6972.062.0742.00312.0378
0.7096 1.52640.136 2.8780.306 2.0960.3085 2.10780.302 2.03620.2356 2.0282
0.81682.7421.791.79931.73421.7926
50 α 0.3058 1.75080.001 5.4050.1 2.560.1023 2.57890.0943 2.36450.0226 2.3963
1.44515.4042.462.47662.27022.3737
0.42 1.63650.003 3.9540.123 2.3310.1247 2.42120.108 1.88940.0336 1.8826
1.21653.9512.2082.29651.78141.849
β 0.8974 2.31280.0001 6.3661.34 2.6981.3499 2.70341.298 2.6781.273 2.6831
1.41546.3661.3581.35351.381.4101
1.0094 2.20090.0001 5.4061.511 2.5871.5126 2.59021.5067 2.57461.4919 2.5774
1.19155.4061.0761.07761.06791.0855
γ 1.2528 1.57780.221 5.0291.838 4.1521.8567 4.16211.7721 4.11351.7663 4.1333
0.3254.8082.3142.30542.34142.367
1.2785 1.55210.306 4.2232.217 4.0122.2498 4.01782.0887 3.98472.0971 3.9982
0.27363.9171.7951.7681.89591.9011
λ 1.2233 1.8750.514 4.0910.385 0.8990.3856 0.90230.3829 1.96930.3742 0.865
0.65173.5770.5140.51670.50150.4908
1.2748 1.82350.559 3.5850.457 0.830.4567 1.78170.4566 0.82440.454 0.8156
0.54863.0260.3730.37510.36780.3617
μ 0.6922 1.51270.121 3.2240.46 2.3840.4603 2.39540.4584 2.33460.4376 2.3406
0.82063.1031.9241.9351.87621.903
0.7571 1.44780.164 2.6450.487 2.1220.4891 2.1310.4815 2.0810.4666 2.082
0.69082.4811.6351.64231.60021.6154
100 α 0.3124 1.47610.001 4.5140.05 1.6250.0517 1.63230.0448 1.5910.0064 1.4982
1.16374.5131.5751.58061.54621.4918
0.4044 1.3840.003 3.4820.118 1.5380.1203 1.5880.1119 1.50040.0404 1.4756
0.97963.4791.421.46771.38851.4351
β 1.415 1.89420.0001 6.320.921 2.3850.9269 2.3880.8982 2.37280.8775 2.3747
0.47926.321.4641.46111.47461.4971
1.4529 1.85630.005 5.0871.4 2.3261.4074 2.32641.3685 2.30951.3358 2.312
0.40345.0820.9260.9190.9410.9762
γ 1.2533 1.71010.333 4.5762.575 3.9822.5849 3.99162.5329 3.94122.5413 3.9615
0.45684.2431.4071.40671.40831.4202
1.2894 1.6740.457 3.8882.792 3.9582.8012 3.96742.7596 3.9192.7688 3.938
0.38463.4311.1661.16621.15941.1692
λ 1.1975 1.50120.516 3.1840.442 0.7010.4425 0.70080.4421 1.5770.441 0.6994
0.30372.6680.2590.25820.25810.2584
1.2215 1.47720.58 2.7090.472 0.6690.4723 1.4510.4713 0.66650.4677 0.6615
0.25572.1290.1970.19730.19510.1939
μ 0.7975 1.43880.137 3.5790.401 2.0940.4032 2.09850.3908 2.07810.3474 2.0788
0.64143.4421.6931.69531.68731.7314
0.8482 1.38810.175 2.590.664 1.9750.6647 1.98550.6592 1.9280.6513 1.9226
0.53992.4151.3111.32091.26881.2713
200 α 0.5583 1.54010.0001 5.4070.249 2.3940.2504 2.40390.2411 2.35440.0914 2.36
0.98185.4072.1452.15352.11332.2686
0.636 1.46250.003 4.4350.309 2.220.3154 2.25370.2854 2.15070.2 2.1515
0.82654.4321.9111.93831.86531.9515
β 1.4217 2.1110.001 6.5111.5 2.3061.5017 2.30711.4917 2.30381.4881 2.3041
0.68946.510.8060.80540.81210.8161
1.4762 2.05650.003 5.4471.622 2.2351.6248 2.23631.6093 2.23121.606 2.2315
0.58035.4440.6130.61150.62190.6255
γ 1.2828 1.65690.316 3.8892.518 3.6532.5368 3.65442.4471 3.64912.4608 3.651
0.37413.5731.1351.11761.20211.1902
1.3124 1.62730.461 3.4092.648 3.6262.6622 3.62822.6182 3.61842.6227 3.6219
0.31492.9480.9780.9661.00020.9992
λ 1.155 1.35660.526 2.9910.471 0.6850.4715 0.68650.4712 1.4760.4707 0.6708
0.20152.4650.2140.2150.20880.2001
1.171 1.34060.593 2.6920.5 0.6640.4997 1.3750.4995 0.66350.4989 0.6603
0.16972.0990.1640.16410.16410.1614
μ 0.8439 1.42740.151 3.4920.521 1.620.5312 1.62390.4891 1.60310.4738 1.5991
0.58353.3411.0991.09271.1141.1253
0.89 1.38130.196 2.5670.58 1.5080.5814 1.52750.5725 1.4820.5499 1.4801
0.49122.3710.9280.94610.90950.9302
400 α 0.6056 1.28390.001 5.0970.258 1.9370.263 1.94010.2417 1.92310.1578 1.9229
0.67835.0961.6791.67711.68151.7651
0.6593 1.23030.003 3.6620.427 1.8410.4336 1.85050.415 1.80480.3641 1.8016
0.5713.6591.4141.4171.38981.4374
β 1.506 2.00580.001 6.2541.691 2.3361.6921 2.33941.6886 2.32391.6881 2.3257
0.49986.2530.6450.64730.63530.6376
1.5455 1.96630.003 5.0511.792 2.2711.7921 2.27181.7908 2.2691.7906 2.2693
0.42085.0480.4790.47970.47830.4786
γ 1.2996 1.60050.312 3.6482.888 3.5862.8926 3.58672.8679 3.58512.8738 3.5857
0.30093.3360.6980.69420.71730.7119
1.3234 1.57670.492 32.92 3.5642.921 3.56482.9142 3.5612.9159 3.5623
0.25332.5080.6440.64380.64680.6464
λ 1.1423 1.29340.55 2.9480.505 0.6550.5052 0.6550.505 1.3370.5047 0.6522
0.15122.3980.150.14980.1490.1475
1.1542 1.28150.635 2.2250.517 0.6380.5169 1.2620.5167 0.63720.5162 0.6363
0.12731.590.1210.12070.12050.1201
μ 0.9521 1.410.165 3.7340.672 1.6560.6736 1.65690.668 1.63370.6597 1.6229
0.4583.5690.9840.98330.96570.9633
0.9883 1.37380.219 2.9130.789 1.5840.7912 1.59150.7815 1.55790.7691 1.5519
0.38552.6940.7950.80020.77640.7827
The first and second lines show the credible HPD interval and the corresponding width of the parameter, respectively. In addition, 95% and 90% interval estimate, respectively.
Table 4. Descriptive statistics of the carbon dioxide emissions data.
Table 4. Descriptive statistics of the carbon dioxide emissions data.
Min.1st Qu.MedianMean3rd Qu.Max
0.59811.65062.16862.60033.33145.6114
Table 5. Estimated values of the APMW-W and the competing models.
Table 5. Estimated values of the APMW-W and the competing models.
Model α ^ β ^ γ ^ λ ^ μ ^
APMW-W9.61765293.29506439.85585157.48545010.5631985
APMW-E9.4615735.9480585.685677-6.868609
TW-D---0.014228450.29302159
Table 6. Kolmogorov–Smirnov test.
Table 6. Kolmogorov–Smirnov test.
ModelKSp-Value
APMW-W0.166590.903
APMW-E0.212420.683
TW-D0.292330.298
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

Emam, W.; Tashkandy, Y. Modeling the Amount of Carbon Dioxide Emissions Application: New Modified Alpha Power Weibull-X Family of Distributions. Symmetry 2023, 15, 366. https://doi.org/10.3390/sym15020366

AMA Style

Emam W, Tashkandy Y. Modeling the Amount of Carbon Dioxide Emissions Application: New Modified Alpha Power Weibull-X Family of Distributions. Symmetry. 2023; 15(2):366. https://doi.org/10.3390/sym15020366

Chicago/Turabian Style

Emam, Walid, and Yusra Tashkandy. 2023. "Modeling the Amount of Carbon Dioxide Emissions Application: New Modified Alpha Power Weibull-X Family of Distributions" Symmetry 15, no. 2: 366. https://doi.org/10.3390/sym15020366

APA Style

Emam, W., & Tashkandy, Y. (2023). Modeling the Amount of Carbon Dioxide Emissions Application: New Modified Alpha Power Weibull-X Family of Distributions. Symmetry, 15(2), 366. https://doi.org/10.3390/sym15020366

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