Next Article in Journal
PFVAE: A Planar Flow-Based Variational Auto-Encoder Prediction Model for Time Series Data
Next Article in Special Issue
A Novel Spatio-Temporal Fully Meshless Method for Parabolic PDEs
Previous Article in Journal
Efficient Malware Classification by Binary Sequences with One-Dimensional Convolutional Neural Networks
Previous Article in Special Issue
A Study on the Experimental Design for the Lifetime Performance Index of Rayleigh Lifetime Distribution under Progressive Type I Interval Censoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls

Department of Statistics, Tamkang University, Taipei 251301, Taiwan
Mathematics 2022, 10(4), 609; https://doi.org/10.3390/math10040609
Submission received: 23 December 2021 / Revised: 31 January 2022 / Accepted: 15 February 2022 / Published: 16 February 2022
(This article belongs to the Special Issue Applications of Partial Differential Equations in Engineering)

Abstract

:
Under heteroscedasticity, we propose one-stage multiple comparison procedures for several treatment groups compared with several control groups in terms of exponential mean lifetimes. The simultaneous confidence intervals including one-sided and two-sided confidence intervals for the difference between the mean lifetime from the i-th treatment group and the mean lifetime from the j-th control group are developed in this research. The required critical values are obtained and tabulated for the practical use of users. The experimenters can use these simultaneous confidence intervals to determine whether the treatment mean lifetimes are better than several controls or worse than several controls under a specified confidence level. Finally, one example of comparing the mean duration of remission using four drugs for treating leukemia is used for the aims of illustrations.

1. Introduction

Our research is related to the field of ranking and selection. For normal distributions, Bechhofer [1] and Gupta [2] are the pioneers in this field. Instead of normal distributions, exponential distributions are widely used to model the lifetimes of products. For exponential lifetime distributions, k(≥2) independent treatment groups and l control groups are considered. Let X i 1 , , X i m be the random sample of size m from the i-th treatment group π i , I = 1, …, k, where π i follows a two-parameter exponential distribution denoted by E ( θ i , σ i ) , I = 1, …, k. The parameters θ 1 , , θ k are the unknown location parameters (also known as the guaranteed time in reliability and engineering) and the parameters σ 1 , , σ k are the unknown and unequal scale parameters. Likewise, let X j 1 * , , X j m * be the random sample of size m from the j-th control group π j * , j = 1, …, l, where π j * follows a two-parameter exponential distribution denoted by E ( η j , ξ j ) , j = 1, …, l. The parameters η 1 , , η l are the unknown location parameters and parameters ξ 1 , , ξ l are the unknown and unequal scale parameters. Furthermore, the mean lifetime for the i-th treatment group is μ i = θ i +   σ i , I = 1, …, k and the mean lifetime for the j-th control group is μ j * = η j +   ξ j , j = 1, …, l. When scale parameters for exponential distributions are unknown and possibly unequal, Lam and Ng [3] proposed the design-oriented two-stage multiple comparisons with the control. Please see Lam and Ng [3] for the detailed introduction for the two-stage procedures. The disadvantage for the two-stage procedure is that the required sample size for the second stage may be large. To remedy this problem, the one-stage sampling procedure addressed above can be considered instead. Based on the one-stage sample from k populations, Wu et al. [4] developed the multiple comparison procedures for exponential location parameters with the control when scale parameters are unknown and unequal. Instead of comparing with the control, Wu [5] proposed one-stage multiple comparisons with the average for exponential location parameters. Maurya [6] investigated the one-stage multiple comparisons with more than one control under heteroscedasticity. Maurya [7] proposed another one-stage multiple comparisons with a control under heteroscedasticity. Instead of comparing the location parameters for exponential distributions, the one-stage multiple comparisons with the average in terms of mean lifetimes are proposed by Wu [8]. Unlike the multiple comparison procedure with the control based on mean lifetimes in Wu [8], the aim and the novelty of this paper is to compare several treatment populations with more than one control populations based on the mean lifetimes. The research methodologies are addressed in Section 2 by using the techniques in Lam’s research ([9,10]). The research results in Wu [8] can be treated as a special case of our proposed method in this paper when the number of control populations l = 1. Therefore, this paper has more general applications compared with the results in Wu [8]. In order to illustrate our proposed simultaneous confidence intervals in Section 2, we use an example of comparing two treatment drugs with two control drugs for treating leukemia in Section 3. In this example, each drug group consists of 20 patients and the dataset for each group passed the fitness test for exponential distribution by the use of the Gini test (Gail and Gastwirth [11]). In Lawless [12], they conducted a likelihood ratio asymptotic chi-squared test for the homogeneity of scale parameters and the results revealed that the scale parameters for four exponential distributions are significantly different. In Section 4, the conclusions related to our research methods are made.

