Next Article in Journal
Finite Element Solution of the Corona Discharge of Wire-Duct Electrostatic Precipitators at High Temperatures—Numerical Computation and Experimental Verification
Next Article in Special Issue
Closed-Form Solutions and Conserved Vectors of a Generalized (3+1)-Dimensional Breaking Soliton Equation of Engineering and Nonlinear Science
Previous Article in Journal
A Multi-Criteria Pen for Drawing Fair Districts: When Democratic and Demographic Fairness Matter
Previous Article in Special Issue
On Generalized Fourier’s and Fick’s Laws in Bio-Convection Flow of Magnetized Burgers’ Nanofluid Utilizing Motile Microorganisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity

Department of Statistics, Tamkang University, Tamsui, New Taipei City 251301, Taiwan
Mathematics 2020, 8(9), 1405; https://doi.org/10.3390/math8091405
Submission received: 21 July 2020 / Revised: 18 August 2020 / Accepted: 18 August 2020 / Published: 21 August 2020
(This article belongs to the Special Issue Applications of Partial Differential Equations in Engineering)

Abstract

:
When the additional sample for the second stage may not be available, one-stage multiple comparisons for exponential median lifetimes with the control under heteroscedasticity including one-sided and two-sided confidence intervals are proposed in this paper since the median is a more robust measure of central tendency compared to the mean. These intervals can be used to identify treatment populations that are better than the control or worse than the control in terms of median lifetimes in agriculture, stock market, pharmaceutical industries. Tables of critical values are obtained for practical use. An example of comparing the survival days for four categories of lung cancer in a standard chemotherapeutic agent is given to demonstrate the proposed procedures.

1. Introduction

In the lifetime test problems, the lifetime of some products follows an exponential distribution (see Lawless [1]). As the lifetime of products possesses a two-parameter exponential distribution, this research focuses on the development of multiple comparison procedures with the control population in terms of median lifetime. We consider k   ( 2 ) independent populations π 1 , , π k , where π i possesses a two-parameter exponential distribution denoted by E ( θ i , σ i ) , i = 1 , , k . The kth population is regarded as the control population and the first k − 1 populations are regarded as the treatment populations. The location parameters θ 1 , , θ k are unknown and usually called the guaranteed time in reliability analysis. The unequal and unknown scale parameters σ 1 , , σ k are regarded as the mean lifetime minus their location parameters θ 1 , , θ k since the mean lifetime for the ith population is μ i = θ i + σ i , i = 1 , , k . Regarding the multiple comparisons with the control population, Ng et al. [2] proposed a procedure in terms of the location parameter under the assumption of equal scale parameters. Under heteroscedasticity (unequal scale parameters), Lam and Ng [3] developed a design-oriented two-stage multiple comparison procedure in terms of location parameters. For the problem of multiple comparisons with the average under heteroscedasticity, Wu and Wu [4] investigated the two-stage procedures in terms of exponential location parameters. However, it may happen that the experimenters are not able to collect the additional sample for the second stage for the two-stage procedure. Because of this reason, Wu et al. [5] propose one-stage multiple comparisons with the average instead. Wu [6] proposes an one-stage multiple comparisons with the control for exponential distributions in terms of mean lifetimes under heteroscedasticity. Instead of doing the multiple comparisons with the control based on mean lifetimes, the median lifetimes should be considered since the median lifetimes are more robust for measuring the central tendency of the exponential lifetime distributions than mean lifetimes. Therefore, we consider the multiple comparisons with the control in terms of median lifetimes instead of mean lifetimes in this study. The median lifetime for the ith population is obtained as δ 0.5 , i = θ i + ln 2   σ i , i = 1 , , k . The one-sided and two-sided confidence interval of the ith median lifetime deviated from the control median lifetime denoted by δ 0.5 , i δ 0.5 , k , i = 1 , , k − 1 are proposed in the next section. All critical values are computed and listed in a table for the application of users. In Section 3, an example of comparing the survival days for four categories of lung cancer is considered for the illustrative aims to illustrate the implementation of proposed methods. Finally, our conclusions are summarized in Section 4.

2. One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity

