Comput Stat (2014) 29:1083–1094 DOI 10.1007/s00180-014-0481-6 ORIGINAL PAPER
Comparing exponential location parameters with several controls under heteroscedasticity A. Malekzadeh · M. Kharrati-Kopaei · S. M. Sadooghi-Alvandi
Received: 9 September 2013 / Accepted: 15 January 2014 / Published online: 2 February 2014 © Springer-Verlag Berlin Heidelberg 2014
Abstract Suppose that random samples are taken from k treatment groups and l control groups, where the observations in each group have a two-parameter exponential distribution. We consider the problem of constructing simultaneous confidence intervals for the differences between location parameters of the treatment groups and the control groups when the scale parameters may be unequal. Using the parametric bootstrap approach, we develop a new multiple comparisons procedure when the scale parameters and sample sizes are possibly unequal. We then present a simulation study in which we compare the performance of our proposed procedure with two other procedures. The results of our simulations indicate that our proposed procedure performs better than other procedures. The usefulness of our proposed procedure is illustrated with an example. Keywords Multiple comparisons · Simultaneous confidence intervals · Parametric bootstrap · Coverage probability · Monte Carlo method · Simulation Mathmetics Subject Classification
62F25 · 62F40
A. Malekzadeh · M. Kharrati-Kopaei (B) · S. M. Sadooghi-Alvandi Department of Statistics, Shiraz University, Shiraz, Iran e-mail:
[email protected];
[email protected] A. Malekzadeh e-mail:
[email protected];
[email protected] S. M. Sadooghi-Alvandi e-mail:
[email protected];
[email protected]
123
1084
A. Malekzadeh et al.
1 Introduction Suppose that random samples are taken from k treatment groups and l control groups, where the observations in each group have a two-parameter exponential distribution. Multiple comparisons of several exponential distributions have applications in engineering, life testing, and medical sciences. In reliability analysis, the location parameter is known as the guaranteed mean life and the scale parameter is known as the mean life. In dose-response analysis, the location and scale parameters are known as the guaranteed and mean effective duration, respectively. For more detailed applications of the exponential distribution, see Johnson and Kotz (1970), Zelen (1966), Lawless and Singhal (1980), and Maurya et al. (2011a); Maurya et al. (2011b). We use the notation Exp (ξ, σ ) to denote the exponential distribution with location parameter ξ ∈ (−∞, +∞), scale parameter σ ∈ (0, +∞), and density function x −ξ 1 I[ξ,∞) (x) , f (x; ξ, σ ) = exp − σ σ where I A (.) is the indicator function of event A. Suppose that the observations in the ith treatment group have a two-parameter exponential distribution Exp ξiT , σiT and the observations in the jth control group have a two-parameter exponential distribution C C Exp ξ j , σ j . We consider the problem of constructing simultaneous confidence intervals (SCIs) for ξiT − ξ Cj , i = 1, . . ., k, j = 1, . . ., l, under heteroscedasticity, i.e., when the scale parameters may be unequal. Singh and Abebe (2009) proposed a multiple comparison procedure when the sample sizes and scale parameters are equal; see also Ng et al. (1993). A procedure for multiple comparisons with a control under heteroscedasticity was first proposed by Lam and Ng (1990). However, their method is based on a two-stage sampling procedure. For single-stage sampling, Wu et al. (2010) proposed multiple comparisons with a control under heteroscedasticity when the sample sizes are equal. Maurya et al. (2011a) extended these procedures to comparisons with more than one control. The procedures proposed by these authors are based on Lam’s (1987, 1988) technique. However, using this technique leads to SCIs that are too conservative; see Wu et al. (2010) and Maurya et al. (2011a). More recently, Karrati-Kopaei et al. (2013) proposed simultaneous fiducial generalized confidence intervals (SFGCIs) for the successive differences of the location parameters of several exponential distributions and indicated how their proposed procedure can be modified to obtain SFGCIs for comparisons of the location parameters with one or several controls. In this paper, we use the parametric bootstrap approach to obtain a new multiple comparisons procedure when the scale parameters and sample sizes are possibly unequal. The parametric bootstrap approach has been successfully used in other contexts, to solve inference problems (mainly testing problems) involving nuisance parameters; see, e.g. Krishnamoorthy et al. (2007), Ma and Tian (2009), Tian et al. (2009), Krishnamoorthy and Lu (2010), Li et al. (2011), Xu et al. (2012), and SadooghiAlvandi and Jafari (2013).
123
Comparing exponential location parameters
1085
The organization of this paper is as follows: In Sect. 2, we present parametric bootstrap SCI (PBSCIs) when the scale parameters and sample sizes are possibly unequal. In Sect. 3, we present a simulation study in which we compare the coverage probability of our proposed procedure with the procedures proposed by Maurya et al. (2011a) and Karrati-Kopaei et al. (2013). The results of our simulation study indicate that our proposed PBSCIs have coverage probabilities close to the nominal confidence coefficient and, consequently, the resulting confidence intervals are typically shorter. In Sect. 4, we illustrate our proposed procedure with an example. 2 Parametric bootstrap simultaneous confidence intervals sample of sizes m i from the ith treatLet X i = X i1 , . . . , Xim i denote a random ment group and Y j = Y j1 , . . . , Y jn j denote a random sample of sizes n j from the jth control group (i = 1, . . ., k, j = 1, . . ., l), and let ξ T = ξ1T , . . . , ξkT , ξ C = C ξ1 , . . . , ξlC , σ T = σ1T , . . . , σkT , and σ C = σ1C , . . . , σlC denote the vectors of the location and scale parameters of the treatment and control groups. Let X i(1) and Y j (1) denote the first order statistics of X i and Y j , respectively, and let
SX i
m
nj
t=1
t=1
i 1 1 = (X it − X i(1) ) and SY j = (Y jt − Y j (1) ) mi − 1 nj − 1
denote the unbiased estimators of σiT and σ jC , respectively. It is well-known that X i(1) ∼ Exp ξiT , σiT /m i and Y j(1) ∼ Exp ξ Cj , σ jC /n j ,
(2.1)
independently of SX i ∼
σiT χ2 (2m i − 2) 2m i −2
and
SY j ∼
σ jC 2n j − 2
2 χ2n , j −2
(2.2)
where χr2 denotes a chi-square distribution with r degrees of freedom. Preliminaries. Suppose that the values of σiT ’s and σ jC ’s were known. Let X C T i(1) − Y j (1) − ξi − ξ j , T∗ = max √ 1≤i≤k,1≤ j≤l γi j where
γi j = Var X i(1) − Y j (1) =
σiT mi
2 +
σ jC nj
(2.3)
2 .
(2.4)
123
1086
A. Malekzadeh et al.
Then from (2.1), it follows T ∗ is a pivotal quantity, and 100(1-α) % two-sided SCIs for ξiT − ξ Cj are √ X i(1) − Y j (1) ± qα∗ γi j , i = 1, . . . , k;
j = 1, . . . , l,
(2.5)
where qα∗ denotes the (1 − α)th quantile of the distribution of T ∗ . (The approximate value of qα∗ may be obtained by Monte Carlo methods.) For the realistic case when the values of σiT ’s and σ jC ’s are unknown, we propose to replace their values with their estimates and base SCIs on the random quantity X C T i(1) − Y j (1) − ξi − ξ j , T = max 1≤i≤k,1≤ j≤l Vi j where Vi j =
(m i − 1)S X2 i m i3
+
(n j − 1)SY2 j
(2.6)
(2.7)
n 3j
is the minimum variance unbiased estimator of γi j . We propose to approximate the quantiles of T using a parametric bootstrap (PB) approach. As noted in the Introduction, the PB approach has been successfully used in various other contexts. Parametric bootstrap SCI (PBSCI). The PB distribution of T is obtained by finding its sampling distribution when the values of the nuisance parameters σiT and σ jC are replaced with the observed values of their estimates sxi and syj . It follows from (2.1) and (2.2) that
X i(1) − ξiT ∼ Exp 0, σiT /m i , Y j(1) − ξ Cj ∼ Exp 0, σ jC /n j ,
and that the PB pivotal variable corresponding to T is TB =
max
1≤i≤k,1≤ j≤l
B Vi Bj , X i(1) − Y jB(1)
(2.8)
where
Vi Bj =
2 (m i − 1) S XBi m i3
+
2 n j − 1 SYBj n 3j
,
with B ∼ sxi /m i Exp (0, 1) and Y jB(1) ∼ s y j /n j Exp(0, 1), X i(1) sy j sxi 2 χ2m χ2 SXBi ∼ and SYBj ∼ . i −2 (2n j − 2) 2n j −2 (2m i − 2)
123
(2.9)
Comparing exponential location parameters
1087
Then the distribution of T B provides the PB approximation to the distribution of T and 100(1 − α) % two-sided PBSCIs for ξiT − ξ Cj are X i(1) − Y j (1) ± qαB Vi j , i = 1, . . . , k;
j = 1, . . . , l
(2.10)
where qαB denotes the (1 − α)th quantile of the distribution of T B . Computational procedure. The value of qαB can be estimated by the following Monte Carlo procedure: (i) Calculate sxi and syj (the observed values of SX i and SY j ) from the B , Y B , S B , S B , using (2.9) [by generating data. (ii) Generate a realization of X i(1) Xi Yj j (1) values from the standard exponential distributions Exp(0, 1)and chi-square distribu2 2 and χ2n ], and calculate the value of T B , as given by (2.8). (iii) Repeat tions χ2m i −2 j −2 step (ii) a large number of times, N , to obtain N realizations of T B . From these N values of T B , obtain the empirical distribution of T B and calculate its empirical (1 − α)th quantile as an estimate of qαB . Remark The PB method can also be used for constructing SCIs for all pair-wise comparisons of the location parameters of several exponential distributions under heteroscedasticity (details are omitted). The following theorem shows that the PBSCIs in (2.10) have correct coverage probability asymptotically. l k Theorem 1 Let n = i=1 m i + j=1 n j and suppose that n → ∞ such that m i /n → u i ∈ (0, 1) and n j /n → v j ∈ (0, 1), for i = 1, . . ., k, j = 1, . . ., l. Then Pr ξiT −ξ Cj ∈ X i(1) −Y j (1) ± qαB Vi j , all i=1, . . . , k, j=1, . . . , l → 1 − α. Proof First, note that Pr
ξiT − ξ Cj ∈ X i(1) − Y j (1) ± qαB Vi j , all i = 1, . . . , k; j = 1, . . . , l ⎧ ⎫ X C T ⎨ ⎬ i(1) − Y j (1) − ξi − ξ j ≤ q B = Pr T ≤ q B . = Pr max α α ⎩1≤i≤k,1≤ j≤l ⎭ Vi j
Now, let X (1) = X 1(1) , . . . , X k(1) and Y(1) = (Y1(1) , . . . , Yl(1) ). Then it is easily shown that X (1) − ξ T D W ∗ , as n → ∞ → n Y (1) − ξ C Z∗ D
→ denotes convergence in distribution and the components of W ∗ = ∗ where ∗ ∗ ∗ ∗ exponential variW1 , . . . , Wk and Z = Z 1 , . . . , Z l are mutually independent C ∗ T ∗ ables, with Wi ∼ Exp 0, σi /u i and Z j ∼ Exp 0, σ j /v j . But it is easily ver-
123
1088
A. Malekzadeh et al.
2 ified, using (2.2), that E(SX i ) = σiT and Var(SX i ) = 2 σiT /(2m i − 2). Therefore, using Chebyshev’s inequality, it follows that S X i converges in probability to σiT (i = 1, . . ., k). Similarly, SY j converges in probability to σ jC ( j = 1, . . . , l). Hence, D
using Slutsky’s theorem (see Ferguson 1996, p. 6), it follows that T → T∞ , as n → ∞, where Wi∗ − Z ∗j T∞ = max 2 . 1≤i≤k,1≤ j≤l T 2 C σi /u i + σ j /v j D
Now we show that T B → T∞ , as n → ∞. First note that SXi P χ2 → σiT (2m i − 2) 2m i −2
and
SY j (2n j − 2)
P
2 χ2n → σ jC , as n → ∞, j −2
P
where → denotes convergence in probability. Hence n
2
Vi Bj
n 2 (m i − 1) B 2 n 2 n j − 1 B 2 SX i + SY j = m i3 n 3j 2 2 P → σiT /u i + σ jC /v j , as n → ∞.
(indepenNow let W1 , . . ., Wk and Z 1 , . . ., Z l be independent exponential variables dent of X i j and Yi j ), with Wi ∼ Exp 0, σiT /u i and Z j ∼ Exp 0, σ jC /v j . Note that SY j SXi B T B X i(1) = u i /σi Wi and Y j(1) = v j /σ jC Z j . mi nj Then B B Vi Bj X i(1) − Y j(1) 1≤i≤k,1≤ j≤l B = max − Y jB(1) n 2 Vi Bj n X i(1) 1≤i≤k,1≤ j≤l
nu S
Sy j nv j i xi B 2 V = max − n W Z i j ij 1≤i≤k,1≤ j≤l mi nj σiT σ jC 2 T 2 C D σi u i + σ j /v j → max , as n → ∞ Wi − Z j 1≤i≤k,1≤ j≤l
TB =
max
D
Thus T B → T∞ as n → ∞. Note that since the limiting distribution of T B is continuous, we have qαB → qα , as n → ∞, where qc is the (1 − α)th quantile of the
123
Comparing exponential location parameters
1089
! distribution of T∞ . Therefore, as n → ∞, Pr T ≤ qαB → Pr {T∞ ≤ qα } = 1 − α. This completes the proof.
3 Simulation results In this section, we present the results of our simulation studies in which we examine the coverage probabilities of our proposed parametric bootstrap procedure and compare its performance with the method proposed by Maurya et al. (2011a) and also with the simultaneous fiducial generalized confidence intervals proposed by Karrati-Kopaei et al. (2013). We use the abbreviations PBSCI, M-SCI, and SFGCI to refer to these three methods. We note that Maurya et al. (2011a) compared M-SCI with the Bonferroni SCIs and concluded that M-SCI method performs better than the Bonferroni method. Following Maurya et al. (2011a), in our simulations we considered the case of k = 3 treatments and l = 2 controls. Note that it is easily verified that the coverage probabilities of these methods do not depend on the values of the location parameters ξ T and ξ C . Therefore, without loss of generality, in our simulations we assumed that ξ T = 0 and ξ C = 0. We first briefly review the methods SFGCI and M-SCI and then present our simulation results. SFGCI. This procedure is based on the concept of fiducial generalized pivotal quantities (FGPQ); see Hannig et al. (2006a, b) and Chang and Huang (2009). Let X i∗ and Y ∗j denote independent copies of X i and Y j , and let Ri j = X i(1) −
SY j ∗ SX i ∗ X i(1) − ξiT − Y j(1) + Y j(1) − ξ Cj . SX ∗i SY ∗j
Then, it is easily verified that Ri j are simultaneous fiducial generalized pivotal quantities for ξiT − ξ Cj , for i = 1, . . ., k and j = 1, . . ., l. Now, let
T
F
X i(1) − Y j (1) − Ri j = max , 1≤i≤k;1≤ j≤l Vi j
where Vi j ’s are the same as in (2.7). Then 100(1−α) % two-sided SFGCIs for ξiT −ξ Cj are X i(1) − Y j (1) ± qαF Vi j , i = 1, . . . , k, j = 1, . . . , l, where qαF is the(1-α)th quantile of the conditional distribution of T F given the observed values of X i ’s and Y j ’s. (The conditional distribution of T F does not depend on any of the parameters and depends on the observed values of X i ’s and Y j ’s only through sxi and syj .) The value of qαF is easily obtained by a Monte Carlo procedure similar to the computational procedure given in Sect. 2. For details of the SFGCI procedure, see Karrati-Kopaei et al. (2013).
123
1090
A. Malekzadeh et al.
Table 1 The coverage probabilities (CP) and average volumes (AV) of 95 % SCIs, when σ T = (1, 1, 1) and σ C = (1, 1)
Method
Sample sizes 10
20
30
CP
AV
CP
AV
CP
AV
PBSCI
0.9625
2.177
0.9596
0.0142
0.9558
0.00092
SFGCI
0.9790
4.764
0.9694
0.0208
0.9628
0.00118
M-SCI
0.9978
42.620
0.9937
0.1181
0.9897
0.00549
M-SCI. In this procedure, it is assumed that all the sample sizes are equal. Let m denote the common sample size; i.e., in our notation, m i = n j = m, for i = 1, . . ., k and j = 1, . . ., l. The two-sided conservative SCIs for ξiT − ξ Cj , i = 1, . . ., k, j = 1, . . ., l, are −1 (1 − α)1/(k+l) , X i(1) − Y j (1) ± d F2,2m−2 where d = max
max S X i /m , max SY j /m
1≤i≤k
1≤ j≤l
−1 and F2,2m−2 (x) is the xth quantile
of an F distribution with 2 and 2m − 2 degrees of freedom. For details, see Maurya et al. (2011a). Simulation procedure. We used the following procedure to estimate the true coverage probabilities: I Generate X i and Y j , random samples of sizes m i and n j from exponential distribu T T C C tions Exp ξi , σi and Exp ξ j , σ j (i = 1, . . ., k, j = 1, . . ., l) and calculate the observed values of X i(1) , SX i , Y j (1) , and SY j . II Calculate the value of qαB using the computational procedure given in Sect. 2, with N = 1,000,000. Also, calculate qαF using a similar computational procedure; see Karrati-Kopaei et al. (2013). III For each method, construct two-sided SCIs and record whether or not all the values of ξiT − ξ Cj , i = 1, . . ., k, j = 1, . . ., l, fall in their corresponding SCIs. IV Repeat steps I–III, a large number of times, M = 10,000. Then, for each method, the fraction of times that all ξiT − ξ Cj (i = 1, . . ., k, j = 1, . . ., l) are in their corresponding SCIs provides an estimate of the true coverage probability. Simulation results I. We first considered the two cases considered by Maurya et al. (2011a), with k = 3, l = 2, and sample sizes equal to m, with m = 10, 20 and 30 (as small, moderate, and large sample sizes). As noted before, without loss of generality, we assumed that ξ T = 0 and ξ C = 0. We also considered two cases for the scale parameters: (i) σ T = (1, 1, 1), σ C = (1, 1) and (ii) σ T = (1, 1.1, 1.3), σ C = (1.2, 1.4). The results of simulations for these two cases with N = 1, 000, 000, M = 10, 000, and 1 − α = 0.95 are shown in Tables 1 and 2, respectively. It is seen from Tables 1 and 2 that, compared to the M-SCI method, the coverage probabilities of our proposed PBSCI method are closer to the nominal confidence
123
Comparing exponential location parameters Table 2 The coverage probabilities (CP) and average volumes (AV) of 95 % SCIs, when σ T = (1, 1.1, 1.3) and σ C = (1.2, 1.4)
Table 3 The coverage probabilities (CP) and average volumes (AV) of 95 % SCIs, when σ T = (1, 1.1, 1.3) and σ C = (0.6, 1.2)
Method
1091 Sample sizes 10 CP
20 AV
30
CP
AV
CP
AV
PBSCI
0.9650
7.413
0.9587
0.0487
0.9549
0.0031
SFGCI
0.9784
24.412
0.9696
0.1096
0.9648
0.0063
M-SCI
0.9976
164.553
0.9934
0.4675
0.9899
0.0229
Method
Sample sizes 10 CP
20 AV
CP
30 AV
CP
AV
PBSCI
0.9548
4.165
0.9535
0.0298
0.9525
0.0020
SFGCI
0.9726
8.301
0.9648
0.0409
0.9567
0.0025
M-SCI
0.9972
99.740
0.9928
0.2943
0.9904
0.0142
coefficient and, consequently, the resulting confidence intervals are typically shorter (as reflected in the average volume). The SFGCI method also performs better than the M-SCI method, but not as well as the PBSCI method. Simulation results II. In the cases considered by Maurya et al. (2011a), the scale parameters are not much different. But if the scale parameters are more unequal, then the differences between the methods become more apparent, as can be seen from Table 3. Simulation results III. In M-SCI method, it is assumed that the sample sizes are equal to m (i.e., m i = n j = m, for i = 1, . . ., k and j = 1, . . ., l); but note that the PBSCI and SFGCI methods can be used with unequal sample sizes. In view of this, we decided to compare the PBSCI and SFGCI methods when the sample sizes are unequal. As before, we took k = 3 and l = 2. The results of simulation for different values of (m 1 , m 2 , m 3 ) and (n 1 , n 2 ), when σ T = (1, 1.1, 1.3) and σ C = (0.6, 1.2), are shown in Table 4. It is seen from Table 4 that, compared to the SFGCI method, the coverage probabilities of our proposed PBSCI method are closer to the nominal confidence level.
4 Illustrative example The data given in Table 5, taken from Maurya et al. (2011a), concern the effectiveness of four drugs (two test drugs and two control drugs) used in the treatment of Leukemia, as measured by the duration of remission time.
123
1092 Table 4 The coverage probabilities (CP) and average volumes (AV) of 95 % SCIs for the PBSCI and SFGCI methods, when σ T = (1, 1.1, 1.3) and σ C = (0.6, 1.2) and the sample sizes are unequal
Table 5 Remission duration by four drugs
A. Malekzadeh et al. Sample sizes
Method PBSCI
SFGCI
(m 1 , m 2 , m 3 )
(n 1 , n 2 )
CP
AV
CP
AV
(10, 10, 10)
(10, 30)
0.9498
2.1984
0.9658
3.5154
(10, 10, 10)
(30, 10)
0.9472
7.8477
0.9573
10.7208
(10, 10, 10)
(30, 30)
0.9464
2.6389
0.9559
3.2741
(30, 10, 10)
(10, 10)
0.9477
2.6611
0.9659
4.8440
(30, 30, 10)
(10, 10)
0.9508
0.9515
0.9653
1.5565
(30, 30, 30)
(10, 10)
0.9486
0.2654
0.9594
0.3821
(10, 10, 30)
(10, 10)
0.9557
1.5337
0.9713
2.9800
(10, 30, 30)
(10, 10)
0.9539
0.6395
0.9670
1.0936
(30, 10, 10)
(30, 10)
0.9431
2.0294
0.9554
2.6781
(30, 10, 10)
(10, 30)
0.9452
0.5022
0.9614
0.8369
(10, 10, 30)
(10, 30)
0.9495
0.3180
0.9645
0.5547
(10, 10, 30)
(30, 10)
0.9504
1.6573
0.9617
2.2796
Test drug 1
Test drug 2
Control drug 1
Control drug 2
1.034
5.115
2.214
4.158
2.344
4.498
4.976
4.025
1.266
4.617
8.154
5.170
1.563
4.651
2.686
11.909
1.169
4.533
2.271
4.912
4.118
4.513
3.139
4.629
1.013
7.641
2.214
3.955
1.509
5.971
4.480
6.735
1.109
12.130
8.847
3.140
1.965
4.699
2.239
12.446
5.136
4.914
3.473
8.777
1.533
17.169
2.761
6.321
1.716
5.497
2.833
3.256 8.250
2.778
11.332
2.381
2.546
18.922
3.548
3.759
2.626
13.712
2.414
5.205
3.413
6.309
2.832
3.071
1.929
10.086
5.551
3.147
2.061
9.293
3.376
9.773
2.951
11.787
2.968
10.218
To obtain the value of qαB and qαF , we used the procedure given in Sect. 2 with N = 1,000,000 runs. The estimated values of qαB and qαF for α = 0.05, 0.025, and 0.01 are shown in Table 6.
123
Comparing exponential location parameters Table 6 The estimated values of qαB and qαF for two-sided PBSCIs and SFGCIs
Table 7 Two-sided PBSCIs, SFGCIs, and M-SCIs for comparing the tests drugs with the control drugs
1093
α
qαB
qαF
0.050
3.65
3.86
0.025
4.40
4.66
0.010
5.39
5.76
Parameters α = 0.050
α = 0.025
α = 0.010
PBSCIs ξ1T − ξ1C (−1.552, −0.851) (−1.623, −0.779) (−1.718, −0.684) (1.352, 3.216) (1.141, 3.427) ξ2T − ξ1C (1.509, 3.059) ξ1T − ξ2C (−2.675, −1.441) (−2.799, −1.317) (−2.967 −1.149) ξ2T − ξ2C (0.501, 2.354) Volume 2.48
(0.313, 2.541)
(0.061, 2.793)
5.20
11.74
SFGCIs ξ1T − ξ1C (−1.571, −0.831) (−1.648, −0.754) (−1.754, −0.648) ξ2T − ξ1C (1.465, 3.103) (1.296, 3.272) (1.062, 3.506) ξ1T − ξ2C (−2.709,−1.407) (−2.844, −1.272) (−3.030, −1.086) ξ2T − ξ2C (0.448, 2.406)
(0.246, 2.608)
(−0.034, 2.888)
Volume
6.56
15.35
3.09
M-SCIs ξ1T − ξ1C (−2.200, −0.202) (−2.384, −0.018) (−2.635, 0.233) ξ2T − ξ1C (1.285, 3.283) (1.101, 3.467) (0.850, 3.718) ξ1T − ξ2C (−3.057, −1.059) (−3.241, −0.875) (−3.492, −0.624) ξ2T − ξ2C (0.428, 2.426) (0.244, 2.610) (−0.007, 2.861) Volume
15.94
31.34
67.66
The corresponding two-sided PBSCIs for comparing the test drugs with the control drugs are shown in Table 7. For comparisons, the SCIs obtained by Maurya et al. (2011a) and FGPQ method are also shown in this table. It is seen that the PBSCI procedure results in shorter intervals (the volumes of MSCIs are about 6 times the volumes of PBSCIs). Note, in particular, that the conclusions based on 99 % PBSCIs are different from the conclusions based on 99 % M-SCIs and SFGCIs: The 99 % PBSCIs indicate that test drug 2 is better than control drugs 1 and 2, and test drug1 is worse than control drugs 1 and 2; whereas the 99 % M-SCIs indicate only that test drug 2 is better than control drug 1, and test drug1 is worse than control drug 2; and the 99 % SFGCIs do not detect any difference between test drug 2 and control drug 2. Therefore, in this example our proposed PBSCI procedure performs better than M-SCI and SFGCI methods. Acknowledgments We are grateful to two reviewers for their valuable comments and suggestions. This work was supported by the Research Council of Shiraz University.
123
1094
A. Malekzadeh et al.
References Chang YP, Huang WT (2009) Simultaneous fiducial generalized confidence intervals for all pairwise comparisons of means. Int J Inf Manag Sci 20:459–467 Ferguson TS (1996) A course in large sample theory. Chapman & Hall, London Hannig J, Iyer H, Patterson P (2006a) Fiducial generalized confidence intervals. J Am Stat Assoc 101:254– 269 Hannig J, Lidong E, Abdel-Karim A, Iyer H (2006b) Simultaneous fiducial generalized confidence intervals for ratios of means of Lognormal distributions. Aus J Stat 35:261–269 Johnson NL, Kotz S (1970) Continuous univariatedistributions. Wiley, New York Karrati-Kopaei M, Malekzadeh A, Sadooghi-Alvandi SM (2013) Simultaneous fiducial generalized confidence intervals for the successive differences of exponential location Parameters under heteroscedasticity. Stat Probab Lett 83:1547–1552 Krishnamoorthy K, Lu F, Mathew T (2007) A parametric bootstrap approach for ANOVA with unequal variances: fixed and random models. Comput Stat Data Anal 51:5731–5742 Krishnamoorthy K, Lu F (2010) A parametric bootstrap solution to the MANOVA under heteroscedasticity. J Stat Comput Simul 80(8):837–887 Lam K (1987) Subset selection of normal population under heteroscedasticity, in Proceedings of the second international advanced seminar/workshop on inference procedures associated with statistical ranking and selection, Sydney, Australia Lam K (1988) An improved two-stage selection procedure. Commun Stat Simul Comput 17(3):995–1005 Lam K, Ng CK (1990) Two-stage procedures for comparing several exponential populations with a control when scale parameters are unknown and unequal. Sequ Anal 9(2):151–164 Lawless JF, Singhal K (1980) Analysis of data from life test experiments under an exponential model. Naval Res Logist Q 27:323–334 Li X, Wang J, Liang H (2011) Comparing of several means: a fiducial based approach. Comput Stat Data Anal 55:1993–2002 Ma CX, Tian LL (2009) A parametric bootstrap approach for testing equality of inverse Guassian means under heterogeneity. Comm Stat Simul Comput 38:1153–1160 Maurya V, Goyal A, Gill AN (2011a) Multiple comparisons with more than one control for exponential location parameters under heteroscedasticity. Comm Stat Simul Comput 40:621–644 Maurya V, Goyal A, Gill AN (2011b) Simultaneous testing for the successive differences of exponential location parameters under heteroscedasticity. Stat Probab Lett 8:1507–1517 Ng CK, Lam K, Chen HJ (1993) Multiple comparisons of exponential location parameters with the best under type II censoring. Am J Math Manag Sci 12:383–402 Sadooghi-Alvandi SM, Jafari AA (2013) A parametric bootstrap approach for one-way ANCOVA with unequal variances. Comm Stat Theory Methods 42:2473–2498 Singh P, Abebe A (2009) Comparing several exponential populations with more than one control. Stat Methods Appl 18:359–374 Tian LL, Ma CX, Vexler A (2009) A parametric bootstrap test for comparing heteroscedastic regression models. Comm Stat Simul Comput 38:1026–1036 Wu S, Lin YP, Yu YR (2010) One-stage multiple comparisons with the control for exponential location parameters under heteroscedasticity. Comput Stat Data Anal 54:1372–1380 Xu LW, Yang FQ, Abula A, Qin S (2012) A parametric bootstrap approach for two-way ANOVA in presence of possible interaction with unequal variances. J Multi Anal 115:172–180 Zelen M (1966) Application of exponential models problems in cancer research. J R Stat Soc 129:368–398
123