2. Comparing with Several Controls for Exponential Mean Lifetimes

For the one-stage sample X i 1 , , X i m of size m from the i-th treatment group, we let Y i = min ( X i 1 , , X i m ) and S i = j = 1 m ( X i j Y i ) / ( m 1 ) , I = 1, …, k. Likewise, for the one-stage sample X j 1 * , , X j m * of size m from the j-th control group, we let Y j * = min ( X j 1 * , , X j m * ) and S j * = t = 1 m ( X j t * Y j * ) / ( m 1 ) , j = 1, …, l.
Before we start to construct the simultaneous confidence interval (SCI) for μ i μ j * , I = 1, …, k, j = 1, …, l, we need to consider the following random variables:
G i = m S i / σ i + m m ( Y i θ i ) / σ i / S i / σ i ,   I = 1 ,   ,   k   and
G j * = m S j * / ξ j + m m ( Y j * η j ) / ξ j / S j * / ξ j ,   j = 1 ,   . . . ,   l .
From Roussas [13], they indicated that Q i = 2 ( m 1 ) S i / σ i ~ χ 2 m 2 2 and Q j * = 2 ( m 1 ) S j * / ξ ~ χ 2 m 2 2 , where χ 2 m 2 2 represents a chi-squared distribution with 2m − 2 degrees of freedom (df). Furthermore,   E i = m ( Y i θ i ) / σ i ~Exp(1) and E j * = m ( Y j * η j ) / ξ j ~Exp(1), where Exp(1) represents a standard exponential distribution. In addition, the random variables E i   ,   Q i   ,   i = 1 , , k and   E j * , Q j * , j = 1 , , l are stochastically independent.
In terms of E i   ,   Q i   ,   i = 1 , , k and   E j * , Q j * , j = 1 , , l , the random variables G i , I = 1, …, k and G j * , j = 1, …, l can be rewritten as follows:
G i = m S i / σ i + m m ( Y i θ i ) / σ i S i / σ i = 2 m ( m 1 ) S i / σ i + 2 ( m 1 ) ( m m ( Y i θ i ) / σ i ) 2 ( m 1 ) S i / σ i = m + 2 ( m 1 ) ( m E i ) Q i
and
G j * = m S j * / ξ j + m m ( Y j * η j ) / ξ j S j * / ξ j = m + 2 ( m 1 ) ( m E j * ) Q j * .
Once the independent random variables E i ~ E x p ( 1 ) and Q i ~ χ 2 m 2 2 are generated, the random variable G i can be thus generated, I = 1, …, k. Likewise, once the independent random variables E i ~ E x p ( 1 ) and Q i ~ χ 2 m 2 2 , are generated, the random variable G j * can be generated, j = 1, …, l.
Now we can propose the one-stage multiple comparison procedures for exponential mean lifetimes with several controls denoted by μ i μ j * , I = 1, …, k, j = 1, …, l as follows:
Theorem 1.
For a given 0 < P < 1 , let c ˜ = max i = 1 , , k ,   j = 1 , , l   ( S i m , S j * m ) ,
G i = m S i / σ i + m m ( Y i θ i ) / σ i / S i / σ i ,   I = 1 ,   ,   k and
G j * = m S j * / ξ j + m m ( Y j * η j ) / ξ j / S j * / ξ j ,   j = 1 ,   ,   l .
We have the following results:
(a) 
If s ˜ U is the 100P-th percentile of the distribution of max ( G j * , G i , G i G j * ,     i = 1 , , k , j = 1 , , l ) , then we have
P ( μ i μ j * Y i + S i Y j * S j * + c ˜ s ˜ U ,   i = 1 , , k , j = 1 , , l ) P .
Thus, ( , Y i + S i Y j * S j * + c ˜ s ˜ U ) is a set of upper confidence intervals for μ i μ j * with confidence coefficient P, I = 1, …, k, j = 1, …, l.
(b) 
If s ˜ L is the 100P-th percentile of the distribution of max ( G i , G j * , G j * G i , i = 1 , , k , j = 1 , , l ) , then we have
P ( μ i μ j * Y i + S i Y j * S j * c ˜ s ˜ L ,   i = 1 , , k , j = 1 , , l ) P .
Thus, ( Y i + S i Y j * S j * c ˜ s ˜ L , ) is a set of lower confidence intervals for μ i μ j * with confidence coefficient P, I = 1, …, k, j = 1, …, l.
(c) 
If s ˜ t is the 100P-th percentile of the distribution of max ( | G i | , | G j * | , | G j * G i | , i = 1 , , k , j = 1 , , l ) , then we have
P ( Y i + S i Y j * S j * c ˜ s ˜ t μ i μ j * Y i + S i Y j * S j * + c ˜ s ˜ t , i = 1 , , k , j = 1 , , l ) P
Thus, ( Y i + S i Y j * S j * ± c ˜ s ˜ t ) is a set of simultaneous two-sided confidence intervals for μ i μ j * with confidence coefficient P, I = 1, …, k, j = 1, …, l.
The technique given in Lam ([9,10]) is described in the following lemma:
Lemma 1.
If X and Y are two random variables, a and b are two positive constants, then [ a X b Y d max ( a , b ) ] [ X d ,   Y d   and   X Y d ] .
The proof of Lemma 1 is given in Lam [1] and it is briefly given in the section of Appendix A. Using the technique in Lemma 1, the proof of Theorem 1 is also given in Appendix A.
Since the critical values s ˜ U , s ˜ L and s ˜ t in Theorem 1 are difficult to obtain, the Monte Carlo simulation method is used to find the critical values and the algorithm to obtain these critical values is given as follows:
Step 1:
Generate k independent random variables E i ~ E x p ( 1 ) and l independent random variables E j * ~ E x p ( 1 ) . Generate another k independent random variables Q i ~ χ 2 ( m 1 ) 2 and l independent random variables Q j * ~ χ 2 ( m 1 ) 2 .
Then, the k + l independent random variables:
G i = m + 2 ( m 1 ) ( m E i ) / Q i   = 1 ,   ,   k and
G j * = m + 2 ( m 1 ) ( m E j * ) / Q j * ,     = 1 ,   ,   l   can   be   generated .
Step 2:
The values of max ( G j * , G i , G i G j * ,     i = 1 , , k , j = 1 , , l ) ,
max ( G i , G j * , G j * G i , , i = 1 , , k , j = 1 , , l )   and
max ( | G i | , | G j * | , | G j * G i | , i = 1 , , k , j = 1 , , l )   can   be   computed .
Step 3:
Repeat Step 1 and Step 2 100,000 times. Then, the critical values s ˜ U , s ˜ L and s ˜ t can be obtained as the 100P-th empirical percentiles of
max ( G j * , G i , G i G j * ,     i = 1 , , k , j = 1 , , l ) ,  
max ( G i , G j * , G j * G i , , i = 1 , , k , j = 1 , , l ) and  
max ( | G i | , | G j * | , | G j * G i | , i = 1 , , k , j = 1 , , l )   respectively .  
The approximate critical values s ˜ U , s ˜ L , and s ˜ t under (k,l) = (2,2), (2,3), (2,4), (3,2), (3,3), (4,2), m = 2(1)10(5)30 and P = 0.90, 0.95, and 0.975 are listed in Table 1. From Table 1, it can be seen that the approximate critical values s ˜ U , s ˜ L , and s ˜ t are increasing while k + l is increasing for any given P and m or while P is increasing for any given k, l, and m. Let L1 be the length of the two-sided confidence intervals for μ i μ j * , I = 1, …, k, j = 1, …, l for the one-stage procedure, then we have L1 = 2 c ˜ s ˜ t . From this equation, we can see that the larger the k + l, the larger the value of s ˜ t and then the larger confidence length of L1 for any given P and m. Likewise, we can also see that when P increases, the confidence length increases for any given k, l, and m.