The exponential distribution is primarily used in reliability applications and this distribution is used to model data with a constant failure rate, see for example, Johnson et al. [7]. The probability density function (pdf) and cumulative distribution function (cdf) for the ith exponentially distributed population are defined as
f ( x ) = 1 σ i exp ( x θ i σ i ) ,   x θ i , θ i > 0 , σ i > 0 ,
and
F ( x ) = 1 exp ( x θ i σ i ) , x θ i , θ i > 0 , σ i > 0 , i = 1 , , k .
Location parameters θ 1 , , θ k and scale parameters σ 1 , , σ k are unknown and possibly unequal. The survival function for the ith exponentially distributed population is S ( x ) = exp ( x θ i σ i ) . The qth quantile of the ith exponential distribution denoted by η q , i can be obtained by solving F ( x ) = 1 exp ( x θ i σ i ) = q and results in δ q , i = θ i ln ( 1 q )   σ i , i = 1 , , k . That is there are at least 1 − q percentage of products having lifetime longer than δ q , i units. Let q = 0.5 and then the median lifetimes can be obtained as δ 0.5 , i = θ i + ln 2   σ i , i = 1 , , k . In other words, there are at least 50% of products having lifetime longer than δ 0.5 , i units.
Lam and Ng [3] proposed two-stage multiple comparison procedures with the control. But for some reasons like lacking budget or encountering experimental difficulties, it is possible that the experimenters cannot collect the additional sample for the second stage. In this case, one-stage procedures should be considered instead. Therefore, we propose one-stage multiple comparison procedures for exponential median lifetimes with the control as follows:
Take a random sample of size m   ( 2 ) from π i denoted by X i 1 , , X i m for the one-stage procedures. Let Y i = min ( X i 1 , , X i m ) and S i = j = 1 m ( X i j Y i ) / ( m 1 ) and let
c * = max i = 1 , , k   S i m .
It is well-known that the complete sufficient statistics for ( θ i , σ i ) are ( Y i , S i ). From Roussas [8], the following three distributional results are observed.
(1) Q i = 2 ( m 1 ) S i / σ i = , i = 1 , , k has a chi-square distribution with 2 m 2 degrees of freedom (df) denoted by χ 2 m 2 2 .
(2) E i = m ( Y i θ i ) / σ i has a standard exponential distribution denoted by E x p ( 1 ) .
(3) E i and Q i are independent.
Using the distribution results of (1) and (2), we can find the uniformly minimum variance unbiased estimator (UMVUE) for ( θ i , σ i ) as ( Y i S i / m , S i ). Furthermore, the UMVUE for the ith median lifetime δ 0.5 , i = θ i + ln 2   σ i is Y i S i / m + ln 2 S i = Y i + ( m ln 2 1 ) S i / m . Then we find the UMVUE for the ith median lifetime deviated from the control median lifetime denoted by δ 0.5 , i δ 0.5 , k as Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m . Based on this estimator, we are going to propose the simultaneous confidence intervals for δ 0.5 , i δ 0.5 , k = θ i + ln 2   σ i θ k + ln 2   σ k , i = 1 , , k 1 in Theorem 1.
For the ith population, we consider the pivotal quantity
   G i = δ 0.5 , i Y i ( m ln 2 1 ) S i / m S i / m = θ i + ln 2   σ i Y i ( m ln 2 1 ) S i / m S i / m , = ( m ln 2 1 ) ( S i / σ i ) + m ln 2 m ( Y i θ i ) / σ i S i / σ i , i = 1 , , k .
Based on these pivotal quantities, we propose the one-stage multiple comparison procedures with the control in terms of median lifetimes in Theorem 1.
Theorem 1.
For a given 0 < P * < 1 , we can find the upper confidence bounds, lower confidence bounds, and two-sided confidence intervals for δ 0.5 , i δ 0.5 , k , i = 1, …, k−1 as follows:
(a) 
If s U * is the 100 P t h percentile of the distribution of U, where U = max ( G k , G i , G i G k , i = 1 , , k 1 ), then the simultaneous P* upper confidence bounds for δ 0.5 , i δ 0.5 , k are Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m + c * s U * , i = 1 , , k 1 .
(b) 
If s L * is the 100 P t h percentile of the distribution of L, where L= max ( G i , G k , G k G i , i = 1 , , k 1 ) , then the simultaneous P* lower confidence bounds for δ 0.5 , i δ 0.5 , k are Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m c * s L * , i = 1 , , k 1 .
(c) 
If s t * is the 100 P t h percentile of the distribution of T, where T = max ( | G i | , | G k | , | G k G i | , i = 1 , , k 1 ) , then the simultaneous P* two-sided confidence intervals for δ 0.5 , i δ 0.5 , k are
( Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m ± c * s t * ) ,   i = 1 , , k 1 .
The technique we use to prove the above Theorem is the following Lemma given in Lam [9,10]:
Lemma 1.
Suppose 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 ] .
Proof of Theorem 1.
 
