Next Article in Journal
LW-YOLO: Lightweight Deep Learning Model for Fast and Precise Defect Detection in Printed Circuit Boards
Next Article in Special Issue
Improving the Robustness of the Theil-Sen Estimator Using a Simple Heuristic-Based Modification
Previous Article in Journal
A Comparative Study of the Mode-Decomposed Characteristics of the Asymmetricity of a Vortex Rope with Flow Rate Variation
Previous Article in Special Issue
Statistical Inference and Application of Asymmetrical Generalized Pareto Distribution Based on Peaks-Over-Threshold Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Inference for the Gamma Zero-Truncated Poisson Distribution with an Application to Real Data

by
Patchanok Srisuradetchai
* and
Ausaina Niyomdecha
Department of Mathematics and Statistics, Faculty of Science and Technology, Thammasat University, Pathum Thani 12120, Thailand
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(4), 417; https://doi.org/10.3390/sym16040417
Submission received: 12 March 2024 / Revised: 28 March 2024 / Accepted: 31 March 2024 / Published: 2 April 2024

Abstract

:
This article presents Bayesian estimation methods applied to the gamma zero-truncated Poisson (GZTP) and the complementary gamma zero-truncated Poisson (CGZTP) distributions, encompassing both one-parameter and two-parameter models. These distributions are notably flexible and useful for modeling lifetime data. In the one-parameter model case, the Jeffreys prior is mathematically derived. The use of informative and noninformative priors, combined with the random walk Metropolis algorithm within a Bayesian framework, generates samples from the posterior distributions. Bayesian estimators’ effectiveness is examined through extensive simulation studies, in comparison with the maximum likelihood method. Results indicate that Bayesian estimators provide more precise parameter estimates, even with smaller sample sizes. Furthermore, the study and comparison of the coverage probabilities (CPs) and average lengths (ALs) of the credible intervals with those from Wald intervals suggest that Bayesian credible intervals typically yield shorter ALs and higher CPs, thereby demonstrating the effectiveness of Bayesian inference in the context of GZTP and CGZTP distributions. Lastly, Bayesian inference is applied to real data.

1. Introduction

The gamma zero-truncated Poisson (GZTP) and complementary gamma zero-truncated Poisson (CGZTP) distributions hold significant importance in statistical modeling, particularly for lifetime data exhibiting nonmonotonic hazard functions. The GZTP distribution provides a flexible model for phenomena where an event is guaranteed to occur, effectively handling datasets where zero counts are inapplicable. This makes it ideal for reliability analyses where the time to first failure is of interest. The CGZTP further extends this utility by modeling the bathtub-shaped hazard function, which is characterized by an initial decrease, followed by a constant rate, and then an increase. Such a capability to fit a bathtub hazard function makes the CGZTP a robust tool for complex survival data, overcoming the limitations of traditional models such as the gamma distribution.
The GZTP distribution, introduced by Niyomdecha et al. [1], is derived from compounding the gamma and zero-truncated Poisson distributions using the minimum function. Consider a collection of N independent and identically distributed random variables X 1 , X 2 , , X N , each following a gamma distribution with a probability density function f ( x ; α , β ) defined as f ( x ; α , β ) = β α x α 1 e β x / Γ ( α ) , x > 0 , where α > 0 is a shape parameter and β > 0 is a rate parameter. The variable N follows a zero-truncated Poisson distribution, given by: P ( N = n ) = e λ λ n / n ! ( 1 e λ ) ,   n = 1 , 2 , , and λ > 0 . Assuming X and N are independent, Z is defined as the minimum of X 1 , X 2 , , X N . The probability density function (pdf) of the GZTP distribution is denoted as
f G Z T P ( z ; λ , α , β ) = ( λ β α z α 1 e λ β z Γ ( α ) ( 1 e λ ) ) exp [ λ ( Γ ( α , β z ) Γ ( α ) ) ] , z > 0 , λ > 0 , α > 0 , β > 0 ,
where Γ ( α , β z ) = β z t α 1 e t d t is the upper incomplete gamma function. The CGZTP distribution utilizes the same compounding principle as the GZTP but employs the maximum function instead [2]. Consider Y to be the maximum of X 1 , X 2 , , X N . The pdf for the CGZTP distribution is then given as follows:
f C G Z T P ( y ; λ , α , β ) = ( λ β α y α 1 e λ β y Γ ( α ) ( 1 e λ ) ) exp [ λ ( 1 Γ ( α , β y ) Γ ( α ) ) ] ,   y > 0 , λ > 0 , α > 0 , β > 0 .
These distributions exhibit flexible shapes of distribution functions, as shown in Figure 1. Previous studies have covered inferential procedures for the parameters of the GZTP and CGZTP distributions. Niyomdecha et al. [1] employed maximum likelihood estimation (MLE) to estimate GZTP parameters and then examine their asymptotic properties, while the MLEs and asymptotic confidence intervals for CGZTP parameters were discussed by Niyomdecha and Srisuradetchai [2]. The MLEs exhibited accurate estimations, and the confidence interval achieved the nominal coverage probability in the case of a large sample size. Several studies have been conducted on compound distributions using Bayesian methods. Xu et al. [3] investigated Bayesian estimators of Exponential-Poisson (EP) parameters by employing general noninformation prior distributions under symmetric and asymmetric loss functions. Yan et al. [4] determined the Bayesian estimators of the parameters in the EP distribution under general entropy, LINEX, and a scaled squared loss function based on type-II censoring. In a study conducted by Pathak et al. [5], the Bayesian estimators of Weibull–Poisson (WP) parameters were obtained by assuming that these parameters follow independently distributed prior distributions.
The loss function is instrumental in measuring how much an estimated parameter value deviates from its true value. The squared error loss function, a symmetric loss function, is commonly used in practice, especially when overestimation and underestimation are equally problematic [6,7]. Elbatal et al. [8] addressed parameter estimation for the Nadarajah–Haghighi distribution with progressive Type-1 censoring, employing the squared error loss function to produce Bayes estimates and credible intervals for maximum posterior density. Eliwa et al. [9] utilized balanced linear exponential and general entropy loss functions to estimate parameters for the new Weibull-Pareto distribution. Similarly, Abdel-Aty et al. [10] applied squared error, LINEX, and general entropy loss functions for future failure times in a joint type-II censored sample from multiple exponential populations.
In many studies, the posterior distributions often become complicated and cannot be simplified into any closed form. The samples were obtained from the posteriors using a Markov Chain Monte Carlo (MCMC) method, such as the Gibbs sampler (see [11]) and the Metropolis–Hastings algorithm (see [12]). Additionally, a summary of the predicted posteriors is provided based on a sample-based approach. The MCMC based on Metropolis–Hastings algorithms was used in [13] to estimate the unknown parameters of the alpha-power Weibull distribution under Type II hybrid censoring. In [14], the parameters of the unit-log-logistic distribution were estimated using a Bayesian approach. Noninformative priors were used, and samples from the joint posterior distribution were obtained using the random walk Metropolis algorithm. El-Sagheer et al. [15] employed Gibbs sampling to estimate the parameters for an asymmetric distribution and various lifetime indices, including reliability and hazard rate functions.
It is widely regarded that the conjugate prior for the Poisson parameter and the gamma rate parameter follow a gamma distribution. However, there is no proper conjugate prior for the gamma shape parameter [16,17]. Several papers explore Bayesian inference for estimating the parameters of the gamma distribution. Naji and Rasheed [18] derived Bayes estimators for the shape and scale parameters of the gamma distribution using the precautionary loss function. They assumed gamma and exponential priors for the shape and scale parameters, respectively, to represent prior information. Moala et al. [19] studied various noninformative priors, including Jeffreys prior, reference prior, maximal data information prior, Tibshirani prior, and a novel prior based on the copula method. Additionally, Pradhan and Kundu [20] assumed that the scale parameter follows a gamma distribution prior, while the shape parameter follows a log-concave distribution prior.
The existing literature has not addressed Bayesian inference on parameters of the GZTP and CGZTP) distributions. While Niyomdecha et al. [1] and Niyomdecha and Srisuradetchai [2] have conducted MLE and Wald’s interval analyses, their findings suggest that the mean square errors of MLEs remain high and that Wald’s interval coverage probabilities are below the nominal level for small sample sizes. This study, therefore, seeks to explore Bayesian inference for the GZTP and CGZTP distributions.
This paper is structured as follows: Section 2 delves into maximum likelihood estimation along with the corresponding interval estimation, which will be compared with the Bayesian credible interval. Section 3 elaborates on the prior and posterior distributions, estimation procedures based on the squared error loss function, and the application of the random walk Metropolis algorithm for simulating posterior samples. Simulation studies, which are conducted for scenarios involving one and two unknown parameters within both GZTP and CGZTP distributions, are presented in Section 4. Section 5 demonstrates two example applications using real data. Finally, the paper concludes with a discussion in Section 6.

