# Inference on Reliability of Stress-Strength Models for Poisson Data

Authors: Alessandro Barbiero
Journal: Journal of Quality and Reliability Engineering
 (2013)

---

## Abstract

Researchers in reliability engineering regularly encounter variables that are discrete in nature, such as the number of events (e.g., failures) occurring in a certain spatial or temporal interval. The methods for analyzing and interpreting such data are often based on asymptotic theory, so that when the sample size is not large, their accuracy is suspect. This paper discusses statistical inference for the reliability of stress-strength models when stress and strength are independent Poisson random variables. The maximum likelihood estimator and the uniformly minimum variance unbiased estimator are here presented and empirically compared in terms of their mean square error; recalling the delta method, confidence intervals based on these point estimators are proposed, and their reliance is investigated through a simulation study, which assesses their performance in terms of coverage rate and average length under several scenarios and for various sample sizes. The study indicates that the two estimators possess similar properties, and the accuracy of these estimators is still satisfactory even when the sample size is small. An application to an engineering experiment is also provided to elucidate the use of the proposed methods.


---

## Body

## 1. Introduction

A stress-strength model, in the simplest terms, considers a unit/system that is subjected to an external stress, modeled by r.v.X, against which the unit sets its own strength, modeled by r.v. Y, in order to properly operate. The probability that the unit withstands the stress is then given by R=P(X<Y), which is usually called reliability.A great deal of work has been done about this topic: most of it deals with the computation of reliability, if the distributions of stress and strength are known, or its estimation under various parametric assumptions onX and Y, when samples from X and Y are available. A complete review is available in [1]. Many applications of the stress-strength model, for its own nature, are related to engineering or military problems, where it is also referred to as a load-strength model [2]. However, there are also natural applications in medicine or psychology, which involve the comparison of two r.v., representing, for example, the effect of a specific drug or treatment administered to two groups (control and test); here, reliability assumes a wider meaning.Almost all of these papers consider continuous distributions forX and Y, since many practical applications of the stress-strength model in engineering fields presuppose continuous quantitative data. A relatively small amount of work is devoted to discrete or categorical data. Data may be discrete by nature, for example, the number of events occurring in a certain spatial or temporal interval; sometimes discrete data are derived from continuous ones by grouping or discretization or censoring, and then, instead of numerical measurements on X and Y, they are presented in a form of ordered categories.Among the r.v. modeling discrete data, the Poisson can be of interest in several practical applications. The Poisson r.v. is often used to model rare events such as the number of claims in automobile insurance, the number of times a website is accessed, the number of calls to a phone operator, the number of words mistyped per page in a book, and so forth [3, 4]. The distribution of the difference between two independent r.v. each having a Poisson distribution has already attracted some attention [5]. Strackee and van der Gon [6] stated that “in a steady state the number of light quanta, emitted or absorbed in a definite time, is distributed according to a Poisson distribution. In view thereof, the physical limit of perceptible contrast in vision can be studied in terms of the difference between two independent variates each following a Poisson distribution”. Irwin [7] studied the case when the two variables X and Y each have the same expected value; Skellam [8] was the first to discuss the problem when E(X)=λ1≠E(Y)=λ2. Strackee and van der Gon [6] gave tables of the approximate values of the cumulative probability P(Y-X<d) for several combinations of the values of the parameters λ1 and λ2. More recently, Karlis and Ntzoufras [9] used the Poisson difference distribution to model the difference in the decayed, missing, and filled teeth index before and after treatment; Karlis and Ntzoufras [10] applied it to model the difference in the number of goals in football games.In this paper, we examine point and interval estimation for the reliability of the stress-strength model with independent Poisson stress and strength. Although the maximum likelihood (ML) and uniformly minimum variance unbiased (UMVU) estimators of reliability have a known analytical expression, their statistical properties cannot be easily derived and thus need to be assessed through a Monte Carlo simulation study. Confidence intervals for reliability based on approximate expression for variance are also presented, and their performances in terms of coverage rate and average width are empirically investigated.The paper is laid out as follows: in Section2 reliability for Poisson stress-strength model and its ML and UMVU estimators are presented and discussed. Section 3 introduces approximate variance estimators and confidence intervals for reliability. Section 4 is devoted to a Monte Carlo (MC) study, which empirically investigates the performance of ML and UMVU estimators, and the corresponding confidence intervals for different combinations of distributional parameters and sample sizes. Section 5 describes an application, and Section 6 gives final remarks.