For (a), we have
P ( δ 0.5 , i δ 0.5 , k Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m + c * s U * ,   i = 1 , , k 1 ) = P ( θ i + ln 2 σ i θ k ln 2 σ k Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m + c * s U * ,   i = 1 , , k 1 ) = P ( ( m ln 2 1 ) S i / m + ln 2 σ i + θ i Y i ( m ln 2 1 ) S k / m + ln 2 σ k + θ k Y k + c * s U * ,   i = 1 , , k 1 ) = P ( S i m σ i S i m ( ( m ln 2 1 ) S i / m + ln 2 σ i + θ i Y i ) σ i S k m σ k S k m ( ( m ln 2 1 ) S k / m + ln 2 σ k + θ k Y k ) σ k + c * s U * ,   i = 1 , , k 1 ) = P ( S k m G k S i m G i c * s U * ,   i = 1 , , k 1 ) E S 1 , , S k P ( S k m G k S i m G i max ( S i m , S k m )   s U * ,   i = 1 , , k 1 ) P ( G k s U * ,   G i s U * ,   G k G i s U * ,   i = 1 , , k 1 ) ( By   using   Lemma 1 ) = P ( G k s U * ,   G i s U * ,   G i G k s U * ,   i = 1 , , k 1 ) = P ( max ( G k , G i , G i G k ,   i = 1 , , k 1 ) s U * ) = P *
It is clear that s U * represents the 100 P t h percentile of the distribution of U and thus the proof is completed.
For (b), we have
P ( δ 0.5 , i δ 0.5 , k Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m c * s L * ,   i = 1 , , k 1 ) = P ( θ i + ln 2 σ i θ k ln 2 σ k Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m c * s L * ,   i = 1 , , k 1 ) = P ( S i m σ i S i m ( ( m ln 2 1 ) S i / m + ln 2 σ i + θ i Y i ) σ i S k m σ k S k m ( ( m ln 2 1 ) S k / m + ln 2 σ k + θ k Y k ) σ k c * s L * ,   i = 1 , , k 1 ) = P ( S i m G i S k m G k c * s L * ,   i = 1 , , k 1 ) E S 1 , , S k P ( S i m G i S k m G k max ( S i m , S k m ) s L * ,   i = 1 , , k 1 ) P ( G i s L * ,   G k s L * ,   G i G k s L * ,   i = 1 , , k 1 ) ( By   using   Lemma   1 ) P ( G i s L * ,   G k s L * ,   G k G i s L * ,   i = 1 , , k 1 ) = P ( max ( G i , G k , G k G i , i = 1 , , k 1 ) s L * ) = P *
It is clear that s L * represents the 100 P t h percentile of the distribution of L and thus the proof is completed.
For (c), combining (a) and (b), we have
P ( Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m c * s t * δ 0.5 , i δ 0.5 , k Y i + ( m ln 2 1 ) S i / m Y k ( m ln 2 1 ) S k / m + c * s t * ,   i = 1 , , k 1 ) = E S 1 , , S k P ( G i s t * ,   G k s t * ,   G k G i s t *         G k s t * ,   G i s t * ,   G i G k s t * ,   i = 1 , , k 1 ) P ( max ( | G i | , | G k | , | G k G i | , i = 1 , , k 1 ) s t * ,   i = 1 , , k 1 ) = P * .
It is clear that s t * represents the 100 P t h percentile of the distribution of T and thus the proof is completed. □
When the lifetime of products follows a two-parameter exponential distribution, this theorem can be used to find the upper confidence bounds and the lower confidence bounds for the parameters of δ 0.5 , i δ 0.5 , k , i = 1, …, k − 1, where δ 0.5 , k represents the median lifetime of the control population. This theorem can also be used to find the two-sided simultaneous confidence intervals for parameters δ 0.5 , i δ 0.5 , k , i = 1, …, k − 1. Based on these estimations, experimenters can identify better-than-the-control, worse-than-the-control, and not-much-different-from-the-control treatment populations in terms of median lifetimes. The real-life example to demonstrate the application of this theorem is given in Section 3.
It is very difficult to derive the p.d.f. or c.d.f. for U, L, and T. Using the above three distributional results (1)~(3), we observe that
G i = ( m ln 2 1 ) ( ν S i / σ i ) + m ν ln 2 m ν ( Y i θ i ) / σ i ν S i / σ i ( m ln 2 1 ) + ν ( m ln 2 E x p ( 1 ) ) χ 2 m 2 2 , i = 1 , , k , ν = 2 m−2. If we can generate independent random variables G i , i = 1 , , k , then we can find the empirical distribution of U, L, and T and the critical values s U * , s L * , and s t * are the empirical 100 P t h percentiles of the distributions of U, L, and T, through Monte-Carlo simulation methods.
The steps to find the critical values of s U * , s L * , and s t * in theorem 1 are enumerated as follows:
Step 1: Generate k independent random variables E i ~ E x p ( 1 ) and another k independent random variables Q i ~ χ 2 m 2 2 and then obtain the k independent random variables G i = ( m ln 2 1 ) + ν ( m ln 2 E i ) Q i ,   i = 1 , , k .
Step 2: Compute U = max ( G k , G i , G i G k , i = 1 , , k 1 ) ;
L = max ( G i , G k , G k G i , i = 1 , , k 1 ) and
T = max ( | G i | , | G k | , | G k G i | , i = 1 , , k 1 ) , i = 1 , , k .
Step 3: Repeat Steps 1,2 for 100,000 times. After sorting, we have U ( 1 ) U ( 100,000 ) ; L ( 1 ) L ( 100,000 ) ; T ( 1 ) T ( 100,000 ) .
Step 4: The critical values are obtained as s U * = U ( [ 100,000 * P * ] + 1 ) ; s L * = L ( [ 100,000 * P * ] + 1 ) ; s t * = T ( [ 100,000 * P * ] + 1 ) , where [x] is the largest integer less than or equal to x.
Remark: 
For unequal initial sample sizes denoted as m i , i = 1 , , k , Theorem 1 can be modified by replacing m by m i , i = 1 , , k .
For the practical use of application, we find the critical values s U * , s L * , and s t * by using the above algorithm under k = 3,4,…,10, m = 2,3, …, 10,15,20,25,30 and P* = 0.90, 0.95 and 0.975. The critical values are listed in the following table. From part (c) of Theorem 1, we observe that the length of the two-sided confidence intervals for δ 0.5 , i δ 0.5 , k is L1 = 2 c * s t * . The larger the critical values, the larger the confidence length when c* is fixed. From Table 1, we observe that the critical value s t * increases when P * increases for fixed k and m or when k increases for fixed P * and m. Therefore, the confidence length L1 increases when P * increases for fixed k and m or when k increases for fixed P * and m.