3. Example

From Wu et al. [4], data consisting of the duration of remission achieved by four drugs for treating leukemia are used to illustrate our proposed simultaneous confidence intervals in Theorem 1. The data are listed in Table 2. We regard drugs 1,2 as the treatment groups π 1 ,   π 2 and drugs 3,4 as the control groups π 1 * ,   π 2 * .
The longer mean duration of remission time ( μ i or μ j * , i = 1 , 2 , j = 1 , 2 ) is desired for this example. Lawless [12] conducted the likelihood ratio asymptotic χ 2 test to show that there is a significant difference among the four scale parameters. Thus, the data analysis one-stage multiple comparison procedures with several controls μ i μ j * , i = 1 , 2 , j = 1 , 2 in Theorem 1 for exponential mean lifetimes under heteroscedasticity can be applied. The required statistics and critical values of s ˜ U , s ˜ L , and s ˜ t for P = 0.90, 0.95, and 0.975 are summarized in Table 3.
Applying parts (a) and (b) of Theorem 1, we can find the upper and lower confidence bounds for μ i μ j * , i = 1 , 2 , j = 1 , 2 with confidence coefficients 0.90, 0.95, and 0.975 given in Table 4. With respect to the upper confidence bounds, drug 1 is selected in a subset of worse than several controls (drug 3 and 4) and drugs 1 and 2 are selected in a subset of worse than the control (drug 3) with the probability of correct selection being at least 0.90, 0.95, and 0.975 in terms of mean lifetimes since their one-sided confidence intervals do not contain zero.
Applying part (c) of Theorem 1, we can find the two-sided confidence bounds for μ i μ j * , i = 1 , 2 , j = 1 , 2 with confidence coefficients 0.90, 0.95, and 0.975 listed in Table 5. For confidence coefficients 0.90, 0.95, and 0.975, we can conclude that drug 1 is worse than several controls (drugs 3 and 4) and drug 2 is worse than drug 4 since both bounds of the two-sided confidence intervals for these two drugs are negative. Furthermore, drug 2 is not significantly different from drug 3 in terms of mean remission times.