## 2. Point Estimators

LetX and Y be independent r.v. modeling stress and strength, respectively, with X~Poisson(λ1) and Y~Poisson(λ2). Then, the reliability R=P(X<Y) of the stress-strength model is given by (see [1])
(1)R=P(X<Y)=∑x=0+∞‍λ1xe-λ1x![1-∑y=0x‍λ2ye-λ2y!]=limk→∞∑x=0k‍λ1xe-λ1x![1-∑y=0x‍λ2ye-λ2y!].The terms of the external sum rapidly converge to zero: reliability can be actually computed taking into account only its first terms. As an example, we compute the reliabilityR when λ1=1 and λ2=4; the partial sums are reported in Table 1: the value of R is already stable at the 7th decimal digit when k=10.Table 1
Partial sums for the computation ofR for a Poisson stress-strength model (λ1=1,λ2=4).
k
=
1
k
=
2
k
=
3
k
=
4
k
=
5
0.6953312
0.8354743
0.87021
0.8758994
0.876558
k
=
6
k
=
7
k
=
8
k
=
9
k
=
10
0.8766146
0.8766183
0.8766185
0.8766186
0.8766186If two simple random samplesx of size n1 and y of size n2 from X and Y, respectively, are available, reliability can be estimated with the ML estimator, obtained by substituting in (1) the maximum likelihood estimators of the unknown parameters λ1 and λ2:
(2)R~=∑x=0+∞‍x-xe-x-x![1-∑y=0x‍y-ye-y-y!].Otherwise, one can use the UMVU estimator [1]:
(3)R^=∑x=0U‍(n1x-x)(n1-1)n1x--xn1n1x-[1-∑y=0x‍(n2y-y)(n2-1)n2y--yn2n2y-],
where U=min(n1x-,n2y--1). Note that formula (3) is represented via a finite sum, whereas formula (2) contains a rapidly converging series. The number of calculations that formula (3) performs depends on the sample means and the sample sizes, which jointly define the number of terms of the external sum; in formula (2) the terms of the external sum rapidly converge to zero, so that it may practically need fewer calculations than (3).In Table2, the values for the UMVU estimator are reported when x-=1 and y-=2;3;4, for different combinations of sample sizes n1 and n2. Note that the values of R^ are very close to the value of R~ even for small sample sizes and get closer as the sample sizes increase. These results are pictorially displayed in Figure 1 for x-=1 and y-=4.Table 2
Values of the UMVU estimator ofR for a Poisson stress-strength model (x-=1,y-=2;3;4) with varying sample sizes.
n
1
n
2
x
-
=
1
,
y
-
=
2
x
-
=
1
,
y
-
=
3
x
-
=
1
,
y
-
=
4
10
10
0.612649
0.786485
0.887889
10
20
0.607790
0.781051
0.883313
10
30
0.606199
0.779271
0.881805
20
10
0.613908
0.785907
0.886621
20
20
0.609078
0.780562
0.882117
20
30
0.607496
0.778811
0.880632
30
10
0.614334
0.785732
0.886212
30
20
0.609512
0.780416
0.881731
30
30
0.607932
0.778674
0.880254
50
50
0.607031
0.777192
0.878786
100
100
0.606364
0.776097
0.877697
R
~
0.605703
0.775015
0.876619Figure 1
Values of the UMVU estimator ofR for a Poisson stress-strength model (x-=1,y-=4) with varying sample sizes. The ML estimator is represented by the horizontal line.Due to the complex expressions involved, the bias ofR~ and the variance of either estimators R~ and R^ cannot be analytically derived; a comparison of their performance (in terms of mean square error) can be carried out through MC simulations.


