Inferences on stress-strength reliability from weighted Lindley distributions
|
|
- Willis Barrett
- 5 years ago
- Views:
Transcription
1 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 R = P (Y < X), when X and Y are two independent weighted Lindley random variables with a common shape parameter. The stress strength parameter cannot be obtained in closed form and it has been expressed in terms of hyper-geometric functions. When the common shape parameter is known, the minimum variance unbiased estimator of R exists. It is quite complicated, and it is not pursued any further here. In this case the maximum likelihood estimators (MLEs) can be obtained explicitly. When the common shape parameter is unknown, the MLEs of the unknown parameters cannot be obtained in explicit forms. The MLEs can be obtained by maximizing the profile log-likelihood function in one dimension. The asymptotic distribution of the MLEs are also obtained, and they have been used to construct the asymptotic confidence interval of R. Bootstrap confidence intervals are also proposed. Monte Carlo simulations are performed to verify the effectiveness of the different methods, and data analysis has been performed for illustrative purposes. Key Words and Phrases: Lindley distribution; Maximum likelihood estimator; Asymptotic distribution; Uniformly minimum variance unbiased estimator; Bootstrap confidence interval. Department of Statistics and Operations Research, Faculty of Science, Kuwait University, P.O. Box 5969 Safat, Kuwait Department of Mathematics and Statistics, Indian Institute of Technology Kanpur, Pin , INDIA. Corresponding author, Phone: , Fax: ; e- mail: kundu@iitk.ac.in. 1
2 2 1 Introduction Lindley (1958) originally proposed a lifetime distribution which has the following probability density function (PDF) f(x; θ) = θ2 1 + θ (1 + x)e θx, x > 0, θ > 0. (1) The distribution function which has the PDF (1) is well known as the Lindley distribution. This distribution was proposed in the context of Bayesian statistics, and it is used as a counter example of fudicial statistics. Ghitany et al. (2008) discussed different properties of the Lindley distribution. It is known that the shape of the PDF of the Lindley distribution is either unimodal or a decreasing function, and the hazard function is always an increasing function. Although, Lindley distribution has some interesting properties, due to presence of only one parameter, it may not very useful in practice. Due to this reason, Ghitany et al. (2011) proposed a two-parameter extension of the Lindley distribution, and named it as the weighted Lindley (WL) distribution. The WL distribution is more flexible than the Lindley distribution, and Lindley distribution can be obtained as a special case of the WL distribution. The WL distribution can be defined as follows. A random variable X is said to have a WL distribution if it has the PDF f(x) = θ c+1 (θ + c) Γ(c) xc 1 (1 + x)e θx, x > 0, c, θ > 0. (2) From now on a WL distribution with the parameters c and θ, and PDF (2) will be denoted by WL(c, θ). The WL distribution has several interesting properties. A brief description of the WL distribution, and some of its new properties are discussed in Section 2. The main aim of this paper is to consider the inference on the stress-strength parameter R = P (Y < X), where X WL(c, θ 1 ), Y WL(c, θ 2 ) and they are independently distributed. Here the notation means follows or has the distribution. The estimation of the
3 3 stress-strength parameter arises quite naturally in the area of system reliability. For example, if X is the strength of a component which is subject to stress Y, then R is a measure of system performance and arises in the context of mechanical reliability of a system. The system fails if and only if at any time the applied stress is greater than its strength. In this paper first we consider the case when all the three parameters are unknown. The MLEs of the unknown parameters cannot be obtained in explicit forms. For a given c, the MLEs of θ 1 and θ 2 can be obtained explicitly, and finally the MLE of c can be obtained by maximizing the profile likelihood function of c. The asymptotic distribution of the MLEs of the unknown parameters and also of R are obtained. The asymptotic distribution has been used to construct asymptotic confidence interval of the stress strength parameter. Since the asymptotic variance of the MLE of R is quite involved, we have proposed to use parametric and non-parametric Bootstrap confidence intervals also. Different methods are compared using Monte Carlo simulations and data analysis has been performed for illustrative purposes. We further consider a special case when the common shape parameter c is known. In this case it is observed that the uniformly minimum variance unbiased estimator (UMVUE) of R exists, although the exact form is very complicated and it is not not pursued any further. When c is known, the MLEs of θ 1 and θ 2 can be obtained in explicit forms. The asymptotic distribution of the MLE of R also can be obtained quite conveniently, and which can be used to construct asymptotic confidence interval of R. It may be mentioned that the estimation of the stress-strength parameter has received considerable attention in the statistical literature starting with the pioneering work of Birnbaum (1956), where the author has provided an interesting connection between the classical Mann-Whitney statistic and the stress-strength model. Since then extensive work has been done on developing the inference procedure on R for different distributions, both from the Bayesian and frequentist points of view. The monograph by Kotz et al. (2003) provided
4 4 an excellent review on this topic till that time. For some of the recent references the readers are referred to the articles by Kundu and Gupta (2005, 2006), Kim and Chung (2006), Krishnomoorthy et al. (2007), Saraçoğlu and Kaya (2007), Jiang and Wong (2008), Raqab et al. (2008), Saraçoğlu et al. (2009), Kundu and Raqab (2009), Saraçoğlu et al. (2012), Al-Mutairi et al. (2013), and see the references cited therein. The rest of the paper is organized as follows. In Section 2, we briefly discuss about the WL distribution, and discuss some of its properties. In Section 3, we discuss the MLE of R and derive its asymptotic properties when all the parameters are unknown. Different confidence intervals are discussed in Section 4. The estimation of R when the common shape parameter c is known is discussed in Section 5. Monte Carlo simulation results are provided in Section 6, and the data analysis has been presented in Section 7. Finally, we conclude the paper in Section 8. 2 Weighted Lindley Distribution A random variable X is said to have a WL distribution with parameter c and θ, if it has the PDF (2). The cumulative distribution function (CDF) and the survival function (SF) take the following forms; F (x) = 1 (θ + c)γ(c, θx) + (θx)c e θx, x > 0, c, θ > 0 (3) (θ + c) Γ(c) S(x) = (θ + c)γ(c, θx) + (θx)c e θx, x > 0, c, θ > 0, (4) (θ + c) Γ(c) where Γ(a, z) = z y a 1 e y dy, a > 0, z 0. It has been observed by Ghitany et al. (2011) that the PDF of the WL distribution is either decreasing, unimodal or decreasing-increasing-decreasing shape depending on the values of c
5 5 and θ. The hazard function of the WL distribution is either bathtub shaped or an increasing function depending on the value of c. It is clear from (3) and (4) that the inversion of the CDF or the SF cannot be performed analytically. The following mixture representation is useful in generating random sample and also deriving some other properties for the WL distribution. The PDF of WL distribution can be written as Here p = θ/(θ + c), and f GA (x; c + j 1) = f(x) = pf GA (x; c, θ) + (1 p)f GA (x; c + 1, θ). (5) θ c+j 1 Γ(c + j 1) xc+j 2 e θx ; x > 0, c, θ > 0, j = 1, 2, is the PDF of the gamma distribution with the shape parameter c+j 1 and scale parameter θ. From now on it will be denoted by Gamma(c+j 1, θ), j = 1, 2. Since the generation from a gamma distribution can be performed quite efficiently, see for example Kundu and Gupta (2007), the generation from a WL distribution becomes quite straight forward. Using (5), it immediately follows that if X WL(c, θ), then the moment generating function (MGF) of X is M X (t) = E(e tx ) = θ ( 1 t ) c + c ( 1 t ) (c+1). (6) θ + c θ θ + c θ The following result will be useful to derive the UMVUE of R, and it may have some independent interest also. Theorem 1: If X 1,..., X n are independent identically distributed WL(c, θ) random variables, then the PDF of Z = X 1, X n n f Z (t) = p n,k (c, θ)f GA (t; nc + n k, θ), k=0 where p n,k (c, θ) = ( ) n θ k c n k k (c + θ). n
6 6 Proof: The proof mainly follows using the MGF of Z. The MGF of Z is M Z (t) = E(e tz ) = = [ θ θ + c ( n n k k=0 ( 1 t ) c + c θ ) ( θ θ + c Now the result follows by using the fact that θ + c ) k ( c θ + c ( 1 t ) ] (c+1) n θ ) n k ( 1 t θ ) (nc+n k). ( 1 t θ ) a is the MGF of Gamma(a, θ). It may be mentioned that the WL distribution belongs to the two-parameter exponential n n family. If x 1,..., x n, is a random sample of size n from WL(c, θ), then ( x i, x i ) is jointly a complete sufficient statistic for (c, θ). We further have the following result regarding likelihood ratio ordering. i=1 i=1 Theorem 2: Let X WL(c 1, θ 1 ) and Y WL(c 2, θ 2 ). If c 1 c 2 and θ 1 θ 2, then X is larger than Y in the sense of likelihood ratio order, written as X lr Y. Proof: From the definition of the likelihood ratio, f(x; c 1, θ 1 ) f(x; c 2, θ 2 ) = (c 2 + θ 2 ) Γ(c 2 ) θ1 c1 1 (c 1 + θ 1 ) Γ(c 1 ) θ c x c 1 c 2 e (θ 1 θ 2 ) x 2 1 increases in x, since c 1 c 2 and θ 1 θ 2. Therefore the result follows. 2 Remarks. (i) If c 1 = c 2, then X lr Y when θ 1 θ 2. It further implies the stochastic ordering also between X and Y. 3 Maximum Likelihood Estimator of R In this section first we provide the expression of R, and then provide the maximum likelihood estimator of R when all the parameters are unknown.
7 7 3.1 Expression of R To derive the expression of R, we will be using the following notation: 0 where Γ(a, x) x s 1 Γ(a + s) β e x dx = 2F s(1 + β) a+s 1 (a + s, 1; s + 1; ), s > 0, β > β 2F 1 (a, b; c; z) = 2 F 1 (b, a; c; z) = n=0 (a) n (b) n (c) n z n n 1 n!, z < 1, (a) 0 = 1, (a) n = (a + k), k=0 is the hyper-geometric function. We assume c 0, 1, 2,... to prevent the denominators from vanishing, see Erdélyi et al. (1954), p. 325, formula (16). Let X WL(c, θ 1 ) and Y WL(c, θ 2 ) be independent random variables. Then R = P (X > Y ) = = = = = = 0 0 P (X > Y Y = y). f Y (y) dy S X (y). f Y (y) dy θ2 c+1 { (θ1 + c) (θ 1 + c)(θ 2 + c) Γ 2 (c) 0 +θ c 1 Γ(c, θ 1 y). (y c 1 + y c )e θ 2y dy 0 (y 2c 1 + y 2c )e (θ 1+θ 2 )y dy } θ2 c+1 { (θ 1 + c) (θ 1 + c)(θ 2 + c) Γ 2 (c) θ1 c+1 Γ(c, t). (θ 1 t c 1 + t c )e (θ 2/θ 1 )t dt 0 + θc 1 Γ(2c) (θ 1 + θ 2 ) (θ 2c θ 2 + 2c) } θ1 c θ2 c+1 Γ(2c) { (θ 1 + c)(θ 1 + θ 2 ) 2F (θ 1 + c)(θ 2 + c)(θ 1 + θ 2 ) 2c+1 Γ 2 1 (2c, 1; c + 1; (c) c + 2c(θ 1 + c) 2F 1 (2c + 1, 1; c + 2; c + 1 θ1 c θ2 c+1 Γ(2c) { θ 1 + θ 2 θ 2 2F (θ 2 + c)(θ 1 + θ 2 ) 2c+1 Γ 2 1 (2c, 1; c + 1; ) (c) c θ 1 + θ 2 θ 2 θ 1 + θ 2 ) θ 2 ) + θ 1 + θ 2 + 2c } θ 1 + θ 2 + 2c θ 2 c F 1 (2c + 1, 1; c + 2; ) + θ 2 + c θ 1 + θ 2 θ 1 + c + 1}
8 8 Remarks. (i) When θ 1 = θ 2, R = 0.5. This, of course, is expected since, in this case, X and Y are i.i.d. and there is an equal chance that X is bigger than Y. 3.2 Maximum Likelihood Estimator of R Suppose x 1, x 2,..., x n1 is a random sample of size n 1 from the WL(c, θ 1 ) and y 1, y 2,..., y n2 is an independent random sample of size n 2 from the WL(c, θ 2 ). The log-likelihood function l l(c, θ 1, θ 2 ) = l(θ) based on the two independent random samples is given by l = n 1 n 2 ln[f X (x i )] + ln[f Y (y j )] i=1 j=1 = n 1 [(c + 1) ln(θ 1 ) ln(θ 1 + c) ln(γ(c))] + ln(1 + x i ) + (c 1) ln(x i ) n 1 θ 1 x + n 2 [(c + 1) ln(θ 2 ) ln(θ 2 + c) ln(γ(c))] + n 1 i=1 n 2 ln(1 + y j ) + (c 1) n 1 i=1 n 2 j=1 j=1 ln(y j ) n 2 θ 2 ȳ, where x and ȳ are the the sample means of x 1,..., x n1 and y 1,..., y n2, respectively. The maximum likelihood estimator (MLE) θ of θ is the solutions of the non-linear equations: l c = n 1 [ ln(θ 1 ) 1 θ 1 + c n 1 n 2 + ln(x i ) + ln(y i ) = 0, i=1 j=1 l c(θ 1 + c + 1) = n 1 n 1 x = 0, θ 1 θ 1 (θ 1 + c) l c(θ 2 + c + 1) = n 2 n 2 ȳ = 0, θ 2 θ 2 (θ 2 + c) ] [ + n 2 ln(θ 2 ) 1 ] (n 1 + n 2 )ψ(c) θ 2 + c where ψ(c) = d ln(γ(c)) is the digamma function. It is clear that the MLEs of the unknown dc parameters cannot be obtained in explicit forms. The MLEs are as follows: θ 1 θ 1 (ĉ) = ĉ( x 1) + ĉ 2 ( x 1) 2 + 4ĉ(ĉ + 1) x 2 x
9 9 θ 2 θ 2 (ĉ) = ĉ(ȳ 1) + ĉ 2 (ȳ 1) 2 + 4ĉ(ĉ + 1)ȳ, 2ȳ where ĉ is the solution of the non-linear equation: 1 n 1 [ ln( θ 1 (ĉ)) θ 1 (ĉ) + ĉ ]+n 1 2[ ln( θ 2 (ĉ)) θ 2 (ĉ) + ĉ ] (n 1 +n 2 )ψ(ĉ)+ ln(x i )+ ln(y i ) = 0. Some iterative procedure needs to be used to solve the above non-linear equation. n 1 i=1 n 2 j=1 Now we provide the elements of the expected Fisher information matrix, which will be needed to construct the asymptotic confidence intervals of the unknown parameters. We will use the following notations. p 1 = lim n 1,n 2 n 1 n 1 + n 2, p 2 = lim n 1,n 2 n 2 n 1 + n 2 and H(θ) = 2 l c 2 2 l θ 1 c 2 l θ 2 c 2 l c θ 1 2 l θ l θ 2 θ 1 2 l c θ 2 2 l θ 1 θ 2 The elements of the symmetric expected Fisher information matrix of (c, θ 1, θ 2 ) is given by In this case 2 l θ 2 2. [ ] I(θ) = ((I ij )) = lim E 1 H(θ). n 1,n 2 n 1 + n 2 I 11 = ψ (c) I 22 = p 1 [ c + 1 θ 2 1 I 33 = p 2 [ c + 1 θ 2 2 p 1 (θ 1 + c) 2 p 2 (θ 2 + c) 2, I 12 = p 1 [ 1 θ 1 + I 13 = p 2 [ 1 θ 2 + I 23 = 0, 1 (θ 1 + c) 2 ], 1 ], (θ 2 + c) 2 1 ], (θ 1 + c) 2 1 (θ 2 + c) 2 ],
10 10 where ψ (c) = d ψ(c) is the trigamma function. Since WL family satisfies all the regularity dc of the consistency and asymptotic normality of the MLEs, see for example Lehmann and Casella (1998), pp , we have the following result. Theorem 3: The asymptotic distribution of the MLE θ of θ satisfies n1 + n 2 ( θ θ) d N 3 (0, I 1 (θ)), where d denotes convergence in distribution and I 1 (θ) is the inverse of the matrix I. The asymptotic variance-covariance matrix of the MLE θ is given by V ar(ĉ) Cov(ĉ, ˆθ 1 ) Cov(ĉ, ˆθ 2 ) 1 I 1 (θ) = Cov(ĉ, n 1 + n 2 ˆθ 1 ) V ar( ˆθ 1 ) Cov( ˆθ 1, ˆθ 2 ). Cov(ĉ, ˆθ 2 ) Cov( ˆθ 1, ˆθ 2 ) V ar( ˆθ 2 ) We need the following integral representation for further development, see for example Andrews (1998), p. 364: 2F 1 (a, b; c; z) = Γ(c) Γ(b)Γ(c b) t b 1 (1 t) c b 1 (1 zt) a dt, zt < 1, c > b > 0, 2F 1 (2c, 1; c + 1; z) = c (1 t) c 1 (1 zt) 2c dt 0 J 1 (c, z) = { } 1 c c 2 F 1 (2c, 1; c + 1; z) ( ) 1 1 t = (1 t) c 1 (1 zt) 2c ln dt 0 (1 zt) 2 1 2F 1 (2c + 1, 1; c + 2; z) = (c + 1) (1 t) c (1 zt) 2c 1 dt 0 J 2 (c, z) = { } 1 c c F 1 (2c + 1, 1; c + 2; z) ( ) 1 1 t = (1 t) c (1 zt) 2c 1 ln dt. 0 (1 zt) 2 In order to establish the asymptotic normality of R, we further define where T is the transpose operation and d(θ) = ( R c, R θ 1, R θ 2 ) T = ( d1, d 2, d 3 ) T,
11 11 d 1 = θ1 c θ2 c+1 Γ(2c) { (θ1 + θ (θ 2 + c)(θ 1 + θ 2 ) 2c+1 Γ 2 2 ) J 1 (c, (c) + 2 θ 2 θ 1 + θ 2 ) + 2c J 2 (c, ) θ 1 + θ 2 θ 2 c F 1 (2c + 1, 1; c + 2; ) + θ 1 θ } 2 θ 1 + θ 2 (θ 1 + c) 2 +R { ( ) θ1 θ 2 ln (θ 1 + θ 2 ) 2 θ 2 + 2ψ(2c) 2ψ(c) 1 θ 2 + c } d 2 = d 2 (θ 1, θ 2 ) = Γ(c ){ 4c 3 + (θ 1 + 1)(θ 1 + θ 2 ) 2 + c 2 (2 + 8θ 1 + 4θ 2 ) + c[5θ1 2 + θ 2 (θ 2 + 2) + θ 1 (4 + 6θ 2 )] } 2 1 2c θ c, π Γ(c)(θ1 + c) 2 (θ 2 + c)(θ 1 + θ 2 ) 2c+2 1 θ2 c 1 d 3 = d 2 (θ 2, θ 1 ). Therefore, using Theorem 3 and delta method, we obtain the asymptotic distribution of R, the MLE of R, as n1 + n 2 ( ˆR R) d N ( 0, d T (θ) I 1 (θ) d(θ) ). (7) From (7) the asymptotic variance of R is obtained as V ar( R) = 1 n 1 + n 2 d T (θ) I 1 (θ) d(θ) = d 2 1 V ar(ĉ) + d 2 2 V ar( θ 1 ) + d 2 3 V ar( θ 2 ) + 2 d 1 d 2 Cov(ĉ, θ 1 ) +2 d 1 d 3 Cov(ĉ, θ 2 ) + 2 d 2 d 3 Cov( θ 1, θ 2 ). Hence, using (7), an asymptotic 100(1 α)% confidence interval for R can be obtained as R z α 2 V ar( ˆR), where z α 2 is the upper α/2 quantile of the standard normal distribution. Since, the expression of asymptotic variance of R is quite involved, we consider different bootstrap confidence intervals of R.
12 12 4 Bootstrap Confidence Intervals In this section we provide different parametric and non-parametric bootstrap confidence intervals of R. It is assumed that we have independent random samples x 1,, x n1 and y 1,, y n2 obtained from WL(c, θ 1 ) and WL(c, θ 2 ), respectively. First we propose to use the following method to generate non-parametric bootstrap samples, as suggested by Efron and Tibshirani (1998), from the given random samples. Algorithm: (Non-parametric bootstrap sampling) Step 1: Generate independent bootstrap samples x 1,, x n 1 and y 1,, y n 2 taken with replacement from the given samples x 1,, x n1 and y 1,, y n2, respectively. Based on the bootstrap samples, compute the MLE (ĉ, θ 1, θ 2) of (c, θ 1, θ 2 ) as well as R = R(ĉ, θ 1, θ 2) of R. Step 2: Repeat Step 1, B times to obtain a set of bootstrap samples of R, say { R j, j = 1,, B}. Using the above bootstrap samples of R we obtain three different bootstrap confidence interval of R. The ordered R j for j = 1, B will be denoted as: R (1) < < R (B). (i) Percentile bootstrap (p-boot) confidence interval: Let R (τ) be the τ percentile of { R j, j = 1, 2,..., B}, i.e. R (τ) is such that 1 B B j=1 I( R j R (τ) ) = τ, 0 < τ < 1, where I( ) is the indicator function.
13 13 A 100(1 α)% p-boot confidence interval of R is given by ( R (α/2), R (1 α/2) ). (ii) Student s t bootstrap (t-boot) confidence interval: Let se ( R) be the sample standard deviation of { R j, j = 1, 2,..., B}, i.e. se ( R) = 1 B B ( R j R ) 2. j=1 Also, let t (τ) be the τ percentile of { R j R se ( R), j = 1, 2,..., B}, i.e. t (τ) is such that 1 B B j=1 I( R j R se ( R) t (τ) ) = τ 0 < τ < 1. A 100(1 α)% t-boot confidence interval of R is given by R ± ˆt (α/2) se ( R). (iii) Bias-corrected and accelerated bootstrap (BC a -boot) confidence interval: Let z (τ) and ẑ 0, respectively, be such that z (τ) = Φ 1 (τ) and ẑ 0 = Φ 1( 1 B B j=1 I( R j R) ), where Φ 1 (.) is the inverse CDF of the standard normal distribution. The value ẑ 0 is called bias-correction. Also, let â = ni=1 ( R (.) R (i) ) 3 6 [ ni=1 ( R (.) R (i) ) 2 ] 3/2 where R (i) is the MLE of R based of (n 1) observations after excluding the ith observation and R ( ) = 1 n ni=1 R (i). The value â is called acceleration factor.
14 14 A 100(1 α)% BC a -boot confidence interval of R is given by ( R (ν 1), R (ν 2) ), where ν 1 = Φ ( ẑ 0 + ẑ 0 + z (α/2) ), ν2 = Φ ( ẑ 0 + z (1 α/2) ) ẑ 1 â(ẑ 0 + z (α/2) 0 +. ) 1 ˆα(ẑ 0 + z (1 α/2) ) Now we provide the following method to generate parametric bootstrap samples from the given random samples, and they can be used to construct different parametric bootstrap confidence intervals. Algorithm: (Parametric bootstrap sampling) Step 1: Compute the MLEs of c, θ 1 and θ 2 from the given random samples, say ĉ, θ 1 and θ 2, respectively. Generate independent bootstrap samples x 1,, x n 1 and y1,, yn 2, from WL(ĉ, θ 1 ) and WL(ĉ, θ 2 ), respectively. Based on the bootstrap samples, compute the MLE (ĉ, θ 1, θ 2) of (c, θ 1, θ 2 ) as well as R = R(ĉ, θ 1, θ 2) of R. Step 2: Repeat Step 1, B times to obtain a set of bootstrap samples of R, say { R j, j = 1,, B}. Using the above bootstrap samples of R we can obtain three different parametric bootstrap confidence interval of R same as the non-parametric ones. It may be mentioned that all the bootstrap confidence intervals can be obtained even in the logit scale also, and we have presented those results in Section 6. In the next section we provide two estimators of R when the common shape parameter c is known.
15 15 5 Estimation of R when c is Known 5.1 Maximum Likelihood Estimator of R Based on the random samples {x 1, x 2,..., x n1 } from WL(c, θ 1 ) and {y 1, y 2,..., y n2 } from the WL(c, θ 2 ), the MLEs of θ 1 and θ 2 can be obtained as ˆθ 1 = c( x 1) + c 2 ( x 1) 2 + 4c(c + 1) x (8) 2 x ˆθ 2 = c(ȳ 1) + c 2 (ȳ 1) 2 + 4c(c + 1)ȳ. (9) 2ȳ The asymptotic distribution of the MLE ( θ 1, θ 2 ) of (θ 1, θ 2 ) satisfies n1 + n 2 { θ 1 θ 1, θ 2 θ 2 } d N 2 (0, A), where A = [ I 1 ] I33 1, and the asymptotic distribution of the MLE R of R satisfies n1 + n 2 ( ˆR R) d N ( ) 0, d 2 I d 3 I33 1. The elements I 22, I 33, d 2 and d 3 are same as defined in Section 3.2. The asymptotic variance of ˆR is given by V ar( ˆR) = 1 ( ) d 2 n 1 + n 2 I d 2 3I The asymptotic 100(1 α)% confidence interval for R is given by R z α 2 V ar( ˆR), where z α 2 is the upper α 2 quantile of the standard normal distribution.
16 Uniformly Minimum Variance Unbiased Estimator of R n If x 1,..., x n is a random sample from WL(c, θ), then x i is a complete sufficient statistic i=1 for θ when c is known. We need the following result for further development. Theorem 4: If X 1,..., X n are i.i.d. Lindley (c, θ) random variables, then the conditional n PDF of X 1 given Z = X i is where i=1 f X1 Z=z(x) = xc 1 (1 + x) Γ(c)A n (z) n 1 k=0 ( ) n 1 c n k 1 C k,n = k Γ(nc + n c k 1) Proof: It mainly follows from the fact and using Theorem 1. C k,n (z x) nc+n c k 2 ; 0 < x < z, and A n (z) = n k=0 f X1 Z=z(x) = f X 1 (x)f X X n (z x) f Z=z (z) ( ) n z nc+n k 1 k Γ(nc + n k). Now we can provide the UMVUE of R in the following Theorem. Theorem 5: Suppose form n 1 x i = u and i=1 n 2 y j = v, then the UMVUE of R has the following j=1 R UMV UE = here 1 A n1 (u)a n2 (v)[γ(c)] 2 I(u, v, a, b) = n 1 1 m=0 min{u,v} u 0 n 1 i=1 y n 2 1 n=0 C m,n1 C n,n2 I(u, v, n 1 c+n 1 m c 2, n 2 c+n 2 n c 2), (xy) c 1 (1 + x)(1 + y)(u x) a (v y) b dxdy. (10) Proof: Let us denote U = X i and V = n 2 j=1 Y j. For known c, since (U, V ) is a complete sufficient statistic for (θ 1, θ 2 ), therefore, the UMVUE of R can be obtained as R UMV UE = E(ϕ(X 1, Y 1 ) U = u, V = v),
17 17 where Therefore, R UMV UE = 1 if x > y ϕ(x, y) = 0 if x y min{u,v} u 0 y f X1 U=u(x)f Y1 V =v(y)dxdy. Now the result follows using Theorem 4. It may be noted that if c is an integer, then a and b in (10) can be taken to be integers, and in that case I(u, v, a, b, ) can be expressed in terms of double summation. For general a and b I(u, v, a, b) = u a+c min{u,v} 0 here B(x; a, b) is the incomplete Beta function, i.e. [ y c 1 (v y) b (1+y) B (1 y ) u ; a + 1, c + ub (1 y )] u ; a + 1, c + 1 dy, B(x; a, b) = x 0 t a 1 (1 t) b 1 dt. Clearly, the UMVUE has a quite complicated form, and it will not be pursued further. 6 Monte Carlo Simulations In this section we present some Monte Carlo simulation results mainly to compare the performances of different methods for different sample sizes and for different parameter values. We mainly investigate the performance of the point and interval estimation of the reliability R = P (X > Y ) based on maximum likelihood procedure when all the parameters are unknown. Specifically, we investigate the bias and mean square error (MSE) of the simulated MLEs. Also, we investigate the coverage probability and the length of the simulated 95% confidence intervals based on maximum likelihood method and different bootstrap procedures presented in Section 4.
18 18 For this purpose, we have generated 3, 000 samples with B = 1, 000 bootstrap samples from each of independent WL(c, θ 1 ) and WL(c, θ 2 ) distributions where (c, θ 1, θ 2 ): (0.75, 2.5, 1), (2, 1, 1) and (4, 2, 5). These parameter values correspond to R = 0.26, 0.5 and 0.91, respectively. We have taken different sample sizes namely (n 1, n 2 ): (20,20), (20,30), (20,50), (30,20), (30,30),(30,50), (50,20), (50,30), (50,50). In Tables 3-5, we report the average biases and the mean squared errors of the estimates of R based on MLEs, parametric bootstrap and non-parametric bootstrap methods. In Tables 6-8, we provide the coverage percentages and the average lengths of the confidence intervals of R based on MLE of R and using different bootstrap methods. Some of the points are quite clear from these experiments. Even for small sample sizes the performances of the MLEs are quite satisfactory in terms of biases and MSEs. It is observed that when n 1 = n 2 = n, and increases, the bias and MSEs decrease. It verifies the consistency property of the MLE of R. For fixed n 1 (n 2 ) as n 2 (n 1 ) increases, the MSEs decrease. Comparing the average confidence lengths and coverage percentages, it is observed that the performance of the confidence intervals of R, based on the asymptotic distribution of MLE is quite satisfactory. It maintains the nominal coverage percentages even for small sample sizes. Among the different bootstrap confidence intervals the biased corrected parametric bootstrap confidence intervals perform the best in terms of the coverage percentages. In most of the cases considered here, the confidence intervals based on MLE performs better than the biased corrected bootstrap confidence intervals in terms of shorter average confidence intervals. 7 Data Analysis In this section we consider two data sets and describe all the details for illustrative purposes. The two data sets were originally reported by Bader and Priest (1982), on failure stresses
19 19 (in GPa) of single carbon fibers of lengths 20 mm and 50 mm, respectively. We present the data sets below. Length 20 mm: X (n 1 = 69): 1.312, 1.314, 1.479, 1.552, 1.700, 1.803, 1.861, 1.865, 1.944, 1.958, 1.966, 1.997, 2.006, 2.021, 2.027, 2.055, 2.063, 2.098, 2.14, 2.179, 2.224, 2.240, 2.253, 2.270, 2.272, 2.274, 2.301, 2.301, 2.359, 2.382, 2.382, 2.426, 2.434, 2.435, 2.478, 2.490, 2.511, 2.514, 2.535, 2.554, 2.566, 2.57, 2.586, 2.629, 2.633, 2.642, 2.648, 2.684, 2.697, 2.726, 2.770, 2.773, 2.800, 2.809, 2.818, 2.821, 2.848, 2.88, 2.954, 3.012, 3.067, 3.084, 3.090, 3.096, 3.128, 3.233, 3.433, 3.585, Length 50 mm: Y (n 2 = 65): 1.339, 1.434, 1.549, 1.574, 1.589, 1.613, 1.746, 1.753, 1.764, 1.807, 1.812, 1.84, 1.852, 1.852, 1.862, 1.864, 1.931, 1.952, 1.974, 2.019, 2.051, 2.055, 2.058, 2.088, 2.125, 2.162, 2.171, 2.172, 2.18, 2.194, 2.211, 2.27, 2.272, 2.28, 2.299, 2.308, 2.335, 2.349, 2.356, 2.386, 2.39, 2.41, 2.43, 2.431, 2.458, 2.471, 2.497, 2.514, 2.558, 2.577, 2.593, 2.601, 2.604, 2.62, 2.633, 2.67, 2.682, 2.699, 2.705, 2.735, 2.785, 3.02, 3.042, 3.116, Before progressing further first we perform some preliminary data analysis. The scaled TTT plots of both X and Y (see Figure 1) indicate that the empirical hazard functions will be increasing functions in both the cases. Therefore, WL can be used to analyze the two data sets. First it is assumed that X WL(c 1, θ 1 ) and Y WL(c 2, θ 2 ). The MLEs of the unknown parameters are as follows: ĉ 1 = , θ 1 = , ĉ 2 = , θ 2 = , and the associated log-likelihood value is L 1 = Second suppose that X WL(c, θ 1 ) and Y WL(c, θ 2 ). Then, the MLEs of the unknown parameters are as follows: ĉ = , θ 1 = , θ 2 = , and the associated log-likelihood value is L 0 = We perform the following testing of hypothesis; H 0 : c 1 = c 2, vs. H 1 : c 1 c 2, and in this case 2(L 0 L 1 ) = Hence, the null hypothesis cannot be rejected.
20 20 Therefore, in this case the assumption of c 1 = c 2 is justified (a) (b) Figure 1: Scaled TTT plots: (a) Length 20 mm data (b) Length 50 mm data. To compute the MLEs of c, θ 1 and θ 2, we use the profile likelihood method. Figure 2 provides the profile log-likelihood of c. It indicates that it has a unique maximum. In Table 1 we provide the MLEs of the parameters c, θ 1, θ 2 as well as the Kolmogorov-Smirnov (K-S) and Anderson-Darling (A-D) goodness-of-fit tests. The table shows that both the tests accept the null hypothesis that each data set is drawn from WL distribution. Such conclusions are also supported by various diagnostic plots in Figures 3 and 4. Based on the MLEs ĉ, ˆθ 1, ˆθ 2, the point estimate of R is and the 95% confidence interval of R is (0.5164, ) with confidence length MLE, parametric and non-parametric bootstrap estimates of R are provided in Table 2.
21 21 Table 1: MLEs and K-S and A-D goodness-of-fit tests. Plane MLEs K-S statistic p-value A-D statistic p-value Length 20 mm ĉ = ˆθ 1 = Length 50 mm ĉ = ˆθ 2 = Profile log likelihood of c c Figure 2: The profile log-likelihood of c for the given data. Table 2: Statistical inference of R for the different estimation methods. Estimation Method ˆR 95% C.I. of R confidence length MLE (0.5286, ) Parametric Bootstrap p-boot (0.5265, ) t-boot (0.5265, ) BC a -boot (0.5218, ) Non-Parametric Bootstrap p-boot (0.5364, ) t-boot (0.5364, ) BC a -boot (0.5328, )
22 22 Empirical and fitted densities QQ plot Density Sample Quantiles Failure Stress (in GPa) Fitted Weighted Lindley Quantiles Empirical and fitted CDFs PP plot CDF Empirical CDF Failure Stress (in GPa) Fitted Weighted Lindley CDF Figure 3: Diagnostic plots for the fitted weighted Lindley model for length 20 mm data.
23 23 Empirical and fitted densities QQ plot Density Sample Quantiles Failure Stress (in GPa) Fitted Weighted Lindley Quantiles Empirical and fitted CDFs PP plot CDF Empirical CDF Failure Stress (in GPa) Fitted Weighted Lindley CDF Figure 4: Diagnostic plots for the fitted weighted Lindley model for length 50 mm data.
24 24 8 Concluding Remarks In this paper we consider the statistical inference of R = P (X > Y ), where X and Y are independent weighted Lindley random variables with a common shape parameter. This probability is a measure of discrimination between two groups and has been studied quite extensively under different conditions by various authors. We investigate Maximum likelihood, p-, t-, and BC a -bootstrapping estimation methods (point and interval) of R and their performances are examined by extensive simulations. It is observed that the maximum likelihood method provides very satisfactory results both for point and interval estimation of R. An example is provided to illustrate these results. It is hoped that our investigation will be useful for researchers dealing with the kind of data considered in this paper. One natural question arises what will happen if the shape parameter c is different in two populations? In this case, the expression of R will be different, although it can be written in terms of hyper-geometric functions. The estimation of R and the construction of confidence interval of R can be performed along the same line. References [1] Al-Mutairi, D.K., Ghitany, M.E. and Kundu, D. (2013). Inference on stress-strength reliability from Lindley distribution, Communications in Statistics - Theory and Methods, (to appear) [2] Andrews, L.C. (1998). Special Functions of Mathematics for Engineers. Second edition. Oxford University Press, Oxford. [3] Bader, M.G. and Priest, A.M. (1982). Statistical aspects of fiber and bundle strength in hybrid composites, In Progress in Science and Engineering Composites, (T. Hayashi, K. Kawata, and S. Umekawa, eds.), vol. ICCM-IV, pp [4] Birnbaum, Z.W. (1956). On a use of Mann-Whitney statistics, In Proceedings of the third Berkley Symposium in Mathematics, Statistics and Probability, vol. 1, [5] Church, J.D. and Harris, B. (1970). The estimation of reliability from stress-strength relationships. Technometrics 12:49 54.
25 25 [6] Efron, B., Tibshirani, R.J. (1998). An Introduction to Bootstrap. New York: Chapman & Hall. [7] Erdélyi, A., Magnus, W., Oberhettinger, F. and Tricomi, F.G. (1954). Tables of Integral Transforms. McGraw-Hill, New York. [8] Ghitany, M.E., Atieh, B., Nadarajah, S. (2008). Lindley distribution and its application, Mathematics and Computers in Simulation 78: [9] Ghitany, M.E., Al-Qalaf, F., Al-Mutairi, and Husain, H.A. S. (2011). A two-parameter weighted Lindley distribution and its applications to survival data, Mathematics and Computers in Simulation 8: [10] Gupta, R.C., Ramakrishnan, S. and Zhou, X. (1999) Point and interval estimation of P(X < Y ): The normal case with common coefficient of variation. Annals of the Institute of Statistical Mathematics 51: [11] Gupta, R.C. and Brown, N. (2001) Reliability studies of the skew-normal distribution and its application to strength-stress models. Communications in Statistics, Theory and Methods 30: [12] Gupta, R.C. and Li, X. (2006) Statistical inference for the common mean of two lognormal distributions and some applications in reliability. Computational Statistics and Data Analysis 50: [13] Gupta, R.C. and Peng, C. (2009) Estimating reliability in proportional odds ratio models. Computational Statistics and Data Analysis 53: [14] Jiang L. and Wong, A.C.M. (2008). A note on inference for P (X < Y ) for right truncated exponentially distributed data, Statistical Papers 49: [15] Kim, C. and Chung, Y. (2006). Bayesian estimation of P (Y < X) from Burr-Type X model containing spurious observations, Statistical Papers 47: [16] Kotz, S., Lumelskii, Y. and Pensky, M. (2003). The Stress-Strength Model and its Generalizations: Theory and Applications. Singapore: World Scientific Press. [17] Krishnamoorthy, K., Mukherjee, S. and Guo, H. (2007). Inference on reliability in twoparameter exponential stress-strength model. Metrika 65: [18] Kundu, D. and Raqab, M.Z. (2009). Estimation of R = P (Y < X) for three-parameter Weibull distribution. Statistics and Probability Letters 79: [19] Kundu, D. and Gupta, R.D. (2007). A convenient way of generating gamma random variables. Computational Statistics and Data Analysis 51: [20] Kundu, D. and Gupta, R.D. (2006). Estimation of R = P [Y < X] for Weibull distributions. IEEE Transactions on Reliability 55: [21] Kundu, D. and Gupta, R.D. (2005). Estimation of P [Y < X] for generalized exponential distribution. Metrika 61:
26 26 [22] Lehmann, L.E. and Casella, G. (1998) Theory of Point Estimation, 2nd Edition, New York: Springer. [23] Lindley, D.V. (1958). Fudicial distributions and Bayes theorem. Journal of the Royal Statistical Society B 20: [24] Raqab, M.Z. and Kundu, D. (2005). Comparison of different estimators of P [Y < X] for a scaled Burr type X distribution. Communications in Statistics-Simulation and Computation 34: [25] Raqab. M.Z., Madi, M.T. and Kundu, D. (2008). Estimation of R = P (Y < X) for the 3-parameter generalized exponential distribution, Communications in Statistics - Theory and Methods 37: [26] Saraçoğlu, B., Kaya, M.F. (2007). Maximum likelihood estimation and confidence intervals of system reliability for Gompertz distribution in stress-strength models, Seluk Journal of Applied Mathematics 8: [27] Saraçoğlu, B., Kaya, M.F., Abd-Elfattah, A.M. (2009). Comparison of estimators for stress-strength reliability in Gompertz Case, Hacettepe Journal of Mathematics and Statistics 38: [28] Saraçoğlu, B., Kinaci, I. and Kundu, D. (2012). On estimation of R = P (Y < X) for exponential distribution under progressive type-ii censoring. Journal of Statistical Computation and Simulation 82: Table 3: Average bias (mean squared error) of different estimators of R = P (X > Y ) when c = 0.75, θ 1 = 2.5 and θ 2 = 1. n 1, n 2 MLE Parametric Boot Non-Parametric Boot 20, (0.0046) (0.0047) (0.0047) 20, (0.0037) (0.0038) (0.0037) 20, (0.0029) (0.0029) (0.0030) 30, (0.0037) (0.0039) (0.0037) 30, (0.0029) (0.0031) (0.0029) 30, (0.0023) (0.0023) (0.0024) 50, (0.0031) (0.0033) (0.0032) 50, (0.0023) (0.0023) (0.0023) 50, (0.0018) (0.0017) (0.0018)
27 27 Table 4: Average bias (mean squared error) of different estimators of R = P (X > Y ) when c = 2, θ 1 = 1 and θ 2 = 1. n 1, n 2 MLE Parametric Boot Non-Parametric Boot 20, (0.0078) (0.0079) (0.0080) 20, (0.0065) (0.0065) (0.0067) 20, (0.0055) (0.0051) (0.0056) 30, (0.0064) (0.0067) (0.0065) 30, (0.0050) (0.0052) (0.0052) 30, (0.0040) (0.0041) (0.0040) 50, (0.0052) (0.0056) (0.0053) 50, (0.0040) (0.0042) (0.0040) 50, (0.0030) (0.0030) (0.0030) Table 5: Average bias (mean squared error) of different estimators of R = P (X > Y ) when c = 4, θ 1 = 2 and θ 2 = 5. n 1, n 2 MLE Parametric Boot Non-Parametric Boot 20, (0.0016) (0.0015) (0.0016) 20, (0.0013) (0.0012) (0.0013) 20, (0.0010) (0.0011) (0.0010) 30, (0.0014) (0.0013) (0.0013) 30, (0.0011) (0.0010) (0.0011) 30, (0.0009) (0.0009) (0.0008) 50, (0.0010) (0.0011) (0.0010) 50, (0.0008) (0.0008) (0.0008) 50, (0.0006) (0.0007) (0.0006)
28 28 Table 6: Coverage probability (average confidence length) of different estimators of R = P (X > Y ) when c = 0.75, θ 1 = 2.5 and θ 2 = 1. n 1, n 2 MLE Parametric Bootstrap Non-parametric Bootstrap p-boot t-boot BC a -boot p-boot t-boot BC a -boot 20, (0.2473) (0.2546) (0.3050) (0.2635) (0.2496) (0.3026) (0.2566) 20, (0.2265) (0.2299) (0.2764) (0.2379) (0.2252) (0.2739) (0.2310) 20, (0.2048) (0.2075) (0.2474) (0.2142) (0.2003) (0.2432) (0.2057) 30, (0.2277) (0.2348) (0.2648) (0.2397) (0.2312) (0.2624) (0.2347) 30, (0.2044) (0.2086) (0.2392) (0.2129) (0.2048) (0.2358) (0.2085) 30, (0.1815) (0.1845) (0.2128) (0.1884) (0.1805) (0.2101) (0.1840) 50, (0.2090) (0.2163) (0.2288) (0.2174) (0.2113) (0.2234) (0.2119) 50, (0.1842) (0.1879) (0.2036) (0.1897) (0.1855) (0.2017) (0.1868) 50, (0.1599) (0.1621) (0.1792) (0.1640) (0.1597) (0.1770) (0.1614)
29 29 Table 7: Coverage probability (average confidence length) of different estimators of R = P (X > Y ) when c = 2, θ 1 = 1 and θ 2 = 1. n 1, n 2 MLE Parametric Bootstrap Non-parametric Bootstrap p-boot t-boot BC a -boot p-boot t-boot BC a -boot 20, (0.3172) (0.3361) (0.3391) (0.3380) (0.3355) (0.3398) (0.3365) 20, (0.2918) (0.3060) (0.3135) (0.3071) (0.3035) (0.3122) (0.3041) 20, (0.2690) (0.2791) (0.2904) (0.2795) (0.2756) (0.2880) (0.2759) 30, (0.2920) (0.3059) (0.3031) (0.3073) (0.3040) (0.3012) (0.3049) 30, (0.2635) (0.2740) (0.2760) (0.2749) (0.2724) (0.2743) (0.2728) 30, (0.2373) (0.2444) (0.2496) (0.2447) (0.2432) (0.2492) (0.2433) 50, (0.2693) (0.2786) (0.2713) (0.2797) (0.2747) (0.2661) (0.2755) 50, (0.2372) (0.2448) (0.2422) (0.2455) (0.2423) (0.2394) (0.2427) 50, (0.2069) (0.2121) (0.2132) (0.2124) (0.2107) (0.2117) (0.2109)
30 30 Table 8: Coverage probability (average confidence length) of different estimators of R = P (X > Y ) when c = 4, θ 1 = 2 and θ 2 = 5. n 1, n 2 MLE Parametric Bootstrap Non-parametric Bootstrap p-boot t-boot BC a -boot p-boot t-boot BC a -boot 20, (0.1641) (0.1473) (0.1305) (0.1731) (0.1426) (0.1254) (0.1646) 20, (0.1490) (0.1368) (0.1226) (0.1543) (0.1346) (0.1200) (0.1507) 20, (0.1325) (0.1262) (0.1155) (0.1363) (0.1235) (0.1124) (0.1330) 30, (0.1485) (0.1344) (0.1190) (0.1546) (0.1305) (0.1147) (0.1478) 30, (0.1331) (0.1235) (0.1100) (0.1377) (0.1208) (0.1071) (0.1337) 30, (0.1179) (0.1113) (0.1007) (0.1200) (0.1102) (0.0994) (0.1183) 50, (0.1327) (0.1227) (0.1094) (0.1370) (0.1191) (0.1054) (0.1319) 50, (0.1170) (0.1097) (0.0980) (0.1203) (0.1077) (0.0955) (0.1175) 50, (0.1020) (0.0981) (0.0884) (0.1047) (0.0960) (0.0862) (0.1023)
inferences 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 informationExact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring
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
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 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 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 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 informationAnalysis 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 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 informationA Quasi Gamma Distribution
08; 3(4): 08-7 ISSN: 456-45 Maths 08; 3(4): 08-7 08 Stats & Maths www.mathsjournal.com Received: 05-03-08 Accepted: 06-04-08 Rama Shanker Department of Statistics, College of Science, Eritrea Institute
More informationOptimal Unbiased Estimates of P {X < Y } for Some Families of Distributions
Metodološki zvezki, Vol., No., 04, -9 Optimal Unbiased Estimates of P {X < } for Some Families of Distributions Marko Obradović, Milan Jovanović, Bojana Milošević Abstract In reliability theory, one of
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 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 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 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 informationGeneralized Exponential Distribution: Existing Results and Some Recent Developments
Generalized Exponential Distribution: Existing Results and Some Recent Developments Rameshwar D. Gupta 1 Debasis Kundu 2 Abstract Mudholkar and Srivastava [25] introduced three-parameter exponentiated
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 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 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 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 informationDifferent methods of estimation for generalized inverse Lindley distribution
Different methods of estimation for generalized inverse Lindley distribution Arbër Qoshja & Fatmir Hoxha Department of Applied Mathematics, Faculty of Natural Science, University of Tirana, Albania,, e-mail:
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 informationAn Extension of the Generalized Exponential Distribution
An Extension of the Generalized Exponential Distribution Debasis Kundu and Rameshwar D. Gupta Abstract The two-parameter generalized exponential distribution has been used recently quite extensively to
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 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 informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2008 Prof. Gesine Reinert 1 Data x = x 1, x 2,..., x n, realisations of random variables X 1, X 2,..., X n with distribution (model)
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 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 informationthe Presence of k Outliers
RESEARCH ARTICLE OPEN ACCESS On the Estimation of the Presence of k Outliers for Weibull Distribution in Amal S. Hassan 1, Elsayed A. Elsherpieny 2, and Rania M. Shalaby 3 1 (Department of Mathematical
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 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 informationOn the estimation of stress strength reliability parameter of inverted exponential distribution
International Journal of Scientific World, 3 (1) (2015) 98-112 www.sciencepubco.com/index.php/ijsw c Science Publishing Corporation doi: 10.14419/ijsw.v3i1.4329 Research Paper On the estimation of stress
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 informationMoments of the Reliability, R = P(Y<X), As a Random Variable
International Journal of Computational Engineering Research Vol, 03 Issue, 8 Moments of the Reliability, R = P(Y
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 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 informationA TEST OF FIT FOR THE GENERALIZED PARETO DISTRIBUTION BASED ON TRANSFORMS
A TEST OF FIT FOR THE GENERALIZED PARETO DISTRIBUTION BASED ON TRANSFORMS Dimitrios Konstantinides, Simos G. Meintanis Department of Statistics and Acturial Science, University of the Aegean, Karlovassi,
More informationEstimation of the Stress-Strength Reliability for the Dagum Distribution
156 Journal of Advanced Statistics, Vol 1, No 3, September 2016 Estimation of the Stress-Strength Reliability for the Dagum Distribution Ber Al-Zahrani Samerah Basloom Department of Statistics, King Abdulaziz
More informationAn Extension of the Freund s Bivariate Distribution to Model Load Sharing Systems
An Extension of the Freund s Bivariate Distribution to Model Load Sharing Systems G Asha Department of Statistics Cochin University of Science and Technology Cochin, Kerala, e-mail:asha@cusat.ac.in Jagathnath
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 informationA THREE-PARAMETER WEIGHTED LINDLEY DISTRIBUTION AND ITS APPLICATIONS TO MODEL SURVIVAL TIME
STATISTICS IN TRANSITION new series, June 07 Vol. 8, No., pp. 9 30, DOI: 0.307/stattrans-06-07 A THREE-PARAMETER WEIGHTED LINDLEY DISTRIBUTION AND ITS APPLICATIONS TO MODEL SURVIVAL TIME Rama Shanker,
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 informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
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 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 informationStatistics Ph.D. Qualifying Exam: Part I October 18, 2003
Statistics Ph.D. Qualifying Exam: Part I October 18, 2003 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. 1 2 3 4 5 6 7 8 9 10 11 12 2. Write your answer
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 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 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 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 informationMaximum Likelihood and Bayes Estimations under Generalized Order Statistics from Generalized Exponential Distribution
Applied Mathematical Sciences, Vol. 6, 2012, no. 49, 2431-2444 Maximum Likelihood and Bayes Estimations under Generalized Order Statistics from Generalized Exponential Distribution Saieed F. Ateya Mathematics
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 informationEstimation of Quantiles
9 Estimation of Quantiles The notion of quantiles was introduced in Section 3.2: recall that a quantile x α for an r.v. X is a constant such that P(X x α )=1 α. (9.1) In this chapter we examine quantiles
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 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 informationInference about Reliability Parameter with Underlying Gamma and Exponential Distribution
Florida International University FIU Digital Commons FIU Electronic Theses and Dissertations University Graduate School 9-3-211 Inference about Reliability Parameter with Underlying Gamma and Exponential
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 informationIntroduction of Shape/Skewness Parameter(s) in a Probability Distribution
Journal of Probability and Statistical Science 7(2), 153-171, Aug. 2009 Introduction of Shape/Skewness Parameter(s) in a Probability Distribution Rameshwar D. Gupta University of New Brunswick Debasis
More informationPROPERTIES AND DATA MODELLING APPLICATIONS OF THE KUMARASWAMY GENERALIZED MARSHALL-OLKIN-G FAMILY OF DISTRIBUTIONS
Journal of Data Science 605-620, DOI: 10.6339/JDS.201807_16(3.0009 PROPERTIES AND DATA MODELLING APPLICATIONS OF THE KUMARASWAMY GENERALIZED MARSHALL-OLKIN-G FAMILY OF DISTRIBUTIONS Subrata Chakraborty
More informationUNIVERSITY OF CALGARY. A New Hybrid Estimation Method for the. Generalized Exponential Distribution. Shan Zhu A THESIS
UNIVERSITY OF CALGARY A New Hybrid Estimation Method for the Generalized Exponential Distribution by Shan Zhu A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
More informationChapter 4 HOMEWORK ASSIGNMENTS. 4.1 Homework #1
Chapter 4 HOMEWORK ASSIGNMENTS These homeworks may be modified as the semester progresses. It is your responsibility to keep up to date with the correctly assigned homeworks. There may be some errors in
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 informationEstimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators
Estimation theory Parametric estimation Properties of estimators Minimum variance estimator Cramer-Rao bound Maximum likelihood estimators Confidence intervals Bayesian estimation 1 Random Variables Let
More informationBeta-Linear Failure Rate Distribution and its Applications
JIRSS (2015) Vol. 14, No. 1, pp 89-105 Beta-Linear Failure Rate Distribution and its Applications A. A. Jafari, E. Mahmoudi Department of Statistics, Yazd University, Yazd, Iran. Abstract. We introduce
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 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 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 informationBy Godase, Shirke, Kashid. Published: 26 April 2017
Electronic Journal of Applied Statistical Analysis EJASA, Electron. J. App. Stat. Anal. http://siba-ese.unisalento.it/index.php/ejasa/index e-issn: 2070-5948 DOI: 10.1285/i20705948v10n1p29 Tolerance intervals
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 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 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 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 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 informationThe Marshall-Olkin Flexible Weibull Extension Distribution
The Marshall-Olkin Flexible Weibull Extension Distribution Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash arxiv:169.8997v1 math.st] 25 Sep 216 Department of Mathematics, Faculty of Science,
More informationBivariate Geometric (Maximum) Generalized Exponential Distribution
Bivariate Geometric (Maximum) Generalized Exponential Distribution Debasis Kundu 1 Abstract In this paper we propose a new five parameter bivariate distribution obtained by taking geometric maximum of
More informationBivariate Sinh-Normal Distribution and A Related Model
Bivariate Sinh-Normal Distribution and A Related Model Debasis Kundu 1 Abstract Sinh-normal distribution is a symmetric distribution with three parameters. In this paper we introduce bivariate sinh-normal
More informationAnalysis of Gamma and Weibull Lifetime Data under a General Censoring Scheme and in the presence of Covariates
Communications in Statistics - Theory and Methods ISSN: 0361-0926 (Print) 1532-415X (Online) Journal homepage: http://www.tandfonline.com/loi/lsta20 Analysis of Gamma and Weibull Lifetime Data under a
More informationPractice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:
Practice Exam 1 1. Losses for an insurance coverage have the following cumulative distribution function: F(0) = 0 F(1,000) = 0.2 F(5,000) = 0.4 F(10,000) = 0.9 F(100,000) = 1 with linear interpolation
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 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 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 informationA New Method for Generating Distributions with an Application to Exponential Distribution
A New Method for Generating Distributions with an Application to Exponential Distribution Abbas Mahdavi & Debasis Kundu Abstract Anewmethodhasbeenproposedtointroduceanextraparametertoafamilyofdistributions
More informationEstimation of P (Y < X) in a Four-Parameter Generalized Gamma Distribution
AUSTRIAN JOURNAL OF STATISTICS Volume 4 202, Number 3, 97 20 Estimation of P Y < X in a Four-Parameter Generalized Gamma Distribution M. Masoom Ali, Manisha Pal 2 and Jungsoo Woo 3 Department of Mathematical
More informationProblem Selected Scores
Statistics Ph.D. Qualifying Exam: Part II November 20, 2010 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. Problem 1 2 3 4 5 6 7 8 9 10 11 12 Selected
More informationOn Five Parameter Beta Lomax Distribution
ISSN 1684-840 Journal of Statistics Volume 0, 01. pp. 10-118 On Five Parameter Beta Lomax Distribution Muhammad Rajab 1, Muhammad Aleem, Tahir Nawaz and Muhammad Daniyal 4 Abstract Lomax (1954) developed
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 informationReview. December 4 th, Review
December 4 th, 2017 Att. Final exam: Course evaluation Friday, 12/14/2018, 10:30am 12:30pm Gore Hall 115 Overview Week 2 Week 4 Week 7 Week 10 Week 12 Chapter 6: Statistics and Sampling Distributions Chapter
More informationUNIVERSITÄT POTSDAM Institut für Mathematik
UNIVERSITÄT POTSDAM Institut für Mathematik Testing the Acceleration Function in Life Time Models Hannelore Liero Matthias Liero Mathematische Statistik und Wahrscheinlichkeitstheorie Universität Potsdam
More informationarxiv: v1 [stat.ap] 31 May 2016
Some estimators of the PMF and CDF of the arxiv:1605.09652v1 [stat.ap] 31 May 2016 Logarithmic Series Distribution Sudhansu S. Maiti, Indrani Mukherjee and Monojit Das Department of Statistics, Visva-Bharati
More informationMaster s Written Examination
Master s Written Examination Option: Statistics and Probability Spring 016 Full points may be obtained for correct answers to eight questions. Each numbered question which may have several parts is worth
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 informationTHE FIVE PARAMETER LINDLEY DISTRIBUTION. Dept. of Statistics and Operations Research, King Saud University Saudi Arabia,
Pak. J. Statist. 2015 Vol. 31(4), 363-384 THE FIVE PARAMETER LINDLEY DISTRIBUTION Abdulhakim Al-Babtain 1, Ahmed M. Eid 2, A-Hadi N. Ahmed 3 and Faton Merovci 4 1 Dept. of Statistics and Operations Research,
More informationJournal of Environmental Statistics
jes Journal of Environmental Statistics February 2010, Volume 1, Issue 3. http://www.jenvstat.org Exponentiated Gumbel Distribution for Estimation of Return Levels of Significant Wave Height Klara Persson
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 information1. Fisher Information
1. Fisher Information Let f(x θ) be a density function with the property that log f(x θ) is differentiable in θ throughout the open p-dimensional parameter set Θ R p ; then the score statistic (or score
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 informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
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 informationPoisson-Size-biased Lindley Distribution
International Journal of Scientific and Research Publications, Volume 4, Issue 1, January 2014 1 Poisson-Size-biased Lindley Distribution Mr. Tanka Raj Adhikari and Prof. R.S. Srivastava Department of
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 information