2. Maximum Likelihood Estimation

Let Z 1 , Z 2 , , Z n be random samples from a GZTP distribution and Y 1 , Y 2 , , Y n be random samples from a CGZTP distribution. The likelihood functions based on the observed random sample of size n will be as follows:
L G Z T P ( λ , α , β | z ) = ( λ e λ 1 e λ ) n ( β n α ( Γ ( α ) ) n ) ( i = 1 n z i ) α 1 exp [ β ( i = 1 n z i ) + λ Γ ( α ) i = 1 n Γ ( α , β z i ) ] ,
L C G Z T P ( λ , α , β | y ) = ( λ 1 e λ ) n ( β n α ( Γ ( α ) ) n ) ( i = 1 n y i ) α 1 exp [ β ( i = 1 n y i ) λ Γ ( α ) i = 1 n Γ ( α , β y i ) ] .
The corresponding log-likelihood function of the GZTP distribution is
l G Z T P ( λ , α , β | z ) = n ( log λ λ log ( 1 e λ ) ) + n α log β n log Γ ( α ) + ( α 1 ) i = 1 n log z i β ( i = 1 n z i ) + λ Γ ( α ) i = 1 n Γ ( α , β z i ) ,
and the log-likelihood function of the CGZTP distribution is
l C G Z T P ( λ , α , β | y ) = n ( log λ log ( 1 e λ ) ) + n α log β n log Γ ( α ) + ( α 1 ) i = 1 n log y i β ( i = 1 n y i ) λ Γ ( α ) i = 1 n Γ ( α , β y i ) .
The maximum likelihood estimators of λ , α and β for the GZTP and CGZTP distributions are obtained by maximizing (5) and (6). This process is accomplished by solving the first derivatives with respect to each parameter of the log-likelihood function. These first derivatives are difficult and complex to solve, making it impossible to find the MLE of λ , α and β analytically. Consequently, numerical techniques such as the simulated annealing method are employed to estimate λ , α and β that maximize the likelihood function.
The MLEs are asymptotically normally distributed with a multivariate normal (MVN) distribution given by
( λ ^ , α ^ , β ^ ) M V N ( ( λ , α , β ) , I 1 ( λ , α , β ) )   as   n ,
where I ( λ , α , β ) is the Fisher information matrix with element I i j = 𝜕 l 2 / 𝜕 θ i 𝜕 θ j ,     i , j = 1 , 2 , 3 and θ = ( λ , α , β ) [21]. The asymptotic variances of MLEs can be obtained from the inverse Fisher information matrix. Then, the corresponding ( 1 α ) 100 % Wald confidence intervals for θ i are given by θ ^ i ± z 1 α / 2 I ^ i i , where I ^ i i is the ii-th element of the inverse of [ I ( θ ^ ) 1 ] , i.e., I ^ i i = [ I ( θ ^ ) 1 ] i i and z 1 α / 2 is the ( 1 α / 2 ) -th quantile of the standard normal [22].

3. Bayesian Estimation

This section presents the formulation of prior distributions for each parameter, acknowledging their independence, and the subsequent derivation of joint posterior distributions.

3.1. Prior and Posterior Distributions

3.1.1. Case 1: α and β Are Unknown

To estimate the parameters for the GZTP or CGZTP distributions when α and β are unknown but λ is known, we assume that α and β have priors p 1 ( ) and p 2 ( ) , which correspond to G a m m a ( a ,   b ) and G a m m a ( c ,   d ) , respectively, and they are independently distributed. The prior distributions for α and β are obtained as follows:
p 1 ( α ; a , b ) = b a α a 1 e b α Γ ( a )   and   p 2 ( β ; c , d ) = d c β c 1 e d β Γ ( c ) ,
where hyperparameters β > 0 ,   c > 0 ,   d > 0 .
Let Z 1 , Z 2 , , Z n be random samples from a GZTP distribution, so the joint posterior distribution given data z = ( z 1 , z 2 , , z n ) is as follows:
p G Z T P ( α , β | z ) = L ( α , β ; z ) p 1 ( α ) p 2 ( β ) 0 0 L ( α , β ; z ) p 1 ( α ) p 2 ( β ) d α d β   L ( α , β ; z ) p 1 ( α ) p 2 ( β ) ( λ e λ 1 e λ ) n ( β n α ( Γ ( α ) ) n ) ( i = 1 n z i ) α 1 exp [ β ( i = 1 n x i ) + λ Γ ( α ) i = 1 n Γ ( α , β z i ) ] × ( b a α a 1 e b α Γ ( a ) ) ( d c β c 1 e d β Γ ( c ) ) ,
and Y 1 , Y 2 , , Y n are the random samples from the CGZTP distribution, and the corresponding joint posterior distribution given data y = ( y 1 , y 2 , , y n ) is given by:
p C G Z T P ( α , β | y ) ( λ 1 e λ ) n ( β n α ( Γ ( α ) ) n ) ( i = 1 n y i ) α 1 exp [ β ( i = 1 n y i ) λ Γ ( α ) i = 1 n Γ ( α , β y i ) ] × ( b a α a 1 e b α Γ ( a ) ) ( d c β c 1 e d β Γ ( c ) ) .
The marginal posterior distributions of α and β have no closed form; as a consequence, the MCMC method is employed to provide Bayesian estimation.

3.1.2. Case 2: λ Is Unknown

Gamma Priors

To estimate parameter λ for the GZTP or CGZTP when α and β are known. Assuming that λ has a prior G a m m a ( m ,   n ) :
p 3 ( λ ) = n m λ m 1 e n λ Γ ( m ) , λ > 0 ,   m > 0 ,   n > 0 .
Let Z 1 , Z 2 , , Z n be random samples from a GZTP distribution. The corresponding posterior distribution, given data z = ( z 1 , z 2 , , z n ) , is as follows:
p G Z T P ( λ | z )   L ( λ ; z ) p 3 ( λ ) = ( λ e λ 1 e λ ) n ( β n α ( Γ ( α ) ) n ) ( i = 1 n z i ) α 1 exp [ β ( i = 1 n z i ) + λ Γ ( α ) i = 1 n Γ ( α , β z i ) ] × ( n m λ m 1 e n λ Γ ( m ) ) .
Furthermore, let Y 1 , Y 2 , , Y n be random samples from a CGZTP distribution. The corresponding posterior distribution given data y = ( y 1 , y 2 , , y n ) is:
p C G Z T P ( λ | y )   L ( λ ; y ) p 3 ( λ )   = ( λ 1 e λ ) n ( β n α ( Γ ( α ) ) n ) ( i = 1 n y i ) α 1 exp [ β ( i = 1 n y i ) λ Γ ( α ) i = 1 n Γ ( α , β y i ) ] × ( n m λ m 1 e n λ Γ ( m ) ) .
Because the posterior distributions above are also complicated to derive analytically, MCMC is used to simulate samples from them.

Jeffreys Prior

Jeffreys prior is proposed as a widely known prior to represent a situation in which there is little information regarding the parameters. The Jeffreys prior for one parameter is proportional to the square root of the expected Fisher information [23]. From the likelihood functions in cases when α and β are known, the associated gradients are:
𝜕 l G Z T P ( λ ; z ) 𝜕 λ = n ( 1 λ 1 e λ 1 e λ ) + i = 1 n Γ ( α , β z i ) Γ ( α ) ,
𝜕 l C G Z T P ( λ ; y ) 𝜕 λ = n ( 1 λ e λ 1 e λ ) i = 1 n Γ ( α , β y i ) Γ ( α ) .
By differentiating (7) and (8), the observed Fisher information values of λ for the GZTP and CGZTP distributions are the same:
I G Z T P ( λ ) = I C G Z T P ( λ ) = n ( 1 + e 2 λ e λ ( λ 2 + 2 ) ) ( e λ 1 ) 2 λ 2
and the expected Fisher information is
J G Z T P ( λ ) = J C G Z T P ( λ ) = E [ n ( 1 + e 2 λ e λ ( λ 2 + 2 ) ) ( e λ 1 ) 2 λ 2 ] = n ( 1 + e 2 λ e λ ( λ 2 + 2 ) ) ( e λ 1 ) 2 λ 2 .
Thus, from (9), the Jeffreys prior for the λ parameter is given by:
p 4 ( λ ) n ( 1 + e 2 λ e λ ( λ 2 + 2 ) ) ( e λ 1 ) 2 λ 2 .