4. Conclusions

From the critical values in Table 1 and the equation of L1 = 2 c ˜ s ˜ t , it is found that the length of SCI is getting wider when we compare more populations for fixed P and m.
In practical applications, we may have several control populations instead of one control population. Applying Lam’s [9,10] technique, we propose a multiple comparison procedure with several controls for exponential mean lifetimes under heteroscedasticity. In addition, the critical values for (k,l) = (2,2), (2,3), (2,4), (3,2), (3,3), (4,2), m = 2(1)10(5)30 and P* = 0.90, 0.95, and 0.975 are provided for the use of this investigation problem in Table 1. Finally, we use one example to illustrate the practical use of our proposed procedures in Theorem 1. The multiple comparison procedure with one control for exponential mean lifetimes in Wu [8] is only a special case of our new procedure in Theorem 1 of this research when l = 1. Therefore, the results in this paper have more general methods in many fields of applications.

Funding

This research was funded by the Ministry of Science and Technology, Taiwan, MOST 108-2118-M-032-001 and MOST 109-2118-M-032-001-MY2 and the APC was funded by MOST 109-2118-M-032-001-MY2.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available in a publicly accessible repository. The data presented in this study are openly available in Wu et al. [4].

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Proof of Lemma 1.
Consider the following three cases:
Case 1: For b a , [ X Y d and Y d ] implies
a X a ( Y d ) b ( Y d ) b Y d max ( a , b ) .
Case 2: For b < a and Y 0 ,   [ X d ] implies
a X a d b Y d max ( a , b ) .
Case 3: For b < a and Y > 0 ,   [ X Y d ] implies
a X a ( Y d ) b Y d max ( a , b ) .
The proof is established. □
Proof of Theorem 1.
For (a), we have
P μ i μ j * Y i + S i Y j * S j * + c ˜ s ˜ U ,   i = 1 , , k , j = 1 , , l
= P θ i + σ i η j ξ j Y i + S i Y j * S j * + c ˜ s ˜ U ,   i = 1 , , k , j = 1 , , l
= P S i + σ i + θ i Y i S j * + η j + ξ j Y j * + c ˜ s ˜ U ,   i = 1 , , k , j = 1 , , l
= P S i m σ i S i m ( S i + σ i + θ i Y i ) σ i S j * m η j S j * m ( S j * + η j + ξ j Y j * ) η j + c ˜ s ˜ U
,   i = 1 , , k , j = 1 , , l
= P S j * m G j * S i m G i c ˜ s ˜ U ,   i = 1 , , k , j = 1 , , l
E S 1 , , S k P S j * m G j * S i m G i max ( S i m , S j * m )   s ˜ U ,   i = 1 , , k , j = 1 , , l
P G j * s ˜ U ,   G i s ˜ U ,   G j * G i s ˜ U ,   i = 1 , , k , j = 1 , , l
= P G j * s ˜ U ,   G i s U * ,   G i G k s ˜ U ,   i = 1 , , k , j = 1 , , l
= P ( max ( G j * , G i , G i G j * ,     i = 1 , , k , j = 1 , , l ) s ˜ U ) = P
Satisfying the above equation, the value of s ˜ U is determined as the 100P-th percentile of the distribution of max ( G j * , G i , G i G j * ,     i = 1 , , k , j = 1 , , l ) and the proof is thus established.
For (b), we have
P ( μ i μ j * Y i + S i Y j * S j * c ˜ s ˜ L ,   i = 1 , , k , j = 1 , , l )
P θ i + σ i η j ξ j Y i + S i Y j * S j * + c ˜ s ˜ L ,   i = 1 , , k , j = 1 , , l
= P S i + σ i + θ i Y i S j * + η j + ξ j Y j * + c ˜ s ˜ L ,   i = 1 , , k , j = 1 , , l
= P S i m σ i S i m ( S i + σ i + θ i Y i ) σ i S j * m η j S j * m ( S j * + η j + ξ j Y j * ) η j + c ˜ s ˜ L
,   i = 1 , , k , j = 1 , , l
= P S i m G i S j * m G j * c ˜ s ˜ L ,   i = 1 , , k , j = 1 , , l
E S 1 , , S k P S i m G i i S j * m G j * max ( S i m , S j * m )   s ˜ L ,   i = 1 , , k , j = 1 , , l
P G i   s ˜ L ,   G j *   s ˜ L ,   G i G j *   s ˜ L ,   i = 1 , , k , j = 1 , , l
P G i   s ˜ L ,   G j *   s ˜ L ,   G j * G i   s ˜ L ,   i = 1 , , k , j = 1 , , l
= P ( max ( G i , G j * , G j * G i , i = 1 , , k , j = 1 , , l ) s ˜ L ) = P
Satisfying the above equation, the value of   s ˜ L is determined as the 100P-th percentile of the distribution of max ( G i , G j * , G j * G i , i = 1 , , k , j = 1 , , l ) and the proof is thus established.
For (c), combining (a) and (b), we have
P ( Y i + S i Y j * S j * c ˜ s ˜ t μ i μ j * Y i + S i Y j * S j * + c ˜ s ˜ t , i = 1 , , k , j = 1 , , l )
= P Y i + S i Y j * S j * c ˜ s ˜ t θ i + σ i η j ξ j Y i + S i Y j * S j * + c ˜ s ˜ t ,
,   i = 1 , , k , j = 1 , , l
= E S 1 , , S k P G i s ˜ t , G j * s ˜ t ,   G j * G i s ˜ t       G j * s ˜ t ,   G i s ˜ t ,   G i G j * s ˜ t ,   i = 1 , , k , j = 1 , , l
P ( max ( | G i | , | G j * | , | G j * G i | , i = 1 , , k , j = 1 , , l ) s ˜ t ) = P
Satisfying the above equation, the value of s ˜ t is determined as the 100P-th percentile of the distribution of max ( | G i | , | G j * | , | G j * G i | , i = 1 , , k , j = 1 , , l ) and the proof is thus established. □

