Exact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring
|
|
- Silvester Stevenson
- 5 years ago
- Views:
Transcription
1 Exact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring A. Ganguly, S. Mitra, D. Samanta, D. Kundu,2 Abstract Epstein [9] introduced the Type-I hybrid censoring scheme as a mixture of Type-I and Type-II censoring schemes. Childs et al. [5] introduced the Type-II hybrid censoring scheme as an alternative to Type-I hybrid censoring scheme, and provided the exact distribution of the maximum likelihood estimator of the mean of a one parameter exponential distribution based on Type-II hybrid censored samples. The associated confidence interval also has been provided. The main aim of this paper is to consider a two-parameter exponential distribution, and to derive the exact distribution of the maximum likelihood estimators of the unknown parameters based on Type-II hybrid censored samples. The marginal distributions and the exact confidence intervals are also provided. The results can be used to derive the exact distribution of the maximum likelihood estimator of the percentile point, and to construct the associated confidence interval. Different methods are compared using extensive simulations and one data analysis has been performed for illustrative purposes. Key Words and Phrases: Hybrid censoring; two-parameter exponential distribution; joint moment generating function; maximum likelihood estimators; bootstrap confidence interval. Department of Mathematics and Statistics, IIT Kanpur, Kanpur, Pin 20806, India. 2 Corresponding author. kundu@iitk.ac.in Introduction In life testing experiments, often the data are censored. Among the different censoring schemes, Type-I and Type-II censoring schemes are the two most popular censoring schemes. In Type-I censoring scheme, the experimental time is fixed, but the number of failures is random, whereas Part of this work has been supported by grants from DST and CSIR, Government of India.
2 2 in Type-II censoring scheme, the experimental time is random but the number of failures is fixed. A mixture of Type-I and Type-II censoring scheme is known as the hybrid censoring scheme and it can be described as follows: Suppose a total of n units is placed on a life testing experiment, and the lifetimes of the sample units are independent and identically distributed i.i.d. random variables. Let the ordered lifetimes of these items are denoted by T :n,..., T n:n respectively. The test is terminated when r, a pre-chosen number, out of n items are failed, or when a pre-determined time T on test has been reached, i.e., the test is terminated at the time point T = min{t r:n, T. This particular censoring scheme was introduced by Epstein [9], and it is popularly known now as the Type-I hybrid censoring scheme. It has been used quite extensively in reliability acceptance test in MIL-STD-78C [2]. Like conventional Type-I censoring scheme, the main disadvantage of Type-I hybrid censoring scheme is that most of the inferential results are obtained under the condition that the number of observed failures is at least one, and moreover, there may be very few failures at the termination point of the experiment. In that case the efficiency of the estimators may be very low. Due to this reason, Childs et al. [5] introduced the Type-II hybrid censoring scheme, as an alternative to the Type-I hybrid censoring scheme, that would terminate the experiment at the time point T = max{t r:n, T. It has the advantage of guaranteeing that at least r failures are observed at the end of the experiment. Childs et al. [5] considered the exact distribution of the maximum likelihood estimator MLE of the mean of an one-parameter exponential distribution based on a Type-II hybrid censored samples. The exact confidence interval of the mean also has been obtained based on the exact distribution. Recently, Banerjee and Kundu [3] considered the statistical inference of the two-parameter Weibull distribution based on Type-II hybrid censored samples. Very recently, Childs et al. [6] obtained the exact distributions of the MLEs of the scale and location parameters of a two-parameter exponential distribution, when the data are Type-I hybrid censored. The purpose of this paper is to consider a Type-II hybrid censoring scheme when the lifetime distributions of the n experimental units are assumed to be i.i.d. two-parameter exponential distribution. It is assumed that the probability density function of the lifetime of an experimental unit,
3 3 for > 0, is ft; µ, = e t µ if t > µ, and 0 otherwise. First we obtain the MLEs of the unknown parameters of µ R and > 0, and provide the joint moment generating function. Based on the joint moment generating function MGF, we obtain the marginal MGFs, and the marginal distribution functions of the MLEs. From the marginal distribution functions, using the same idea as in Chen and Bhattacharya [4] the confidence interval of is obtained by solving non-linear equations. Since the confidence interval based on the MLEs is quite difficult to implement, we propose to use different bootstrap confidence intervals, whose implementation are quite straight forward. Rest of the paper is organized as follows. In Section 2, we present the maximum likelihood estimators MLEs of the unknown parameters, and the joint and marginal moment generating functions of the MLEs are derived in the same section. Different confidence intervals are proposed in Section 3. Simulation results and the analysis of a data set are presented in Section 4. Finally we conclude the paper in Section 5. Proof of the theorems are provided in the appendix. 2 Maximum Likelihood Estimators and Their Marginal Distributions Let N be the number of units failed before the time T, then the likelihood of the observed data is given by L µ, Data = n! e P N n N! N i= t i:n µ n NT µ if t r:n < T n! e P r n r! r i= t i:n µ n rtr:n µ if t r:n > T. 2 From 2, it is clear that MLE of does not exist when r = and N = 0. So for r =, the conditional MLEs, conditioning on the event {N is given by µ = T :n, { = N T i:n n T :n + n N T. N i=2 3
4 4 For r 2, MLEs of the unknown parameters exist for all values of N and they are given by µ = T :n, = { N T i:n n T :n + n N T N i=2 { r T i:n n T :n + n r + T r:n r i=2 if T r:n < T if T r:n > T. 4 Next we consider the distributions of the MLEs. When r =, MLE of does not exist if N = 0. So for r =, we consider the conditional distribution of conditioning on the event {N. For r 2, MLE of exists for all values of N. In order to find the distribution of µ and, we first find the joint moment generating function MGF of and µ, and then inverting it to get the distribution of the MLEs. Theorem 2.. For r =, the conditional joint MGF of, µ conditioning on the event {N exists for all < ω < and < ω 2 < and is given by E[e ω b +ω 2 bµ N ] = q n n + c ij i=2 + c nj c 0 e µ 0ω +µ ω 2 + ω λ 0 ν 0 d 0 e µ ijω +µ ω 2 e T ω 2 + ω λ 0 ν 0 + ω λ i αi + ω λ ij d ij e µ njω +µ ω 2 ω αn + ω j d nj e µ ω 2 q n e T ω 2 ω αn ν n e T ω 2 ω λ i αi + ω λ ij e T ω 2 ω αn + ω j.
5 5 For r 2, the joint MGF of and µ exists for ω < r and ω 2 < + n r E[e ω b +ω 2 bµ ] ω and is given by = q n e T ω 2 e ω α µ ω 2 q n e T ω 2 + λ ν n ω αn ν n n e µ ij ω +µ ω 2 + c ij i= ω αi λ i + ω λ ij n e T ω 2 d ij ν i= j ω αi λ i + ω λ ij e µ nj ω +µ ω 2 + c nj ω αn + ω when µ < T and for µ T where q = e T µ j E[e ω b +ω 2 bµ ] = n j T µ r µ ij = n j T µ i r α i = i e T ω 2 d nj ν j ω αn + ω e µ ω 2 ω α, λ ν n j for j = 0,, i, i =,, r for j = 0,, i, i = r,, n. for i =,, r for i = r,, n 5 6 ν i = i + r λ i = λ ij = i r j + n j i j + n j c ij = i j n i d ij = i j n i i q n j j + i q n j + for i = 0,, n for i =,, r for i = r,, n for j = 0,, i, i =,, r for j = 0,, i, i = r,, n for j = 0,, i, i =,, n for j = 0,, i, i =,, n 7
6 6 Remark. Note that as T, the joint MGF of and µ at ω, ω 2 reduces to e µ ω 2 ω αn ν n which is the joint MGF of and µ in case of complete sample. That means as T, 2n / is distributed as χ 2 random variable with 2n 2 degrees of freedom, n µ µ/ is a standard exponential random variable, and they are independently distributed. Remark 2. For r 2, if T is smaller than the left end point of the support µ, i.e., T µ, the joint MGF of and µ at ω, ω 2, as given in 5, becomes the joint MGF of and µ in case of ordinary Type-II censoring scheme. In this case 2r / is distributed as χ 2 random variable with 2r 2 degrees of freedom, n µ µ/ is a standard exponential random variable, and they are independently distributed. Remark 3. It is obvious from the expression of the joint MGF of and µ that the marginal distributions belong to the generalized mixtures of the well known distributions. From the MGF of, PDF of can be obtained by using the inversion technique as suggested by Chen and Bhattacharya [4], and it is as follows: Theorem 2.2. For r =, conditional PDF of conditioned on {N for < t < is f b t = g t; α n, + q n [c 0 g t + µ 0 ;, λ 0 d 0 g t;, λ 0 n + {c ij g 2 t µ ij ; α i, λ i, λ ij d ij g 2 t; α i, λ i, λ ij i=2 ] + {c nj g 2 t µ nj ; α n,, j d nj g 2 t; α n,, j. 8 For r 2, the PDF of for < t < is g 3 t f b t = g t; α, λ if µ < T if µ T 9
7 7 where, g t; α, λ = g 2 t; α, λ, λ 2 = α λ α Γ α e λt t α if t 0, 0 otherwise α p k g t; α k, λ + p k g t,, λ 2 for t R, p k = λ 2 λ k. λ +λ 2 λ +λ 2 g 3 t = q n g t; α, λ + q n g t; α n, n + {c ij g 2 t µ ij ; α i, λ i, λ ij d ij g 2 t; α i, λ i, λ ij i= + {c nj g 2 t µ nj ; α n,, j d nj g 2 t; α n,, j. Note that g 3 t depends on µ ij, α i, λ i, and λ ij ; i, j =,..., n, however, for brevity we do not write it explicitly. Since the integration of density function over the whole range is one, we have, from the Theorem 2.2, the following identity. n c ij d ij + c nj d nj = 0. i= Theorem 2.3. When < t <, for r =, conditional PDF of µ conditioning on {N is given by f bµ t = q n [c 0 g t µ;, ν 0 d 0 g t T ;, ν 0 + g t µ;, ν n n q n g t T ;, ν n + {c ij g t µ;, d ij g t T ;, i=2 ] + {c nj g t µ;, d nj g t T ;,,
8 8 and for r 2, PDF of µ is where f bµ t = g 4 t µ if µ < T g t µ;, ν n if µ T, n g 4 t = g t;, ν n + {c ij g t;, d ij g t + µ T ;, i= + {c nj g t;, d nj g t + µ T ;,. Note that g 4 t depends on ν i ; i =,..., n, however, for brevity we do not write it explicitly. From PDFs of, the corresponding moments can be easily obtained. The first two moments of are as follows. For r =, E[ ] = + A µ, + q n B µ,, E[ 2 ] = C µ, + D µ, + q n E µ,. For r 2, E[ ] = + A 2 µ, + B µ, if µ < T if µ T, r E[ 2 ] = C 2 µ, + D 2 µ, + E µ, if µ < T 2 if µ T, r
9 9 where A µ, = [ n q n i=2 j n c ij d ij { i 2 r c 0 d 0 n qn n p jk i k i ], i 2 p jk n j ij + B µ, = C µ, = n c ij µ ij, i= [ n i= j n c ij d ij D µ, = 2 q n [ n i= { i 2 +2n c 0 d 0 qn n +n c 0 µ 0 ], j n i k p jk i k i + 2 i 2 2 { i 2 c ij µ ij ] q n, p jk i k i p jk i 2 p jk n j i 2 j + 2 n j ij + E µ, = n c ij µ 2 ij, i= A 2 µ, = q n n r n r + c ij d ij C 2 µ, = + i= n i=r j n r c ij d ij i= + n i=r j n qn r qn, n c ij d ij { r 2 c ij d ij { r 2 { i 2 p jk r k r p jk i k i r 2 p jk i 2 p jk r k p jk r 2 r k + 2 r 2 p jk { i 2 p jk i k i 2 i k + 2 i 2 p jk n j rj + n j ij + n j 2 r 2 j + 2, n j 2 i 2 j + 2
10 0 { r r 2 r k D 2 µ, = 2 c ij µ ij p jk r i= +2 n i=r j n { i 2 c ij µ ij p jk i k i r 2 p jk i 2 n j rj + p jk n j ij +. The CDFs of and µ can be easily obtained from their respective PDFs. The CDFs are as follows. When < t <, the CDF of for r = is F b t = Γ t, α n + q n [ c 0 e λ 0 µ 0 t d 0 e λ 0 t and for r 2, the CDF is F b t = n + {c ij Gt; µ ij, α i, λ i, λ ij d ij G t; 0, α i, λ i, λ ij i=2 ] + {c nj Gt; µ nj, α n,.j d nj G t; 0, α n,, j, q n Γ λ t, α + q n Γ t, α n n + {c ij Gt; µ ij, α i, λ i, λ ij d ij G t; 0, α i, λ i, λ ij i=2 n 2 + {c nj Gt; µ nj, α n,.j d nj G t; 0, α n,, j if µ < T Γ λ t, α if µ T, where, x denotes the max{0, x, and 0 if t 0 Γ t; α = t Γα e x x α dx if t > 0, Gt; µ, α, λ, λ 2 = 0 α p k Γ λ 2 t; if t 0 α p k Γ λ t; α k if t 0. The CDF of µ for r = is F bµ t = 0 if t < µ q n [ e n t µ] if µ t < T if t T,
11 and for r 2, the same is F bµ t = { 0 if t < µ e n t µ if t µ. Note that the p th quantile of the exponential distribution given in is η p = µ + a p, where a p = ln p. Therefore the MLE of η p can be found by replacing µ and by their respective MLEs. The moment generating function of the MLE of η p, say η p can be obtained easily from the joint MGF of µ and in Theorem Confidence Intervals: Note that the distribution of µ is same as that of the lowest order statistic and this has been extensively studied in the literature, hence we do not pursue it here. In this section, we present different methods for construction of confidence intervals CIs for the unknown parameter. From the Theorem 2.2, we can find the approximate CI for. However, as PDF of is quite complicated, we also present the bootstrap CI for the same parameter. 3. Approximate Confidence interval: From CDF of, approximate CI can be found, based on the assumption that CDF of is a strictly decreasing function of. Several authors including Chen and Bhattacharyya [4], Gupta and Kundu [0], Kundu and Basu [], Childs et al. [5], and Balakrishnan et al. [2] used this method to find the CI of. Though it is difficult to verify the assumption, an extensive numerical study supports the monotonicity assumption. Suppose obs is the estimate of. Then a two-sided 00 α% approximate CI, say L, U, for can be constructed by solving the equations F L obs = α 2 and F U obs = α 2 for L the lower bound of and U the upper bound of replacing µ by its MLE. However, they are nonlinear equations, and we need to solve them by some numerical procedure, e.g., bisection method or Newton-Raphson method.
12 2 3.2 Bootstrap Confidence Interval: The exact CIs presented in the previous section are computationally quite complicated, specially when sample size is large. So we consider the bootstrap CIs. Here we consider two types of bootstrap CIs, viz., percentile bootstrap CI and Bias adjusted percentile BCa interval; see Efron and Tibshirani [8] for details Bootstrap sample: Step. Given T, r, n, and the original Type-II sample, µ and are obtained from 3 or 4. Step 2. Based on T, r, n, µ, and, a random sample of size n from Uniform0, distribution is generated, and order them to get U :n,..., U n:n. Step 3. Let t i:n = µ log U i:n Step 4. If t r:n < T, then find N such that t N :n < T t N +:n, and set N = Now, {t :n,..., t N :n is the bootstrap sample. { N if t r:n < T r if t r:n T. Step 5. Based on n, T, r, and the bootstrap sample, µ and are obtained form 3 or 4. Step 6. Step -5 are repeated B times, and s are ordered in ascending order to obtain the bootstrap sample { [], [2],..., [B] Percentile bootstrap CI: A two-sided 00 α% bootstrap confidence interval for is [ α 2 B], [ α 2 B], where, [x] denotes the largest integer less than or equal to x.
13 Bias adjusted percentile BCa interval: A two-sided 00 α% BCa bootstrap confidence interval for is [α B], [α 2 B], where, { ẑ 0 + z α/2 α = Φ ẑ 0 + aẑ 0 + z α/2 and { ẑ 0 + z α/2 α 2 = Φ ẑ 0 +. aẑ 0 + z α/2 Here Φ is the CDF of the standard normal distribution, z α is the upper α-point of standard normal distribution and ẑ 0 = Φ { # of [j] < B, j =,..., B. A estimate of the acceleration a is â = N i= { N 6 i= [ i] 3 [ i] 2 3/2, where i is the MLE of based on the original sample with the i th observation deleted, and = N N i= i. 4 Simulation Results and Data Analysis In this section the results of Monte Carlo simulation is presented to study the performance of the inference procedures described in the Sections 3. We choose the value of the location parameter µ to be zero and different values for the scale parameter, viz.,.00, 2.00, 3.00, 4.00 and We also take n = 20, r = 6 and different choices for T, viz.,.5, 2.5, 3.5. The coverage percentage of different confidence interval discussed in the section 3 are calculated based on the 0, 000 Monte Carlo simulations and B = 000. These values are presented in the Tables 2, 3 and 4. It is clear from the Table 2 that the approximate method of constructing confidence interval is always maintaining its coverage percentage to its pre-fixed nominal level. Among the bootstrap
14 4 methods for constructing confidence interval, adjusted percentile bootstrap method is better then percentile bootstrap method with respect to the coverage percentage. From the Table 3, we can see that the coverage percentage of the percentile bootstrap method is quite lower than its pre-fixed nominal level, while the same of the adjusted bootstrap method is somewhat close to its nominal level see Table 4. Though approximate method of construction of confidence interval is always better with respect to coverage percentage and average length, it is very complicated to calculate specially when n is large. So we suggest to use the BCa bootstrap CI for the large values of n. Next, we consider a data set to illustrate the procedures described in the previous sections. Here we consider the data provided in Bain []. Suppose a sample of 20 items are put on the test and the test is terminated after 50 hours. There are 3 failures within this time and the failures times are 3, 9, 23, 26, 37, 38, 4, 45,58, 84, 90, 09, and 38. To make an artificial hybrid Type-II data from this data set one can take any r less than or equal to 3. Here we take r = 2. Under the assumption that lifetime has two-parameters exponential distribution, MLEs of µ and are 3 and respectively. Different types of confidence intervals are reported in the Table. 5 Conclusion In this paper, we have considered the Hybrid Type-II censoring scheme, when life times have twoparameters exponential distribution, under frequentist approach. We have found the MLEs for both the parameters. We have considered different methods for construction of confidence interval. We have seen that the approximate and BCa bootstrap methods of construction of confidence interval are quite good. So we recommend to use BCa bootstrap method of construction of confidence interval specially when n is large. Acknowledgements The authors would like to thank the referees and the editor-in-chief Professor N. Balakrishnan for many valuable suggestions which had helped us to improve the previous version of the manuscript.
15 5 A Appendix We need the the following lemmas to prove the theorems. Lemma A.. Let T :n <... < T n:n be the order statistics of a random sample of size n from a continuous distribution with PDF fx and the corresponding CDF is F. Let T be a pre-fixed number such that F T is positive and N denote the number of order statistics less than or equal to T. The conditional joint PDF of T :n,..., T N:n conditioned on the event N = i is given by n! i f t, t 2,..., t i N = i = f t j { F T n i n i!p [N = i] j= where F is the distribution function of f. if t <... < t i < T, Proof: See [5]. Lemma A.2. Let T :n <... < T n:n be the order statistics of a random sample of size n from a continuous distribution with PDF fx and the corresponding CDF is F. Let T be a pre-fixed number such that F T is positive and N denote the number of order statistics less than or equal to T. Let r {, 2,..., n be a pre-fixed integer. Then PDF of T :n, T 2:n,..., T r:n conditioned on the event N = 0 is given by f t, t 2,..., t r N = 0 = n! n r!p [N = 0] r f t j { F t r n r j= if T < t <... < t r <. For i =, 2,..., r, PDF of T :n, T 2:n,..., T r:n conditioned on the event N = i is given by n! r f t, t 2,..., t r N = i = f t j { F t r n r n r!p [N = i] Proof: See [5]. j= if t <... < t i < T < t i+ <... < t r <. Lemma A.3. Let X be a Gammaα, random variable, Y be a standard Exponential random variable, and they are independently distributed. The PDF of X is Γα xα e x for x > 0 f x; α, λ = 0 otherwise.
16 6 The PDF of Y is e y for y > 0 f 2 y = 0 otherwise. Then for any arbitrary constant A, λ and λ 2 the MGF of A + λ X + λ 2 Y is given by M A+λ X+λ 2 Y ω = e ωa λ ω α λ 2 ω. This MGF exists if ω { λ 2, λ if λ 2 < 0, min{ λ, λ2 if λ 2 0. Proof: This lemma can be proved using the joint distribution of X, Y, and therefore the proof is omitted. Lemma A.4. Let X be a Gammaα, λ with α integer and Y be a Exponentialλ 2 random variable and they be independently distributed. Then the PDF of X Y is given by α α g 2 t; α, λ, λ 2 = p k g t; α k, λ + p k g t,, λ 2 for t R, where p k = λ 2 λ +λ 2 λ λ +λ 2 k. Proof: See [6]. Proof of Theorem 2.: Note that the number of unit failed before time T, N, is a non-negative random variable with the probability mass function PMF P [N = i] = n i e T µ i e n i T µ if i = 0,,..., n. 0 Case-I: r = The joint conditional MGF of, µ conditioned on the event {N can be written as E[e ω b +ω 2 bµ N ] = n E[e ω +ω b 2 bµ N = i] P [N = i N ]. i=
17 7 Using 3 and Lemma A., E[e ω b +ω 2 bµ N = ] = = e T µ e n T ω + µ T µ e ω n T t+ω 2 t t µ dt {e n µ ω +µ ω 2 µ e T n T ω +T ω 2 e T µ + n ω. Using the above expression and 0 E[e ω +ω b 2 bµ N = ] P [N = N ] [ = q n e µ 0 ω +µ ω 2 ] c 0 + n ω d e T ω 2 0, 2 + n ω where q, c 0, d 0 and µ 0 are as in 6 and 7. For i = 2, 3,..., n, using 3 and Lemma A., we have E[e ω b +ω 2 bµ N = i] = i! i i e T µ T T µ t... T t i e ω i { P i j=2 t j n t +n it+ω 2 t P i j= t j µ dt i... dt = T µ +T ω 2 T i! T e i i... i e T µ µ t = e i T µ +T ω 2 e T µ i T t i e + n i ω ω 2t T ω i P i j=2 t j T dt i... dt i i j j + i ω e j+ + n j Using the above expression and 0, one will get for i = 2,, n ω ω 2T µ i. + n j ω ij+ ω j+ 2 E[e ω b +ω 2 bµ N = i] P [N = i N ] T µ = e n +T ω 2 T µ n e n i j i i j + i ω e j+ + n j ω i ω 2T µ i + n j ij+ ω j+ ω 2.
18 8 So using the above expression and 2 in, we have E[e ωb N ] = q n n + c ij i=2 + c nj c 0 e µ 0ω +µ ω 2 + ω λ 0 ν 0 d 0 e µ ijω +µ ω 2 e T ω 2 + ω λ 0 ν 0 + ω λ i αi + ω λ ij d ij e µ njω +µ ω 2 ω αn + ω j d nj e µ ω 2 q n e T ω 2 ω αn ν n e T ω 2 ω λ i αi + ω λ ij e T ω 2 ω αn + ω j. where q, α i, λ i, µ ij, λ ij, c ij and d ij are as in 6 and 7. Now using Lemmas A.3 and A.4, one will get 8. Case-II: r 2 In this case, the joint MGF of and µ can be expressed as E[e ω b +ω 2 bµ ] = n E[e ω +ω b 2 bµ N = i] P [N = i]. 3 i=0 Now using 4 and Lemma A.2, for ω < r and ω 2 < + n r ω = E[e ω +ω b 2 bµ N = 0] n! r n r!p [N = 0] max{t,µ t... t r e ω r { P r j=2 t j n t +n r+t r+ω 2 t P r j= tj µ n r tr µ dtr... dt = e max{t,µ µ n +ω 2 max{t,µ P [N = 0] r ω r n ω 2. Hence for ω < r and ω 2 < + n r ω E[e ω +ω b 2 bµ max{t,µ µ n +ω N = 0] P [N = 0] = e 2 max{t,µ r ω r n ω 2. 4
19 9 Now for i =, 2,..., r and ω < r, using 4 and Lemma A.2 E[e ω b +ω 2 bµ N = i] = = n! n r!p [N = i] r T µ T... t i T... nµ n! e n r!p [N = i] r { T T T... µ t {... T t i+ it n i r T ω +ω 2 T t r e ω r { P r j=2 t j n t +n r+t r+ω 2 t P r j= tj µ t i e + n r ω ω 2t T ω r P i j=2 t j T dt i... dt t r e ω r P r j=i+ t j n r+ ω r t r dt r... dt i+ n r tr µ dtr... dt T µ = e n +T ω 2 P [N = i] n i j i Hence for i =, 2,..., r and ω < r, E[e ω b +ω 2 bµ N = i] P [N = i] T µ n +T ω = e 2 i j n i i j + i j + Now for i = r, r +,..., n, using 4 and Lemma A. E[e ω b +ω 2 bµ N = i] r ω r ω e j+ + n j ω r ω 2T µ r. + n j ω rj+ ω j+ 2 e j+ + n j ω r ω 2T µ r. 5 + n j ω rj+ ω j+ 2 = n! n i! i P [N = i] T T T µ t... T µ = e n +T ω 2 P [N = i] t i e ω i { P i j=2 t j n t +n it+ω 2 t P i j= n i j i Hence for i = r, r +,..., n, E[e ω b +ω 2 bµ N = i] P [N = i] T µ n +T ω = e 2 i j n i i j + i j + i ω i ω tj µ n i T µ dti... dt e j+ + n j ω i ω 2T µ i + n j ij+ ω j+ ω 2 e j+ + n j ω i ω 2T µ i + n j ij+ ω j+ ω 2..
20 20 So using the above expression, 4 and 5 in 3, we have for ω < r and ω 2 < n E[e ω b +ω 2 bµ ] = q n e T ω 2 e ω α µ ω 2 q n e T ω 2 + λ ν n ω αn ν n n e µ ij ω +µ ω 2 + c ij i= ω αi λ i + ω e µ nj ω +µ ω 2 + c nj ω αn + ω λ ij j n e T ω 2 d ij ν i= j ω αi λ i + ω e T ω 2 d nj ν j ω αn + ω j λ ij where q, µ ij, α i, λ i, λ ij, c ij and d ij are as in 6 and 7. Now using Lemmas A.3, A.4, one will get 9., Table : Lower and upper limits of different confidence intervals for Approx. CI Per Boot. CI BCa Boot. CI Level LL UL LL UL LL UL 90% % % Table 2: Coverage percentage and average length of approximate confidence intervals based on 0000 simulations with µ = 0, n = 20 and r = 6. 90% C.I. 95% C.I. 99% C.I. T Cov. Per. Ave. Len. Cov. Per. Ave. Len. Cov. Per. Ave. Len
21 2 Table 3: Coverage percentage and average length of bootstrap confidence intervals based on 0000 simulations with µ = 0, n = 20, r = 6 and B = % C.I. 95% C.I. 99% C.I. T Cov. Per. Ave. Len. Cov. Per. Ave. Len. Cov. Per. Ave. Len Table 4: Coverage percentage and average length of BCa bootstrap confidence intervals based on 0000 simulations with µ = 0, n = 20, r = 6 and B = % C.I. 95% C.I. 99% C.I. T Cov. Per. Ave. Len. Cov. Per. Ave. Len. Cov. Per. Ave. Len
22 22 References [] Bain, L. J. 978, Statistical analysis of reliability and life testing models, Marcel and Dekker Inc., New York. [2] Balakrishnan, N., Xie, Q. and Kundu, D. 2009, Exact inference for a simple step-stress model from the exponential distribution under time constraint, Annals of the Institute of Statistical Mathematics, vol. 6, [3] Banerjee, A. and Kundu, D. 2008, Inference based on Type-II hybrid censored data from a Weibull distribution, IEEE Transaction on Reliability, vol. 57, [4] Chen, S. and Bhattacharya, G.K. 988, Exact confidence bounds for an exponential parameter under hybrid censoring, Communications in Statistics - Theory and Methods, vol. 6, [5] Childs, A., Chandrasekar, B., Balakrishnan, N. and Kundu, D. 2003, Exact likelihood inference based on Type-I and Type-II hybrid censored samples from the exponential distribution, Annals of the Institute of Statistical Mathematics, vol. 55, [6] Childs, A., Balakrishnan, N. and Chandrasekar, B. 20, Exact distribution of the MLEs of the parameters and of the quantiles of two-parameter exponential distribution under hybrid censoring, Statistics, DOI: 0.080/ [7] Draper, N. and Guttman, I. 987, Bayesian analysis of hybrid life tests with exponential failure times, Annals of the Institute of Statistical Mathematics, vol. 39, [8] Efron, B. and Tibshirani 998, An introduction to the bootstrap, Boca Ration: Chapman and Hall/CRC CRC Press. [9] Epstein, B. 954, Truncated life-test in exponential case, Annals of Mathematical Statistics, vol. 25, [0] Gupta R. D. and Kundu D. 998, Hybrid censoring scheme with exponential failure distribution, Communication in Statistics-Theory and Methods, vol. 27,
23 23 [] Kundu D. and Basu A. 2000, Analysis of incomplete data in presence of competing risks, Journal of Statistical Planing and Inference, vol. 87, [2] MIL-STD-78-C977, Reliability design qualification and production acceptance tests: exponential distribution, Washington, D.C., U.S. Government Printing Office.
Analysis of Type-II Progressively Hybrid Censored Data
Analysis of Type-II Progressively Hybrid Censored Data Debasis Kundu & Avijit Joarder Abstract The mixture of Type-I and Type-II censoring schemes, called the hybrid censoring scheme is quite common in
More informationBayesian Analysis of Simple Step-stress Model under Weibull Lifetimes
Bayesian Analysis of Simple Step-stress Model under Weibull Lifetimes A. Ganguly 1, D. Kundu 2,3, S. Mitra 2 Abstract Step-stress model is becoming quite popular in recent times for analyzing lifetime
More informationOn Simple Step-stress Model for Two-Parameter Exponential Distribution
On Simple Step-stress Model for Two-Parameter Exponential Distribution S. Mitra, A. Ganguly, D. Samanta, D. Kundu,2 Abstract In this paper, we consider the simple step-stress model for a two-parameter
More informationSTEP STRESS TESTS AND SOME EXACT INFERENTIAL RESULTS N. BALAKRISHNAN. McMaster University Hamilton, Ontario, Canada. p.
p. 1/6 STEP STRESS TESTS AND SOME EXACT INFERENTIAL RESULTS N. BALAKRISHNAN bala@mcmaster.ca McMaster University Hamilton, Ontario, Canada p. 2/6 In collaboration with Debasis Kundu, IIT, Kapur, India
More informationHybrid Censoring; An Introduction 2
Hybrid Censoring; An Introduction 2 Debasis Kundu Department of Mathematics & Statistics Indian Institute of Technology Kanpur 23-rd November, 2010 2 This is a joint work with N. Balakrishnan Debasis Kundu
More informationA New Two Sample Type-II Progressive Censoring Scheme
A New Two Sample Type-II Progressive Censoring Scheme arxiv:609.05805v [stat.me] 9 Sep 206 Shuvashree Mondal, Debasis Kundu Abstract Progressive censoring scheme has received considerable attention in
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationPoint and Interval Estimation of Weibull Parameters Based on Joint Progressively Censored Data
Point and Interval Estimation of Weibull Parameters Based on Joint Progressively Censored Data Shuvashree Mondal and Debasis Kundu Abstract Analysis of progressively censored data has received considerable
More informationinferences on stress-strength reliability from lindley distributions
inferences on stress-strength reliability from lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter R = P (Y
More informationOptimum Hybrid Censoring Scheme using Cost Function Approach
Optimum Hybrid Censoring Scheme using Cost Function Approach Ritwik Bhattacharya 1, Biswabrata Pradhan 1, Anup Dewanji 2 1 SQC and OR Unit, Indian Statistical Institute, 203, B. T. Road, Kolkata, PIN-
More informationAnalysis of incomplete data in presence of competing risks
Journal of Statistical Planning and Inference 87 (2000) 221 239 www.elsevier.com/locate/jspi Analysis of incomplete data in presence of competing risks Debasis Kundu a;, Sankarshan Basu b a Department
More informationHybrid Censoring Scheme: An Introduction
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline 1 2 3 4 5 Outline 1 2 3 4 5 What is? Lifetime data analysis is used to analyze data in which the time
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by:[indian Institute of Technology] On: 10 December 2007 Access Details: [subscription number 768117315] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales
More informationPoint and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples
90 IEEE TRANSACTIONS ON RELIABILITY, VOL. 52, NO. 1, MARCH 2003 Point and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples N. Balakrishnan, N. Kannan, C. T.
More informationThe comparative studies on reliability for Rayleigh models
Journal of the Korean Data & Information Science Society 018, 9, 533 545 http://dx.doi.org/10.7465/jkdi.018.9..533 한국데이터정보과학회지 The comparative studies on reliability for Rayleigh models Ji Eun Oh 1 Joong
More informationReliability analysis under constant-stress partially accelerated life tests using hybrid censored data from Weibull distribution
Hacettepe Journal of Mathematics and Statistics Volume 45 (1) (2016), 181 193 Reliability analysis under constant-stress partially accelerated life tests using hybrid censored data from Weibull distribution
More informationBayes Estimation and Prediction of the Two-Parameter Gamma Distribution
Bayes Estimation and Prediction of the Two-Parameter Gamma Distribution Biswabrata Pradhan & Debasis Kundu Abstract In this article the Bayes estimates of two-parameter gamma distribution is considered.
More informationOn Weighted Exponential Distribution and its Length Biased Version
On Weighted Exponential Distribution and its Length Biased Version Suchismita Das 1 and Debasis Kundu 2 Abstract In this paper we consider the weighted exponential distribution proposed by Gupta and Kundu
More informationExponentiated Rayleigh Distribution: A Bayes Study Using MCMC Approach Based on Unified Hybrid Censored Data
Exponentiated Rayleigh Distribution: A Bayes Study Using MCMC Approach Based on Unified Hybrid Censored Data ABSTRACT M. G. M. Ghazal 1, H. M. Hasaballah 2 1 Mathematics Department, Faculty of Science,
More informationOn Generalized Progressive Hybrid Censoring in presence of competing risks
On Generalized Progressive Hybrid Censoring in presence of competing risks Arnab Koley & Debasis Kundu Abstract The progressive Type-II hybrid censoring scheme introduced by Kundu and Joarder (Computational
More informationSTATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University
STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment of the Requirements for the Degree
More informationEstimation in an Exponentiated Half Logistic Distribution under Progressively Type-II Censoring
Communications of the Korean Statistical Society 2011, Vol. 18, No. 5, 657 666 DOI: http://dx.doi.org/10.5351/ckss.2011.18.5.657 Estimation in an Exponentiated Half Logistic Distribution under Progressively
More informationINFERENCE FOR A GAMMA STEP-STRESS MODEL UNDER CENSORING
INFERENCE FOR A GAMMA STEP-STRESS MODEL UNDER CENSORING INFERENCE FOR A GAMMA STEP-STRESS MODEL UNDER CENSORING BY LAILA ALKHALFAN, M.Sc. McMaster University a thesis submitted to the school of graduate
More informationBayesian Inference and Life Testing Plan for the Weibull Distribution in Presence of Progressive Censoring
Bayesian Inference and Life Testing Plan for the Weibull Distribution in Presence of Progressive Censoring Debasis KUNDU Department of Mathematics and Statistics Indian Institute of Technology Kanpur Pin
More informationOn Sarhan-Balakrishnan Bivariate Distribution
J. Stat. Appl. Pro. 1, No. 3, 163-17 (212) 163 Journal of Statistics Applications & Probability An International Journal c 212 NSP On Sarhan-Balakrishnan Bivariate Distribution D. Kundu 1, A. Sarhan 2
More informationOn the Comparison of Fisher Information of the Weibull and GE Distributions
On the Comparison of Fisher Information of the Weibull and GE Distributions Rameshwar D. Gupta Debasis Kundu Abstract In this paper we consider the Fisher information matrices of the generalized exponential
More informationConditional independence of blocked ordered data
Conditional independence of blocked ordered data G. Iliopoulos 1 and N. Balakrishnan 2 Abstract In this paper, we prove that blocks of ordered data formed by some conditioning events are mutually independent.
More informationBAYESIAN ESTIMATION OF THE EXPONENTI- ATED GAMMA PARAMETER AND RELIABILITY FUNCTION UNDER ASYMMETRIC LOSS FUNC- TION
REVSTAT Statistical Journal Volume 9, Number 3, November 211, 247 26 BAYESIAN ESTIMATION OF THE EXPONENTI- ATED GAMMA PARAMETER AND RELIABILITY FUNCTION UNDER ASYMMETRIC LOSS FUNC- TION Authors: Sanjay
More information(Received August 7, 2001; revised March 18, 2002)
Ann. Inst. Statist. Math. Vol. 55, No. 2, 319-330 (2003) (~)2003 The Institute of Statistical Mathematics EXACT LIKELIHOOD INFERENCE BASED ON TYPE-I AND TYPE-II HYBRID CENSORED SAMPLES FROM THE EXPONENTIAL
More informationPrediction for future failures in Weibull Distribution under hybrid censoring
Prediction for future failures in Weibull Distribution under hybrid censoring A. Asgharzadeh a, R. Valiollahi b, D. Kundu c a Department of Statistics, School of Mathematical Sciences, University of Mazandaran,
More informationConstant Stress Partially Accelerated Life Test Design for Inverted Weibull Distribution with Type-I Censoring
Algorithms Research 013, (): 43-49 DOI: 10.593/j.algorithms.01300.0 Constant Stress Partially Accelerated Life Test Design for Mustafa Kamal *, Shazia Zarrin, Arif-Ul-Islam Department of Statistics & Operations
More informationEstimation of P (X > Y ) for Weibull distribution based on hybrid censored samples
Estimation of P (X > Y ) for Weibull distribution based on hybrid censored samples A. Asgharzadeh a, M. Kazemi a, D. Kundu b a Department of Statistics, Faculty of Mathematical Sciences, University of
More informationBayesian Analysis for Partially Complete Time and Type of Failure Data
Bayesian Analysis for Partially Complete Time and Type of Failure Data Debasis Kundu Abstract In this paper we consider the Bayesian analysis of competing risks data, when the data are partially complete
More informationAn Analysis of Record Statistics based on an Exponentiated Gumbel Model
Communications for Statistical Applications and Methods 2013, Vol. 20, No. 5, 405 416 DOI: http://dx.doi.org/10.5351/csam.2013.20.5.405 An Analysis of Record Statistics based on an Exponentiated Gumbel
More informationExact Linear Likelihood Inference for Laplace
Exact Linear Likelihood Inference for Laplace Prof. N. Balakrishnan McMaster University, Hamilton, Canada bala@mcmaster.ca p. 1/52 Pierre-Simon Laplace 1749 1827 p. 2/52 Laplace s Biography Born: On March
More informationInferences on stress-strength reliability from weighted Lindley distributions
Inferences on stress-strength reliability from weighted Lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter
More informationWeighted Marshall-Olkin Bivariate Exponential Distribution
Weighted Marshall-Olkin Bivariate Exponential Distribution Ahad Jamalizadeh & Debasis Kundu Abstract Recently Gupta and Kundu [9] introduced a new class of weighted exponential distributions, and it can
More informationEstimation for Mean and Standard Deviation of Normal Distribution under Type II Censoring
Communications for Statistical Applications and Methods 2014, Vol. 21, No. 6, 529 538 DOI: http://dx.doi.org/10.5351/csam.2014.21.6.529 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation for Mean
More informationEstimation of Stress-Strength Reliability Using Record Ranked Set Sampling Scheme from the Exponential Distribution
Filomat 9:5 015, 1149 116 DOI 10.98/FIL1505149S Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Estimation of Stress-Strength eliability
More informationParameter Estimation of Power Lomax Distribution Based on Type-II Progressively Hybrid Censoring Scheme
Applied Mathematical Sciences, Vol. 12, 2018, no. 18, 879-891 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.8691 Parameter Estimation of Power Lomax Distribution Based on Type-II Progressively
More informationBurr Type X Distribution: Revisited
Burr Type X Distribution: Revisited Mohammad Z. Raqab 1 Debasis Kundu Abstract In this paper, we consider the two-parameter Burr-Type X distribution. We observe several interesting properties of this distribution.
More informationEstimation of Parameters of the Weibull Distribution Based on Progressively Censored Data
International Mathematical Forum, 2, 2007, no. 41, 2031-2043 Estimation of Parameters of the Weibull Distribution Based on Progressively Censored Data K. S. Sultan 1 Department of Statistics Operations
More informationHacettepe Journal of Mathematics and Statistics Volume 45 (5) (2016), Abstract
Hacettepe Journal of Mathematics and Statistics Volume 45 (5) (2016), 1605 1620 Comparing of some estimation methods for parameters of the Marshall-Olkin generalized exponential distribution under progressive
More informationSome Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution
Journal of Probability and Statistical Science 14(), 11-4, Aug 016 Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution Teerawat Simmachan
More informationConfidence Intervals in Ridge Regression using Jackknife and Bootstrap Methods
Chapter 4 Confidence Intervals in Ridge Regression using Jackknife and Bootstrap Methods 4.1 Introduction It is now explicable that ridge regression estimator (here we take ordinary ridge estimator (ORE)
More informationESTIMATION AND PREDICTION FOR TYPE-II HYBRID CENSORED DATA FOLLOW FLEXIBLE WEIBULL DISTRIBUTION
STATISTICA, anno LXXVII, n. 4, 2017 ESTIMATION AND PREDICTION FOR TYPE-II HYBRID CENSORED DATA FOLLOW FLEXIBLE WEIBULL DISTRIBUTION Vikas Kumar Sharma 1 Department of Mathematics, Institute of Infrastructure
More informationAnalysis of Type-II Hybrid Censored Competing Risks Data
Analysis of Type-II Hybrid Censored Competing Risks Data Arnab Koley, Debasis Kundu,2 & Ayon Ganguly 3 Abstract Kundu and Gupta [22] provided the analysis of Type-I hybrid censored competing risks data,
More informationSTAT 6385 Survey of Nonparametric Statistics. Order Statistics, EDF and Censoring
STAT 6385 Survey of Nonparametric Statistics Order Statistics, EDF and Censoring Quantile Function A quantile (or a percentile) of a distribution is that value of X such that a specific percentage of the
More informationOptimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests
International Journal of Performability Engineering, Vol., No., January 24, pp.3-4. RAMS Consultants Printed in India Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests N. CHANDRA *, MASHROOR
More informationStatistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS040) p.4828 Statistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions
More informationON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION
ON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION ZHENLINYANGandRONNIET.C.LEE Department of Statistics and Applied Probability, National University of Singapore, 3 Science Drive 2, Singapore
More informationINFERENCE FOR BIRNBAUM-SAUNDERS, LAPLACE AND SOME RELATED DISTRIBUTIONS UNDER CENSORED DATA
INFERENCE FOR BIRNBAUM-SAUNDERS, LAPLACE AND SOME RELATED DISTRIBUTIONS UNDER CENSORED DATA INFERENCE FOR BIRNBAUM-SAUNDERS, LAPLACE AND SOME RELATED DISTRIBUTIONS UNDER CENSORED DATA By Xiaojun Zhu A
More informationContinuous Distributions
A normal distribution and other density functions involving exponential forms play the most important role in probability and statistics. They are related in a certain way, as summarized in a diagram later
More informationDiscriminating Between the Bivariate Generalized Exponential and Bivariate Weibull Distributions
Discriminating Between the Bivariate Generalized Exponential and Bivariate Weibull Distributions Arabin Kumar Dey & Debasis Kundu Abstract Recently Kundu and Gupta ( Bivariate generalized exponential distribution,
More informationAnalysis of Middle Censored Data with Exponential Lifetime Distributions
Analysis of Middle Censored Data with Exponential Lifetime Distributions Srikanth K. Iyer S. Rao Jammalamadaka Debasis Kundu Abstract Recently Jammalamadaka and Mangalam (2003) introduced a general censoring
More informationDistributions of Functions of Random Variables. 5.1 Functions of One Random Variable
Distributions of Functions of Random Variables 5.1 Functions of One Random Variable 5.2 Transformations of Two Random Variables 5.3 Several Random Variables 5.4 The Moment-Generating Function Technique
More informationOn the Hazard Function of Birnbaum Saunders Distribution and Associated Inference
On the Hazard Function of Birnbaum Saunders Distribution and Associated Inference Debasis Kundu, Nandini Kannan & N. Balakrishnan Abstract In this paper, we discuss the shape of the hazard function of
More informationEstimation for inverse Gaussian Distribution Under First-failure Progressive Hybird Censored Samples
Filomat 31:18 (217), 5743 5752 https://doi.org/1.2298/fil1718743j Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Estimation for
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationEstimation of Stress-Strength Reliability for Kumaraswamy Exponential Distribution Based on Upper Record Values
International Journal of Contemporary Mathematical Sciences Vol. 12, 2017, no. 2, 59-71 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijcms.2017.7210 Estimation of Stress-Strength Reliability for
More informationSTAT 512 sp 2018 Summary Sheet
STAT 5 sp 08 Summary Sheet Karl B. Gregory Spring 08. Transformations of a random variable Let X be a rv with support X and let g be a function mapping X to Y with inverse mapping g (A = {x X : g(x A}
More informationA Reliability Sampling Plan to ensure Percentiles through Weibull Poisson Distribution
Volume 117 No. 13 2017, 155-163 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu A Reliability Sampling Plan to ensure Percentiles through Weibull
More information3 Continuous Random Variables
Jinguo Lian Math437 Notes January 15, 016 3 Continuous Random Variables Remember that discrete random variables can take only a countable number of possible values. On the other hand, a continuous random
More informationAn Extended Weighted Exponential Distribution
Journal of Modern Applied Statistical Methods Volume 16 Issue 1 Article 17 5-1-017 An Extended Weighted Exponential Distribution Abbas Mahdavi Department of Statistics, Faculty of Mathematical Sciences,
More informationA Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints
Noname manuscript No. (will be inserted by the editor) A Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints Mai Zhou Yifan Yang Received: date / Accepted: date Abstract In this note
More informationProblem 1 (20) Log-normal. f(x) Cauchy
ORF 245. Rigollet Date: 11/21/2008 Problem 1 (20) f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.0 0.2 0.4 0.6 0.8 4 2 0 2 4 Normal (with mean -1) 4 2 0 2 4 Negative-exponential x x f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.5
More informationBayesian Life Test Planning for the Weibull Distribution with Given Shape Parameter
Statistics Preprints Statistics 10-8-2002 Bayesian Life Test Planning for the Weibull Distribution with Given Shape Parameter Yao Zhang Iowa State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu
More informationStatistics 3858 : Maximum Likelihood Estimators
Statistics 3858 : Maximum Likelihood Estimators 1 Method of Maximum Likelihood In this method we construct the so called likelihood function, that is L(θ) = L(θ; X 1, X 2,..., X n ) = f n (X 1, X 2,...,
More informationAnalysis of Progressive Type-II Censoring. in the Weibull Model for Competing Risks Data. with Binomial Removals
Applied Mathematical Sciences, Vol. 5, 2011, no. 22, 1073-1087 Analysis of Progressive Type-II Censoring in the Weibull Model for Competing Risks Data with Binomial Removals Reza Hashemi and Leila Amiri
More informationCOMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH *
Jordan Journal of Mathematics and Statistics JJMS 4,, pp. 6-78 COMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH * Sanku Dey ABSTRACT:
More informationINVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION
Pak. J. Statist. 2017 Vol. 33(1), 37-61 INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION A. M. Abd AL-Fattah, A.A. EL-Helbawy G.R. AL-Dayian Statistics Department, Faculty of Commerce, AL-Azhar
More informationTHE MODIFIED EXPONENTIAL DISTRIBUTION WITH APPLICATIONS. Department of Statistics, Malayer University, Iran
Pak. J. Statist. 2017 Vol. 33(5), 383-398 THE MODIFIED EXPONENTIAL DISTRIBUTION WITH APPLICATIONS Mahdi Rasekhi 1, Morad Alizadeh 2, Emrah Altun 3, G.G. Hamedani 4 Ahmed Z. Afify 5 and Munir Ahmad 6 1
More informationReview Quiz. 1. Prove that in a one-dimensional canonical exponential family, the complete and sufficient statistic achieves the
Review Quiz 1. Prove that in a one-dimensional canonical exponential family, the complete and sufficient statistic achieves the Cramér Rao lower bound (CRLB). That is, if where { } and are scalars, then
More informationSupporting Information for Estimating restricted mean. treatment effects with stacked survival models
Supporting Information for Estimating restricted mean treatment effects with stacked survival models Andrew Wey, David Vock, John Connett, and Kyle Rudser Section 1 presents several extensions to the simulation
More informationInference on reliability in two-parameter exponential stress strength model
Metrika DOI 10.1007/s00184-006-0074-7 Inference on reliability in two-parameter exponential stress strength model K. Krishnamoorthy Shubhabrata Mukherjee Huizhen Guo Received: 19 January 2005 Springer-Verlag
More informationEstimation for generalized half logistic distribution based on records
Journal of the Korean Data & Information Science Society 202, 236, 249 257 http://dx.doi.org/0.7465/jkdi.202.23.6.249 한국데이터정보과학회지 Estimation for generalized half logistic distribution based on records
More informationDistribution Theory. Comparison Between Two Quantiles: The Normal and Exponential Cases
Communications in Statistics Simulation and Computation, 34: 43 5, 005 Copyright Taylor & Francis, Inc. ISSN: 0361-0918 print/153-4141 online DOI: 10.1081/SAC-00055639 Distribution Theory Comparison Between
More informationConfidence Estimation Methods for Neural Networks: A Practical Comparison
, 6-8 000, Confidence Estimation Methods for : A Practical Comparison G. Papadopoulos, P.J. Edwards, A.F. Murray Department of Electronics and Electrical Engineering, University of Edinburgh Abstract.
More informationDouble Bootstrap Confidence Interval Estimates with Censored and Truncated Data
Journal of Modern Applied Statistical Methods Volume 13 Issue 2 Article 22 11-2014 Double Bootstrap Confidence Interval Estimates with Censored and Truncated Data Jayanthi Arasan University Putra Malaysia,
More informationGoodness-of-fit tests for randomly censored Weibull distributions with estimated parameters
Communications for Statistical Applications and Methods 2017, Vol. 24, No. 5, 519 531 https://doi.org/10.5351/csam.2017.24.5.519 Print ISSN 2287-7843 / Online ISSN 2383-4757 Goodness-of-fit tests for randomly
More informationA Convenient Way of Generating Gamma Random Variables Using Generalized Exponential Distribution
A Convenient Way of Generating Gamma Random Variables Using Generalized Exponential Distribution Debasis Kundu & Rameshwar D. Gupta 2 Abstract In this paper we propose a very convenient way to generate
More informationParameters Estimation for a Linear Exponential Distribution Based on Grouped Data
International Mathematical Forum, 3, 2008, no. 33, 1643-1654 Parameters Estimation for a Linear Exponential Distribution Based on Grouped Data A. Al-khedhairi Department of Statistics and O.R. Faculty
More informationThe Nonparametric Bootstrap
The Nonparametric Bootstrap The nonparametric bootstrap may involve inferences about a parameter, but we use a nonparametric procedure in approximating the parametric distribution using the ECDF. We use
More informationEstimation of parametric functions in Downton s bivariate exponential distribution
Estimation of parametric functions in Downton s bivariate exponential distribution George Iliopoulos Department of Mathematics University of the Aegean 83200 Karlovasi, Samos, Greece e-mail: geh@aegean.gr
More informationHazard Function, Failure Rate, and A Rule of Thumb for Calculating Empirical Hazard Function of Continuous-Time Failure Data
Hazard Function, Failure Rate, and A Rule of Thumb for Calculating Empirical Hazard Function of Continuous-Time Failure Data Feng-feng Li,2, Gang Xie,2, Yong Sun,2, Lin Ma,2 CRC for Infrastructure and
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationBAYESIAN PREDICTION OF WEIBULL DISTRIBUTION BASED ON FIXED AND RANDOM SAMPLE SIZE. A. H. Abd Ellah
Serdica Math. J. 35 (2009, 29 46 BAYESIAN PREDICTION OF WEIBULL DISTRIBUTION BASED ON FIXED AND RANDOM SAMPLE SIZE A. H. Abd Ellah Communicated by S. T. Rachev Abstract. We consider the problem of predictive
More informationESTIMATOR IN BURR XII DISTRIBUTION
Journal of Reliability and Statistical Studies; ISSN (Print): 0974-804, (Online): 9-5666 Vol. 0, Issue (07): 7-6 ON THE VARIANCE OF P ( Y < X) ESTIMATOR IN BURR XII DISTRIBUTION M. Khorashadizadeh*, S.
More informationExact Solutions for a Two-Parameter Rayleigh Distribution
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 11 (2017), pp. 8039 8051 Research India Publications http://www.ripublication.com/gjpam.htm Exact Solutions for a Two-Parameter
More informationBest linear unbiased and invariant reconstructors for the past records
Best linear unbiased and invariant reconstructors for the past records B. Khatib and Jafar Ahmadi Department of Statistics, Ordered and Spatial Data Center of Excellence, Ferdowsi University of Mashhad,
More informationNoninformative Priors for the Ratio of the Scale Parameters in the Inverted Exponential Distributions
Communications for Statistical Applications and Methods 03, Vol. 0, No. 5, 387 394 DOI: http://dx.doi.org/0.535/csam.03.0.5.387 Noninformative Priors for the Ratio of the Scale Parameters in the Inverted
More informationCOMPETING RISKS WEIBULL MODEL: PARAMETER ESTIMATES AND THEIR ACCURACY
Annales Univ Sci Budapest, Sect Comp 45 2016) 45 55 COMPETING RISKS WEIBULL MODEL: PARAMETER ESTIMATES AND THEIR ACCURACY Ágnes M Kovács Budapest, Hungary) Howard M Taylor Newark, DE, USA) Communicated
More information4 Resampling Methods: The Bootstrap
4 Resampling Methods: The Bootstrap Situation: Let x 1, x 2,..., x n be a SRS of size n taken from a distribution that is unknown. Let θ be a parameter of interest associated with this distribution and
More informationSpring 2012 Math 541B Exam 1
Spring 2012 Math 541B Exam 1 1. A sample of size n is drawn without replacement from an urn containing N balls, m of which are red and N m are black; the balls are otherwise indistinguishable. Let X denote
More informationThe Relationship Between Confidence Intervals for Failure Probabilities and Life Time Quantiles
Statistics Preprints Statistics 2008 The Relationship Between Confidence Intervals for Failure Probabilities and Life Time Quantiles Yili Hong Iowa State University, yili_hong@hotmail.com William Q. Meeker
More informationA comparison of inverse transform and composition methods of data simulation from the Lindley distribution
Communications for Statistical Applications and Methods 2016, Vol. 23, No. 6, 517 529 http://dx.doi.org/10.5351/csam.2016.23.6.517 Print ISSN 2287-7843 / Online ISSN 2383-4757 A comparison of inverse transform
More informationStatistical Inference Using Progressively Type-II Censored Data with Random Scheme
International Mathematical Forum, 3, 28, no. 35, 1713-1725 Statistical Inference Using Progressively Type-II Censored Data with Random Scheme Ammar M. Sarhan 1 and A. Abuammoh Department of Statistics
More informationUnivariate and Bivariate Geometric Discrete Generalized Exponential Distributions
Univariate and Bivariate Geometric Discrete Generalized Exponential Distributions Debasis Kundu 1 & Vahid Nekoukhou 2 Abstract Marshall and Olkin (1997, Biometrika, 84, 641-652) introduced a very powerful
More informationLoad-strength Dynamic Interaction Principle and Failure Rate Model
International Journal of Performability Engineering Vol. 6, No. 3, May 21, pp. 25-214. RAMS Consultants Printed in India Load-strength Dynamic Interaction Principle and Failure Rate Model LIYANG XIE and
More information