3.2. Point and Interval Estimations

This section explores the process of obtaining Bayesian estimates and constructing credible intervals for unknown parameters of the GZTP and CGZTP distributions. The squared error loss function is a symmetric function, defined by L ( θ , θ ^ ) = ( θ - θ ^ ) 2 , where θ ^ is an estimate of θ [24,25]. For example, for given data z = ( z 1 , z 2 , , z n ) , under the squared error loss function, the Bayesian estimator of θ is θ ^ B G Z T P = E ( θ | z ) .
Bayesian interval estimates for θ are also calculated based on the posterior distribution p ( θ | z ) . They are referred to as credible intervals to differentiate them from confidence intervals. For a given value of γ ( 0 , 1 ) , a γ 100 % credible interval is determined by values l and u that satisfy
l u p ( θ | z ) d θ = γ ,
where γ is called the credible level of the credible interval [ l ,   u ] .
For GZTP, the Bayesian estimates of the unknown parameter(s) are given by:
  • Case 1: α and β are unknown:
    g ^ B G Z T P ( α , β ) = E ( g ( α , β ) | z ) = 0 0 g ( α , β ) p G Z T P ( α , β | z ) d α d β ,
  • Case 2: λ is unknown:
    λ ^ B G Z T P = E ( λ | z ) = 0 λ p G Z T P ( λ | z ) d λ .
From CGZTP for the given data y 1 , y 2 , , y n , the Bayes estimates of the unknown parameter(s) are given by
  • Case 1: α and β are unknown:
    g ^ B C G Z T P ( α , β ) = E ( g ( α , β ) | y ) = 0 0 g ( α , β ) p C G Z T P ( α , β | y ) d α d β ,
  • Case 2: λ is unknown:
    λ ^ B C G Z T P = E ( λ | y ) = 0 λ p C G Z T P ( λ | y ) d λ .
Due to the complexity involved in constructing explicit forms of Bayes estimates for both cases, the random walk Metropolis algorithm, which is a variant of the MCMC methods, will be utilized to derive the Bayes estimates of the unknown parameters.

3.3. Random Walk Metropolis Algorithm

In this paper, random walk Metropolis is implemented. The random walk Metropolis (RWM) algorithm, a subset of the Metropolis–Hastings algorithms, is favored in Bayesian computation for its conceptual simplicity and operational ease. It is particularly advantageous when the posterior distribution is unknown or complex, offering a straightforward mechanism to generate sample values for parameter estimation. The strength of RWM lies in its local exploration capability, allowing it to meticulously probe the parameter space using a symmetric proposal density, which simplifies the acceptance criteria. This simplicity facilitates easier tuning and implementation, often requiring only the adjustment of the proposal distribution’s scale to balance the acceptance rate and the chain’s mixing [26,27]. Concurrently, credible intervals for unknown parameters are generated during this procedure. Since the density of the posterior distribution is proportional to the product of the likelihood and the density of the prior distribution, we use L ( α , β ; z ) p 1 ( α ) p 2 ( β ) or L ( α , β ; y ) p 1 ( α ) p 2 ( β ) as the target density for generating random samples from the joint posterior distribution of α and β . The RWM algorithm for generating random samples from the joint posterior distributions of α and β is shown as follows:
  • Choose starting values of θ 0 = [ α 0 β 0 ] and define σ 1 .
  • At step i , we draw [ ε 1 ε 2 ] i M V N ( [ 0 0 ] , σ 1 2 [ 1 0 0 1 ] ) and draw a new value
    [ α ˜ β ˜ ] i = [ α β ] i 1 + [ ε 1 ε 2 ] i
  • The candidate [ α ˜ β ˜ ] i will be accepted with a probability given by the Metropolis ratio:
    r ( [ α β ] i 1 , [ α ˜ β ˜ ] i ) = min { L G Z T P ( α ˜ , β ˜ ; z ) p 1 ( α ˜ ) p 2 ( β ˜ ) L G Z T P ( α , β ; z ) p 1 ( α ) p 2 ( β ) , 1 }   for   GZTP ,
    and   r ( [ α β ] i 1 , [ α ˜ β ˜ ] i ) = min { L C G Z T P ( α ˜ , β ˜ ; y ) p 1 ( α ˜ ) p 2 ( β ˜ ) L C G Z T P ( α , β ; y ) p 1 ( α ) p 2 ( β ) , 1 }   for   CGZTP .
  • Repeat steps 2–3 M times to obtain samples ( α 1 , β 1 ) , ( α 2 , β 2 ) , , ( α M , β M ) and discard the first N values of the chain for burn-in.
  • The Bayesian estimates of parameters α and β are computed as
    α ^ B = j = N + 1 M α j M N   and   β ^ B = j = N + 1 M β j M N .
  • To compute the credible intervals of α and β , order α i and β i , then γ 100 % credible interval of α and β can be given, respectively, by
    [ α ( ( M N ) ( 1 γ ) / 2 ) ,   α ( ( M N ) ( 1 + γ ) / 2 ) ]   and   [ β ( ( M N ) ( 1 γ ) / 2 ) ,   β ( ( M N ) ( 1 + γ ) / 2 ) ] ,
    where α ( . ) and β ( . ) are the ordered statistics from the MCMC samples for α and β , after discarding the burn-in, and is the floor function.
For the unknown λ, L ( λ ; z ) p ( λ ) or L ( λ ; y ) p ( λ ) is considered the target density for generating random samples from the posterior distribution of λ . The RWM algorithm for generating random samples from the posterior distribution of λ is given below:
  • Choose starting values of λ 0 and define σ 2 .
  • At step i , we draw ε i N ( 0 , σ 2 2 ) and draw a new value λ ˜ i = λ i 1 + ε i .
  • The candidate λ ˜ i will be accepted with a probability given by the Metropolis ratio:
    r ( λ i 1 , λ ˜ i ) = min ( L G Z T P ( λ ˜ ; z ) p 3 ( λ ) L G Z T P ( λ ; z ) p 3 ( λ ) , 1 )   for   GZTP ,
    and
    r ( λ i 1 , λ ˜ i ) = min ( L C G Z T P ( λ ˜ ; y ) p 3 ( λ ) L C G Z T P ( λ ; y ) p 3 ( λ ) , 1 )   for   CGZTP .
  • Repeat steps 2–3 M times to obtain samples λ 1 , λ 2 , , λ M and remove the first N values of the chain for burn-in.
  • The Bayesian estimates of parameters λ is computed by λ ^ B = j = N + 1 M λ j / ( M N ) .
  • To compute the credible intervals of λ , order λ i , then γ 100 % credible interval of λ can be given by [ λ ( [ ( M N ) ( 1 γ ) / 2 ] ) ,   λ ( ( M N ) ( 1 + γ ) / 2 ) ] , where λ ( . ) is the ordered statistics.

4. Simulation Study