References

  1. Bechhofer, R.E. A single sample multiple decision procedure for ranking means of normal populations with known variances. Ann. Math. Stat. 1954, 25, 16–39. [Google Scholar] [CrossRef]
  2. Gupta, S.S. On a Decision Rule for a Problem in Ranking Means. Doctoral Dissertation, Mimeograph Series No. 150. Institute of Statistics, University of North Carolina, Chapel Hill, NC, USA, 1956. [Google Scholar]
  3. Lam, K.; Ng, C.K. Two-stage procedures for comparing several exponential populations with a control when the scale parameters are unknown and unequal. Seq. Anal. 1990, 9, 151–164. [Google Scholar] [CrossRef]
  4. Wu, S.F.; Lin, Y.P.; Yu, Y.R. One-stage multiple comparisons with the control for exponential location parameters under heteroscedasticity. Comput. Stat. Data Anal. 2010, 54, 1372–1380. [Google Scholar] [CrossRef]
  5. Wu, S.F. New one-stage multiple comparisons procedures with the average for exponential location parameters under heteroscedasticity. J. Stat. Comput. Simul. 2016, 86, 2740–2748. [Google Scholar] [CrossRef]
  6. Maurya, V.; Goyal, A.; Gill, A.N. Multiple comparisons with more than one control for exponential location parameters under heteroscedasticity. Commun. Stat. Simul. Comput. 2011, 40, 621–644. [Google Scholar] [CrossRef]
  7. Maurya, V.; Gill, A.N.; Singh, P. Multiple comparisons with a control for exponential location parameters under heteroscedasticity. J. Appl. Stat. 2013, 40, 1817–1830. [Google Scholar] [CrossRef]
  8. Wu, S.F. One stage multiple comparisons of k-1 treatment mean lifetimes with the control for exponential distributions under heteroscedasticity. Commun. Stat. Simul. Comput. 2018, 47, 2968–2978. [Google Scholar] [CrossRef]
  9. Lam, K. Subset Selection of Normal Populations under Heteroscedasticity. IPASRAS-II: Proceedings and Discussions of the Second International Conference on Inference Procedures Associated with Statistal Ranking and Selection on the Frontiers of Modern Statistical Inference Procedures, II; ACM: New York, NY, USA, 1992; pp. 307–344. [Google Scholar]
  10. Lam, K. An improved two-stage selection procedure. Commun. Stat. Simul. Comput. 1988, 17, 995–1006. [Google Scholar] [CrossRef]
  11. Gail, M.H.; Gastwirth, J.L. A scale-free goodness-of-fit Test for the Exponential Distribution Based on the Gini Statistic. J. R. Stat. Soc. Ser. B 1978, 40, 350–357. [Google Scholar] [CrossRef]
  12. Lawless, J.F. Statistical Models and Methods for Lifetime Data; Wiley: New York, NY, USA, 2003. [Google Scholar]
  13. Roussas, G.G. A Course in Mathematical Statistics; Academic Press: San Diego, CA, USA, 1997. [Google Scholar]