## 3. Variance Estimators and Confidence Intervals

Whereas the exact value of the variance or the mean square error of either estimator introduced in Section2 is almost impracticable to derive, an approximate value can be easily supplied recalling the delta method [11]. For the ML estimator R~(x-,y-), since X- and Y- are independent estimators of λ1 and λ2, the variance of R~ can be approximated as
(4)V(R~)≈(∂R~∂x-)2|x-=λ1,y-=λ2V(x-)V(R~)+(∂R~∂y-)2|x-=λ1,y-=λ2V(y-),
with V(x-)=λ1/n1, V(y-)=λ2/n2 and
(5)∂R~∂x-=∑x=0∞‍e-x-x-x-1x!(x-x-)[1-∑y=0x‍y-ye-y-y!],∂R~∂y-=-∑x=0∞‍x-xe-x-x!∑y=0x‍e-y-y-y-1y!(y-y-).An analogous approximation can be carried out for the variance of the UMVU estimator; remembering thatx!=Γ(x+1), R^ can be rewritten as
(6)R^=∑x=0U‍Γ(n1x-+1)Γ(x+1)Γ(n1x--x+1)(n1-1)n1x--xn1n1x-×[1-∑y=0x‍Γ(n2y-+1)Γ(y+1)Γ(n2y--y+1)(n2-1)n2y--yn2n2y-],
and then the two first-order partial derivatives are given by
(7)∂R^∂x-=∑x=0U[1-∑y=0x‍(n2y-y)(n2-1)n2y--yn2n2y-]×1(n1-1)xn1(n1x-x)(n1-1n1)n1x-×[Γ′(n1x-+1)Γ(n1x-+1)-Γ′(n1x--x+1)Γ(n1x--x+1)+log(n1-1n1)],∂R^∂y-=-∑x=0U∑y=0x(n1x-x)(n1-1)n1x--xn1n1x-×1(n2-1)yn2(n2y-y)(n2-1n2)n2y-×[Γ′(n2y-+1)Γ(n2y-+1)-Γ′(n2y--y+1)Γ(n2y--y+1)+log(n2-1n2)].The approximate variances ofR~ and R^ derived through the delta method can be estimated substituting in (4) the sample means to the unknown parameters and thus getting
(8)v(R~)=(∂R~∂x-)2x-n1+(∂R~∂y-)2y-n2
and an analogous result for v(R^).The Gamma functionΓ(x) and its first derivative, Γ′(x)=∫0+∞‍tx-1e-tlogtdt, involved in the partial derivatives of R^, have to be numerically computed. In the R software environment [12] this task is easily performed through the gamma and digamma functions, the latter providing the ratio Γ′(x)/Γ(x).Once one has computedv(R~), an approximate (1-α)·100% confidence interval for R can be built, recalling the asymptotic normality of R~:
(9)(R~+zα/2v(R~),R~+z1-α/2v(R~)),
and in an analogous way for R^. Since R is bounded in [0,1], special care has to be given when R^ is close to one (close to zero) and/or sample sizes are small: the upper bound may exceed one (the lower bound may fall below zero), and then the CI in (9) will be modified as follows:
(10)(max(0,R~+zα/2v(R~)),5min(1,R~+z1-α/2v(R~))).More sophisticated asymptotic confidence intervals forR can be built recalling some normalizing transformations, such as logit and arcsine [13].


## 4. Simulation Study

