Journal of Statistical Computation and Simulation. The Log-Exponentiated Generalized Gamma Regression Model for Censored Data
|
|
- Ashlee Small
- 6 years ago
- Views:
Transcription
1 The Log-Exponentiated Generalized Gamma Regression Model for Censored Data Journal: Journal of Statistical Computation and Simulation Manuscript ID: GSCS-00-0.R Manuscript Type: Original Paper Date Submitted by the Author: 0-Feb-0 Complete List of Authors: Ortega, Edwin; Universidade de São Paulo, Ciências Exatas Cordeiro, Gauss; UFRPE de Pascoa, Marcelino; Universidade de São Paulo, Ciências Exatas Couto, Epamonondas; Universidade de São Paulo, Ciências Extas Areas of Interest: CENSORED DATA, REGRESSION ANALYSIS, RELIABILITY, SURVIVAL ANALYSIS Note: The following files were submitted by the author for peer review, but cannot be converted to PDF. You must view these files (e.g. movies) online. Epa_Edwin_Gauss_final_0_0_0.tex
2 Page of Journal of Statistical Computation and Simulation The Log-Exponentiated Generalized Gamma Regression Model for Censored Data Edwin M. M. Ortega Departamento de Ciências Exatas, ESALQ - USP Gauss M. Cordeiro Departamento de Estatística, UFPE Marcelino A. R. Pascoa Departamento de Ciências Exatas, ESALQ - USP Epaminondas V. Couto Departamento de Ciências Exatas, ESALQ - USP Abstract For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution (Cordeiro et al., 00) that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub and unimodal shaped hazard functions. A second advantage is that it includes classical distributions reported in lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely-known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data. Keywords: Censored data; Exponentiated generalized gamma distribution; Log-gamma generalized regression; Survival function. Introduction Standard lifetime distributions usually present very strong restrictions to produce bathtub curves, and thus appear to be unappropriate for interpreting data with this characteristic. The gamma distribution is the most popular model for analyzing skewed data. The generalized gamma Address for correspondence: Departamento de Ciências Exatas, ESALQ/USP, Av. Pádua Dias - Caixa Postal, -00, Piracicaba - São Paulo - Brazil. edwin@esalq.usp.br
3 Page of distribution (GG) was introduced by Stacy () and includes as special sub-models: the exponential, Weibull, gamma and Rayleigh distributions, among others. The GG distribution is suitable for modeling data with hazard rate function of different forms: increasing, decreasing, bathtub and unimodal, which makes it useful for estimating individual hazard functions as well as both relative hazards and relative times (Cox, 00). Recently, the GG distribution has been used in several research areas such as engineering, hydrology and survival anlaysis. Ortega et al. (00) discussed influence diagnostics in GG regression models, Nadarajah and Gupta (00) used this distribution with application to drought data, Cox et al. (00) presented a parametric survival analysis and taxonomy of GG hazard functions and Ali et al. (00) derived the exact distributions of the product X X and the quotient X /X, when X and X are independent GG random variables providing applications of their results to drought data from Nebraska. Further, Gomes et al. (00) focused on the parameter estimation, Ortega et al. (00) compared three types of residuals based on the deviance component in GG regression models under censored observations, Cox (00) discussed and compared the F -generalized family with the GG model, Almpanidis and Kotropoulos (00) presented a text-independent automatic phone segmentation algorithm based on the GG distribution and Nadarajah (00a) analyzed some incorrect references with respect to the use of this distribution in electrical and electronic engineering. More recently, Barkauskas et al. (00) modeled the noise part of a spectrum as an autoregressive moving average (ARMA) model with innovations following the GG distribution, Malhotra et al. (00) provided a unified analysis for wireless system over generalized fading channels that is modeled by a two parameter GG model and Xie and Liu (00) analyzed three-moment auto conversion parametrization based on the GG distribution. Further, Ortega et al. (00) proposed a modified GG regression model to allow the possibility that long-term survivors may be presented in the data and Cordeiro et al. (00) proposed the exponentiated generalized gamma (EGG) distribution. This distribution due to its flexibility in accommodating many forms of the risk function seems to be an important model that can be used in a variety of problems in survival analysis. In the last decade, new classes of distributions for modeling survival data based on extensions of the Weibull distribution were developed. Mudholkar et al. () introduced the exponentiated Weibull (EW) distribution, Xie and Lai () presented the additive Weibull distribution, Lai et al. (00) proposed the modified Weibull (MW) distribution and Carrasco et al. (00) defined the generalized modified Weibull (GMW) distribution. Furthermore, the main motivation for the use of the EGG distribution is that it contains as special sub-models several distributions such as the generalized gamma (GG), EW, exponentiated exponential (EE) (Gupta and Kundu, ) and generalized Rayleigh (GR) (Kundu and Raqab, 00) distributions. The EGG distribution can model four types of the failure rate function (i.e. increasing, decreasing, unimodal and bathtub) depending on the values of its parameters. It is also suitable for testing goodness-of-fit of some special sub-models, such as the GG, EW, Weibull and GR distributions. Different forms of regression models have been proposed in survival analysis. Among them, the location-scale regression model (Lawless, 00) is distinguished since it is frequently used in clinical trials. In this article, we propose an extension of the EGG distribution and obtain some
4 Page of Journal of Statistical Computation and Simulation of its structural properties. Further, we propose a location-scale regression model based on this distribution, referred to as the log-exponentiated generalized gamma (LEGG) regression model, which is a feasible alternative for modeling the four existing types of failure rate functions. The article is organized as follows. In Section, we define an extended version of the EGG distribution. In Section, we consider its moments, generating function, mean deviations, reliability, order statistics and their moments. In Section, we provide a simulation study. In Section, we define the LEGG distribution and derive an expansion for its moments. In Section, we propose a LEGG regression model for censored data. We consider the method of maximum likelihood to estimate the model parameters and derive the observed information matrix. In Section, we give two applications using well-known data sets to demonstrate the applicability of the proposed regression model. Section ends with some concluding remarks. The Exponentiated Generalized Gamma Distribution Cordeiro et al. (00) proposed the EGG distribution with four parameters α > 0, τ > 0, k > 0 and λ > 0 to extend the GG distribution (Stacy, ) that should be able to fit various types of data. The probability density function (pdf) of the EGG distribution has the form f(t) = λτ ( t ) τk [ ( t τ ] exp αγ(k) α α) { [ ( t τ ] γ k, α) } λ, t > 0, () where Γ( ) is the gamma function, γ(k, x) = x 0 wk e w dw is the incomplete gamma function and γ (k, x) = γ(k, x)/γ(k) is the incomplete gamma function ratio. The function γ (k, x) is simply the cumulative distribution function (cdf) of a standard gamma distribution with shape parameter k. In (), α is a scale parameter and τ, k and λ are shape parameters. The Weibull and GG distributions arise as special sub-models of () when λ = k = and λ =, respectively. The EGG distribution approaches the log-normal distribution when λ = and k. If T is a random variable with density (), we write T EGG(α, τ, k, λ). The survival and hazard rate functions corresponding to () are and [ ( t τ ] S(t) = F (t) = {γ k, α) } λ h(t) = αγ(k)( λτ t ) τk [ ( t τ ] exp α α) { [ ( t τ ] γ k, α) } λ { ( ( t ) τ )] } λ [γ k,, () α respectively. The EGG distribution has a hazard rate function that involves a complicated function, although it can be easily computed numerically. Moreover, it is quite flexible for modeling survival data. ()
5 Page of We consider a generalized form of the EGG density function () given by f(t) = λ τ ( t ) τk [ ( t τ ] exp αγ(k) α α) { [ ( t τ ] γ k, α) } λ, t > 0, () where τ is not zero and the other parameters are positive. If τ > 0, () is the EGG distribution and if λ = and τ > 0 it becomes the GG distribution. The special case λ = k = and τ > 0 gives the Weibull distribution. For λ = k = and τ < 0, it yields to the reciprocal Weibull (or inverse Weibull) distribution. Other special sub-models of the EGG distribution are discussed by Cordeiro et al. (00). Plots of the EGG density function for selected values of τ > 0 and τ < 0 are given in Figure. f(t) (a) (b) (c) t α=; τ=0.; k=; λ= α=; τ=; k=; λ= α=; τ=; k=; λ= α=; τ=; k=; λ= α=; τ=; k=; λ= f(t) t α=; τ= 0.; k=; λ= α=; τ= ; k=; λ= α=; τ= ; k=; λ= α=; τ= ; k=; λ= α=; τ= ; k=; λ= f(t) t α=.; τ=.; k=.; λ=0. α=; τ=.; k=; λ=0. α=.; τ= ; k=0.; λ=. α=.; τ= 0.; k=0.; λ=0 α=; τ=; k=; λ=0 α=.; τ= 0.; k=; λ=0. Figure : The EGG density function: (a) For some values of τ > 0. (b) For some values of τ < 0. (c) For some values of τ > 0 and τ < 0. The cdf of the EGG distribution can be defined by { [ γ k, ( t ]} α )τ λ F (t) = { [ γ k, ( t ]} α )τ λ if τ > 0, if τ < 0. ()
6 Page of Journal of Statistical Computation and Simulation The survival function corresponding to () is S(t) = { γ [ k, ( t α )τ ]} λ { γ [ k, ( t α )τ ]} λ if τ > 0, if τ < 0. Let g α,τ,k (t) be the density function of the GG(α, τ, k) distribution (Stacy and Mihram, ) given by g α,τ,k (t) = τ ( ) t τk ( ) ] t τ exp[, t > 0. () αγ(k) α α Theorem. If T EGG(α, τ, k, λ), we have the representation f(t) = m,i=0 w(m, i, k, λ) g α,τ,k(m+)+i (t), t > 0, () where g α,τ,k(m+)+i (t) is the GG(α, τ, k(m + )i) density function defined by () and the weighted coefficients w(m, i, k, λ) are given by w(m, i, k, λ) = m,i=0 () λ s m (λ) c m,i Γ[k(m + ) + i] Γ(k) m+, () where the quantity s m (λ) is given by () and the constants c m,i can be determined from the recurrence equation (0). Clearly, the coefficients satisfy m,i=0 w(m, i, k, λ) =, Some mathematical properties of the EGG distribution can follow directly from those properties of the GG distribution, since equation () is expressed in terms of a linear combination of GG densities. For example, the ordinary, inverse and factorial moments and moment generating function (mgf) of the EGG distribution can be derived directly from those quantities of the GG distribution. The proof of the theorem is in Appendix B. General Properties. Moments Here, we give two different expansions for the moments of the EGG distribution. Let µ r = E(T r ) be the rth ordinary moment of the EGG distribution. First, we obtain an infinite sum representation for µ r from equation (). The rth moment of the GG(α, τ, k) distribution is µ r,gg = αr Γ(k + r/τ). Γ(k)
7 Page of Equation () then immediately gives µ r = λ α r m,i=0 s m (λ) c m,i Γ([k(m + ) + i] + r/τ) Γ(k) m+. (0) Equation (0) has the inconvenient that depends on the constants c m,i that can be calculated recursively from (0). Theorem. If T EGG(α, τ, k, λ), the rth moment of T reduces to µ r = λ αr sgn(τ) Γ(k) j=0 m=0 where A(k, r/τ, m) = k m Γ ( r/τ + k(m + ) ) and j ( ) ( ) λ j ( ) j+m I j m (k, r τ, m ), () I (k, r ) τ, m = A(k, r/τ, m) F (m) ( ) r/τ + k(m + ); k,, k; k +,, k + ;,,. () A Hence, as an alternative representation to (0), the rth moment of the EGG distribution follows from both equations () and () as an infinite weighted sum of the Lauricella functions of type A. The proof of the theorem is in Appendix B.. Moment Generating Function Theorem. If T EGG(α, τ, k, λ), the mgf of T reduces to M(t) = m,i=0 where w(m, i, k, λ) is defined by () and sgn(τ) M α,τ,k(m+)+i (s) = Γ(k(m + ) + i) The proof of the theorem is in Appendix B.. Deviations w(m, i, k, λ) M α,τ,k(m+)+i (s), () ( m ) (sα) m Γ τ + k(m + ) + i. m! m=0 If T has the EGG density function f(t), the deviations about the mean δ = E( T µ ) and about the median δ = E( T m ) can be calculated from the relations δ = µ F (µ ) I(µ ) and δ = µ I(m), ()
8 Page of Journal of Statistical Computation and Simulation where F (a) is easily obtained from (), m is the solution of F (m) = / and I(s) = s 0 t f(t)dt. This integral can be determined from () as where I(s) = m,i=0 w(m, i, k, λ) J(α, τ, [k(m + ) + i], s), () J(α, τ, k, s) = s 0 t g α,τ,k (t)dt. We can obtain from the density of the GG(α, τ, k) distribution by setting x = t/α J(α, τ, k, s) = α τ Γ(k) s/α 0 x τk exp( x τ )dx. The substitution w = x τ yields J(α, τ, k, s) in terms of the incomplete gamma function J(α, τ, k, s) = α sgn(τ) Γ(k) (s/α) τ Hence, inserting the last result into () gives I(s) = m,i=0 0 w k+τ exp( w)dw = γ(k + τ, (s/α) τ ). w(m, i, k, λ)γ(k(m + ) + i + τ, (s/α) τ ). () The mean deviations of the EGG distribution can be obtained from equations () and ().. Reliability In the context of reliability, the stress-strength model describes the life of a component which has a random strength X that is subjected to a random stress X. The component fails at the instant that the stress applied to it exceeds the strength, and the component will function satisfactorily whenever X > X. For X and X independent random variables having a common EGG distribution, the reliability R = R = P r(x < X ) = 0 f(t) F (t)dt, where f(t) and F (t) are calculated from () and (), respectively, can be written explicitly as follows:
9 Page of For τ > 0, R = [ τ w(m, i, k, λ) ( ) ] t τ ( ) t τ[k(m+)+i] exp αγ[k(m + ) + i] m,i=0 0 α α ( [ ( ) ]} t τ λ ) {γ k, dt. α Setting x = ( t α) τ in the last equation gives R = m,i=0 sgn(τ) w(m, i, k, λ) Γ[k(m + ) + i] By equation () in Appendix A, we obtain where and R = j m,i,j=0 m =0 0 x k(m+)+i exp( x) γ (k, x) λ dx. v (m, m, i, j, k, λ) I[k(m + ), i, m ], () v (m, m, i, j, k, λ) = ( )j+m ( ) ( ) sgn(τ) λ s m (λ) c m,i λ j Γ(k) m+, j I[k(m + ), i, m ] = 0 m x k(m+)+i exp( x) γ (k, x) m dx. Defining C = [k(m + )] m Γ ( i + k(m + m + ) ) and using the Lauricella function of type A (see Appendix B), this integral can be expressed as (see Nadarajah, 00b, equation ()) I [k(m + ), i, m ] = C F (m ( ) i + k(m + m + ); k(m + ),..., k(m + ); For τ < 0, R = A k(m + ) +,..., k(m + ) + ;,..., ). () [ τ w(m, i, k, λ) ( ) ] t τ ( ) t τ[k(m+)+i] exp α Γ[k(m + ) + i] m,i=0 0 α α [ ( ) ]} t τ λ {γ k, dt. α
10 Page 0 of Journal of Statistical Computation and Simulation Setting x = ( t τ α) in the last equation yields sgn(τ) w(m, i, k, λ) R = Γ[k(m + )] + i m,i=0 Again, by equation (), we obtain where R = j m,i,j=0 m =0 0 x k(m+)+i exp( x) γ (k, x) λ dx. v (m, m, i, j, k, τ, λ) I [k(m + ), i, m ], () v (m, m, i, j, k, τ, λ) = ( )j+m ( )( ) sgn(τ) λ j, Γ[k(m + ) + i] j w(m, i, k, λ) is just defined after () and I [k(m + ), i, m ] is given by ().. Order Statistics The density function f i:n (t) of the ith order statistic, for i =,..., n, from random variables T,..., T n having density (), can be written as n i f i:n (t) = B(i, n i + ) f(t) ( ) n i ( ) l F (t) i+l, (0) l l=0 where f(t) and F (t) are the pdf and cdf of the EGG distribution, respectively, and B(, ) denotes the beta function. Let f α,τ,k,λ (t) be the density function of the EGG(α, τ, k, λ) distribution. Plugging () and () in the last equation and after some algebra, we can write: For τ > 0, n i ( ) l ( ) n i f i:n (t) = f B(i, n i + ) (i + l) l α,τ,k,λ(i+l) (t). () l=0 Equation () gives the density function of the ith order statistic as a finite linear combination of EGG densities. Hence, the moments of the order statistics can be calculated directly from () using equations () and (). The rth moment of the ith order statistic, say µ r(i:n), reduces to µ r(i:n) = λ αr sgn(τ) n i B(i, n i + ) j l=0 j=0 m=0 where I(k, r/τ, m) comes from (). ( ) l+j+m Γ(k) m ( )( )( ) n i λ(i + l) j I l j m (k, r τ, m ), ()
11 Page of For τ < 0, f i:n (t) = ( ) λ τ t τk ( ) [ t τ ( ) ]} t τ λ exp[ ]{γ k, αγ(k)b(i, n i + ) α α α n i ( ( [ ( ) n i t τ ]} ) λ i+l )( ) l {γ k,. l α l=0 Using the binomial expansion in the last expression, we have n i f i:n (t) = B(i, n i + ) l=0 i+l j =0 ( ) l+j ( )( ) n i i + l f ( + j ) l α,τ,k,λ(+j )(t). () Equation () shows that the density function of the ith order statistic is a finite linear combination of EGG densities and then its moments can be determined directly from () using equations () and (). We obtain µ r(i:n) = Simulation study λ α r sgn(τ) n i B(i, n i + ) l=0 j =0 ( )( ) λ( + j ) j I j m i+l j j=0 m=0 j ( ) l+j+j +m Γ(k) ( )( ) n i i + l (k, r τ, m ). () We perform some Monte Carlo simulation studies to assess on the finite sample behavior of the maximum likelihood estimators (MLEs) of α, τ, k and λ. The results were obtained from 000 Monte Carlo replications and the simulations were carried out using the software R. In each replication, a random sample of size n is drawn from the EGG(α, τ, k, λ) distribution and the parameters were estimated by maximum likelihood. The EGG random variable was generated using the inversion method. The true parameter values used in the data generating processes are α = 0., τ = 0., k =. and λ =.0. The mean estimates of the four model parameters and the corresponding root mean squared errors (RMSEs) for sample sizes n = 0, n = 00, n = 00 and n = 00 are listed in Table. The figures in Table show that the biases and RMSEs of MLEs of α, τ, k and λ decay toward zero when the sample size increases, as expected. There is a small sample bias in the estimation of the model parameters. Future research should be conducted to obtain bias corrections for these estimators. 0 l j
12 Page of Journal of Statistical Computation and Simulation Table : Mean estimates and RMSEs of α, τ, k and λ n Parameter Mean RMSE 0 α τ k.0. λ α τ k.. λ α τ k.0.0 λ.. 00 α τ k..0 λ.. The Log-Exponentiated Generalized Gamma Distribution Henceforth, let T be a random variable having the EGG density function () and Y = log(t ). It is easy to verify that the density function of Y obtained by replacing k = q, τ = (σ k) and α = exp[µ τ log(k)] reduces to f(y) = λ q (q ) q { σγ(q exp q ( y µ ) [ ( y µ )]} q exp q ) σ σ [ { ( y µ )}]} λ, {γ q, q exp q () σ where < y <, < µ <, σ > 0, λ > 0 and q is different from zero. We consider an extended form including the case q = 0 (Lawless, 00). Thus, the density of Y can be written as [ f(y) = {γ λ q (q ) q σγ(q ) { exp q ( y µ σ q, q exp { q ( y µ σ [ λ πσ exp ( y µ σ ) q exp [ q ( y µ )] } )} ]} λ ) ]Φ (λ )( y µ ) σ σ if q 0, if q = 0, ()
13 Page of where Φ( ) is the standard normal cumulative distribution. The sub-model defined by q = 0 is exactly the skew normal distribution, whereas if λ = and q =, we obtain the extreme value distribution. Further, the sub-model λ = corresponds to the log-gamma generalized distribution and, if in addition q =, it reduces to the log-inverse Weibull distribution. We refer to equation () as the LEGG distribution, say Y LEGG(µ, σ, q, λ), where µ R is the location parameter, σ > 0 is the scale parameter and q and λ are shape parameters. So, if T EGG(α, τ, k, λ) then Y = log(t ) LEGG(µ, σ, q, λ). The plots of the density function () for selected values of λ, µ = 0 and σ = for q < 0, q > 0 and q = 0 are given in Figure. These plots show that the LEGG distribution could be very flexible for modeling its kurtosis. The corresponding survival function is f(y) (a) (b) (c) 0 y q= 0.; λ=0 q= 0.; λ= q= ; λ= q= ; λ=. q=.; λ=. f(y) q=0.; λ=0. q=0.; λ= q=.; λ= q=.; λ=. q=; λ= 0 y f(y) λ=0. λ= λ= λ=. λ= 0 Figure : The LEGG density curves: (a) For some values of q < 0. (b) For some values of q > 0. (c) For some values of q = 0. q > 0 ( S(y) = F (y) = P (Y > y) = P (µ + σz > y) = P Z > y µ ) = P (Z > z) σ λq(q ) q { [ λ du; = Γ(q exp q u q exp(qu) }{γ q, q exp(qu)]} ) z y
14 Page of Journal of Statistical Computation and Simulation q < 0 ( S(y) = F (y) = P (Y > y) = P (µ + σz > y) = P Z > y µ ) = P (Z > z) σ λq(q ) q { [ λ du. = Γ(q exp q u q exp(qu) }{γ q, q exp(qu)]} ) z These integrals have explicit forms { ( γ q, q exp [ q( y µ σ )])} λ { S(y) = ( γ q, q exp [ q( y µ )])} λ Φ λ( y µ ) σ σ if q > 0, if q < 0, if q = 0. Theorem. If Y LEGG(µ, σ, q, λ), then the rth moment is given by i. For q 0, µ r = m,i=0 l=0 where Γ(p) = Γ(p) p. ii. For q = 0, r ( ) [ λ sgn(q)s m (λ) c m,i r σ Γ (q ) m+ l q log( q ) + µ µ r = λ p=0 j=0 r ( r s p (λ) j ] r l [ Γ(q (m + ) + i)] σl q (), () ) σ j µ r j τ j,p, () where s p (λ) is given by (), and when j + p l is even, we obtain p ( ) ( p j + p l + τ j,p = j/ π (p+/) l π l Γ l F (p l) A l=0 (j+p l) even ( j + p l + ) ;,, ;,, ;,, ). (0) Inserting (0) in () yields the rth moment of the LEGG distribution when q = 0. The proof of Theorem is given in Appendix B. The skewness and kurtosis measures can be calculated from the ordinary moments using wellknown relationships. Plots of the skewness and kurtosis for selected values of µ and σ versus λ are shown in Figures and, respectively.
15 Page of Skewness 0 0 µ=0.0 µ=0. µ=0. µ=0. (a) 0 0 λ Kurtosis (b) µ=0.0 µ=0. µ=0. µ= Figure : Skewness and kurtosis of the LEGG distribution for some values of µ with σ = 0., λ = and q =.. Skewness 0 σ= σ= σ=0 σ= (a) 0 0 λ Kurtosis σ= σ= σ=0 σ= (b) λ 0 0 Figure : Skewness and kurtosis of the LEGG distribution for some values of σ with µ = 0.00, λ = and q = 0.. The Log-Exponentiated Generalized Gamma Regression Model In many practical applications, the lifetimes are affected by explanatory variables such as the cholesterol level, blood pressure, weight and many others. Parametric models to estimate univariate survival functions for censored data regression problems are widely used. A parametric model that provides a good fit to lifetime data tends to yield more precise estimates of the quantities of interest. If Y LEGG(µ, σ, q, λ), we define the standardized random variable Z = (Y µ)/σ λ
16 Page of Journal of Statistical Computation and Simulation having density function f(z) = λ q Γ(q ) (q ) q exp { q z q exp(qz) }{ [ γ q, q exp(qz) ]} λ if q 0, λ π exp ( ) z Φ (λ ) (z) if q = 0. We write Z LEGG(0,, q, λ). Now, we propose a linear location-scale regression model linking the response variable y i and the explanatory variable vector x T i = (x i,..., x ip ) by () y i = x T i β + σz i, i =,..., n, () where the random error z i has density function (), β = (β,..., β p ) T, σ > 0, λ > 0 and < q < are unknown parameters. The parameter µ i = x T i β is the location of y i. The location parameter vector µ = (µ,..., µ n ) T is is given by a linear model µ = Xβ, where X = (x,..., x n ) T is a known model matrix. The LEGG model () opens new possibilities for fitted many different types of data. It contains as special sub-models the following well-known regression models. For λ = q =, we obtain the classical Weibull regression model (see, Lawless, 00). If σ = and σ = 0., in addition to λ = q =, the LEGG regression model reduces to the exponential and Rayleigh regression models, respectively. For q =, we obtain the log-exponentiated Weibull regression model introduced (Mudholkar et al., ). See, also, Cancho et al. (), Ortega et al. (00), Cancho et al. (00) and Hashimoto et al. (00). If σ =, in addition to q =, the LEGG regression model reduces to the log-exponentiated exponential regression model. If σ = 0., in addition to q =, the LEGG model reduces to the log-generalized Rayleigh regression model. For λ =, we obtain the log-gamma generalized (LGG) regression model (Lawless, 00; Ortega et al., 00). More recently, the LGG distribution has been used in several research areas. See, for example, Ortega et al. (00, 00). For q =, it gives the log-generalized inverse Weibull regression model. If λ =, in addition to q =, we obtain the log-inverse Weibull regression model (Gusmão et al., 00). Finally, for q = 0, it yields the log-exponentiated normal regression model.. Maximum Likelihood Estimation Consider a sample (y, x ),..., (y n, x n ) of n independent observations, where each random response is defined by y i = min{log(t i ), log(c i )}. We assume non-informative censoring such that the observed lifetimes and censoring times are independent. Let F and C be the sets of individuals for which y i is the log-lifetime or log-censoring, respectively. Conventional likelihood estimation techniques can be applied here. The log-likelihood function for the vector of parameters θ = (λ, σ, q, β T ) T from model () has the form l(θ) = l i (θ) + l (c) i (θ), where l i (θ) = log[f(y i )], l (c) i (θ) = log[s(y i )], f(y i ) is the density () and S(y i ) is the survival function () of Y i. The total i F i C
17 Page of log-likelihood function for θ reduces to { log λq(q ) q exp [ q z σγ(q ) i q exp(qz i ) ]{ [ γ q, q exp(qz i ) ]} } λ i F + { log { γ [q, q exp(qz i )] } } λ, if q > 0, i C { l(θ) = log λq(q ) q exp [ q z σγ(q ) i q exp(qz i ) ]{ [ γ q, q exp(qz i ) ]} } λ i F +λ { } log γ [q, q exp(qz i )], if q < 0, i C } log{ λ πσ exp( z i )Φ(λ ) (z i ) + log [ Φ λ (z i ) ], if q = 0, i F () where z i = (y i x T i β)/σ. The MLE θ of the vector of unknown parameters can be calculated by maximizing the log-likelihood (). The maximization of () follows the same two steps for obtaining the MLE of θ for the uncensored case. In general, it is reasonable to expect that the shape parameter q belongs to the interval [, ]. We fixed, in the first step of the iterative process, different values of q in this interval. Then, we obtain the MLEs λ(q), σ(q) and β(q) and the maximized log-likelihood function L max (q) is determined. We use, in this step, a subroutine NLMixed of SAS. In the second step, the log-likelihood L max (q) is maximized, and then q is obtained. The MLEs of λ, σ and β are, respectively, given by λ = λ( q), σ = σ( q) and β = β( q). Initial values for β and σ are obtained by fitting the Weibull regression model with λ = and q =. The fit of the LEGG model yields the estimated survival function for y i (ẑ i = (y i x T i ˆβ)/ˆσ) given by { γ (ˆq, ˆq exp[ˆq(ẑ i )] )}ˆλ, if q > 0, Ŝ(y i ; ˆλ, ˆσ, ˆq, β T { ) = γ (ˆq, ˆq exp[ˆq(ẑ i )] )}ˆλ, if q < 0, i C Φˆλ ( ẑ i ), if q = 0. Under conditions that are fulfilled for the parameter vector θ in the interior of the parameter space but not on the boundary, the asymptotic distribution of n( θ θ) is multivariate normal N p+ (0, K(θ) ), where K(θ) is the information matrix. The asymptotic covariance matrix K(θ) of θ can be approximated by the inverse of the (p+) (p+) observed information matrix L(θ). The elements of the observed information matrix L(θ), namely L λλ, L λσ, L λβj, L σσ, L σβj and L βj β s for j, s =,..., p, are given in Appendix C. The approximate multivariate normal distribution N p+ (0, L(θ) ) for θ can be used in the classical way to construct approximate confidence regions for some parameters in θ. ()
18 Page of Journal of Statistical Computation and Simulation We can construct likelihood ratio (LR) statistics for comparing some special sub-models with the LEGG model. We consider the partition θ = (θ T, θ T ) T, where θ is a subset of parameters of interest and θ is a subset of remaining parameters. The LR statistic for testing the null hypothesis H 0 : θ = θ (0) versus the alternative hypothesis H : θ θ (0) is w = {l( θ) l( θ)}, where θ and θ are the estimates under the null and alternative hypotheses, respectively. The statistic w is asymptotically (as n ) distributed as χ k, where k is the dimension of the subset of parameters θ of interest. Applications In this section, we provide two applications to show the usefulness of the proposed regression model. Example : The diabetic retinopathy study We consider a data set analyzed by Huster et al. (), Liang et al. () and Wada and Hotta (000). Patients with diabetic retinopathy in both eyes and 0/00 or better visual acuity for both eyes were eligible for the study. One eye was randomly selected for the treatment and the other was observed without treatment. The patients were followed for two consecutively completed month follow-ups and the endpoint was the occurrence of visual acuity less than /00. We choose only the treatment time. A 0% sample of the high-risk patients defined by diabetic retinopathy criteria was taken for the data set (n = ) and the percentage of censored observations was.%. The variables involved in the study are: t i - failure time for the treatment (in min); cens i - censoring indicator (0=censoring, =lifetime observed); x i - age (0 = patient is an adult diabetic, = patient is a juvenile diabetic). We adopt the model y i = β 0 + β x i + σz i, where the random variable y i has the LEGG distribution () for i =,,...,. The MLEs of the model parameters are calculated using the NLMixed procedure in SAS. In order to estimate λ, σ and β of the LEGG regression model, we take some values for q. We choose the value ˆq =. that maximizes the likelihood function over several values of q (, ). Hence, this value is assumed for the MLE of the parameter q. Iterative maximization of the logarithm of the likelihood function () starts with initial values for β and σ taken from the fitted log-weibull regression model (with λ = ). The MLEs (approximate standard errors and p-values in parentheses) are: ˆλ =.(0.), ˆσ =.(0.), ˆβ 0 =.00(.)(0.0) and ˆβ = 0.0(0.)(0.0). The explanatory variable x is marginally significant for the model at the significance level of 0%. In order to assess if the model is appropriate, the empirical survival function and the estimated survival function () from the fitted LEGG regression model are plotted in Figure a. In fact, the LEGG regression model provides a good fit for these data.
19 Page of Example : Voltage data Lawless (00) reports an experiment in which specimens of solid epoxy electrical-insulation were studied in an accelerated voltage life test. The sample size was n = 0 and the percentage of censored observations was 0%. We considered three levels of voltage.,.0 and.. The variables involved in the study are: t i - failure time for epoxy insulation specimens (in min); cens i - censoring indicator (0=censoring, =lifetime observed); x i - voltage (kv). Now, we consider the model y i = β 0 + β x i + σz i, where the random variable y i follows the LEGG distribution () for i =,,..., 0. We choose the value of q that maximizes the likelihood function over selected values of q (, ) yielding ˆq = 0.. Hence, this value is assumed for the MLE of the parameter q. Iterative maximization of the logarithm of the likelihood function () starts with initial values for β and σ taken from the fitted log-weibull regression model (with λ = ). The MLEs (approximate standard errors and p-values in parentheses) are: ˆλ =.0(0.), ˆσ = 0.0(0.), ˆβ 0 =.0(.)(0.000) and ˆβ = 0.(0.0)(0.00). The fitted LEGG regression model shows that x is significant at (%) and that there is a significant difference between the voltages.,.0 and. for the survival times. In order to assess if the model is appropriate, the empirical survival function and estimated survival functions of the LEGG regression model are plotted in Figure b for different voltage levels. We conclude that the LEGG regression model provides a good fit to these data. Concluding Remarks We introduce a generalized form of the exponentiated generalized gamma (EGG) distribution whose hazard rate function accommodates the four types of shape forms, i.e. increasing, decreasing, bathtub and unimodal. We derive expansions for its moments, moment generating function, mean deviations, reliability, order statistics and their moments. Further, we define the log-exponentiated generalized gamma (LEGG) distribution and propose a LEGG regression model very suitable for modeling censored and uncensored lifetime data. The new regression model serves as a good alternative for lifetime data analysis, since we can adopt goodness of fit tests for several widely known regression models as special sub-models. We demonstrate in two applications to real data that the LEGG model can produce better fits than the usual models. Acknowledgment: The authors are grateful to the Editor and an anonymous referee for very useful comments and suggestions. This work was supported by CNPq and CAPES.
20 Page 0 of Journal of Statistical Computation and Simulation S(y x) (a) Kaplan Meier LEGG Regression model (x=0) LEGG Regression model (x=) 0 y S(y x) (b) Kaplan Meier LEGG Regression model (x=.) LEGG Regression model (x=.0) LEGG Regression model (x=.) Figure : Estimated survival function by fitting the LEGG regression model and the empirical survival (a) For each level of the diabetic retinopathy study. (b) For each level of the voltage data. Appendix A We derive an expansion for γ (k, x) λ for any λ > 0 real non-integer. By simple binomial expansion since 0 < γ (k, x) <, we can write γ (k, x) λ = j=0 m=0 j ( ) ( ) λ j ( ) j+m γ (k, x) m, () j m which always converges. We can substitute j j=0 m=0 for m=0 j=m to obtain γ (k, x) λ = m=0 s m (λ) γ (k, x) m, () y
21 Page of where s m (λ) = ( ) ( ) λ j ( ) j+m. () j m j=m A power series expansion for the incomplete gamma function ratio is given by xk γ (k, x) = Γ(k) i=0 ( x) i (k + i)i!. () We use an equation in Section 0. of Gradshteyn and Ryzhik (000) for a power series raised to a positive integer m ( ) m a i x i = c m,i x i, () i=0 where the coefficients c m,i (for i =,,...) satisfy the recurrence relation c m,i = (ia 0 ) i=0 i (mp i + p) a p c m,i p (0) p= and c m,0 = a m 0. Here, c m,i can be calculated from c m,0,..., c m,i and also be expressed explicitly as a function of a 0,..., a i, although it is not necessary for programming numerically our expansions using any software with numerical facilities. By equation (), we obtain xkm γ (k, x) m = Γ(k) m c m,i x i, () whose coefficients c m,i are determined from (0) with a p = ( ) p /(k + p)p!. Combining () and (), we can rewrite () as γ (k, x) λ = m,i=0 i=0 s m (λ) c m,i Γ(k) m xkm+i. () Appendix B Proof of Theorem : From () and (), the EGG(α, τ, k, λ) density function can be written as [ f(t) = λ τ ( ) ] t τ α Γ(k) exp ( ) s m (λ) c m,i t τ[k(m+)+i] α Γ(k) m, t > 0. α m,i=0 0
22 Page of Journal of Statistical Computation and Simulation We can express f(t) as a linear combination given by f(t) = m,i=0 w(m, i, k, λ) g α,τ,k(m+)+i (t), t > 0, where g α,τ,k(m+)+i (t) denotes the GG(α, τ, k(m + )i) density function defined by () and the weighted coefficients w(m, i, k, λ) are w(m, i, k, λ) = m,i=0 λ s m (λ) c m,i Γ[k(m + ) + i] Γ(k) m+. Proof of Theorem : The rth moment of the EGG distribution comes from () as µ r = λ τ αr Γ(k) Setting x = ( t α) τ, we have 0 µ r = λ αr sgn(τ) Γ(k) Equation () for γ (k, x) λ leads to γ (k, x) λ = ( ) t τk+r ( ) ] [ t τ ( ) ]} t τ λ exp[ {γ k, dt. α α α 0 j=0 m=0 x k+ r τ exp( x) γ (k, x) λ dx. () j ( ) ( ) λ j ( ) j+m γ (k, x) m. j m Inserting the last equation into () and interchanging terms gives where µ r = λ αr sgn(τ) Γ(k) I (k, r ) τ, m = 0 j=0 m=0 j ( ) ( ) λ j ( ) j+m I (k, r ) j m τ, m, 0 x k+ r τ exp( x) γ (k, x) m dx. We calculate the last integral using the series expansion () for the incomplete gamma function I (k, r ) [ ] m τ, m = x k+ r τ x k ( x) i exp( x) dx. Γ(k) (k + i)i! i=0
23 Page of This integral can be determined from equations () and () of Nadarajah (00b) in terms of the Lauricella function of type A (Exton, ; Aarts, 000) defined by F (n) A (a; b,..., b n ; c,..., c n ; x,..., x n ) = (a) m +...+m n (b ) m (b n ) mn x m... x mn n (c ) m (c n ) mn m!... m n!, m =0 m n =0 where (a) i = a(a + ) (a + i ) is the ascending factorial defined by (with the convention that (a) 0 = ). Numerical routines for the direct computation of the Lauricella function of type A are available, see Exton () and Mathematica (Trott, 00). We obtain I (k, r ) τ, m = k m Γ ( r/τ + k(m + ) ) F (m) ( ) r/τ + k(m + ); k,..., k; k +,..., k + ;,...,. A Proof of Theorem : Suppose that T is a random variable having a GG(α, τ, k) density function (). First, we provide a closed form expression for the mgf of T using the Wright function (Wright, ). Setting u = t/α, we have M α,τ,k (s) = τ Γ(k) 0 exp(sαu) u τk exp( u τ )du. Expanding the exponential in Taylor series and using 0 u τk+m exp( u τ )du = τ Γ(k + m/τ), we obtain M α,τ,k (s) = sgn(τ) ( m ) (sα) m Γ Γ(k) τ + k. () m! m=0 Equation () holds for any τ different from zero. However, if τ >, we can simplify it by using the Wright generalized hypergeometric function (Wright, ) defined by [ ( ) ( ) ] α, A,, αp, A p p pψ q ( ) ( β, B,, βq, B q ) ; x j= = Γ(α j + A j m) x m q j= Γ(β j + B j m) m!. () This function exists if + q j= B j p j= A j > 0. By combining () and (), we have M α,τ,k (s) = [ ] (, /τ) Γ(k) Ψ 0 ; sα. () m=0 Finally, the mgf of the EGG distribution follows from () as M(t) = m,i=0 w(m, i, k, λ) M α,τ,k(m+)+i (s),
24 Page of Journal of Statistical Computation and Simulation where M α,τ,k(m+)+i (s) is easily obtained from equations () or (). Proof of Theorem : The rth moment of the LEGG distribution can be obtained from () For q 0, we have µ r = E(Y r ) = y r λ q ( q ) q { ( ) [ ( )]} y µ y µ σγ(q exp q q exp q ) σ σ { [ ( )]}) y µ λ (γ q, q exp q dy. σ Setting x = q exp [ q ( y µ )] σ in the last equation yields µ r = λ sgn(q) { } σ r Γ (q [log(x) + log( q )] + µ x q ( exp( x) γ q, x ) λ dx. ) q 0 Using expansion () for γ (q, x) λ leads to γ (q, x) λ = m,i=0 s m (λ) c m,i Γ (q ) m x q m+i. Inserting the last equation into µ r and interchanging terms, we obtain µ r = λ sgn(q) Γ (q ) 0 { } σ r [log(x) + log( q )] + µ x q exp( x) q The binomial expansion in { σ q [log(x) + log( q )] + µ } r gives µ r = σl q m,i=0 l=0 0 m,i=0 r ( ) [ λ sgn(q)s m (λ) c m,i r σ Γ (q ) m+ l q log( q ) + µ log(x) exp( x) x q (m+)+i dx. s m (λ) c m,i Γ (q ) m x q m+i. The last integral is given in Prudnikov et al. (, vol, Section.., integral ) and then µ r = r m,i=0 l=0 where Γ(p) = Γ(p) p. λ sgn(q)s m (λ)c m,i Γ (q ) m+ ] r l ( ) [ ] r σ r l [ σl l q log( q ) + µ Γ(q q (m + ) + i)],
25 Page of For q = 0, we obtain µ r = E(Y r ) = λ πσ [ y r exp ( y µ ) ] Setting y = µ + σz and using () for Φ (λ ) (z), we have µ r = λ p=0 j=0 σ r ( ) r s p (λ) σ j µ r j τ j,p, j Φ (λ )( y µ σ ) dy. where s p (λ) is given by (). We define the integral for j and p non-negative integers τ j,p = z j φ(z)φ p (z)dz, where φ(z) is the standard normal density. The standard normal cdf can be written as Φ(x) = { ( )} x + erf, x R. By the binomial expansion and interchanging terms, we obtain τ j,p = j π p ( ) p ( ) x p l x j exp( x /) erf dx. l l=0 Using the series expansion for the error function erf(.) erf(x) = ( ) m x m+ π (m + )m!, m=0 we can solve the last integral by equations ()-() of Nadarajah (00b). When j + p l is even, we obtain τ j,p = j/ π (p+/) p F (p l) A l=0 (j+p l) even ( j + p l + ( ) ( ) p j + p l + l π l Γ l ;,, ;,, ;,, ). Appendix C: Matrix of second derivatives L(θ)
26 Page of Journal of Statistical Computation and Simulation Here, we provide the formulas to obtain the second-order partial derivatives of the log-likelihood function. After some algebraic manipulations, we obtain { L α,α = nτk α τ n (λ )τ n vi k α (v i + τv i ) exp( v i) α τkvk i exp( v i) Γ(k) γ i= i= (k, v i ) γ (k, v i ) } + τv(k+) i exp( v i ) + τ [vi k] [exp( v i )] γ (k, v i ) Γ(k) [γ (k, v i )], { L α,τ = nk α + n (λ ) n vi k [v i + v i log(v i )] exp( v i) + vk i exp( v i) log(vi k) α αγ(k) γ i= i= (k, v i ) γ (k, v i ) } v(k+) i exp( v i ) log(v i ) [vi k] [exp( v i )] log(v i ) γ (k, v i ) Γ(k) [γ (k, v i )], L α,k = nτ α L α,λ = τ αγ(k) τ(λ ) αγ(k) { } n vi k log(v i ) exp( v i ) γ (k, v i ) [ γ (k, v i )] k Γ(k)[γ (k, v i )], i= n vi kexp( v i), L τ,λ = γ (k, v i ) τγ(k) i= n vi kexp( v i) log(v i ), γ (k, v i ) { L τ,τ = n τ n τ v i [log(v i )] (λ ) n vi k + exp( v i) log(v i ) k+ τ Γ(k) γ i= i= (k, v i ) } (v i) k+ exp( v i )[log(v i )] [vk i ] [exp( v i )] [log(v i )] γ (k, v i ) Γ(k)[γ (k, v i )], L τ,k = τ n (λ ) log(v i ) + τγ(k) i= L k,k = nλψ (λ ) (k) + Γ(k) L k,λ = nψ(k) + Γ(k) i= { } n vi k log(v i ) exp( v i ) log(v i ) γ (k, v i ) [ γ (k, v i )] k Γ(k)[γ (k, v i )], i= { } n [ γ (k, v i )] kk {[ γ (k, v i )] k } γ (k, v i ) Γ(k)[γ (k, v i )], i= { } n [ γ (k, v i )] k, L λ,λ = n γ (k, v i ) λ, i=
27 Page of where v i = ( t i ) vi τ, [ γ (k, v i )] k = w k exp( w) log(w)dw, α 0 [ γ (k, v i )] kk = vi 0 w k exp( w)[log(w)] dw, ψ(.) is the digamma function and ψ (.) is the polygamma function. References Aarts, R.M. (000). Lauricella functions, From MathWorld - A Wolfram Web Resource, created by Eric W. Weisstein. Ali, M.M., Woo, J. and Nadarajah, S. (00). Generalized gamma variables with drought application. Journal of the Korean Statistical Society,, -. Almpanidis, G. and Kotropoulos, C. (00). Phonemic segmentation using the generalized gamma distribution and small sample Bayesian information criterion. Speech Communication, 0, -. Barkauskas, D.A., Kronewitter, S.R., Lebrilla, C.B. and Rocke, D.M. (00). Analysis of MALDI FT-ICR mass spectrometry data: A time series approach. Analytica Chimica Acta,, 0-. Cancho, V.G., Bolfarine, H. and Achcar, J.A. (). A Bayesian Analysis for the Exponentiated- Weibull Distribution. Journal of Applied Statistics,, -. Cancho, V.G., Ortega, E.M.M. and Bolfarine, H. (00). The Exponentiated-Weibull Regression Models with Cure Rate. Journal of Applied Probability and Statistics. In press. Carrasco, J.M.F., Ortega, E.M.M. and Cordeiro, G.M. (00). A generalized modified Weibull distribution for lifetime modeling. Computational Statistics and Data Analysis,, 0-. Cordeiro, G.M., Ortega, E.M.M. and Silva, G.O. (00). The exponentiated generalized gamma distribution with application to lifetime data. Journal of Statistical Computation and Simulation. DOI: 0.00/0000. Cox, C. (00). The generalized F distribution: an umbrella for parametric survival analysis. Statistics in Medicine,, 0-.
28 Page of Journal of Statistical Computation and Simulation Cox, C., Chu, H., Schneider, M.F. and Muñoz, A. (00). Parametric survival analysis and taxonomy of hazard functions for the generalized gamma distribution. Statistics in Medicine,, -. Exton, H. (). Handbook of Hypergeometric Integrals: Theory, Applications, Tables, Computer Programs. New York, Halsted Press. Gomes, O., Combes, C. and Dussauchoy, A. (00). Parameter estimation of the generalized gamma distribution. Mathematics and Computers in Simulation,, -. Gradshteyn, I.S. and Ryzhik, I.M. (000). Table of Integrals, Series and Products, th edition. Edited by A. Jeffrey and D. Zwillinger. Gupta, R.D. and Kundu, D. (). Generalized exponential distributions. Australian and New Zealand Journal of Statistics,, -. Gusmão, F.R.S., Ortega, E.M.M. and Cordeiro, G.M. (00). The generalized inverse Weibull distribution. Statistical Papers. DOI 0.00/s Hashimoto, E.M., Ortega, E.M.M., Cancho, V.G. and Cordeiro, G.M. (00). The log-exponentiated Weibull regression model for interval-censored data. Computational Statistics and Data Analysis,, 0-0. Huster, W.J., Brookmeyer, R. and Self, S.G. (). Modelling Paired Survival Data with Covariates. Biometrics,, -. Kundu, D. and Raqab, M.Z. (00). Generalized Rayleigh distribution: different methods of estimation. Computational Statistics and Data Analysis,, -00. Lai, C.D., Xie, M. and Murthy, D.N.P. (00). A modified Weibull distribution. Transactions on Reliability,, -. Liang, K., Self, S.G. and Chang, Y. (). Modelling marginal hazards in multivariate failure time data. Journal of the Royal Statistical Society. Series B (Methodological),, -. Lawless, J.F. (00). Statistical Models and Methods for Lifetime Data. New York, John Wiley. Malhotra, J., Sharma, A.K. and Kaler, R.S. (00). On the performance analysis of wireless receiver using generalized-gamma fading model. Annals of Telecommunications,, -. Mudholkar, G.S., Srivastava, D.K. and Freimer, M. (). The exponentiated Weibull family: A reanalysis of the bus-motor-failure data. Technometrics,, -. Nadarajah, S. (00a). On the use of the generalized gamma distribution. International Journal of Electronic,, 0-0.
29 Page of Nadarajah, S. (00b). Explicit expressions for moments of order statistics. Statistics and Probability Letters,, -0. Nadarajah, S. and Gupta, A.K. (00). A generalized gamma distribution with application to drought data. Mathematics and Computers in Simulation,, -. Ortega, E.M.M., Bolfarine, H. and Paula, G.A. (00). Influence diagnostics in generalized loggamma regression models. Computational Statistics and Data Analysis,, -. Ortega, E.M.M., Cancho, V.G. and Bolfarine, H. (00). Influence diagnostics in exponentiated -Weibull regression models with censored data. Statistics and Operations Research Transactions, 0, -. Ortega, E.M.M., Paula G.A. and Bolfarine, H. (00). Deviance residuals in generalized loggamma regression models with censored observations. Journal of Statistical Computation and Simulation,, -. Ortega, E.M.M., Cancho, V.G. and Paula, G.A. (00). Generalized log-gamma regression models with cure fraction. Lifetime Data Analysis., -0. Prudnikov, A.P., Brychkov, Y.A. and Marichev, O.I. (). Integrals and Series, volume. Amsterdam, Gordon and Breach Science Publishers. Stacy, E.W. (). A generalization of the gamma distribution. Annals of Mathematical Statistics,, -. Stacy, E.W. and Mihram, G.A. (). Parameter estimation for a generalized gamma distribution. Technometris,, -. Trott, M. (00). The Mathematica Guidebook for Symbolics. With DVD ROM (Windows, Macintosh and UNIX). New York, Springer. Xie, M. and Lai, C.D. (). Reliability analysis using an additive Weibull model with bathtubshaped failure rate function. Reliability Engineering and System Safety,, -. Xie, X. and Liu, X. (00). Analytical three-moment autoconversion parameterization based on generalized gamma distribution. Journal of Geophysical Research,, D0, doi: 0.0/00JD0. Wada, C.Y. and Hotta, L.K. (000). Restricted alternatives tests in a bivariate exponential model with covariates. Communication in Statistics - Theory and Methods,, -0. Wright, E.M. (). The asymptotic expansion of the generalized Bessel function. Proceedings of the London Mathematical Society,, 0.
30 Page 0 of Journal of Statistical Computation and Simulation xmm (00 x 00 DPI)
31 Page of xmm (00 x 00 DPI)
32 Page of Journal of Statistical Computation and Simulation xmm (00 x 00 DPI)
33 Page of xmm (00 x 00 DPI)
34 Page of Journal of Statistical Computation and Simulation xmm (00 x 00 DPI)
35 Page of xmm (00 x 00 DPI)
36 Page of Journal of Statistical Computation and Simulation xmm (00 x 00 DPI)
37 Page of xmm (00 x 00 DPI)
38 Page of Journal of Statistical Computation and Simulation xmm (00 x 00 DPI)
39 Page of xmm (00 x 00 DPI)
40 Page 0 of Journal of Statistical Computation and Simulation xmm (00 x 00 DPI)
41 Page of xmm (00 x 00 DPI)
The 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 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 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 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 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 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 informationUniversidade de São Paulo, Brazil c Departamento de Estatística, Universidade Federal de. Pernambuco, Brazil. Available online: 17 Aug 2011
This article was downloaded by: [USP University of Sao Paulo] On: 17 August 2011, At: 02:46 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationThe Gamma Modified Weibull Distribution
Chilean Journal of Statistics Vol. 6, No. 1, April 2015, 37 48 Distribution theory Research Paper The Gamma Modified Weibull Distribution Gauss M. Cordeiro, William D. Aristizábal, Dora M. Suárez and Sébastien
More informationThe Kumaraswamy Generalized Half-Normal Distribution for Skewed Positive Data
Journal of Data Science 10(01, 195-4 The Kumaraswamy Generalized Half-Normal Distribution for Skewed Positive Data Gauss M. Cordeiro 1, Rodrigo R. Pescim and Edwin M. M. Ortega 1 UFPE and ESALQ - USP Abstract:
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 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 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 informationOn the product of independent Generalized Gamma random variables
On the product of independent Generalized Gamma random variables Filipe J. Marques Universidade Nova de Lisboa and Centro de Matemática e Aplicações, Portugal Abstract The product of independent Generalized
More informationINFERENCE FOR MULTIPLE LINEAR REGRESSION MODEL WITH EXTENDED SKEW NORMAL ERRORS
Pak. J. Statist. 2016 Vol. 32(2), 81-96 INFERENCE FOR MULTIPLE LINEAR REGRESSION MODEL WITH EXTENDED SKEW NORMAL ERRORS A.A. Alhamide 1, K. Ibrahim 1 M.T. Alodat 2 1 Statistics Program, School of Mathematical
More informationTHE ODD LINDLEY BURR XII DISTRIBUTION WITH APPLICATIONS
Pak. J. Statist. 2018 Vol. 34(1), 15-32 THE ODD LINDLEY BURR XII DISTRIBUTION WITH APPLICATIONS T.H.M. Abouelmagd 1&3, Saeed Al-mualim 1&2, Ahmed Z. Afify 3 Munir Ahmad 4 and Hazem Al-Mofleh 5 1 Management
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 informationINVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION
Pak. J. Statist. 2017 Vol. 33(1), 37-61 INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION A. M. Abd AL-Fattah, A.A. EL-Helbawy G.R. AL-Dayian Statistics Department, Faculty of Commerce, AL-Azhar
More informationThe Distributions of Sums, Products and Ratios of Inverted Bivariate Beta Distribution 1
Applied Mathematical Sciences, Vol. 2, 28, no. 48, 2377-2391 The Distributions of Sums, Products and Ratios of Inverted Bivariate Beta Distribution 1 A. S. Al-Ruzaiza and Awad El-Gohary 2 Department of
More informationSome Statistical Properties of Exponentiated Weighted Weibull Distribution
Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 4 Issue 2 Version. Year 24 Type : Double Blind Peer Reviewed International Research Journal Publisher: Global Journals
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 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 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 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 informationTransmuted Exponentiated Gumbel Distribution (TEGD) and its Application to Water Quality Data
Transmuted Exponentiated Gumbel Distribution (TEGD) and its Application to Water Quality Data Deepshikha Deka Department of Statistics North Eastern Hill University, Shillong, India deepshikha.deka11@gmail.com
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 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 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 informationTruncated Life Test Sampling Plan Under Odd-Weibull Distribution
International Journal of Mathematics Trends and Technology ( IJMTT ) Volume 9 Number 2 - July Truncated Life Test Sampling Plan Under Odd-Weibull Distribution G.Laxshmimageshpraba 1, Dr.S.Muthulakshmi
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 informationThe Kumaraswamy-Burr Type III Distribution: Properties and Estimation
British Journal of Mathematics & Computer Science 14(2): 1-21, 2016, Article no.bjmcs.19958 ISSN: 2231-0851 SCIENCEDOMAIN international www.sciencedomain.org The Kumaraswamy-Burr Type III Distribution:
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 informationEXPLICIT EXPRESSIONS FOR MOMENTS OF χ 2 ORDER STATISTICS
Bulletin of the Institute of Mathematics Academia Sinica (New Series) Vol. 3 (28), No. 3, pp. 433-444 EXPLICIT EXPRESSIONS FOR MOMENTS OF χ 2 ORDER STATISTICS BY SARALEES NADARAJAH Abstract Explicit closed
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 informationTwo Weighted Distributions Generated by Exponential Distribution
Journal of Mathematical Extension Vol. 9, No. 1, (2015), 1-12 ISSN: 1735-8299 URL: http://www.ijmex.com Two Weighted Distributions Generated by Exponential Distribution A. Mahdavi Vali-e-Asr University
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 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 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 informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Estimation of the primary hazard ratio in the presence of a secondary covariate with non-proportional hazards An undergraduate honors thesis submitted to the Department
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 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 informationSurvival Analysis. Stat 526. April 13, 2018
Survival Analysis Stat 526 April 13, 2018 1 Functions of Survival Time Let T be the survival time for a subject Then P [T < 0] = 0 and T is a continuous random variable The Survival function is defined
More informationSurvival Distributions, Hazard Functions, Cumulative Hazards
BIO 244: Unit 1 Survival Distributions, Hazard Functions, Cumulative Hazards 1.1 Definitions: The goals of this unit are to introduce notation, discuss ways of probabilistically describing the distribution
More informationDistribution Fitting (Censored Data)
Distribution Fitting (Censored Data) Summary... 1 Data Input... 2 Analysis Summary... 3 Analysis Options... 4 Goodness-of-Fit Tests... 6 Frequency Histogram... 8 Comparison of Alternative Distributions...
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 informationA GENERALIZED CLASS OF EXPONENTIATED MODI ED WEIBULL DISTRIBUTION WITH APPLICATIONS
Journal of Data Science 14(2016), 585-614 A GENERALIZED CLASS OF EXPONENTIATED MODI ED WEIBULL DISTRIBUTION WITH APPLICATIONS Shusen Pu 1 ; Broderick O. Oluyede 2, Yuqi Qiu 3 and Daniel Linder 4 Abstract:
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 informationThe Gamma Generalized Pareto Distribution with Applications in Survival Analysis
International Journal of Statistics and Probability; Vol. 6, No. 3; May 217 ISSN 1927-732 E-ISSN 1927-74 Published by Canadian Center of Science and Education The Gamma Generalized Pareto Distribution
More informationLogistic regression model for survival time analysis using time-varying coefficients
Logistic regression model for survival time analysis using time-varying coefficients Accepted in American Journal of Mathematical and Management Sciences, 2016 Kenichi SATOH ksatoh@hiroshima-u.ac.jp Research
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 informationA Bivariate Weibull Regression Model
c Heldermann Verlag Economic Quality Control ISSN 0940-5151 Vol 20 (2005), No. 1, 1 A Bivariate Weibull Regression Model David D. Hanagal Abstract: In this paper, we propose a new bivariate Weibull regression
More informationUniversity of Lisbon, Portugal
Development and comparative study of two near-exact approximations to the distribution of the product of an odd number of independent Beta random variables Luís M. Grilo a,, Carlos A. Coelho b, a Dep.
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 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 informationA Marshall-Olkin Gamma Distribution and Process
CHAPTER 3 A Marshall-Olkin Gamma Distribution and Process 3.1 Introduction Gamma distribution is a widely used distribution in many fields such as lifetime data analysis, reliability, hydrology, medicine,
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 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 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 informationChapter 17. Failure-Time Regression Analysis. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 17 Failure-Time Regression Analysis William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
More informationGOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS
Statistica Sinica 20 (2010), 441-453 GOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS Antai Wang Georgetown University Medical Center Abstract: In this paper, we propose two tests for parametric models
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS330 / MAS83 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-0 8 Parametric models 8. Introduction In the last few sections (the KM
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 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 informationFINITE MIXTURES OF LOGNORMAL AND GAMMA DISTRIBUTIONS
The 7 th International Days of Statistics and Economics, Prague, September 9-, 03 FINITE MIXTURES OF LOGNORMAL AND GAMMA DISTRIBUTIONS Ivana Malá Abstract In the contribution the finite mixtures of distributions
More informationClosed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
The American Statistician ISSN: 3-135 (Print) 1537-2731 (Online) Journal homepage: http://www.tandfonline.com/loi/utas2 Closed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS3301 / MAS8311 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-10 1 13 The Cox proportional hazards model 13.1 Introduction In the
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 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 informationComputation of Signal to Noise Ratios
MATCH Communications in Mathematical in Computer Chemistry MATCH Commun. Math. Comput. Chem. 57 7) 15-11 ISS 34-653 Computation of Signal to oise Ratios Saralees adarajah 1 & Samuel Kotz Received May,
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 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 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 informationTests of independence for censored bivariate failure time data
Tests of independence for censored bivariate failure time data Abstract Bivariate failure time data is widely used in survival analysis, for example, in twins study. This article presents a class of χ
More informationThe Zografos-Balakrishnan odd log-logistic family of distributions: Properties and Applications
The Zografos-Balakrishnan odd log-logistic family of distributions: Properties and Applications Gauss M. Cordeiro Morad Alizadeh Edwin M.M.Ortega and Luis H. Valdivieso Serrano Abstract We study some mathematical
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 informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationSTAT 6350 Analysis of Lifetime Data. Failure-time Regression Analysis
STAT 6350 Analysis of Lifetime Data Failure-time Regression Analysis Explanatory Variables for Failure Times Usually explanatory variables explain/predict why some units fail quickly and some units survive
More informationContents 1. Contents
Contents 1 Contents 6 Distributions of Functions of Random Variables 2 6.1 Transformation of Discrete r.v.s............. 3 6.2 Method of Distribution Functions............. 6 6.3 Method of Transformations................
More informationSurvival Analysis Math 434 Fall 2011
Survival Analysis Math 434 Fall 2011 Part IV: Chap. 8,9.2,9.3,11: Semiparametric Proportional Hazards Regression Jimin Ding Math Dept. www.math.wustl.edu/ jmding/math434/fall09/index.html Basic Model Setup
More informationTHE BETA WEIBULL-G FAMILY OF DISTRIBUTIONS: THEORY, CHARACTERIZATIONS AND APPLICATIONS
Pak. J. Statist. 2017 Vol. 33(2), 95-116 THE BETA WEIBULL-G FAMILY OF DISTRIBUTIONS: THEORY, CHARACTERIZATIONS AND APPLICATIONS Haitham M. Yousof 1, Mahdi Rasekhi 2, Ahmed Z. Afify 1, Indranil Ghosh 3,
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 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 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 Transmuted Generalized Gamma Distribution: Properties and Application
Journal of Data Science 13(2015), 187-206 The Transmuted Generalized Gamma Distribution: Properties and Application Sadraque E.F. Lucena 1, Ana Herm ınia A. Silva 2, Gauss M. Cordeiro 3 Universidade Federal
More informationLogistic-Modified Weibull Distribution and Parameter Estimation
International Journal of Contemporary Mathematical Sciences Vol. 13, 2018, no. 1, 11-23 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijcms.2018.71135 Logistic-Modified Weibull Distribution and
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 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 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 informationReliability Engineering I
Happiness is taking the reliability final exam. Reliability Engineering I ENM/MSC 565 Review for the Final Exam Vital Statistics What R&M concepts covered in the course When Monday April 29 from 4:30 6:00
More informatione 4β e 4β + e β ˆβ =0.765
SIMPLE EXAMPLE COX-REGRESSION i Y i x i δ i 1 5 12 0 2 10 10 1 3 40 3 0 4 80 5 0 5 120 3 1 6 400 4 1 7 600 1 0 Model: z(t x) =z 0 (t) exp{βx} Partial likelihood: L(β) = e 10β e 10β + e 3β + e 5β + e 3β
More informationFirst Year Examination Department of Statistics, University of Florida
First Year Examination Department of Statistics, University of Florida August 20, 2009, 8:00 am - 2:00 noon Instructions:. You have four hours to answer questions in this examination. 2. You must show
More informationModified slash Birnbaum-Saunders distribution
Modified slash Birnbaum-Saunders distribution Jimmy Reyes, Filidor Vilca, Diego I. Gallardo and Héctor W. Gómez Abstract In this paper, we introduce an extension for the Birnbaum-Saunders (BS) distribution
More informationSubject CS1 Actuarial Statistics 1 Core Principles
Institute of Actuaries of India Subject CS1 Actuarial Statistics 1 Core Principles For 2019 Examinations Aim The aim of the Actuarial Statistics 1 subject is to provide a grounding in mathematical and
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 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 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 informationAn Extended Weighted Exponential Distribution
Journal of Modern Applied Statistical Methods Volume 16 Issue 1 Article 17 5-1-017 An Extended Weighted Exponential Distribution Abbas Mahdavi Department of Statistics, Faculty of Mathematical Sciences,
More 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 informationInternational Journal of Physical Sciences
Vol. 9(4), pp. 71-78, 28 February, 2014 DOI: 10.5897/IJPS2013.4078 ISSN 1992-1950 Copyright 2014 Author(s) retain the copyright of this article http://www.academicjournals.org/ijps International Journal
More information