Abstract
For the hierarchical inverse gamma and inverse gamma model, we calculate the Bayes posterior estimator of the rate parameter of the inverse gamma distribution under Stein's loss function which penalizes gross overestimation and gross underestimation equally and the corresponding Posterior Expected Stein's Loss (PESL). We also obtain the Bayes posterior estimator of the rate parameter under the squared error loss and the corresponding PESL. Moreover, we obtain the empirical Bayes estimators of the rate parameter of the inverse gamma distribution with a conjugate inverse gamma prior by the moment and MLE methods. In numerical simulations, we have illustrated four aspects: Consistency, goodness-of-fit, comparison, and marginal densities. The numerical results indicate that the MLEs are better than the moment estimators when estimating the hyperparameters. Finally, the model could potentially be used to fit right skewed data, not left-skewed data.
The authors would like to thank the Editor, the Associate Editor, and the referees for their constructive comments, which have led to a substantial improvement of the article.
Disclosure statement
No potential conflict of interest was reported by the author(s).