The simulation study aims at empirically comparing the performance of the ML and UMVU estimators, in terms of bias and mean square error, and the confidence intervals based on them, in terms of the coverage rate and average length. Since the approximation of the variance derived through the delta method (4) holds for large samples, we will investigate to what extent it still holds for small and moderate sample sizes, and how it affects inferential results. In this MC study, the value of the parameter λ1 of the Poisson distribution for stress X is first set equal to a “reference” value, 1, and the parameter λ2 of the Poisson distribution modeling strength is allowed to vary in order to obtain four different levels of reliability R, namely, 0.5, 0.6, 0.7, 0.8, and 0.9. Note that a value of λ2=1.547 is needed in order to get R=P(X<Y)=0.5 while λ2=λ1=1 lead only to P(X<Y)+0.5·P(X=Y)=0.5. Then, λ1 is set equal to greater values (namely 2, 5, and 10), and λ2 is allowed to vary in order to ensure the five values of reliability R above. The corresponding values of λ2 for each combination of R and λ1 values are reported in Table 3.Table 3
Parameter values for theX and Y distributions and corresponding reliability explored in the simulation study of a Poisson stress-strength model.
λ
1
λ
2
λ
1
λ
2
λ
1
λ
2
λ
1
λ
2
R
1
1.547
2
2.522
5
5.509
10
10.504
0.5
1
1.973
2
3.089
5
6.361
10
11.683
0.6
1
2.497
2
3.764
5
7.342
10
13.015
0.7
1
3.2
2
4.646
5
8.584
10
14.666
0.8
1
4.338
2
6.032
5
10.47
10
17.121
0.9For each couple(λ1,λ2), a huge number (S=2,000) of samples x of size n1 and y of size n2 are drawn from X~Poisson(λ1) and Y~Poisson(λ2) independently. Different and unequal sample sizes are here considered (all the nine possible combinations between the values n1=10,20,50, and n2=10,20,50). The ML and UMVU estimators are computed on each sample, their approximate variances are calculated, and the corresponding 95% confidence intervals for R are built. Some measures of performance for these estimators are supplied. In more detail, the MC root mean square error and the percentage relative bias of the ML estimator are provided:
(11)RMSEMC(R~)=1S∑s=1S‍(R~(s)-R)2RBMC(R~)=((1/S)∑s=1SR~(s)-R)R·100%,
where R~(s) denotes the value of R~ for the sth sample. Analogous indexes are derived for R^, whose bias is null, and for which we then expect the MC relative bias to be close to zero.Regarding estimating the variance, the true varianceV(R~) is approximated by its MC mean:
(12)V(R~)≈VMC(R~)=EMC[(R~-R~-)2]
with R~-=∑s=1S‍R~(s)/S, and then the MC relative bias and RMSE of v(R~) are calculated the same way as for R~.The MC coverage rate of the CIs is simply defined as follows:(13)1S∑s=1S‍I[R~(s)+zα/2v(R~(s))≤R≤R~(s)+z1-α/2v(R~(s))],
where I(E) is the indicator function, taking value 1 if E is true, 0 otherwise. The length of the confidence interval is then equal to 2z1-α/2v(R~(s)). The same performance indexes are derived for R^.The simulation results forλ1=1 are reported in Table 4 (RB and RMSE for ML and UMVU point estimators), Table 5 (RB and RMSE for variance estimators), and Table 6 (coverage rate and average length of confidence intervals).Table 4
Simulation results: bias and root mean square error of ML and UMVU estimators.
R
=
0.9
R
=
0.8
R
=
0.7
R
=
0.6
R
=
0.5
(
n
1
,
n
2
)
=
(
10,10
)
RB(R~)
-
1.214
-
1.539
-
1.539
-
1.216
-
0.932
RB(R^)
-
0.111
-
0.219
-
0.275
-
0.189
-
0.310
RMSE(R~)
0.061
0.091
0.109
0.119
0.124
RMSE(R^)
0.059
0.091
0.111
0.123
0.129
(
n
1
,
n
2
)
=
(
10,20
)
RB(R~)
-
0.753
-
0.906
-
0.796
-
0.492
-
0.073
RB(R^)
-
0.079
-
0.182
-
0.221
-
0.207
-
0.204
RMSE(R~)
0.049
0.075
0.091
0.100
0.103
RMSE(R^)
0.048
0.075
0.093
0.102
0.106
(
n
1
,
n
2
)
=
(
10,50
)
RB(R~)
-
0.386
-
0.349
-
0.138
0.197
0.673
RB(R^)
0.029
0.017
0.023
0.040
0.093
RMSE(R~)
0.041
0.063
0.077
0.084
0.086
RMSE(R^)
0.040
0.064
0.078
0.086
0.088
(
n
1
,
n
2
)
=
(
20,10
)
RB(R~)
-
1.151
-
1.576
-
1.676
-
1.665
-
1.598
RB(R^)
-
0.154
-
0.297
-
0.335
-
0.409
-
0.552
RMSE(R~)
0.055
0.082
0.098
0.109
0.113
RMSE(R^)
0.052
0.082
0.100
0.111
0.116
(
n
1
,
n
2
)
=
(
20,20
)
RB(R~)
-
0.557
-
0.725
-
0.716
-
0.564
-
0.259
RB(R^)
0.001
-
0.050
-
0.068
-
0.040
0.055
RMSE(R~)
0.040
0.062
0.077
0.084
0.088
RMSE(R^)
0.039
0.062
0.078
0.085
0.089
(
n
1
,
n
2
)
=
(
20,50
)
RB(R~)
-
0.280
-
0.313
-
0.235
-
0.060
0.121
RB(R^)
0.017
-
0.004
-
0.010
0.013
-
0.020
RMSE(R~)
0.031
0.049
0.060
0.066
0.068
RMSE(R^)
0.031
0.049
0.061
0.067
0.069
(
n
1
,
n
2
)
=
(
50,10
)
RB(R~)
-
1.036
-
1.438
-
1.704
-
1.622
-
1.841
RB(R^)
-
0.106
-
0.184
-
0.316
-
0.217
-
0.528
RMSE(R~)
0.050
0.076
0.093
0.101
0.106
RMSE(R^)
0.048
0.075
0.093
0.102
0.108
(
n
1
,
n
2
)
=
(
50,20
)
RB(R~)
-
0.403
-
0.508
-
0.450
-
0.060
-
0.330
RB(R^)
0.084
0.134
0.242
0.013
0.254
RMSE(R~)
0.034
0.054
0.066
0.066
0.077
RMSE(R^)
0.034
0.054
0.067
0.067
0.078
(
n
1
,
n
2
)
=
(
50,50
)
RB(R~)
-
0.284
-
0.412
-
0.480
-
0.448
-
0.427
RB(R^)
-
0.059
-
0.138
-
0.218
-
0.238
-
0.304
RMSE(R~)
0.024
0.038
0.047
0.053
0.054
RMSE(R^)
0.024
0.038
0.048
0.053
0.055Table 5
Simulation results: bias and root mean square error of variance estimators.
R
=
0.9
R
=
0.8
R
=
0.7
R
=
0.6
R
=
0.5
(
n
1
,
n
2
)
=
(
10,10
)
RB(v(R~))
0.22
−2.64
−2.37
−3.28
−5.23
RB(v(R^))
1.02
−2.22
−2.09
−3.17
−5.15
RMSE(v(R~))·104
25.9
36.1
34.0
25.3
19.8
RMSE(v(R^))·104
27.5
40.4
39.1
29.5
22.9
(
n
1
,
n
2
)
=
(
10,20
)
RB(v(R~))
3.99
−0.07
−3.34
−2.72
−4.11
RB(v(R^))
4.67
0.29
−3.17
−2.63
−4.08
RMSE(v(R~))·104
15.5
22.7
21.9
15.8
10.7
RMSE(v(R^))·104
16.1
24.7
24.1
17.6
11.9
(
n
1
,
n
2
)
=
(
10,50
)
RB(v(R~))
6.22
2.33
0.70
0.88
−0.55
RB(v(R^))
6.66
2.56
0.82
0.93
−0.54
RMSE(v(R~))·104
10.2
15.6
15.3
11.5
7.20
RMSE(v(R^))·104
10.5
16.6
16.4
12.4
7.71
(
n
1
,
n
2
)
=
(
20,10
)
RB(v(R~))
−2.80
−3.56
−2.55
−3.67
−5.45
RB(v(R^))
−2.58
−3.43
−2.47
−3.57
−5.37
RMSE(v(R~))·104
18.4
26.3
24.8
18.3
14.6
RMSE(v(R^))·104
19.1
28.5
27.6
20.7
16.1
(
n
1
,
n
2
)
=
(
20,20
)
RB(v(R~))
4.76
2.51
−0.60
1.38
−0.72
RB(v(R^))
5.09
2.67
−0.54
1.40
−0.71
RMSE(v(R~))·104
8.84
13.5
12.9
8.84
5.03
RMSE(v(R^))·104
9.04
14.3
13.9
9.60
5.45
(
n
1
,
n
2
)
=
(
20,50
)
RB(v(R~))
7.39
4.45
2.23
2.81
2.20
RB(v(R^))
7.60
4.55
2.28
2.84
2.21
RMSE(v(R~))·104
4.67
7.41
7.24
5.13
2.68
RMSE(v(R^))·104
4.74
7.69
7.58
5.39
2.81
(
n
1
,
n
2
)
=
(
50,10
)
RB(v(R~))
−3.24
−4.60
−5.48
−3.24
−5.42
RB(v(R^))
−3.08
−4.59
−5.52
−3.27
−5.45
RMSE(v(R~))·104
14.3
21.2
20.8
15.0
11.5
RMSE(v(R^))·104
14.5
22.5
22.6
16.6
12.4
(
n
1
,
n
2
)
=
(
50,20
)
RB(v(R~))
1.79
−0.23
1.15
2.81
0.47
RB(v(R^))
1.93
−0.16
1.18
2.84
0.48
RMSE(v(R~))·104
5.64
9.00
8.78
5.13
3.21
RMSE(v(R^))·104
5.68
9.34
9.25
5.39
3.40
(
n
1
,
n
2
)
=
(
50,50
)
RB(v(R~))
8.02
7.02
6.14
5.35
5.68
RB(v(R^))
8.08
7.05
6.15
5.36
5.69
RMSE(v(R~))·104
2.19
3.62
3.63
2.64
1.88
RMSE(v(R^))·104
2.21
3.69
3.73
2.72
1.92Table 6
Simulation results: coverage rate and average length of confidence intervals based on ML and UMVU estimators.
R
=
0.9
R
=
0.8
R
=
0.7
R
=
0.6
R
=
0.5
(
n
1
,
n
2
)
=
(
10,10
)
cov(R~)
0.905
0.916
0.917
0.918
0.916
cov(R^)
0.870
0.893
0.913
0.916
0.914
length(R~)
0.220
0.338
0.413
0.456
0.471
length(R^)
0.214
0.342
0.425
0.472
0.490
(
n
1
,
n
2
)
=
(
10,20
)
cov(R~)
0.912
0.922
0.922
0.924
0.925
cov(R^)
0.893
0.914
0.915
0.925
0.924
length(R~)
0.184
0.284
0.348
0.383
0.394
length(R^)
0.181
0.287
0.355
0.394
0.406
(
n
1
,
n
2
)
=
(
10,50
)
cov(R~)
0.911
0.919
0.924
0.926
0.927
cov(R^)
0.899
0.918
0.921
0.926
0.926
length(R~)
0.158
0.246
0.300
0.330
0.337
length(R^)
0.156
0.248
0.305
0.337
0.345
(
n
1
,
n
2
)
=
(
20,10
)
cov(R~)
0.925
0.930
0.928
0.929
0.922
cov(R^)
0.894
0.914
0.924
0.926
0.922
length(R~)
0.197
0.305
0.374
0.415
0.431
length(R^)
0.191
0.306
0.381
0.426
0.443
(
n
1
,
n
2
)
=
(
20,20
)
cov(R~)
0.927
0.936
0.935
0.938
0.933
cov(R^)
0.913
0.927
0.927
0.937
0.933
length(R~)
0.154
0.241
0.297
0.330
0.342
length(R^)
0.151
0.242
0.301
0.336
0.348
(
n
1
,
n
2
)
=
(
20,50
)
cov(R~)
0.937
0.941
0.942
0.944
0.938
cov(R^)
0.925
0.933
0.940
0.944
0.938
length(R~)
0.123
0.194
0.238
0.263
0.271
length(R^)
0.122
0.195
0.241
0.267
0.275
(
n
1
,
n
2
)
=
(
50,10
)
cov(R~)
0.926
0.928
0.927
0.927
0.921
cov(R^)
0.892
0.916
0.917
0.928
0.924
length(R~)
0.180
0.281
0.347
0.387
0.403
length(R^)
0.174
0.280
0.351
0.394
0.412
(
n
1
,
n
2
)
=
(
50,20
)
cov(R~)
0.935
0.935
0.943
0.944
0.946
cov(R^)
0.917
0.932
0.938
0.944
0.945
length(R~)
0.132
0.209
0.260
0.263
0.304
length(R^)
0.130
0.209
0.262
0.267
0.308
(
n
1
,
n
2
)
=
(
50,50
)
cov(R~)
0.946
0.953
0.948
0.954
0.948
cov(R^)
0.941
0.949
0.947
0.955
0.953
length(R~)
0.097
0.153
0.190
0.212
0.219
length(R^)
0.096
0.154
0.191
0.213
0.221The simulation results show that the ML estimatorR~ always presents a very small bias even for small samples: in absolute value, the MC percentage relative bias is always smaller than 1.841% for all the scenarios considered (whereas the maximum absolute MC percentage relative bias for R^, which is theoretically unbiased, is 0.552%). In 42 scenarios out of 45, R~ underestimates R. Regarding the RMSE, the ML estimator performs better than UMVU in 27 cases out of 45, worse in 7 cases, and in 11 cases the RMSE is equal at the third decimal digit. However, under each scenario, even for smaller sample sizes, the values of RMSE for the ML and UMVUE estimators are very close. The ML outperforms the UMVU estimator as the value of R gets close to 0.5; their performances tend to be alike as n1 and n2 increase. For both estimators, for fixed sample sizes, the RMSE increases as R decreases; for a fixed R, the RMSE increases, as the sample sizes decrease (as expected). Figure 2 displays the MC distribution of the ML and UMVU estimators in the case R=0.7, for three values of sample size; it highlights their very similar behaviour.MC distribution ofR~ (MLE) and R^ (UMVUE) when λ1=1, λ2=2.497 (R=0.7).
(a)
(
n
1
=
n
2
=
10)
(b)
(n1=n2=20)
(c)
(n1=n2=50).Regarding the approximate variance estimators, surprisingly their performance is good even for the moderate sample sizes considered in this study; the percentage relative bias, in absolute value, is never greater than 8%: the worst performance occurs forn1=n2=50 and R=0.9. Indeed, when both sample sizes equal 50, the RB is greater than for small sample sizes, whereas one would expect that the RB decreases in absolute value when sample sizes increase. The results of further simulations not reported here show that the RB actually decreases to zero for n1=n2=100. For both estimators, the rate of underestimates is almost equal to the rate of overestimates. Under each scenario, and especially when R=0.5, the value of RB of the variance estimator v(R^) is quite close to the corresponding value of the RB of v(R~), whereas the RMSE of v(R~) is smaller than the RMSE of v(R^) in each of the 45 cases considered. The RB of the two approximate variance estimators does not present a clear trend in terms of R; while their RMSEs, for each of the couples (n1,n2) here explored, seem to present a maximum near R=0.8 and a minimum for R=0.5.The confidence intervals built upon the point estimators and these variance estimators present coverage that is always greater than 87% for UMVU and 90.5% for ML: the lowest value is obtained forn1=n2=10 and R=0.9. They attain the nominal level (95%) for n1=n2=50; in 31 and 21 cases out of 45, respectively, the coverage rate of the ML and UMVU interval estimators is greater than or equal to 92.5%. Overall, the CIs present better coverage when R is close to 0.5. In fact, in this case, the distributions of R^ and R~ tend to be symmetrical and are more finely approximated by the normal distribution; then, the confidence intervals (9), which assume an underlying normal distribution, show a better performance. The CIs based on the ML estimator almost always show a coverage rate greater than those based on the UMVU estimator; moreover, the latter are always a bit wider, unless when R=0.9. This feature tends to be negligible when the sample sizes are increased. As one would expect, the average length decreases as sample sizes increase, for fixed R, and as R increases, for fixed sample sizes.The results forλ1>1, which are not reported here for the sake of brevity, confirm the previous findings. Even if the study is obviously not exhaustive, since only several scenarios have been covered, nevertheless these general features can be outlined.


