“…Via the "coda" package (by Plummer et al [20]) in R 4.2.2 programming software, to obtain the Bayes point estimates along with their HPD interval estimates of the same unknown parameters, we simulated 12,000 MCMC samples and ignored the first 2000 iterations as burn-in. According to the prior mean and prior variance criteria, two sets called Prior-I and -II of the hyperparameters (a 1 , a 2 , b 1 , b 2 ) were considered as (2.5, 7.5, 5, 5) and (5,15,10,10), respectively. Specifically, the average point estimates (APEs) of δ, θ, R(t), or h(t) (say Ω) were given by…”