Table 1. Critical values of s ˜ U , s ˜ L , and s ˜ t for P = 0.90, 0.95, and 0.975.
Table 1. Critical values of s ˜ U , s ˜ L , and s ˜ t for P = 0.90, 0.95, and 0.975.
klmP = 0.90P = 0.95P = 0.975
s ˜ U s ˜ L s ˜ t s ˜ U s ˜ L s ˜ t s ˜ U s ˜ L s ˜ t
22224.3224.3649.1750.1250.52100.37100.68102.56202.85
311.1311.1316.8417.0617.0625.1425.3525.2636.73
49.519.5213.2313.3413.3818.0718.1018.1924.08
59.189.2012.2712.3312.3916.0116.0316.1020.43
69.229.2212.0212.1112.1115.3115.3515.3619.02
79.409.4212.0612.1412.1315.0415.0815.0918.33
89.659.6412.2212.3012.2815.0715.1215.1018.16
99.909.9312.4412.4812.5315.1915.2215.2318.09
1010.2110.2312.7212.7812.7915.3915.4215.4318.24
1511.6711.6814.2314.2914.2816.8516.8716.8719.47
2013.0213.0315.7015.7515.7518.3518.3618.3920.98
2514.2814.3017.1017.1317.1719.8419.8419.8822.46
3015.4815.4218.4118.4918.4421.2721.3321.2624.01
23226.3035.0861.2453.6571.91124.81107.23145.94252.28
311.6514.0919.3017.7021.3028.6726.0731.4241.94
49.9611.5714.8313.8416.0120.0818.6721.5426.65
59.6110.9013.5512.8214.4417.5216.5418.5222.18
69.6810.8113.2012.5713.9516.6115.8117.4220.47
79.9010.8913.1512.6313.7816.2715.6216.9219.65
810.1411.0913.2812.7813.8716.1815.5816.8019.33
910.4311.3413.4813.0114.0316.2915.7516.8519.31
1010.7311.6513.7513.2914.2916.4815.9417.0119.38
1512.3113.1015.2814.9115.7517.9117.4918.3620.55
2013.7614.5116.8016.4717.2119.4619.0919.8522.05
2515.1315.8618.2917.9518.7121.0120.6421.4023.68
3016.3917.1019.6719.3820.0522.5422.2322.8825.34
24228.1446.0073.8556.7394.80150.80114.03190.63303.81
312.0816.6921.5218.1925.0131.7426.6736.6246.28
410.2513.1716.0714.1518.0521.6618.9624.1028.58
59.9012.2514.5913.1216.0418.7516.8420.4623.57
69.9811.9914.0812.8615.2717.6016.0818.9721.59
710.1911.9913.9712.9215.0117.1515.8818.2920.70
810.4612.1614.0713.0815.0417.0915.9018.1320.32
910.7612.4114.2813.3615.1817.1516.0518.1220.18
1011.1012.6514.5213.6715.3617.3116.3518.1720.24
1512.7614.1316.0515.3316.7818.7017.8819.4221.38
2014.3215.5917.6317.0118.2820.2919.6320.9022.92
2515.7016.9619.1318.5119.7621.8721.2222.4524.51
3017.0418.2520.5720.0121.1923.4122.8124.0026.16
32235.2026.4461.7872.3754.14126.61147.01108.88255.82
314.1111.7319.3421.3617.7628.7531.5026.1941.94
411.549.9514.7915.9613.8220.0121.3918.6726.57
510.909.6213.5614.4412.8217.5418.5716.5222.21
610.819.6913.2113.9512.5816.6517.4815.8320.50
710.889.8713.1513.7912.6116.2516.9315.6019.68
811.0710.1213.2513.8412.7516.1716.7815.5719.28
911.3610.4313.5014.0913.0216.3216.9315.7119.29
1011.6210.7613.7514.2713.3216.5017.0315.9819.37
1513.0912.3015.2515.7114.8717.8918.3517.4620.53
2014.5213.7816.8217.2516.4819.5019.9119.1322.13
2515.8515.1418.2818.6717.9821.0221.3720.7023.69
3017.0816.4019.6620.0419.3822.5122.8422.1825.23
33236.8937.0871.4275.8476.37149.80154.64151.29307.21
314.7514.6520.6221.6221.8430.3332.0332.5845.30
411.9311.9615.1416.5016.3220.7521.6922.0727.54
511.3411.3313.7214.9114.9217.9818.9819.2322.82
611.1811.2113.1914.4114.3816.7517.8117.8020.63
711.3511.3913.0914.2314.3316.2317.4117.3119.67
811.6011.5913.1314.2614.3416.0617.3417.3419.34
911.8911.8413.2714.6314.6316.2417.4317.3719.13
1012.1712.1813.4814.7814.8216.2117.6417.5119.15
1513.7013.8114.7916.3516.3317.4419.0818.8620.20
2015.2315.3316.2017.9517.9818.8720.5920.6721.55
2516.7416.7117.6019.3919.5420.1922.1822.2323.09
3018.0618.0718.8620.9320.8921.7023.6323.7724.37
42246.0628.2173.9394.2657.19150.66189.64114.78301.99
316.6312.0821.4624.9918.1631.7236.6326.6246.31
413.1910.2316.0718.0114.1121.6024.0018.9328.49
512.269.9014.5916.0513.0818.7720.5016.8323.66
611.989.9514.0915.3012.8617.6418.9916.1021.62
712.0210.1914.0115.0612.9217.2118.3615.9220.74
812.1710.4714.0915.0513.1217.1018.1315.9520.35
912.3710.7814.2515.1413.3817.1318.0616.1120.19
1012.6611.1214.5415.3913.6817.3218.2016.3520.24
1514.1212.7516.0416.7615.3118.6719.3917.9021.32
2015.5714.2817.6218.2916.9820.2720.9319.5722.89
2516.9615.6919.1419.7918.5321.8922.5221.2224.55
3018.2317.0120.5321.1619.9623.3923.9922.7726.09
Table 2. Data for the duration of remission by four drugs.
Table 2. Data for the duration of remission by four drugs.
Drug   1   ( μ 1 ) Drug   2   ( μ 2 ) Drug   3   ( μ 1 * ) Drug   4   ( μ 2 * )
1.0342.2144.1585.115
2.3444.9764.0254.498
1.2668.1545.1704.617
1.5632.68611.9094.651
1.1692.2714.9124.533
4.1183.1394.6294.513
1.0132.2143.9557.641
1.5094.4806.7355.971
1.1098.8473.14012.130
1.9652.23912.4464.699
5.1363.4738.7774.914
1.5332.7616.32117.169
1.7162.8333.2565.497
2.7782.3818.25011.332
2.5463.5483.75918.922
2.6262.4145.20513.712
3.4132.8323.0716.309
1.9295.5513.14710.086
2.0613.3769.7739.293
2.9512.96810.21811.787
Table 3. The required statistics and critical values.
Table 3. The required statistics and critical values.
StatisticsDrug 1Drug 2StatisticsDrug 3Drug 4
Yi1.0132.214 Y j * 3.0714.498
Si1.2381.530 S j * 3.2334.075
c ˜ 0.204
Y 1 Y 1 * +
S 1 S 1 *
Y 2 Y 1 * +
S 2 S 1 *
−4.053−2.560
Y 1 Y 2 * +
S 1 S 2 *
Y 2 Y 2 * +
S 2 S 2 *
−6.322−4.829
P s ˜ U s ˜ L s ˜ t
0.90013.0213.0315.70
0.95015.7515.7518.35
0.97518.3618.3920.98
Table 4. The upper confidence bounds and lower confidence bounds under P = 0.90, 0.95, 0.975.
Table 4. The upper confidence bounds and lower confidence bounds under P = 0.90, 0.95, 0.975.
Parameter ( , Y i + S i Y j * S j * + c ˜ s ˜ U ) ,   ( Y i + S i Y j * S j * c ˜ s ˜ L , )
90%95%97.5%
1. μ 1 μ 1 * (− ,−1.397), (−6.711, )(− ,−0.840), (−7.266, )(− ,−0.320), (−7.805, )
2. μ 1 μ 2 * (− ,−3.666), (−8.980, )(− ,−3.109), (−9.535, )(− ,−2.589), (−10.07, )
3. μ 2 μ 1 * (− ,0.096), (−5.218, )(− ,0.653), (−5.773, )(− ,1.173), (−6.312, )
4. μ 2 μ 2 * (− ,−2.173), (−1.616, )(− ,−1.616), (−8.042, )(− ,−1.096), (−8.581, )
Table 5. The two-sided confidence intervals under P = 0.90, 0.95, 0.975.
Table 5. The two-sided confidence intervals under P = 0.90, 0.95, 0.975.
Parameter ( Y i + S i Y j * S j * c ˜ s ˜ t ,   Y i + S i Y j * S j * + c ˜ s ˜ t )
90%95%97.5%
1. μ 1 μ 1 * (−7.256, −0.850)(−7.800, −0.310)(−8.333, 0.227)
2. μ 1 μ 2 * (−9.525, −3.119)(−10.07, −2.579)(−10.60, −2.042)
3. μ 2 μ 1 * (−5.763, 0.643)(−6.303, 1.183)(−6.840, 1.720)
4. μ 2 μ 2 * (−8.032, −1.626)(−8.572, −1.086)(−9.109, −0.549)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, S.-F. Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls. Mathematics 2022, 10, 609. https://doi.org/10.3390/math10040609

AMA Style

Wu S-F. Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls. Mathematics. 2022; 10(4):609. https://doi.org/10.3390/math10040609

Chicago/Turabian Style

Wu, Shu-Fei. 2022. "Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls" Mathematics 10, no. 4: 609. https://doi.org/10.3390/math10040609

APA Style

Wu, S. -F. (2022). Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls. Mathematics, 10(4), 609. https://doi.org/10.3390/math10040609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop