This paper investigates the estimation of the stress–strength reliability parameter R=P(Y<X) when both stress and strength follow independent Garhy distributions under progressive Type-II censoring schemes. A closed-form expression for R is explicitly derived, enabling effective and precise calculation without numerical integration. The Garhy distribution, a flexible one-parameter lifetime model with an increasing hazard function, is confirmed by full-scale goodness-of-fit diagnostics. A Bayesian estimation model is trained on non-informative priors (normal and extended Jeffreys priors) under squared error loss. The posterior expectations are analytically intractable; we adopt two complementary methods of computation: (i) Markov Chain Monte Carlo (MCMC) using the Metropolis–Hastings algorithm and (ii) the Tierney–Kadane (TK) approximation, which provides extremely precise analytical estimates with significantly reduced computational burden. Monte Carlo simulations are large-scale and compare the proposed estimators under different censoring schemes, sample sizes, and parameter configurations in terms of bias and mean squared error (MSE). The methodology is further applied to a real medical dataset comprising kidney dialysis patient survival times, demonstrating its practical relevance in clinical reliability assessment. Results consistently indicate that Bayesian methods, particularly with the extended Jeffreys prior, outperform classical MLEs in terms of stability and accuracy, especially under heavy censoring. Moreover, the TK approximation yields estimates virtually identical to MCMC while requiring only a fraction of the computational effort. We further extend the TK framework to approximate the posterior variance of R and the expected log-likelihood, providing a fully analytical alternative to MCMC for comprehensive Bayesian inference.
Alenezy et al. (Thu,) studied this question.