## 5. An Example of Application

In this section, we apply the inferential techniques presented in Section3 to a real dataset. The application is based on the data from an engineering experiment discussed in [3], carried out in an electric company, under several experimental conditions (called “runs”), corresponding to different combinations of 8 factors. The blackening experiment was conducted in a three-layer oven; when each run was completed, 30 masks from each layer in the oven were collected to examine the number of defects in each mask. The total number of defects in the 30 masks from the upper layer for each experimental run is observed (see Table 7). We focus on runs 1 and 2, where 317 and 184 defects, respectively, are counted.Table 7
Data for the blackening experiment in [3].
Run
Factor
Number of defects
A
B
C
D
E
F
G
H
1
1
1
1
1
1
1
1
1
317
2
1
1
2
2
2
2
2
2
184
3
1
1
3
3
3
3
3
3
528
4
1
2
1
2
1
2
3
3
163
5
1
2
2
3
2
3
1
1
96
6
1
2
3
1
3
1
2
2
300
7
1
3
1
3
2
1
2
3
177
8
1
3
2
1
3
2
3
1
182
9
1
3
3
2
1
3
1
2
75
10
2
1
1
3
3
2
1
2
146
11
2
1
2
1
1
3
2
3
135
12
2
1
3
2
2
1
3
1
232
13
2
2
1
2
3
3
2
1
543
14
2
2
2
3
1
1
3
2
101
15
2
2
3
1
2
2
1
3
282
16
2
3
1
1
2
3
3
2
90
17
2
3
2
2
3
1
1
3
288
18
2
3
3
3
1
2
2
1
554Since the number of defects in a mask is either zero or a positive integer, the appropriate distribution is the Poisson. Denoting withY and X the variables modeling this number for run 1 and run 2, respectively, we are interested in determining a point estimate and an interval estimator for the probability that the number of defects in run 1 is smaller than in run 2, that is, P(X<Y). Since the sample size is 30 for both variables, then x-=6.13 and y-=10.57; then, supposing that X and Y follow a Poisson distribution, the ML and UMVU estimators and their corresponding approximate variances can be computed according to (2), (3), and (8); the associated confidence intervals can be estimated recalling (10).The results are presented in Table8 and show the closeness between the two approaches. All the confidence intervals for P(X<Y) always exclude 0.5, thus meaning that the difference in sample means testifies to the statistical dominance of Y on X: the number of defects under run 1 is stochastically larger than the number of defects under run 2.Table 8
Results for the application: ML and UMVUE point estimates and asymptotic confidence intervals(L,U) for R.
Point est.
90% CI
95% CI
99% CI
L
U
L
U
L
U
ML
0.834
0.759
0.909
0.745
0.923
0.717
0.951
UMVU
0.838
0.763
0.913
0.749
0.927
0.721
0.955


## 6. Conclusions

In this paper, point and interval estimators for the reliability of a Poisson stress-strength model are presented, discussed, and empirically compared through a Monte Carlo simulation study. The results show that the maximum likelihood and uniformly minimum variance unbiased estimators possess similar sampling properties, and the first is slightly preferable to the second in terms of dispersion around the true value of reliability. Moreover, although the variance estimators proposed here are approximate (i.e., biased), being based on the delta method for asymptotically normal r.v., the empirical results emphasize that the estimators' bias is small even for moderate sample sizes, and these estimators can be usefully employed to build approximate confidence intervals, whose coverage is shown to be overall close to the fixed nominal level, especially when the value of reliability is close to0.5. However, when R is close to 1 (or, symmetrically, 0), the intervals can show a poorer performance; then caution is needed when constructing a confidence interval based on a point estimate close to 1 (0). In this case, one can resort, for example, to some variance-stabilizing transformation of the estimate.



---

*Source: 530530-2013-06-05.xml*