The simulation study encompasses various sample sizes and hyperparameter values. Specifically, sample sizes n = 15 ,   25 ,   50 , and 100 are examined. Table 1 presents the hyperparameter values for informative prior distributions. The means of the prior distributions, which have small and large variances, are equal to the true values of the unknown parameters, α and β . For example, for the case of α = 2 and β = 0.5 , with the hyperparameters of Prior 1 ( a = 1 , b = 0.5 ), the variance of α equals 4, and with the hyperparameters of Prior 2 ( a = 2 , b = 1 ), the variance of α equals 2. Thus, the variance of Prior 1 is considered “High” compared to that of Prior 2. Both prior distributions of α have the same mean, 2, which equals to the true value.
Values of α and β are selected to create a variety of distribution shapes, as shown in Figure 1. Additionally, Table 2 details the prior distributions of the parameter λ . The shapes of gamma prior distributions with different hyperparameters, as presented in Table 1 and Table 2 are illustrated in Figure 2.
The RWM algorithm, as described in Section 3.3, is executed for 10,000 iterations with a burn-in period of 1000. In both panels of Figure 3, the examples of the trace plots for α and β suggest that the Markov chains have reached a stationary distribution, evidenced by the dense and fuzzy appearance of the plots, which indicates good mixing of the chains. The variability observed within each plot is consistent with the stochastic nature expected from RWM sampling, and there are no discernible trends or drifts to suggest nonconvergence.
Further examination of the pair plots shown in Figure 4 reveals that, despite the initial starting points being far from the true values, the pairs of samples drawn from the RWM algorithm progress toward a densely clustered area. This dense cluster signifies the region of high probability density within the posterior distribution, illustrating the algorithm’s ability to converge to the region of interest. For unknown λ , the examples of the trace plots are shown in Figure 5.
Monte Carlo simulations are performed to compare the performances of the Bayes estimators with those of the classical estimators. Point estimates for parameters are averaged over 1000 iterations. Subsequently, the mean squared errors (MSEs) of the parameter estimates are calculated. The coverage probabilities (CPs) of 95% Wald confidence intervals and credible intervals and their average lengths (ALs) are determined.
Table 3 presents the detailed MLEs and MSEs obtained from simulated data sets from GZTP where the values of α and β are unknown. As sample sizes increase, estimates become more accurate and the MSE values decrease. For example, for the parameter set (0.5, 2, 0.5), the MLE of α decreases from 2.3637 with an MSE of 0.8896 at sample size 15 to 2.0644 with an MSE of 0.0760 at sample size 100. It is also noticed that the MSEs of α ^ and Bayes estimate α ^ B increase as α increases given that β is fixed. For example, in a case where λ = 1 , β = 1, and the sample size is 15, when α is 1, the MSE for Prior 2 is 0.0566; whereas when α increases to 2, the MSE for the same prior and sample size is 0.2314. Similarly, when α is fixed, the MSEs of β ^ and Bayes estimate β ^ B increase as β increases. When the sample size is small, MLEs tend to have higher MSE compared to Bayesian estimates. However, between the two priors, Prior 2 has the lowest MSE. For instance, for the parameter set (0.5, 2, 0.5) and a sample size of 25, Prior 2 produces an estimate for α with an MSE of 0.1854, which is lower than the MSE of 0.2777 for Prior 1, indicating that the more informative Prior 2 generally results in a lower MSE.
Form Table 4, the 95% credible intervals and Wald confidence intervals for α and β are presented. As the sample size increases, the CPs tend to approach the nominal coverage probability of 0.95, while the ALs decrease. Typically, CPs are generally above 0.95, despite the small sample size of 15. Moreover, the Prior 2 tends to have the smallest ALs with the same sample size.
Figure 6 graphically summarizes the average of estimates, MSEs, CPs, and ALs for a selected case of GZTP. In the first row, which shows estimates of α and β, there are two line graphs, one for each parameter. These results display the average estimates obtained through MLE and two Bayesian methods with different priors (Prior 1 and Prior 2). Prior 2 yields the most accuracy, followed by Prior 1 and MLE. However, as the sample size increases, the estimates from all methods converge, suggesting that larger sample sizes lead to more accurate estimations. From the second row, the bar charts show that the precision of the estimation methods improves with larger samples. Bayesian estimates, particularly those with Prior 2, tend to have lower MSEs than MLEs. In the third row, the CPs from Prior 2 tend to be higher than the others, especially in the small sample size. As the sample size increases, all the methods tend to produce about the same CP. From the last row, the ALs generally decrease as the sample size increases, showing that the intervals become narrower and, thus, more precise with larger samples. The Bayesian estimates with Prior 2 consistently show the shortest ALs across all sample sizes for both parameters.
From simulated CGZTP datasets where the values of α and β are unknown, the conclusions are consistent with those from GZTP. The detailed results are summarized in Appendix A, Table A1 and Table A2 which provide the average Bayesian estimates, MSEs, CPs, and ALs of parameters. As sample sizes increase, estimates become more accurate, and the MSE values decrease. It is observed that, while holding β constant, the MSEs of α ^ and α ^ B increase as α increases. Likewise, as β increases, the MSEs of β ^ and Bayes estimate β ^ B increase when α remains constant. Applying Prior 2 results in the lowest MSE values for α and β . Figure 7 graphically summarizes the averages of estimates, MSEs, CPs, and ALs for α and β of the two-parameter CGZTP distribution with λ = 1 ,   α = 1 , and β = 1 of CGZTP.
For parameter λ of the GZTP distribution, the MSE values decrease as sample sizes increase, as shown in Figure 8. Both the MLE and the Jeffreys prior estimates result in poor estimations with large MSEs. As expected, the performance of these methods improves as the sample size increases. However, Prior 2 yields the lowest MSE, even with a sample size as small as 15. In all situations, the CPs either approach the desired coverage probability or exceed 0.95, and the ALs decrease as sample sizes increase. For additional scenarios, Table A3 and Table A4 in Appendix A provide a summary of all estimates, MSEs, CPs, and ALs. Similarly, the findings for λ of the CGZTP distribution are in line with those for the GZTP distribution. The MSEs of Bayesian estimates under informative priors are lower than those for the MLE and the Jeffreys prior. The CPs of credible intervals from informative priors tend to be higher than those from Wald intervals and the credible interval from the Jeffreys prior. Figure 9 presents estimates, MSEs, CPs, and ALs for three methods in a specific case, while results for other scenarios are compiled in Table A5 and Table A6.
.

5. Application on Real Data

5.1. March Precipitation

The data in Table 5 represent the amount of precipitation (in inches) that fell in March in Minneapolis/St. Paul, including 30 consecutive measurements. These data were first discussed by Hinkley [28]. The MLE and the Bayesian summaries for α and β are presented in Table 6. The MLE of the parameters are α ^ = 3.0597 and β ^ = 1.7214 with the corresponding standard errors as 0.7207 and 0.4656, respectively. The 95% confidence intervals for α ^ and β ^ are (1.6347, 4.4847) and (0.8086, 2.6341), respectively. Using the RWM with gamma priors and assuming λ = 0.38, the Bayesian estimates for α and β are α ^ B = 3.0289 and β ^ B = 1.7005 , respectively. The 95% credible intervals for α and β are given, respectively, as (1.8884, 4.3697) and (0.9805, 2.6144). A histogram with the fitted GZTP curve (the orange line) and Bayesian estimates is illustrated in Figure 10. The pair plot of RWM depicted in Figure 11 illustrates that the estimates are converging to a stationary state and distributing around the posterior means, with cluster points appearing elliptical, suggesting a correlation between these two parameters. This pattern implies that Wald’s confidence intervals may not be optimal.
Notably, the Bayesian estimates are very close to the MLE values; however, the Bayesian credible intervals are consistently shorter than those from Wald’s confidence intervals. The Kolmogorov–Smirnov (K-S) test and the corresponding p-values indicate that both methods are effective, with Bayesian estimates exhibiting a slight performance edge over the MLEs. Note that a higher p-value indicates a better fit model.

5.2. Remission Time of Bladder Cancer Patients

The dataset consists of the number of months that 128 patients with bladder cancer spent in remission, as reported by Lee and Wang [29]. From Table 6, the MLEs of the parameters are α ^ = 1.1719 and β ^ = 0.1258 with the corresponding 95% confidence intervals for α ^ and β ^ ‚ as (0.9147, 1.4290) and (0.0919, 0.1599), respectively. Assuming that λ = 0.0238 with the unknown parameters being α and β, the Bayesian estimates for α and β are α ^ B = 1.1697 and β ^ B = 0.1246 , respectively. The 95% credible intervals for the parameters α and β are given, respectively, by (0.9283, 1.4277) and (0.0921, 0.1596). Moreover, the K-S tests suggest that both methods can be used to model the data at a significance level of 0.05. The histogram with the fitted GZTP curve (the orange line) is illustrated in Figure 12, and the pair plot of RWM is depicted in Figure 13. Initially, there appears to be a wide spread of values indicating that the Markov chain is exploring the parameter space. As the iterations progress, the points seem to converge towards a narrower region of the plot, which suggests that the parameters are settling into a region that could represent the mode of the posterior distribution.

