Let X and Y follow two independent Burr type XII distributions and δ=P(X<Y) . If X is the stress that is applied to a certain component and Y is the strength to sustain the stress, then δ is called the stress–strength parameter. In this study, The Bayes estimator of δ is investigated based on a progressively first failure-censored sample. Because of computation complexity and no closed form for the estimator as well as posterior distributions, the Markov Chain Monte Carlo procedure using the Metropolis–Hastings algorithm via Gibbs sampling is built to collect a random sample of δ via the joint distribution of the progressively first failure-censored sample and random parameters and the empirical distribution of this collected sample is used to estimate the posterior distribution of δ . Then, the Bayes estimates of δ using the square error, absolute error, and linear exponential error loss functions are obtained and the credible interval of δ is constructed using the empirical distribution. An intensive simulation study is conducted to investigate the performance of these three types of Bayes estimates and the coverage probabilities and average lengths of the credible interval of δ . Moreover, the performance of the Bayes estimates is compared with the maximum likelihood estimates. The Internet of Things and a numerical example about the miles-to-failure of vehicle components for reliability evaluation are provided for application purposes.
