Abstract
This article develops a procedure to obtain highly accurate confidence interval estimates for the stress-strength reliability R = P(X > Y) where X and Y are data from independent normal distributions of unknown means and variances. Our method is based on third-order likelihood analysis and is compared to the conventional first-order likelihood ratio procedure as well as the approximate methods of Reiser and Guttman (Citation1986) and Guo and Krishnamoorthy (Citation2004). The use of our proposed method is illustrated by an empirical example and its superior accuracy in terms of coverage probability and error rate are examined through Monte Carlo simulation studies.
Notes
Note that the Reiser and Guttman (Citation1986) and the Guo and Krishnamoorthy (Citation2004) methods use s2x and s2y in their estimation of σ2x and σ2y instead of and
like the MLE.
In contrast to the Reiser & Guttman (Citation1986) and Guo & Krishnamoorthy (Citation2004) methods, whose rates of convergence are unknown.
Practically, the confidence interval I1 − αψ is obtained by numerically solving for ψ in the inequality |r*(ψ)| ⩽ zα/2 using a sufficient number of grid points of ψ chosen in an appropriate range.