6. Conclusions and Discussion

Both point and interval estimation have been studied within Bayesian frameworks. For point estimation, informative gamma priors with both low and high variance, as well as Jeffreys prior, are employed, and the results are compared with those obtained from MLE in terms of MSEs. Since the posterior distributions for GZTP and CGZTP do not have closed forms, the RWM algorithm is utilized to generate posterior samples. Furthermore, the Bayesian credible intervals are compared to Wald’s intervals in terms of coverage probability and average length.
Bayesian estimates using informative priors are obviously superior to the MLE and Bayesian estimates with Jeffreys priors in terms of MSEs. Among the informative priors having the mean equal to the true parameter, the one with a low variance yields a slightly lower MSE compared to that with high variance. In detail, when α is fixed, the MSEs of β ^ and Bayes estimate β ^ B increase as β increases. Similarly, the MSEs of α ^ and Bayesian estimate α ^ B increase as α increases given that β is fixed. In the case of an unknown λ, where the Jeffreys prior can be mathematically derived, the corresponding Bayesian estimates are slightly better than the MLE. However, as the sample size increases, the discrepancy among the MSEs obtained from all methods tends to decrease.
For interval estimation, the Bayesian credible intervals tend to be more conservative, with coverage probabilities exceeding the nominal level of 0.95, particularly for small sample sizes. The ALs of the credible intervals are notably shorter than those of the Wald confidence intervals. It is worth noting that credible intervals can achieve greater coverage with shorter interval lengths because Priors 1 and 2 were deliberately chosen so that the expected value of the prior equals the true parameter value. These informative priors have a substantial impact, often outweighing the data in their influence on the posterior. As the sample size increases, the influence of the data begins to outweigh that of the prior. In such cases, the lengths of the credible intervals tend to converge towards those of the frequentist confidence intervals, such as Wald’s intervals, and the high coverage probabilities adjust closer to the expected levels under the true confidence level.

Author Contributions

Conceptualization, P.S.; methodology, P.S. and A.N.; software, A.N.; validation, P.S. and A.N.; formal analysis, P.S.; investigation, P.S.; resources, A.N.; data curation, A.N.; writing—original draft preparation, P.S. and A.N.; writing—review and editing, P.S.; visualization, A.N.; supervision, P.S.; project administration, P.S.; funding acquisition, P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by Thammasat University Research Fund, Contract No TUFT 71/2566.

Data Availability Statement

The data presented in this study are openly available in reference number [28,29].

Acknowledgments

We express our gratitude to Thammasat University for their financial assistance. Additionally, we extend our appreciation to the editor and reviewers for their insightful remarks and suggestions, which have greatly enhanced the quality of the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. MLE and Bayesian estimates and mean squared errors of α and β for the CGZTP distributions.
Table A1. MLE and Bayesian estimates and mean squared errors of α and β for the CGZTP distributions.
( λ ,   α ,   β ) n ParameterMLEPrior 1Prior 2
EstMSEEstMSEEstMSE
(0.5, 2, 0.5)15 α 2.40641.07452.23000.54262.14700.2972
β 0.61550.07760.56640.03990.53880.0200
25 α 2.27910.55932.13280.29872.09220.1994
β 0.57700.03920.54170.02340.53050.0151
50 α 2.12260.19492.07040.15922.05810.1306
β 0.53410.01300.52320.01120.52000.0091
100 α 2.05050.08382.05980.08332.05510.0756
β 0.51410.00580.51730.00580.51600.0052
(0.5, 2, 1)15 α 2.46321.34652.21940.57222.11800.2863
β 1.24790.37941.12090.16031.07080.0780
25 α 2.27310.55342.15890.36282.11030.2345
β 1.14620.14691.10140.11041.07430.0710
50 α 2.14150.22062.08980.16482.07670.1345
β 1.07780.06121.04900.04641.04180.0374
100 α 2.04380.08492.04760.08872.04360.0805
β 1.02170.02311.02640.02461.02430.0222
(1, 1, 1)15 α 1.23640.30331.13290.17441.05450.0755
β 1.26130.36501.16010.21481.07590.0872
25 α 1.14460.14461.09260.10161.05350.0617
β 1.16390.16621.10750.11651.07540.0749
50 α 1.06560.05661.04910.04271.04810.0365
β 1.06450.0581.05490.05151.05210.044
100 α 1.02140.02131.01890.01891.01740.018
β 1.02940.02521.02240.02271.02220.0208
(1, 2, 1)15 α 2.50581.53432.21720.61292.12470.2988
β 1.25010.36531.11320.15461.07870.0791
25 α 2.30690.66662.13330.36362.08950.2332
β 1.15460.16371.07040.09141.04900.0582
50 α 2.10530.21202.08230.18802.06830.1528
β 1.05200.05111.0440.04461.03730.036
100 α 2.06380.09112.03130.08612.02810.0779
β 1.03220.02201.01280.02131.01140.0192
Table A2. Coverage probabilities and average lengths of intervals for α and β of the CGZTP distributions.
Table A2. Coverage probabilities and average lengths of intervals for α and β of the CGZTP distributions.
( λ ,   α ,   β ) n ParameterMLEPrior 1Prior 2
CPALCPALCPAL
(0.5, 2, 0.5)15 α 0.97003.38580.96772.88860.98542.4848
β 0.96300.90290.97230.77340.98690.6534
25 α 0.96202.47680.97202.19800.98402.0040
β 0.95600.65640.97400.58950.98200.5340
50 α 0.95601.62490.95301.54120.96601.4722
β 0.96700.43040.96000.41220.97100.3924
100 α 0.95301.10810.95101.09540.95401.0699
β 0.96300.29330.94700.29060.94700.2832
(0.5, 2, 1)15 α 0.96403.48340.96692.87760.98312.4506
β 0.96001.83630.96231.52920.98381.2982
25 α 0.96102.46930.96902.22730.97902.0198
β 0.96501.30410.96101.19970.97501.0820
50 α 0.95801.64020.96901.55630.97301.4867
β 0.95700.86850.96200.82410.97200.7857
100 α 0.95601.10420.94301.09100.94801.0680
β 0.96400.58310.93300.57800.94100.5644
(1, 1, 1)15 α 0.96301.78000.96461.53400.98621.2472
β 0.95801.88640.96381.65100.98771.3453
25 α 0.96701.26790.94901.15380.97101.0215
β 0.97101.35290.96001.23310.96701.1009
50 α 0.94800.83100.95800.79540.96600.7598
β 0.96400.87910.94600.84690.95300.8078
100 α 0.96200.56130.94600.54810.95100.5337
β 0.96400.60280.94800.58740.95600.5721
(1, 2, 1)15 α 0.96203.70030.96852.98720.98312.5573
β 0.96701.81370.97001.49380.97921.2863
25 α 0.95802.63260.94402.28040.96602.0762
β 0.95801.29870.94701.13630.96301.0340
50 α 0.95101.69380.94601.60690.95401.5427
β 0.95000.83790.95200.80080.96100.7682
100 α 0.95901.17220.94601.12870.94601.0986
β 0.95800.58110.94700.55720.94500.5438
Table A3. MLE and Bayesian estimates and mean squared errors for λ of the GZTP distributions.
Table A3. MLE and Bayesian estimates and mean squared errors for λ of the GZTP distributions.
( λ ,   α ,   β ) n MLEPrior 1Prior 2Jeffreys
EstMSEEstMSEEstMSEEstMSE
(0.5, 2, 0.5)150.97830.68020.43260.05200.47850.01771.13440.7413
250.80880.41880.43960.05180.48980.01350.92370.3766
500.66050.19440.46930.04920.50340.01040.71840.1554
1000.56370.09610.47670.04820.50190.00730.61550.0809
(0.5, 2, 1)150.97060.68770.44180.04300.48060.01711.10380.6507
250.77410.36250.45710.04600.48620.01420.92350.3736
500.66460.19740.47480.04770.49670.00990.73180.1698
1000.5580.09970.47700.04990.50160.00730.63970.0964
(1, 1, 1)151.27120.67590.92820.21820.97560.07041.42790.6366
251.12560.45310.88530.19700.98110.07361.26350.4066
501.03750.23030.87100.16110.99100.08441.10830.1933
1001.00070.13010.91200.11110.99970.06291.04110.1204
(1, 2, 1)151.27100.67810.91790.21930.98570.06321.40520.6096
251.15460.41630.89260.19990.96700.07821.25360.3643
501.02150.22250.89360.16050.95680.07261.12490.2155
1001.00900.12550.90930.10020.95120.06391.03580.1125
Table A4. Coverage probabilities and average lengths of intervals for λ of the GZTP distributions.
Table A4. Coverage probabilities and average lengths of intervals for λ of the GZTP distributions.
( λ ,   α ,   β ) n MLEPrior 1Prior 2Jeffreys
CPALCPALCPALCPAL
(0.5, 2, 0.5)150.97302.75711.00001.47011.00001.03810.97502.5537
250.96702.15220.99801.32251.00001.00420.96102.0079
500.97101.57700.99301.11841.00000.91390.97301.4764
1000.96701.17730.98500.95050.99900.80640.97101.1332
(0.5, 2, 1)150.97502.74801.00001.44521.00001.03760.96602.5227
250.97502.12121.00001.33971.00000.99450.96102.0065
500.96901.58070.99801.14701.00000.90770.96501.4882
1000.96301.17160.97900.96430.99900.80900.95901.1481
(1, 1, 1)150.98603.03480.99602.27820.99901.79530.97002.9185
250.97602.40570.98201.96910.99901.66060.96202.3558
500.97901.82040.95401.61830.98901.43160.97501.7843
1000.95701.36050.95701.30580.97801.15280.95401.3459
(1, 2, 1)150.98203.03530.99702.26021.00001.81190.97002.8891
250.97802.44140.97401.97520.99501.64440.97302.3516
500.97801.81430.96101.63800.98501.40570.96301.7875
1000.95801.36420.94801.30700.98201.15190.94501.3473
Table A5. MLE and Bayesian estimates and mean squared errors of λ for the CGZTP distributions.
Table A5. MLE and Bayesian estimates and mean squared errors of λ for the CGZTP distributions.
( λ ,   α ,   β ) n MLEPrior 1Prior 2Jeffreys
EstMSEEstMSEEstMSEEstMSE
(0.5, 2, 0.5)150.94920.66560.43310.04220.48680.01871.12090.6926
250.75980.34790.45180.05050.49080.01640.91100.3663
500.64180.18430.47910.05070.49690.00980.69940.1555
1000.54520.09140.48750.04130.49800.00660.59540.0810
(0.5, 2, 1)150.93350.66150.49330.04560.49320.01861.10180.6776
250.77840.38770.46590.04590.49270.01560.90680.3476
500.64260.18630.45350.04610.49440.01070.71640.1577
1000.55180.09500.44160.04490.50250.00710.61430.0864
(1, 1, 1)151.27110.78900.88920.20520.95250.07271.40250.6383
251.11780.40820.89080.19290.97370.07381.19950.3194
501.06210.24680.87590.15590.97000.08221.10130.2088
1000.98950.11850.90030.10700.99870.06161.03070.1115
(1, 2, 1)151.29150.75550.92630.21300.95320.08101.41410.6101
251.11820.42590.89490.19240.96690.07891.22230.3408
501.02820.21430.89250.15640.98140.07141.09490.1891
1001.01170.12880.91690.11200.98490.06291.04250.1101
Table A6. Coverage probabilities and average lengths of intervals of for the CGZTP distributions.
Table A6. Coverage probabilities and average lengths of intervals of for the CGZTP distributions.
( λ ,   α ,   β ) n MLEPrior 1Prior 2Jeffreys
CPALCPALCPALCPAL
(0.5, 2, 0.5)150.97002.72261.00001.47141.00001.03110.96302.5391
250.97802.11141.00001.34361.00000.99440.96801.9886
500.97101.56390.99701.13721.00000.91370.96801.4510
1000.98101.16680.98500.95910.99900.81500.96901.1202
(0.5, 2, 1)150.96802.70631.00001.48291.00001.04010.95902.5048
250.97302.12231.00001.31631.00000.98910.96201.9820
500.96901.56360.99901.14411.00000.91640.96701.4768
1000.98001.16880.97700.96450.99900.82240.96101.1323
(1, 1, 1)150.97203.02140.99302.21461.00001.82790.97502.8808
250.98602.41100.98201.97610.99601.64790.97302.2986
500.96601.83040.95901.62580.99701.42000.96601.7716
1000.95801.36070.94701.30190.97001.15510.96101.3438
(1, 2, 1)150.97503.04530.99602.27511.00001.81000.97302.8994
250.97502.40820.98901.98541.00001.65750.96702.3193
500.98101.82340.96601.64160.99401.41340.96601.7773
1000.94801.36340.94701.30400.98501.15420.95301.3458