3. Example

Referring to Maurya et al. [11], the example of survival days of patients with inoperable lung cancer who were subjected to a standard chemotherapeutic agent is used to illustrate our proposed multiple comparison procedures with the control in Theorem 1. The patients are divided into four categories based on the histological type of their tumor: squamous, small, adeno, and large. The data are a part of a larger data set collected by the Veterans Administrative Lung Cancer Study Group in the United States. The survival days of 9 patients for four kinds of lung cancer are listed in Table 2:
Maurya et al. [11] had indicated that the data in the four categories may be assumed to be drawn from the two-parameter exponential distributions E ( θ i , σ i ) , i = 1 , , k . We regard the first category as the control population.
The required statistics and critical values for P * = 0.90 ,   0.95 and 0.975 are summarized in Table 3.
The upper confidence bounds and the lower confidence bounds for δ 0.5 , i δ 0.5 , 1 , i = 2,3,4 under confidence coefficients 0.90, 0.95, and 0.975 are listed in Table 4 using parts (a) and (b) in Theorem 1. Since all upper bounds are positive, no categories are selected in a subset of all treatment populations which are worse than the control population (Category 1 lung cancer). Since only Category 4 has positive lower bound for all confidence coefficients, we conclude that only Category 4 is selected in a subset of all treatment populations which are better than the control population with the probability of correct selection being at least 0.90, 0.95, and 0.975.
The two-sided confidence intervals for δ 0.5 , i δ 0.5 , 1 , i = 2,3,4 with confidence coefficients 0.90, 0.95, and 0.975 are computed using part (c) in Theorem 1 and the results are listed in Table 5. For confidence levels of 0.90 and 0.95, only the confidence interval for Category 4 does not contain zero and the lower limit is positive. We conclude that only Category 4 has median survival days better than Category 1. For confidence level of 0.975, no categories are identified to have median survival days greater than Category 1 in terms of median survival days.

