UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS
|
|
- Arline Powers
- 5 years ago
- Views:
Transcription
1 ISSN UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS THE POISSON-EXPONENTIAL SURVIVAL DISTRIBUTION FOR LIFETIME DATA Vicente G. Cancho Francisco Louzada-Neto Gladys D.C. Barriga RELATÓRIO TÉCNICO DEPARTAMENTO DE ESTATÍSTICA TEORIA E MÉTODO SÉRIE A Maio/2010 nº 217
2 The Poisson-Exponential Survival Distribution For Lifetime Data Vicente G. Cancho a, Francisco Louzada-Neto b and Gladys D.C. Barriga b a Universidade de São Paulo b Universidade de Federal de São Carlos Abstract In this paper we proposed a new two-parameters lifetime distribution with increasing failure rate. The new distribution arises on a latent complementary risk problem base. The properties of the proposed distribution are discussed, including a formal prove of its probability density function and explicit algebraic formulaes for its reliability and failure rate functions, quantiles and moments, including the mean and variance. A simple EM-type algorithm for iteratively computing maximum likelihood estimates is presented. The Fisher information matrix is derived analytically in order to obtaining the asymptotic covariance matrix. The methodology is illustrated on a real data set. Keywords: Complementary Risks, EM Algorithm, Exponential Distribution, Poisson Distribution, Survival Analysis. 1 Introduction The exponential distribution (ED) provides a simple, elegant and close form solution to many problems in lifetime testing and reliability studies. However, the ED does not provide a reasonable parametric fit for some practical applications where the underlying hazard rates are nonconstant, presenting monotone shapes. In recent years, in order to overcame such problem, new classes of models were introduced based on modifications of the ED. Gupta & Kundu (1999) proposed a generalized ED. This extended family can accommodate data with increasing and decreasing failure rate function. Kus (2007) proposed another modification of the ED with decreasing failure rate function. While Barreto-Souza & Cribari-Neto (2009) generalizes the distribution proposed by Kus (2007) by including a power parameter in his distribution. In this paper, we propose a new distribution family based on the ED with increasing failure rate function. Its genesis is stated on a complementary risk problem base (Basu & Klein, 1982) in presence of latent risks, in the sense that there is no information about which factor was 1
3 responsible for the component failure and only the maximum lifetime value among all risks is observed. The paper is organized as follows. In Section 2, we introduce the new distribution and present its properties. Section 3 outlines an EM-type algorithm for maximum likelihood estimation. In Section 4 the methodology is illustrated in a real data set. Some final comments are presented in Section 5. 2 The Poisson-Exponential Distribution Let Y be a nonnegative random variable denoting the lifetime of an component in some population. The random variable Y is said to have a Poisson-Exponential distribution (PED) with parameters λ > 0 and θ > 0 if its probability density function is given by f(y) = θλe λy θe λy 1 e θ, y > 0. (1) The parameter λ controls the scale of the distribution while the parameter θ controls its shape. As θ approaches zero, the PED converges to an ED with parameter λ. Figure 1 (top panel) shows the probability density function (1) for θ = 0.1, 1, 2, 4, 8, which is decreasing if 0 < θ < 1 and unimodal for θ 1. The modal value λe 1 is obtained at y = log θ/λ. Density θ=0.1 θ=1.0 θ=2.0 θ=4.0 θ=8.0 Hazard function θ=0.1 θ=1.0 θ=2.0 θ=4.0 θ= y y Figure 1: Top panel: probability density function of the PE distribution. Bottom panel: failure rate function of the PED. We fixed λ = 1. The reliability function of a PED random varible is given by S(y) = 1 e θe λy 1 e θ, y > 0. (2) 2
4 From (1) and (2) it is easy verify that the failure rate function is given by h(y) = θλ exp( λy θe λy) 1 exp( θe λy, y > 0. (3) ) The failure rate function (3) is increasing. Figure 1 (bottom panel) shows some failure rate function shapes for same fixed values of θ. The initial and long-term failure rate function values are both finite and are given by h(0) = λθ(e θ 1) and h( ) = λ. We simulate the PED considering Q(u) = F 1 (u) = log ( ) θ log(u(e θ 1)) λ 1, where u has the uniform U(0, 1) distribution and F(y) = 1 S(y) is distribution function of Y. 2.1 Genesis Complementary risks problems (Basu & Klein, 1982) arise in several areas, such as public health, actuarial science, biomedical studies, demography and industrial reliability. In the classical Complementary risks scenarios the lifetime associated with a particular risk is not observable, rather we observe only the maximum lifetime value among all risks. Simplistically, in reliability, we observe only the maximum component lifetime of a parallel system. That is, the observable quantities for each component are the maximum lifetime value to failure among all risks, and the cause of failure. Full statistical procedures and extensive literature are available to deal with Complementary risks problems and interested readers can refer to Lawless (2003), Crowder et al. (1991) and Cox & Oakes (1984). A difficulty arises if the risks are latent in the sense that there is no information about which factor was responsible for the component failure, which can be often observed in field data. We call these latent Complementary risks data. On many occasions this information is not available or it is impossible that the true cause of failure is specified by an expert. In reliability, the components can be totally destroyed in the experiment. Further, the true cause of failure can be masked from our view. In modular systems, the need to keep a system running means that a module that contains many components can be replaced without the identification of the exact failing component. Goetghebeur & Ryan (1995) addressed the problem of assessing covariate effects based on a semi-parametric proportional hazards structure for each failure type when the failure type is unknown for some individuals. Reiser et al. (1995) considered statistical procedures for analysing masked data, but their procedure can not be applied when all observations have an unknown cause of failure. Lu & Tsiatis (2001) presents a multiple imputation method for estimating regression coefficients for risks modeling with missing cause of failure. A comparison of two partial likelihood approaches for risks modeling with missing cause of failure is presented in Lu & Tsiatis (2005). 3
5 Our model can be derived as follows. Let M be a random variable denoting the number of Complementary risks related to the occurrence of an event of interest. Further, assume that M has a zero truncated Poisson distribution with probability mass function given by P(M = m) = e θ θ m m!(1 e θ, m = 1, 2..., θ > 0. (4) ) Let T j (j = 1, 2,...) denote the time-to-event due to the j-th Complementary risks, hereafter lifetime. Given M = m, the random variables T j, j = 1,, m are assumed to be independent and identically distributed according to a ED with a common distribution function with probability density function given by f(t; λ) = λ exp( λt), t > 0, λ > 0. (5) In the latent Complementary risks scenario, the number of causes M and the lifetime T j associated with a particular cause are not observable (latent variables), but only the maximum lifetime Y among all causes is usually observed. So, the component lifetime is defined as Y = max (T 1, T 2...,T M ). (6) The following result shows that the random variable Y have probability density function given by (1). Proposition 2.1 If the random variable Y is defined as in (6), then, considering (4) and (5), Y is distributed according to a PED with probability density function given by (1). Proof 2.1 The conditional density function of (6) given M = m is given by f(y M = m, λ) = mλ[1 e λt ] m 1 e λt, t > 0, m = 1,... Them, the marginal probability density function of Y is given by f(y) = mλ[1 e λt ] m 1 e λt θm e θ m!(1 e θ ) m=1 = θλe θ λt [θ(1 e λt )] m 1 1 e θ (m 1)! This completes the proof. m=1 = θλe λy θe λy 1 e θ. The PED parameters have a direct interpretation in terms of Complementary risks. The θ represents the mean of the number of Complementary risks, while λ denotes the failure rate. In comparison with the Kus (2007) formulation we follows an opposite way, since he defines the component lifetime as Y = min(t 1, T 2...,T M ), while we are considering Y = max (T 1, T 2..., T M as stated in (6). 4
6 2.2 Moments Some of the most important features and characteristics of a distribution can be studied through its moments, such mean, variance, tending, dispersion skewness and kurtosis. A general expression for rth ordinary moment µ r = E(Y r ) of the PED can be obtained analytically considering the generalized hypergeometric function denoted by F p,q (a, b, θ) and defined as F p,q (a, b, θ) = j=0 θ j p Γ(a i + j)γ(a i ) 1 Γ(j + 1) q Γ(b i + j)γ(b i ) 1, (7) where a = [a 1,...,a p ], p is the number of operands of a, b = [b 1,...,b q ], and q is the number of operands of b. We then formulate the following result. Proposition 2.2 Suppose Y denotes a random variable from a PED with parameters λ > 0 and θ > 0 and probability density function given by (1), then µ r = θγ(r + 1) λ r (1 e θ ) F r+1,r+1([1,...,1], [2,...,2], θ), (8) where F p,q (a, b, θ) is the generalized hypergeometric function. The proof of 2.2 is obtained by direct integration and it is then omitted. Considering (8), the mean and variance of the of the PED are given, respectively, by 3 Inference θ E(Y ) = λ(1 e θ ) F 2,2([1, 1], [2, 2], θ), [ θ V ar(y ) = λ 2 (1 e θ F 3,3 ([1, 1, 1], [2, 2, 2], θ) ) ] θ (1 e θ ) F 2,2([1, 2 1], [2, 2], θ). 3.1 Estimation by maximum likelihood The log-likelihood function based on an observed sample of size n, y = (y 1,...,y n ), from a PED is given by l(θ) = n log(θλ) λ y i θ e λy i n log(1 e θ ). (9) The MLEs of θ = (θ, λ) can be directly obtained by maximizing the log-likelihood function 5
7 (9) or, alternatively, by finding the solution for the following two nonlinear equations, l(θ) θ l(θ) λ = n θ n e θ 1 e λy i = 0, (10) = n λ y i + θ y i e λy i = 0. (11) The equation (11) can be solved exactly for θ, i.e., ˆθ = n y i n λ 1 n y ie ˆλy i, (12) conditional upon the value of ˆλ, where ˆθ and ˆλ are the MLEs for θ and λ, respectively. Therefore, it is sufficient to solve the equation (10) iteratively in order to find the MLE for λ. 3.2 An EM-algorithm In this subsection we develop an EM-algorithm (Dempster et al., 1977) for obtain the MLEs of the PED parameters. Let y = (y 1,...,y n ) and m = (m 1,...,m n ). It follows that the complete log-likelihood function associated with (y,m) is given by l c (θ y,m) = log(f(y i, m i θ)), (13) where f(y i, m i θ) = m i θ m i λ[1 e λy i ] m i 1 e θ λy i /Γ(m i + 1)(1 e θ ), m i = 1, 2,..., y i > 0. Considering m i = E[M i θ = θ, y i ] we can direct to the implementation of the M-step, which consists in maximizing the expected complete data function or the Q function over θ, given by Q(θ θ (k) ) = E[l c (θ) y, θ (k) ] = c + n log λ + ( m i 1)log(1 e λy i ) λ y i + m i log θ n log(e θ 1), where c is a constant that is independent of θ, m i = θ(1 e λy i ) and θ (k) is an updated value of θ. Then, the E-step and the M-step of the algorithm are given by: E-step: Given a current estimate θ (k), compute m (k) i 6
8 M-step: Update θ (k) by maximizing Q(θ θ (k) ) over θ, which leads to the following nice expressions θ (k+1) = λ (k+1) = 3.3 Interval Estimation m (k) i (1 e θ (k+1) ), n n ( y i n [ ( m (k) i 1)y i e λ (k+1) y i 1 e λ (k+1) y i Large sample inference for the parameters can be based, in principle, on the MLEs and their estimated standard errors. Following (Cox & Hinkley, 1974), under suitable regularity conditions, it can be shown that ) ( n ( θ θ N 2 0,I 1 (θ) ), ]). where I(θ) is Fisher information matrix, i.e, ( ) I(θ) = E 2 l(θ) θ 2 ) E ( 2 l(θ) θ λ E E ( 2 l(θ) θ λ ) ( ) 2 l(θ) λ 2, where the second derivatives are given by 2 l(θ)/ θ 2 = n/θ 2 + ne θ /(e θ 1) 2, 2 l(θ)/ θ λ = n y ie λy i and 2 l(θ)/ λ 2 = nλ 2 θ n y2 i e λy i. Them, the elements of the Fisher Information matrix are derived as, ( 2 ) l(θ) nθ neθ E θ 2 = 2 (e θ 1) 2, ( 2 ) l(θ) nθ E = θ λ 4λ(1 e θ ) F 2,2([2, 2], [3, 3], θ), ( 2 ) l(θ) nθ 2 E λ 2 = 4λ 2 (1 e θ ) F 3,3([2, 2, 2], [3, 3, 3], θ) + n λ 2. 4 Application In this section we reanalyze the data set extracted from (Lawless, 2003). The lifetimes are the number of million revolutions before failure for each one of the 23 ball bearings on an endurance test of deep groove ball bearings. Firstly, in order to identify the shape of a lifetime data failure rate function we shall consider a graphical method based on the TTT plot (Aarset, 1987). In its empirical version the TTT plot is given by G(r/n) = [( r Y i:n) (n r)y r:n ]/( r Y i:n), 7
9 where r = 1,...,n and Y i:n represent the order statistics of the sample. It has been shown that the failure rate function is increasing (decreasing) if the TTT plot is concave (convex). Although, the TTT plot is only a sufficient condition, not a necessary one for indicating the failure rate function shape, it is used here as a crude indicative of its shape. Figure 2 (top panel) shows the TTT plot for the considered data, which is concave indicating an increasing failure rate function, which can be properly accommodated by a PED. G(r n) Survival Function K M Poisson Exponential Exponential r/n Time Figure 2: Top panel: Empirical scaled TTT-Transform for the data. Bottom panel: Kaplan- Meier curve with estimated RFs of the PED and ED. Then, the PED was fitted to the data. In order to obtain the MLEs of the parameters we used the EM-algorithm, considering as initial guess θ = 1 and λ = 1/ȳ (a moment estimator for λ from an independent identically distributed exponential observation). Of course, this choice is not foolproof and it is advisable to run the method several times from different starting values. Table 1 presents the MLEs and the corresponding 95% confidence intervals of the PED parameters, which were based on the Fisher information matrix at the MLEs, given by, [ ] I( θ) = For comparison of nested models, which is the case when comparing the PED with the ED, we can compute the maximum values of the unrestricted and restricted log-likelihoods to obtain the LRS. For testing H 0 : θ = 0 versus H 1 : θ > 0 we consider the LRS, w n = 2(l PE l E ), where l ED and l PED are the log-likelihoods for the model under the restricted hipothesis H 0 and under the unrestricted hipothesis H 1 under a sample of size n. Taking into account that the test is 8
10 Table 1: MLEs and their corresponding 95% confidence intervals. Parameter Dsitribution θ λ l( ) PED ( ; ) (0.0241; ) ED ( ; ) performed in the boundary of the parameter space, following Maller & Zhou (1995), the LRS, w n, is assumed to be asymptotically distributed as a symmetric mixture of a chi-squared distribution with one degree of freedom and a point-mass at zero. Then, lim n P(w n c) = 1/2 + 1/2 P(χ 2 1 c), where P(χ2 1 c) denotes a random variable with a chi-square distribution with one degree of freedom. Large positive values of w n give favourable evidence to the full model. The last column of Table 1 presents the log-likelihood values for both models. Thus, w n is equals to , with a p-value < , which is a strong evidence in favour of the PED. Also, we compare the PED and ED by inspection of the Akaike s information criterion (AIC), 2l( θ) + 2q, and Schawarz s Bayesian information criterion (BIC), 2l( θ) + q log(n), where q is the number of parameters in the model and n is size sample. The preferred model is the one with the smaller value on each criterion. The PED overcome the corresponding ED in the two considered criterions. The estimated statistics AIC and BIC for the PED are equal to and , respectively. While, the estimated statistics AIC and BIC for the ED are equal to and , respectively. These results are corroborated by the plot in the bottom panel of Figure 2, which compares the estimated rate function obtained via ED and PED fitting on the empirical Kaplan-Meier rate function base. 5 Concluding remarks In this paper we proposed the PED, as a possible extension to the well known ED. We provide a mathematical treatment of this distribution including a formal prove of its probability density function and explicit algebraic formulae for its rate function and failure rate function, quantiles and moments (particularly, for the mean and variance). The parameters of the PED has a direct interpretation in terms of Complementary risks. We discussed MLE, providing an EM-algorithm. The PED allows a straightforwardly nested hypothesis testing procedure for comparison with its ED particular case. The practical relevance and applicability of the PED were demonstrated in a real data set. 9
11 Acknowledgments: The research of Francisco Louzada-Neto is funded by the Brazilian organization CNPq. References Aarset, M. V. (1987). How to identify a bathtub hazard rate. IEEE Transactions on Reliability, 2(1), Barreto-Souza, W. & Cribari-Neto, F. (2009). A generalization of the exponential-poisson distribution. Statistics and Probability Letters, 79, Basu, A. & Klein, J. (1982). Some recent development in competing risks theory. In Survival Analysis, Edited by Crowley, J. and Johnson, R.A., Hayward: IMS, 1(1), Cox, D. & Oakes, D. (1984). Analysis of Survival Data. Chapman and Hall, London. Cox, D. R. & Hinkley, D. V. (1974). Theoretical statistics. Chapman and Hall, London. Crowder, M., Kimber, A., Smith, R. & Sweeting, T. (1991). Statistical Analysis of Reliability Data. Chapman and Hall, London. Dempster, A., Laird, N. & Rubin, D. (1977). Maximum likelihood from incomplete data via the EM algorithm. Journal of the Royal Statistical Society, Series B,, 39, Goetghebeur, E. & Ryan, L. (1995). A modified log rank test for competing risks with missing failure type. Biometrika, 77(2), Gupta, R. & Kundu, D. (1999). Generalized exponential distributions. Australian and New Zealand Journal of Statistics, 41(2), Kus, C. (2007). A new lifetime distribution distributions. Computational Statistics & Data analysis, 11, Lawless, J. F. (2003). Statistical Models and Methods for Lifetime Data. Wiley, New York, NY, second edition. Lu, K. & Tsiatis, A. A. (2001). Multiple imputation methods for estimating regression coefficients in the competing risks model with missing cause of failure. Biometrics, 57(4), Lu, K. & Tsiatis, A. A. (2005). Comparison between two partial likelihood approaches for the competing risks model with missing cause of failure. Lifetime Data Analysis, 11(1),
12 Maller, R. A. & Zhou, S. (1995). Testing for the presence of immune or cured individuals in censored survival data. Biometrics, 51(4), Reiser, B., Guttman, I., Lin, D., Guess, M. & Usher, J. (1995). Bayesian inference for masked system lifetime data. Applied Statistics, 44(1),
13 PUBLICAÇÕES RODRIGUES, J.; CASTRO, M.; BLAKRISHNAN, N.; GARIBAY, V.; Destructive weighted Poisson cure rate models Novembro/2009 Nº 210. COBRE, J.; LOUZADA-NETO, F., PERDONÁ, G.; A Bayesian Analysis for the Generalized Negative Binomial Weibull Cure Fraction Survival Model: Estimating the Lymph Nodes Metastasis Rates Janeiro/2010 Nº 211. DINIZ, C. A. R. ; LOUZADA-NETO, F.; MORITA, L. H. M.; The Multiplicative Heteroscedastic Von Bertalan_y Model Fevereiro/2010 Nº 212. DINIZ, C. A. R. ; MORITA, L. H. M, LOUZADA-NETO, F.; Heteroscedastic Von Bertalanffy Growth Model and an Application to a Kubbard female chicken corporeal weight growth data Fevereiro/2010 Nº 213. FURLAN, C. P. R.; DINIZ, C. A. R.; FRANCO, M. A. P.; Estimation of Lag Length in Distributed Lag Models: A Comparative Study Março/2010 Nº 214. RODRIGUES, J., CANCHO, V. G., CASTRO, M., BALAKRISHNAN, N., A Bayesian destructive weighted Poisson cure rate model and an application to a cutaneous melanoma data, Março/2010 Nº CORDEIRO, G. M.; RODRIGUES, J.; CASTRO, M. The exponential COM-Poisson distribution, Abril/2010 Nº 216 Mais informações sobre publicações anteriores ao ano de 2009 podem ser obtidas via dfln@ufscar.br
UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS
ISSN 0104-0499 UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS THE COMPLEMENTARY EXPONENTIAL GEOMETRIC DISTRIBUTION: MODEL, PROPERTIES AND A COMPARISON WITH ITS COUNTERPART
More informationUNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS
ISSN 0104-0499 UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS A COMPOUND CLASS OF EXPONENTIAL AND POWER SERIES DISTRIBUTIONS WITH INCREASING FAILURE RATE José Flores Delgado
More informationThe complementary exponential power series distribution
Brazilian Journal of Probability and Statistics 2013, Vol. 27, No. 4, 565 584 DOI: 10.1214/11-BJPS182 Brazilian Statistical Association, 2013 The complementary exponential power series distribution José
More informationUNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS
ISSN 14-499 UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS THE EXPONENTIATED EXPONENTIAL-GEOMETRIC DISTRIBUTION: A DISTRIBUTION WITH DECREASING, INCREASING AND UNIMODAL HAZARD
More informationUNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS
ISSN 0104-0499 UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS THE EXPONENTIAL NATIVE BINOMIAL DISTRIBUTION Francisco Louzada-Neto Patrick Borges RELATÓRIO TÉCNICO DARTAMENTO
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 informationThe Poisson-Weibull Regression Model
Chilean Journal of Statistics Vol. 8, No. 1, April 2017, 25-51 Research Paper The Poisson-Weibull Regression Model Valdemiro Piedade Vigas 1, Giovana Oliveira Silva 2,, and Francisco Louzada 3 1 Instituto
More informationA COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky
A COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky Empirical likelihood with right censored data were studied by Thomas and Grunkmier (1975), Li (1995),
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 informationBayesian reference analysis for the Poisson-exponential lifetime distribution
Chilean Journal of Statistics Vol. 4, No. 1, April 213, 99 113 Bayesian Statistics Research Paper Bayesian reference analysis for the Poisson-exponential lifetime distribution Vera L.D. Tomazella 1, Vicente
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 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 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 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 informationTHE WEIBULL GENERALIZED FLEXIBLE WEIBULL EXTENSION DISTRIBUTION
Journal of Data Science 14(2016), 453-478 THE WEIBULL GENERALIZED FLEXIBLE WEIBULL EXTENSION DISTRIBUTION Abdelfattah Mustafa, Beih S. El-Desouky, Shamsan AL-Garash Department of Mathematics, Faculty of
More informationStatistical Analysis of Competing Risks With Missing Causes of Failure
Proceedings 59th ISI World Statistics Congress, 25-3 August 213, Hong Kong (Session STS9) p.1223 Statistical Analysis of Competing Risks With Missing Causes of Failure Isha Dewan 1,3 and Uttara V. Naik-Nimbalkar
More informationModel Selection for Semiparametric Bayesian Models with Application to Overdispersion
Proceedings 59th ISI World Statistics Congress, 25-30 August 2013, Hong Kong (Session CPS020) p.3863 Model Selection for Semiparametric Bayesian Models with Application to Overdispersion Jinfang Wang and
More informationBivariate Weibull-power series class of distributions
Bivariate Weibull-power series class of distributions Saralees Nadarajah and Rasool Roozegar EM algorithm, Maximum likelihood estimation, Power series distri- Keywords: bution. Abstract We point out that
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 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 informationComputational Statistics and Data Analysis
Computational Statistics and Data Analysis 53 (2009) 4482 4489 Contents lists available at ScienceDirect Computational Statistics and Data Analysis journal homepage: www.elsevier.com/locate/csda A log-extended
More informationMultistate Modeling and Applications
Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)
More informationA COMPARISON OF COMPOUND POISSON CLASS DISTRIBUTIONS
The International Conference on Trends and Perspectives in Linear Statistical Inference A COMPARISON OF COMPOUND POISSON CLASS DISTRIBUTIONS Dr. Deniz INAN Oykum Esra ASKIN (PHD.CANDIDATE) Dr. Ali Hakan
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 informationInference for the dependent competing risks model with masked causes of
Inference for the dependent competing risks model with masked causes of failure Radu V. Craiu and Benjamin Reiser Abstract. The competing risks model is useful in settings in which individuals/units may
More informationMax-Erlang and Min-Erlang power series distributions as two new families of lifetime distribution
BULETINUL ACADEMIEI DE ŞTIINŢE A REPUBLICII MOLDOVA. MATEMATICA Number 275, 2014, Pages 60 73 ISSN 1024 7696 Max-Erlang Min-Erlang power series distributions as two new families of lifetime distribution
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 informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More informationMISCELLANEOUS TOPICS RELATED TO LIKELIHOOD. Copyright c 2012 (Iowa State University) Statistics / 30
MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD Copyright c 2012 (Iowa State University) Statistics 511 1 / 30 INFORMATION CRITERIA Akaike s Information criterion is given by AIC = 2l(ˆθ) + 2k, where l(ˆθ)
More informationThe Log-Beta Generalized Half-Normal Regression Model
Marquette University e-publications@marquette Mathematics, Statistics and Computer Science Faculty Research and Publications Mathematics, Statistics and Computer Science, Department of 1-1-13 The Log-Beta
More informationUNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS
ISSN 0104-0499 UNIVERSIDADE FEDERAL DE SÃO CARLOS CENTRO DE CIÊNCIAS EXATAS E TECNOLÓGICAS THE FGM BIVARIATE LIFETIME COPULA MODEL: A BAYESIAN APPROACH Francisco Louzada-Neto Vicente G. Cancho Adriano
More informationChapter 2 Inference on Mean Residual Life-Overview
Chapter 2 Inference on Mean Residual Life-Overview Statistical inference based on the remaining lifetimes would be intuitively more appealing than the popular hazard function defined as the risk of immediate
More information11 Survival Analysis and Empirical Likelihood
11 Survival Analysis and Empirical Likelihood The first paper of empirical likelihood is actually about confidence intervals with the Kaplan-Meier estimator (Thomas and Grunkmeier 1979), i.e. deals with
More informationThe Lindley-Poisson distribution in lifetime analysis and its properties
The Lindley-Poisson distribution in lifetime analysis and its properties Wenhao Gui 1, Shangli Zhang 2 and Xinman Lu 3 1 Department of Mathematics and Statistics, University of Minnesota Duluth, Duluth,
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 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 informationKey Words: survival analysis; bathtub hazard; accelerated failure time (AFT) regression; power-law distribution.
POWER-LAW ADJUSTED SURVIVAL MODELS William J. Reed Department of Mathematics & Statistics University of Victoria PO Box 3060 STN CSC Victoria, B.C. Canada V8W 3R4 reed@math.uvic.ca Key Words: survival
More informationLikelihood Construction, Inference for Parametric Survival Distributions
Week 1 Likelihood Construction, Inference for Parametric Survival Distributions In this section we obtain the likelihood function for noninformatively rightcensored survival data and indicate how to make
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 informationPart III. Hypothesis Testing. III.1. Log-rank Test for Right-censored Failure Time Data
1 Part III. Hypothesis Testing III.1. Log-rank Test for Right-censored Failure Time Data Consider a survival study consisting of n independent subjects from p different populations with survival functions
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 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 informationAFT Models and Empirical Likelihood
AFT Models and Empirical Likelihood Mai Zhou Department of Statistics, University of Kentucky Collaborators: Gang Li (UCLA); A. Bathke; M. Kim (Kentucky) Accelerated Failure Time (AFT) models: Y = log(t
More informationEstimating the parameters of hidden binomial trials by the EM algorithm
Hacettepe Journal of Mathematics and Statistics Volume 43 (5) (2014), 885 890 Estimating the parameters of hidden binomial trials by the EM algorithm Degang Zhu Received 02 : 09 : 2013 : Accepted 02 :
More informationThe Compound Family of Generalized Inverse Weibull Power Series Distributions
British Journal of Applied Science & Technology 14(3): 1-18 2016 Article no.bjast.23215 ISSN: 2231-0843 NLM ID: 101664541 SCIENCEDOMAIN international www.sciencedomain.org The Compound Family of Generalized
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 informationADJUSTED PROFILE LIKELIHOODS FOR THE WEIBULL SHAPE PARAMETER
ADJUSTED PROFILE LIKELIHOODS FOR THE WEIBULL SHAPE PARAMETER SILVIA L.P. FERRARI Departamento de Estatística, IME, Universidade de São Paulo Caixa Postal 66281, São Paulo/SP, 05311 970, Brazil email: sferrari@ime.usp.br
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
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 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 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 informationROBUSTNESS OF TWO-PHASE REGRESSION TESTS
REVSTAT Statistical Journal Volume 3, Number 1, June 2005, 1 18 ROBUSTNESS OF TWO-PHASE REGRESSION TESTS Authors: Carlos A.R. Diniz Departamento de Estatística, Universidade Federal de São Carlos, São
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 informationSemi-parametric Inference for Cure Rate Models 1
Semi-parametric Inference for Cure Rate Models 1 Fotios S. Milienos jointly with N. Balakrishnan, M.V. Koutras and S. Pal University of Toronto, 2015 1 This research is supported by a Marie Curie International
More informationTypical Survival Data Arising From a Clinical Trial. Censoring. The Survivor Function. Mathematical Definitions Introduction
Outline CHL 5225H Advanced Statistical Methods for Clinical Trials: Survival Analysis Prof. Kevin E. Thorpe Defining Survival Data Mathematical Definitions Non-parametric Estimates of Survival Comparing
More informationInference based on the em algorithm for the competing risk model with masked causes of failure
Inference based on the em algorithm for the competing risk model with masked causes of failure By RADU V. CRAIU Department of Statistics, University of Toronto, 100 St. George Street, Toronto Ontario,
More informationOptimization Methods II. EM algorithms.
Aula 7. Optimization Methods II. 0 Optimization Methods II. EM algorithms. Anatoli Iambartsev IME-USP Aula 7. Optimization Methods II. 1 [RC] Missing-data models. Demarginalization. The term EM algorithms
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 informationA unified view on lifetime distributions arising from selection mechanisms
A unified view on lifetime distributions arising from selection mechanisms Josemar Rodrigues a, N. Balakrishnan b, Gauss M. Cordeiro c and Mário de Castro d a Universidade Federal de São Carlos, Departamento
More informationLecture 3. Truncation, length-bias and prevalence sampling
Lecture 3. Truncation, length-bias and prevalence sampling 3.1 Prevalent sampling Statistical techniques for truncated data have been integrated into survival analysis in last two decades. Truncation in
More informationComputational Statistics and Data Analysis. Estimation for the three-parameter lognormal distribution based on progressively censored data
Computational Statistics and Data Analysis 53 (9) 358 359 Contents lists available at ScienceDirect Computational Statistics and Data Analysis journal homepage: www.elsevier.com/locate/csda stimation for
More informationInterval Estimation for Parameters of a Bivariate Time Varying Covariate Model
Pertanika J. Sci. & Technol. 17 (2): 313 323 (2009) ISSN: 0128-7680 Universiti Putra Malaysia Press Interval Estimation for Parameters of a Bivariate Time Varying Covariate Model Jayanthi Arasan Department
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 informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
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 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 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 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 informationLecture 22 Survival Analysis: An Introduction
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 22 Survival Analysis: An Introduction There is considerable interest among economists in models of durations, which
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 informationEmpirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm
Empirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm Mai Zhou 1 University of Kentucky, Lexington, KY 40506 USA Summary. Empirical likelihood ratio method (Thomas and Grunkmier
More informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume VI: Statistics PART A Edited by Emlyn Lloyd University of Lancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester
More informationMultilevel Statistical Models: 3 rd edition, 2003 Contents
Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction
More informationTwo-stage Adaptive Randomization for Delayed Response in Clinical Trials
Two-stage Adaptive Randomization for Delayed Response in Clinical Trials Guosheng Yin Department of Statistics and Actuarial Science The University of Hong Kong Joint work with J. Xu PSI and RSS Journal
More informationStatistical Estimation
Statistical Estimation Use data and a model. The plug-in estimators are based on the simple principle of applying the defining functional to the ECDF. Other methods of estimation: minimize residuals from
More informationSurvival Analysis I (CHL5209H)
Survival Analysis Dalla Lana School of Public Health University of Toronto olli.saarela@utoronto.ca January 7, 2015 31-1 Literature Clayton D & Hills M (1993): Statistical Models in Epidemiology. Not really
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 24 Paper 153 A Note on Empirical Likelihood Inference of Residual Life Regression Ying Qing Chen Yichuan
More informationFrailty Models and Copulas: Similarities and Differences
Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt
More informationSome classes of statistical distributions. Properties and Applications
DOI: 10.2478/auom-2018-0002 An. Şt. Univ. Ovidius Constanţa Vol. 26(1),2018, 43 68 Some classes of statistical distributions. Properties and Applications Irina Băncescu Abstract We propose a new method
More informationMeei Pyng Ng 1 and Ray Watson 1
Aust N Z J Stat 444), 2002, 467 478 DEALING WITH TIES IN FAILURE TIME DATA Meei Pyng Ng 1 and Ray Watson 1 University of Melbourne Summary In dealing with ties in failure time data the mechanism by which
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 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 informationContents. Preface to Second Edition Preface to First Edition Abbreviations PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1
Contents Preface to Second Edition Preface to First Edition Abbreviations xv xvii xix PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1 1 The Role of Statistical Methods in Modern Industry and Services
More informationISI Web of Knowledge (Articles )
ISI Web of Knowledge (Articles 1 -- 18) Record 1 of 18 Title: Estimation and prediction from gamma distribution based on record values Author(s): Sultan, KS; Al-Dayian, GR; Mohammad, HH Source: COMPUTATIONAL
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 informationA Quantile Implementation of the EM Algorithm and Applications to Parameter Estimation with Interval Data
A Quantile Implementation of the EM Algorithm and Applications to Parameter Estimation with Interval Data Chanseok Park The author is with the Department of Mathematical Sciences, Clemson University, Clemson,
More informationSemi-Competing Risks on A Trivariate Weibull Survival Model
Semi-Competing Risks on A Trivariate Weibull Survival Model Cheng K. Lee Department of Targeting Modeling Insight & Innovation Marketing Division Wachovia Corporation Charlotte NC 28244 Jenq-Daw Lee Graduate
More informationSTAT331. Cox s Proportional Hazards Model
STAT331 Cox s Proportional Hazards Model In this unit we introduce Cox s proportional hazards (Cox s PH) model, give a heuristic development of the partial likelihood function, and discuss adaptations
More informationST495: Survival Analysis: Maximum likelihood
ST495: Survival Analysis: Maximum likelihood Eric B. Laber Department of Statistics, North Carolina State University February 11, 2014 Everything is deception: seeking the minimum of illusion, keeping
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationThe Log-generalized inverse Weibull Regression Model
The Log-generalized inverse Weibull Regression Model Felipe R. S. de Gusmão Universidade Federal Rural de Pernambuco Cintia M. L. Ferreira Universidade Federal Rural de Pernambuco Sílvio F. A. X. Júnior
More informationQuantile Regression for Residual Life and Empirical Likelihood
Quantile Regression for Residual Life and Empirical Likelihood Mai Zhou email: mai@ms.uky.edu Department of Statistics, University of Kentucky, Lexington, KY 40506-0027, USA Jong-Hyeon Jeong email: jeong@nsabp.pitt.edu
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 compound class of Poisson and lifetime distributions
J. Stat. Appl. Pro. 1, No. 1, 45-51 (2012) 2012 NSP Journal of Statistics Applications & Probability --- An International Journal @ 2012 NSP Natural Sciences Publishing Cor. A compound class of Poisson
More informationLocal Influence and Residual Analysis in Heteroscedastic Symmetrical Linear Models
Local Influence and Residual Analysis in Heteroscedastic Symmetrical Linear Models Francisco José A. Cysneiros 1 1 Departamento de Estatística - CCEN, Universidade Federal de Pernambuco, Recife - PE 5079-50
More informationIntroduction to Statistical Analysis
Introduction to Statistical Analysis Changyu Shen Richard A. and Susan F. Smith Center for Outcomes Research in Cardiology Beth Israel Deaconess Medical Center Harvard Medical School Objectives Descriptive
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 informationExercises. (a) Prove that m(t) =
Exercises 1. Lack of memory. Verify that the exponential distribution has the lack of memory property, that is, if T is exponentially distributed with parameter λ > then so is T t given that T > t for
More informationHypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations
Hypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations Takeshi Emura and Hisayuki Tsukuma Abstract For testing the regression parameter in multivariate
More informationNew Bayesian methods for model comparison
Back to the future New Bayesian methods for model comparison Murray Aitkin murray.aitkin@unimelb.edu.au Department of Mathematics and Statistics The University of Melbourne Australia Bayesian Model Comparison
More information