References

  1. Niyomdecha, A.; Srisuradetchai, P.; Tulyanitikul, B. Gamma Zero-Truncated Poisson Distribution with the Minimum Compounded Function. Thail. Stat. Thail. 2023, 21, 863–886. [Google Scholar]
  2. Niyomdecha, A.; Srisuradetchai, P. Complementary Gamma Zero-Truncated Poisson Distribution and Its Application. Mathematics 2023, 11, 2584. [Google Scholar] [CrossRef]
  3. Xu, B.; Guo, Y.; Zhu, N. The Parameter Bayesian Estimation of Two-Parameter Exponential-Poisson Distribution and Its Optimal Property. J. Interdiscip. Math. 2016, 19, 697–707. [Google Scholar] [CrossRef]
  4. Yan, W.-A.; Shi, Y.-M.; Liu, Y. Bayesian Estination of Exponential Poisson Distribution under Different Loss Function. Fire Control Command Control 2012, 2, 124–126,131. [Google Scholar]
  5. Pathak, A.; Kumar, M.; Singh, S.K.; Singh, U. Bayesian Inference: Weibull Poisson Model for Censored Data Using the Expectation–Maximization Algorithm and Its Application to Bladder Cancer Data. J. Appl. Stat. 2022, 49, 926–948. [Google Scholar] [CrossRef] [PubMed]
  6. Davidson-Pilon, C. Bayesian Methods for Hackers: Probabilistic Programming and Bayesian Inference; Addison-Wesley Educational: Boston, MA, USA, 2015; pp. 127–130. [Google Scholar]
  7. Cai, Y.; Gui, W. Classical and Bayesian Inference for a Progressive First-Failure Censored Left-Truncated Normal Distribution. Symmetry 2021, 13, 490. [Google Scholar] [CrossRef]
  8. Elbatal, I.; Alotaibi, N.; Alyami, S.A.; Elgarhy, M.; El-Saeed, A.R. Bayesian and Non-Bayesian Estimation of the Nadaraj ah–Haghighi Distribution: Using Progressive Type-1 Censoring Scheme. Mathematics 2022, 10, 760. [Google Scholar] [CrossRef]
  9. Eliwa, M.S.; EL-Sagheer, R.M.; El-Essawy, S.H.; Almohaimeed, B.; Alshammari, F.S.; El-Morshedy, M. General Entropy with Bayes Techniques under Lindley and MCMC for Estimating the New Weibull–Pareto Parameters: Theory and Application. Symmetry 2022, 14, 2395. [Google Scholar] [CrossRef]
  10. Abdel-Aty, Y.; Kayid, M.; Alomani, G. Generalized Bayes Prediction Study Based on Joint Type-II Censoring. Axioms 2023, 12, 716. [Google Scholar] [CrossRef]
  11. Smith, A.F.M.; Roberts, G.O. Bayesian Computation via the Gibbs Sampler and Related Markov Chain Monte Carlo Methods. J. R. Stat. Soc. 1993, 55, 3–23. [Google Scholar] [CrossRef]
  12. Brooks, S. Markov Chain Monte Carlo Method and Its Application. J. Royal. Statistical. Soc. D 1998, 47, 69–100. [Google Scholar] [CrossRef]
  13. Almetwally, E.M.; Alotaibi, R.; Rezk, H. Estimation and Prediction for Alpha-Power Weibull Distribution Based on Hybrid Censoring. Symmetry 2023, 15, 1687. [Google Scholar] [CrossRef]
  14. Shakhatreh, M.K.; Aljarrah, M.A. Bayesian Analysis of Unit Log-Logistic Distribution Using Non-Informative Priors. Mathematics 2023, 11, 4947. [Google Scholar] [CrossRef]
  15. EL-Sagheer, R.M.; Almuqrin, M.A.; El-Morshedy, M.; Eliwa, M.S.; Eissa, F.H.; Abdo, D.A. Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model. Symmetry 2022, 14, 956. [Google Scholar] [CrossRef]
  16. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis, 3rd ed.; Chapman & Hall/CRC: Philadelphia, PA, USA, 2013; pp. 43–44. [Google Scholar]
  17. Fink, D. A Compendium of Conjugate Priors. Available online: https://courses.physics.ucsd.edu/2018/Fall/physics210b/REFERENCES/conjugate_priors.pdf (accessed on 25 December 2023).
  18. Naji, L.F.; Rasheed, H.A. Bayesian Estimation for Two Parameters of Gamma Distribution under Precautionary Loss Function. Ibn AL-Haitham J. Pure Appl. Sci. 2019, 32, 187–196. [Google Scholar] [CrossRef]
  19. Moala, F.A.; Ramos, P.L.; Achcar, J.A. Bayesian Inference for Two-Parameter Gamma Distribution Assuming Different Non-Informative Priors. Rev. Colomb. Estad. 2013, 36, 319–336. [Google Scholar]
  20. Pradhan, B.; Kundu, D. Bayes Estimation and Prediction of the Two-Parameter Gamma Distribution. J. Stat. Comput. Simul. 2011, 81, 1187–1198. [Google Scholar] [CrossRef]
  21. Anastasiou, A.; Gaunt, R.E. Wasserstein Distance Error Bounds for the Multivariate Normal Approximation of the Maximum Likelihood Estimator. Electron. J. Stat. 2021, 15, 5758–5810. [Google Scholar] [CrossRef]
  22. Srisuradetchai, P.; Niyomdecha, A.; Phaphan, W. Wald Intervals via Profile Likelihood for the Mean of the Inverse Gaussian Distribution. Symmetry 2024, 16, 93. [Google Scholar] [CrossRef]
  23. Jeffreys, H. Theory of Probability; Oxford University Press: London, UK, 1983; ISBN 9780198531937. [Google Scholar]
  24. Panahi, H.; Moradi, N. Estimation of the Inverted Exponentiated Rayleigh Distribution Based on Adaptive Type II Progressive Hybrid Censored Sample. J. Comput. Appl. Math. 2020, 364, 112345. [Google Scholar] [CrossRef]
  25. Dutta, S.; Kayal, S. Inference of a Competing Risks Model with Partially Observed Failure Causes under Improved Adaptive Type-II Progressive Censoring. Proc. Inst. Mech. Eng. O J. Risk Reliab. 2023, 237, 765–780. [Google Scholar] [CrossRef]
  26. Junnumtuam, S.; Niwitpong, S.-A.; Niwitpong, S. Bayesian Computation for the Parameters of a Zero-Inflated Cosine Geometric Distribution with Application to COVID-19 Pandemic Data. Comput. Model. Eng. Sci. 2023, 135, 1229–1254. [Google Scholar] [CrossRef]
  27. Bédard, M. On the Robustness of Optimal Scaling for Random Walk Metropolis Algorithms. Ph.D. Thesis, University of Toronto Libraries, Toronto, ON, Canada, 2006. [Google Scholar]
  28. Hinkley, D. On Quick Choice of Power Transformation. J. R. Stat. Soc. Ser. C Appl. Stat. 1977, 26, 67–69. [Google Scholar] [CrossRef]
  29. Lee, E.T.; Wang, J.W. Statistical Methods for Survival Data Analysis, 4th ed.; John Wiley & Sons: Nashville, TN, USA, 2013; ISBN 9781118095027. [Google Scholar]