4. Conclusions

We analyze the impact of confidence levels P* and number of population k on the confidence length in this paper. Instead of doing multiple comparisons with the control for exponential mean lifetimes, we propose multiple comparison procedures with the control in terms of median lifetimes since the measurement of median lifetimes are more robust than mean lifetimes for the measurement of central tendency for exponential lifetime distribution. For the illustrative aim, we give a real life example to illustrate how to find the upper bounds, lower bounds, and two-sided confidence intervals for our parameters related to median lifetimes.

Funding

This research was funded by [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.

Acknowledgments

The author wish to thank an associate editor and referees for their careful reading and valuable suggestions so that the article is more readable and applicable. The author’s research was supported by Ministry of Science and Technology MOST 108-2118-M-032-001- and MOST 109-2118-M-032 -001 -MY2 in Taiwan, ROC.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lawless, J.F. Statistical Models and Methods for Lifetime Data; Wiley: New York, NY, USA, 2003. [Google Scholar]
  2. Ng, C.; Lam, K.; Chen, H. Multiple Comparison of Exponential Location Parameters with the Best under Type II Censoring. Am. J. Math. Manag. Sci. 1992, 12, 383–402. [Google Scholar] [CrossRef]
  3. Lam, K.; Nag, C. 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.; Wu, C.-C. Two stage multiple comparisons with the average for exponential location parameters under heteroscedasticity. J. Stat. Plan. Inference 2005, 134, 392–408. [Google Scholar] [CrossRef]
  5. Wu, S.-F. One stage multiple comparisons with the average for exponential location parameters under heteroscedasticity. Comput. Stat. Data Anal. 2013, 68, 352–360. [Google Scholar] [CrossRef]
  6. 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]
  7. Johnson, N.L.; Kotz, S.; Balakrishnan, N. Continuous Univariate Distributions; Wiley: New York, NY, USA, 1994. [Google Scholar]
  8. Roussas, G.G. A Course in Mathematical Statistics; Academic Press: San Diego, CA, USA, 1997. [Google Scholar]
  9. Lam, K. Subset selection of normal populations under heteroscedasticity. In 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. Maurya, V.; Goyal, A.; Gill, A.N. Simultaneous testing for the successive differences of exponential location parameters under heteroscedasticity. Stat. Probab. Lett. 2011, 81, 1507–1517. [Google Scholar] [CrossRef]
Table 1. Approximate critical values of s U * , s L * , and s t * .
Table 1. Approximate critical values of s U * , s L * , and s t * .
kmP* = 0.90P* = 0.95P* = 0.975
s U * s L * s t * s U * s L * s t * s U * s L * s t *
3214.7811.2826.2630.4023.1053.5061.5246.58107.45
37.315.719.9211.168.7914.7816.5513.0021.59
46.325.038.008.897.1610.9412.039.7814.56
56.124.947.528.286.779.8510.748.8712.60
66.125.007.408.086.729.4810.278.6111.82
76.225.147.458.076.819.3910.088.6211.52
86.365.327.578.146.979.4310.068.7111.42
96.525.477.728.307.129.5310.138.8611.46
106.705.677.908.447.319.6610.249.0111.53
157.556.538.839.338.2910.6011.1110.0212.37
208.367.389.7710.259.2611.5912.0511.0613.37
259.128.1310.6211.0810.1212.5212.9712.0114.36
309.818.8311.4311.9010.9513.4413.8712.9515.34
4221.0513.9434.9642.8728.4870.7785.8257.02142.45
39.176.3311.7213.879.5817.3620.3714.0625.14
47.645.449.1410.567.6312.3714.1010.2916.34
57.245.278.449.587.1310.9312.299.2713.86
67.185.338.269.287.0710.4511.638.9912.93
77.235.478.279.207.1410.3011.348.9512.51
87.345.658.359.247.3010.2711.259.0512.34
97.485.848.499.327.4910.3611.239.2212.32
107.646.038.669.467.6910.4911.329.4012.38
158.536.989.6010.338.7211.3812.1110.4513.19
209.407.8710.5911.289.7512.4313.0911.5714.23
2510.208.6811.4912.1610.6813.3814.0112.5815.22
3010.959.4512.3513.0011.5614.3214.9313.5516.22
5227.2116.5043.6455.7833.6188.55112.9166.81178.20
310.756.8813.2716.0910.3519.5323.5215.0728.26
48.695.7710.0811.867.9913.4915.8010.6717.71
58.125.569.2110.657.4311.8413.609.5914.91
67.965.618.9410.197.3611.2312.679.3013.80
77.975.728.8810.027.4010.9712.269.2013.30
88.075.878.9410.027.5310.9412.129.3013.07
98.196.089.0710.077.7510.9812.069.4813.00
108.376.289.2510.227.9511.1012.129.6713.04
159.247.2810.1911.059.0511.9812.8510.7813.78
2010.138.2211.1712.0110.0813.0113.8411.8914.82
2510.969.1012.1212.9211.0914.0214.7812.9915.84
3011.749.8812.9913.7811.9814.9715.7113.9816.85
6233.1419.1551.7667.6438.55105.07136.1976.71211.14
312.177.3814.6618.1510.9821.5526.4715.8631.09
49.586.0610.8913.008.3214.5717.2611.0819.08
58.855.789.8311.517.6812.6314.619.8415.86
68.625.789.4910.947.5511.8713.509.5014.54
78.575.939.4010.697.6011.5513.009.4313.92
88.656.099.4410.647.7711.4612.799.5313.66
98.766.279.5410.697.9311.5012.759.6813.58
108.906.479.6910.788.1411.5812.749.8413.56
159.787.5210.6411.619.2712.4713.4410.9914.28
2010.708.4811.6512.5810.3613.4914.4112.1615.30
2511.579.3812.6313.5211.3714.5315.4013.2816.37
3012.3510.2113.5014.3612.3215.4616.2814.3117.35
7239.3121.7860.4980.0644.01122.99162.0687.46249.19
313.477.8215.9320.0111.5523.4229.2416.6333.86
410.356.3111.6113.978.6315.4618.4811.4320.25
59.465.9910.3912.247.9113.2515.4510.0916.59
69.165.969.9511.537.7212.3914.209.6715.10
79.086.099.8311.247.7912.0213.619.6214.41
89.136.249.8611.187.9211.9313.389.6914.16
99.236.439.9411.208.0911.9313.279.8113.99
109.366.6310.0811.288.2912.0113.2710.0214.02
1510.247.7011.0212.089.4612.8513.9111.1814.67
2011.148.7012.0213.0210.5513.8714.8412.3815.67
2512.019.6313.0013.9611.6114.8915.8213.5016.73
3012.8410.4713.9014.8512.5615.8616.7814.5217.75
8245.5024.2768.9792.7948.74140.29186.2097.45282.34
314.608.2217.1121.6312.1225.0331.4217.3435.96
411.086.5212.2814.928.8716.3019.6811.6921.34
510.026.1510.8812.928.0813.8516.2810.2517.30
69.636.1210.3712.077.8812.8814.879.8415.69
79.536.2310.2211.777.9212.4914.209.7414.95
89.556.3810.2111.648.0512.3113.879.8114.57
99.656.5710.3111.678.2312.3413.809.9514.46
109.756.7710.4111.698.4412.3513.7010.1514.40
1510.627.8411.3512.479.6113.1914.3411.3715.04
2011.548.8812.3513.4110.7514.2015.2612.5216.03
2512.439.8113.3314.3611.7915.2316.2213.6717.06
3013.2710.6814.2715.2812.7816.2417.1914.7618.12
9251.5526.6777.45105.2953.72158.46213.80108.06319.87
315.708.6018.2323.2412.6226.6633.7318.0538.33
411.706.7412.8515.769.1217.0820.7611.9722.33
510.516.2911.3213.488.2214.3716.9410.4217.92
610.086.2510.7612.608.0013.3215.429.9416.18
79.936.3410.5712.198.0312.8614.669.8615.37
89.926.5110.5512.058.1812.6814.309.9614.96
99.976.6810.5912.018.3312.6214.1510.0614.78
1010.106.8910.7312.088.5412.7014.1210.2814.75
1510.957.9811.6412.829.7313.4814.6811.4815.34
2011.889.0112.6513.7610.8814.4915.5912.6816.30
2512.779.9813.6314.7111.9615.5316.5813.8617.39
3013.6310.8514.5715.6412.9516.5217.5414.9218.41
10257.4029.3486.01117.1459.19175.47237.09118.10352.95
316.708.9519.2424.5913.0428.0235.5918.6340.30
412.286.9213.3916.469.3417.7921.6612.2423.22
510.946.4511.7214.008.4014.8717.5810.6118.54
610.456.3511.1113.028.1313.7315.9310.0916.68
710.276.4310.8812.598.1413.2015.109.9615.78
810.236.6010.8112.388.2712.9914.6910.0515.31
910.336.7910.9012.418.4512.9814.5510.1815.14
1010.436.9811.0112.418.6513.0014.5110.3715.09
1511.238.0911.8813.119.8513.7314.9911.5915.62
2012.179.1412.8814.0410.9914.7415.8912.7916.57
2513.0810.1013.8915.0212.0815.7816.8813.9517.65
3013.9511.0214.8415.9413.0916.7917.8315.0618.66
Table 2. Survival times for four categories of lung cancer.
Table 2. Survival times for four categories of lung cancer.
CategorymSurvival Times
1Squamous97210811101004282511
2Small9301323162118202731
3Adeno98923511713212162395
4Large9177162553200156182143105103
Table 3. The required statistics and critical values.
Table 3. The required statistics and critical values.
StatisticsCategory 1Category 2Category 3Category 4
Yi8133103
Si48.37510.25078.265106.750
c * 11.861
Yi − Y15−595
( m ln 2 1 ) (Si − S1)/m−22.19017.39733.976
YiY1 + ( m ln 2 1 ) (Si − S1)/m−17.19012.397128.976
P* s U * s L * s t *
0.9007.485.848.49
0.9509.327.4910.36
0.97511.239.2212.32
Table 4. The 90%, 95%, and 97.5% upper bounds (the number before comma) and lower bounds (the number after comma) for three categories compared with the control category (Category 1).
Table 4. The 90%, 95%, and 97.5% upper bounds (the number before comma) and lower bounds (the number after comma) for three categories compared with the control category (Category 1).
Y i Y 4 + ( m ln 2 1 ) ( S i S 4 ) / m + c * s U * , Y i Y 4 + ( m ln 2 1 ) ( S i S 4 ) / m c * s L *
90%95%97.5%
δ 0.5 , 2 δ 0.5 , 1 71.530, −86.45893.355, −106.029116.009, −126.548
δ 0.5 , 3 δ 0.5 , 1 101.117, −56.871122.942, −76.442145.596, −96.961
δ 0.5 , 4 δ 0.5 , 1 217.696, 59.708239.521, 40.137262.175, 19.618
Table 5. The 90%, 95% and 97.5% two-sided confidence intervals for three categories compared with the control category (Category 1).
Table 5. The 90%, 95% and 97.5% two-sided confidence intervals for three categories compared with the control category (Category 1).
Y i Y 4 + ( m ln 2 1 ) ( S i S 4 ) / m + c * s t * , Y i Y 4 + ( m ln 2 1 ) ( S i S 4 ) / m c * s t *
90%95%97.5%
δ 0.5 , 2 δ 0.5 , 1 −117.89,83.51−140.070, 105.690−163.318, 128.937
δ 0.5 , 3 δ 0.5 , 1 −88.303, 113.097−110.483, 135.277−133.730, 158.525
δ 0.5 , 4 δ 0.5 , 1 28.276, 229.6766.096, 251.856−17.151, 275.104

Share and Cite

MDPI and ACS Style

Wu, S.-F. One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity. Mathematics 2020, 8, 1405. https://doi.org/10.3390/math8091405

AMA Style

Wu S-F. One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity. Mathematics. 2020; 8(9):1405. https://doi.org/10.3390/math8091405

Chicago/Turabian Style

Wu, Shu-Fei. 2020. "One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity" Mathematics 8, no. 9: 1405. https://doi.org/10.3390/math8091405

APA Style

Wu, S. -F. (2020). One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity. Mathematics, 8(9), 1405. https://doi.org/10.3390/math8091405

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