Figure 1. Probability density functions of GZTP (left) and CGZTP (Right) for various parameter sets θ = ( λ , α , β ) .
Figure 1. Probability density functions of GZTP (left) and CGZTP (Right) for various parameter sets θ = ( λ , α , β ) .
Symmetry 16 00417 g001
Figure 2. Gamma prior distributions with different hyperparameters.
Figure 2. Gamma prior distributions with different hyperparameters.
Symmetry 16 00417 g002
Figure 3. The trace plots of α i and β i from the joint posterior distributions corresponding to: (a) Prior 3 for GZTP with λ = 1 ,   α = 2 ,   β = 1 , and n = 50 , (b) Prior 2 for CGZTP with λ = 0.5 ,   α = 2 ,   β = 0.5 , and n = 25 .
Figure 3. The trace plots of α i and β i from the joint posterior distributions corresponding to: (a) Prior 3 for GZTP with λ = 1 ,   α = 2 ,   β = 1 , and n = 50 , (b) Prior 2 for CGZTP with λ = 0.5 ,   α = 2 ,   β = 0.5 , and n = 25 .
Symmetry 16 00417 g003
Figure 4. The pair plots of α i and β i from the joint posterior distributions corresponding to: (a) Prior 3 for GZTP with λ = 1 ,   α = 2 ,   β = 1 , and n = 50 , (b) Prior 2 for CGZTP with λ = 0.5 ,   α = 2 ,   β = 0.5 , and n = 25 .
Figure 4. The pair plots of α i and β i from the joint posterior distributions corresponding to: (a) Prior 3 for GZTP with λ = 1 ,   α = 2 ,   β = 1 , and n = 50 , (b) Prior 2 for CGZTP with λ = 0.5 ,   α = 2 ,   β = 0.5 , and n = 25 .
Symmetry 16 00417 g004
Figure 5. The trace plots of λ i from posterior distributions corresponding to: (a) Prior 2 for GZTP with λ = 1 ,   α = 1 ,   β = 1 , and n = 50 , (b) Prior 1 for CGZTP with λ = 0.5 ,   α = 2 ,   β = 1 , and n = 50 .
Figure 5. The trace plots of λ i from posterior distributions corresponding to: (a) Prior 2 for GZTP with λ = 1 ,   α = 1 ,   β = 1 , and n = 50 , (b) Prior 1 for CGZTP with λ = 0.5 ,   α = 2 ,   β = 1 , and n = 50 .
Symmetry 16 00417 g005
Figure 6. The average of estimates, MSEs, CPs, and ALs for α and β of the two-parameter GZTP distribution with λ = 0.5 ,   α = 2 and β = 0.5 . Column (a) presents results for the parameter α, and column (b) for the parameter β.
Figure 6. The average of estimates, MSEs, CPs, and ALs for α and β of the two-parameter GZTP distribution with λ = 0.5 ,   α = 2 and β = 0.5 . Column (a) presents results for the parameter α, and column (b) for the parameter β.
Symmetry 16 00417 g006
Figure 7. The average of estimates, MSEs, CPs, and Als for α and β of the two-parameter CGZTP distribution with λ = 1 ,   α = 1 ,   β = 1 . Column (a) presents results for the parameter α, and column (b) for the parameter β.
Figure 7. The average of estimates, MSEs, CPs, and Als for α and β of the two-parameter CGZTP distribution with λ = 1 ,   α = 1 ,   β = 1 . Column (a) presents results for the parameter α, and column (b) for the parameter β.
Symmetry 16 00417 g007
Figure 8. The averages of estimates, MSEs, CP, and ALs for λ of the one-parameter GZTP distribution with λ = 0.5 ,   α = 2 , and β = 1 .
Figure 8. The averages of estimates, MSEs, CP, and ALs for λ of the one-parameter GZTP distribution with λ = 0.5 ,   α = 2 , and β = 1 .
Symmetry 16 00417 g008
Figure 9. The average of estimates, MSEs, CPs, and ALs for λ of the one-parameter CGZTP distribution with λ = 1 ,   α = 1 , and β = 1 .
Figure 9. The average of estimates, MSEs, CPs, and ALs for λ of the one-parameter CGZTP distribution with λ = 1 ,   α = 1 , and β = 1 .
Symmetry 16 00417 g009
Figure 10. Histogram with fitted GZTP distribution for March precipitation data.
Figure 10. Histogram with fitted GZTP distribution for March precipitation data.
Symmetry 16 00417 g010
Figure 11. The pair plot of α and β from the RWM algorithm for March precipitation data.
Figure 11. The pair plot of α and β from the RWM algorithm for March precipitation data.
Symmetry 16 00417 g011
Figure 12. Histogram with fitted CGZTP distribution for remission time data.
Figure 12. Histogram with fitted CGZTP distribution for remission time data.
Symmetry 16 00417 g012
Figure 13. The pair plot of α and β from the RWM algorithm for remission time of bladder cancer patients.
Figure 13. The pair plot of α and β from the RWM algorithm for remission time of bladder cancer patients.
Symmetry 16 00417 g013
Table 1. The prior distributions of parameters α and β for the GZTP and CGZTP distributions.
Table 1. The prior distributions of parameters α and β for the GZTP and CGZTP distributions.
α G a m m a ( a , b ) ,   β G a m m a ( c , d ) Informative
Priors
Hyperparameter Values
True ValuesVariances a b c d
α β α β
20.54 (High)0.5 (High)Prior 110.50.51
2 (Low)0.125 (Low)Prior 22124
214 (High)2 (High)Prior 110.50.50.5
2 (Low)0.5 (Low)Prior 22122
112 (High)2 (High)Prior 10.50.50.50.5
0.5 (Low)0.5 (Low)Prior 22222
Table 2. The prior distributions of parameter λ for the GZTP and CGZTP distributions.
Table 2. The prior distributions of parameter λ for the GZTP and CGZTP distributions.
λ G a m m a ( m , n ) Informative
Priors
Hyperparameter Values
True ValuesVariances m n
0.50.25 (High)Prior 112
0.083 (Low)Prior 236
11 (High)Prior 111
0.333 (Low)Prior 233
Table 3. MLE and Bayesian estimates and mean squared errors of α and β for the GZTP distributions.
Table 3. MLE and Bayesian estimates and mean squared errors of α and β for the GZTP distributions.
( λ ,   α ,   β ) n ParameterMLEPrior 1Prior 2
EstMSEEstMSEEstMSE
(0.5, 2, 0.5)15 α 2.36370.88962.19420.49562.14590.2622
β 0.62400.08870.57340.05450.55500.0252
25 α 2.22840.44362.13280.27772.09480.1854
β 0.57700.04530.54700.02820.53380.0179
50 α 2.10620.16682.09260.15162.08000.1235
β 0.53700.01650.52990.01440.52560.0115
100 α 2.06440.07602.05080.07012.04700.0639
β 0.52180.00720.51810.00650.51680.0059
(0.5, 2, 1)15 α 2.41561.07622.17800.48742.11730.2686
β 1.27440.44851.14080.20681.08210.0921
25 α 2.24380.47682.1530.29412.11230.1933
β 1.15010.1811.09770.11511.07010.0724
50 α 2.10560.16392.07620.15172.06350.1230
β 1.06620.06461.05210.05931.04360.0472
100 α 2.06010.08462.03910.07012.03630.0638
β 1.04040.03171.02410.02631.02210.0237
(1, 1, 1)15 α 1.17370.18481.09360.11191.06730.0566
β 1.33080.57071.19280.29311.14080.1377
25 α 1.10250.09461.06390.06261.04280.0398
β 1.18970.26541.13930.17801.09490.1005
50 α 1.05460.04101.03590.03061.02970.0246
β 1.08310.08251.07430.07731.06130.0583
100 α 1.02220.01571.01330.01261.0120.0114
β 1.03350.03471.02830.03201.02550.0283
(1, 2, 1)15 α 2.41311.10792.15650.41292.11100.2314
β 1.29950.52991.13950.19841.10400.1045
25 α 2.21860.35952.09830.27052.06330.1757
β 1.14920.16421.08890.12771.06150.0782
50 α 2.11230.16172.09420.14122.08060.1149
β 1.07620.07321.06620.06091.05620.0480
100 α 2.05940.06842.02330.06582.02070.0596
β 1.04050.02931.01910.02771.01720.0248
Table 4. Coverage probabilities and average lengths of intervals for α and β of the GZTP distributions.
Table 4. Coverage probabilities and average lengths of intervals for α and β of the GZTP distributions.
( λ ,   α ,   β ) n ParametersMLEPrior 1Prior 2
CPALCPALCPAL
(0.5, 2, 0.5)15 α 0.97303.06060.96852.64280.98462.3063
β 0.97700.96180.96850.82510.98850.7076
25 α 0.95702.22540.96902.02940.98401.8566
β 0.96200.69190.96400.62730.97100.5675
50 α 0.95501.47930.95001.43530.95801.3681
β 0.95600.45700.94800.43970.95600.4172
100 α 0.95501.02350.94601.00190.95000.9777
β 0.95700.31450.94600.30770.95500.3004
(0.5, 2, 1)15 α 0.96603.13360.95772.61520.97002.2779
β 0.96601.96180.96081.63860.98311.3854
25 α 0.96102.24230.97002.05190.98101.8675
β 0.96101.37870.96001.25650.97501.1329
50 α 0.95901.47890.95801.42530.96401.3588
β 0.95100.90740.94800.87430.96400.8314
100 α 0.95001.02100.94200.99760.95200.9753
β 0.95300.62710.94800.60980.94900.5949
(1, 1, 1)15 α 0.97701.37580.96541.22050.98311.0506
β 0.96902.31370.97081.95580.98311.6251
25 α 0.95800.99160.96000.92960.97000.8376
β 0.96301.61890.95501.49390.97201.3081
50 α 0.95100.69560.94200.64270.94700.6119
β 0.95500.99730.95701.02090.96600.9580
100 α 0.94400.47480.95600.44500.96100.4338
β 0.94500.67600.95300.70080.96200.6791
(1, 2, 1)15 α 0.96503.04600.97002.51650.98462.1990
β 0.96202.07390.97151.69880.98921.4562
25 α 0.96902.15010.96301.93750.97401.7639
β 0.97101.43210.96001.29590.97501.1647
50 α 0.95601.43930.94901.39460.95501.3289
β 0.95800.95200.95600.91930.96400.8700
100 α 0.95800.98900.94900.95750.95200.9363
β 0.95600.65210.95500.63010.96100.6140
Table 5. March precipitation data.
Table 5. March precipitation data.
0.771.740.811.201.951.200.471.433.372.20
3.003.091.512.100.521.621.310.320.590.81
2.811.871.181.354.752.480.961.890.902.05
Table 6. Bayesian and Maximum Likelihood Estimates with standard errors, 95% Wald CIs and Bayesian credible intervals, and goodness-of-fit testing results for two datasets.
Table 6. Bayesian and Maximum Likelihood Estimates with standard errors, 95% Wald CIs and Bayesian credible intervals, and goodness-of-fit testing results for two datasets.
DatasetMethodParameterEstimateStandard
Error
IntervalK-Sp-Value
March precipitation Bayes α 3.02890.6281(1.8884, 4.3697)0.06230.9998
β 1.70550.4089(0.9805, 2.6144)
MLE α 3.05970.7270(1.6347, 4.4847)0.06020.9999
β 1.72140.4656(0.8086, 2.6341)
Remission TimeBayes α 1.16970.1292(0.9283, 1.4277)0.07130.5341
β 0.12460.0175(0.0921, 0.1596)
MLE α 1.17190.1312(0.9147, 1.4290)0.07390.4866
β 0.12580.0174(0.0919, 0.1599)
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

Srisuradetchai, P.; Niyomdecha, A. Bayesian Inference for the Gamma Zero-Truncated Poisson Distribution with an Application to Real Data. Symmetry 2024, 16, 417. https://doi.org/10.3390/sym16040417

AMA Style

Srisuradetchai P, Niyomdecha A. Bayesian Inference for the Gamma Zero-Truncated Poisson Distribution with an Application to Real Data. Symmetry. 2024; 16(4):417. https://doi.org/10.3390/sym16040417

Chicago/Turabian Style

Srisuradetchai, Patchanok, and Ausaina Niyomdecha. 2024. "Bayesian Inference for the Gamma Zero-Truncated Poisson Distribution with an Application to Real Data" Symmetry 16, no. 4: 417. https://doi.org/10.3390/sym16040417

APA Style

Srisuradetchai, P., & Niyomdecha, A. (2024). Bayesian Inference for the Gamma Zero-Truncated Poisson Distribution with an Application to Real Data. Symmetry, 16(4), 417. https://doi.org/10.3390/sym16040417

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