Semi-parametric Inference for Cure Rate Models 1
|
|
- Aubrie Horn
- 6 years ago
- Views:
Transcription
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, This research is supported by a Marie Curie International Outgoing Fellowship within the 7th European Community Framework Programme.
2 Outline 1 Cure rate models 2 Piecewise linear approximation 3 Simulation study 4 Illustrative example 5 Conclusions
3 Introduction The primary aim of this work is the study of survival times or more generally, the times till the occurrence of an event. This event may be: the occurrence or recurrence of a disease; the return to prison or rearrest of a released prisoner; the death of a patient; etc
4 Introduction The primary aim of this work is the study of survival times or more generally, the times till the occurrence of an event. This event may be: the occurrence or recurrence of a disease; the return to prison or rearrest of a released prisoner; the death of a patient; etc Using the models found in traditional survival analysis, every item will experience the event of interest, at some point.
5 Introduction The primary aim of this work is the study of survival times or more generally, the times till the occurrence of an event. This event may be: the occurrence or recurrence of a disease; the return to prison or rearrest of a released prisoner; the death of a patient; etc Using the models found in traditional survival analysis, every item will experience the event of interest, at some point. population Traditional survival analysis the event of interest... after some time...
6 Introduction Cure rate models allow for a proportion of items which will never experience the event of interest.
7 Introduction Cure rate models allow for a proportion of items which will never experience the event of interest. Cure rate models the event of interest (susceptibles) population... after some time... cured items (non-susceptibles) These items are called cured or non-susceptibles or long-term survivors or immune or immortals etc.
8 Introduction The cure rate models have been extensively studied during the last decades. We can find more than 200 publication the last years and manymany applications... The reason is the great improvement of treatments for a range of diseases, but also to the fact that in many social phenomena a number of individuals are not susceptible to the event of interest. Applications may be found in biomedical studies, criminology, finance, industrial reliability etc (see the monographs by Maller and Zhou, 1996; Ibrahim, Chen and Sinha, 2005).
9 Binary cure rate model Boag (1949) and Berkson and Gage (1952) are the first works (to the best of my knowledge) which took into account the existence of two subgroups in the population: the susceptibles and the non -susceptibles (cured).
10 Binary cure rate model Boag (1949) and Berkson and Gage (1952) are the first works (to the best of my knowledge) which took into account the existence of two subgroups in the population: the susceptibles and the non -susceptibles (cured). Based on this assumption, the population survival function of the timeto-event T, is given by the following mixture model S P (t)=p(t > t)=p 0 + (1 p 0 )S(t) where p 0 is the probability of a patient to be cured, and S(t) = P(T > t susceptibles) is the survival function of susceptibles. Note that lim t S P (t)=p 0.
11 Binary cure rate model What problems do we have? What is the distribution of susceptibles S(t) = P(T > t susceptibles)? How the covariates can be incorporated into the model (in p 0 or/and S(t) = P(T > t susceptibles))? How to estimate the parameters of the model?
12 Binary cure rate model S P (t)=p(t > t)=p 0 + (1 p 0 )S(t) What problems do we have? What is the distribution of susceptibles S(t) = P(T > t susceptibles)? How the covariates can be incorporated into the model (in p 0 or/and S(t) = P(T > t susceptibles))? How to estimate the parameters of the model?
13 Binary cure rate model S P (t)=p(t > t)=p 0 + (1 p 0 )S(t) What problems do we have? What is the distribution of susceptibles S(t) = P(T > t susceptibles)? How the covariates can be incorporated into the model (in p 0 or/and S(t) = P(T > t susceptibles))? How to estimate the parameters of the model?
14 Binary cure rate model S P (t)=p(t > t)=p 0 + (1 p 0 )S(t) What problems do we have? What is the distribution of susceptibles S(t) = P(T > t susceptibles)? How the covariates can be incorporated into the model (in p 0 or/and S(t) = P(T > t susceptibles))? How to estimate the parameters of the model?
15 Binary cure rate model S P (t)=p(t > t)=p 0 + (1 p 0 )S(t) What problems do we have? What is the distribution of susceptibles S(t) = P(T > t susceptibles)? How the covariates can be incorporated into the model (in p 0 or/and S(t) = P(T > t susceptibles))? How to estimate the parameters of the model? Parametric, semi-parametric and non-parametric approaches can be found in the relevant literature.
16 Binary cure rate model Assumptions for the distribution of susceptibles: Exponential distribution (e.g. Berkson and Gage, 1952; Mould, 1973; Ghitany and Maller, 1992; Ghitany, 1993; Farewell, 1977a; Balakrishnan and Pal, 2013c) S(t)=P(T > t susceptibles)=exp( λt),λ>0; Weibull distribution (e.g. Farewell, 1977b, 1982, 1986; Bentzen et al.,1989; Struthers and Farewell, 1989; Balakrishnan and Pal, 2013b,c) S(t)=exp( λt) α,λ,α>0;
17 Binary cure rate model Assumptions for the distribution of susceptibles: Cox s proportional hazard model (e.g. Kuk and Chen, 1992; Sy and Taylor, 2000, 2001; Peng and Dear, 2000; Fang, Li and Sun, 2005; Zhao et al., 2014) h(t)=h 0 (t)exp(x β),λ>0; Kaplan-Meyer estimator (e.g. Taylor, 1995; Laska and Meisner, 1992; Maller and Zhou, 1992); Generalized F distribution (e.g. Peng, Dear and Denham, 1998); Log-normal distribution (Balakrishnan and Pal, 2013a) and others.
18 Binary cure rate model The estimation depends on the available data (censoring, truncation etc): MLE and EM algorithm; A marginal likelihood estimation approach (especially in Cox s PH model); An iterative least squares method (Berkson and Gage, 1952).
19 Binary cure rate model The estimation depends on the available data (censoring, truncation etc): MLE and EM algorithm; A marginal likelihood estimation approach (especially in Cox s PH model); An iterative least squares method (Berkson and Gage, 1952). Logistic regression model for p 0 (logistic-binary mixture model), i.e. p 0 = p 0 (x;β)= 1 1+exp(x β).
20 Competing cause scenario A number of competing causes left alive after a treatment, say M; let W 1,W 2,...,W M be the random variable for the time-to-event due to ith (i=1,2,...,m) competing cause. The random variables W 1,W 2,...,W M are assumed i.i.d. with cdf F(t)= 1 S(t) and are also independent of M.
21 Competing cause scenario A number of competing causes left alive after a treatment, say M; let W 1,W 2,...,W M be the random variable for the time-to-event due to ith (i=1,2,...,m) competing cause. The random variables W 1,W 2,...,W M are assumed i.i.d. with cdf F(t)= 1 S(t) and are also independent of M. Therefore, the population time-to-event T is given by with P(W 0 = )=1. T = min{w 0,W 1,...,W M },
22 Competing cause scenario Then its survival function has the following form S P (t)=p(t > t)=p(w 0 > t,w 1 > t,...,w M > t) = P(M = 0)+P(W 1 > t,w 2 > t,...,w M > t,m 1) = P(M = 0)+ S(t) i P(M = i) i=1
23 Competing cause scenario Then its survival function has the following form S P (t)=p(t > t)=p(w 0 > t,w 1 > t,...,w M > t) = P(M = 0)+P(W 1 > t,w 2 > t,...,w M > t,m 1) = P(M = 0)+ S(t) i P(M = i) i=1 For example, let M follow a Poisson distribution, with mean θ (0, ); then, S P (t)=exp( θ)+ S(t) i exp( θ) θi i! = exp( θf(t)) i=1 e.g. Tsodikov, Ibrahim and Yakovlev (2003) and Ibrahim et al. (2005).
24 Competing cause scenario S P (t)=p(m = 0)+ S(t) i P(M = i) i=1 Again, the next questions must be answered: what is the distribution of M (the number of competing causes); what is the distribution of W i (time-to-event due to ith competing cause); how the covariates can be incorporated into the model; how to deal with the estimation of the parameters.
25 Competing cause scenario Assumptions for the distribution of the number of competing causes M: Poisson distribution (e.g. Yakovlev et al. 1993; Cantor and Shuster, 1992; Yakovlev, Cantor and Shuster, 1994; Chen, Ibrahim and Sinha, 1999; Hashimoto, Ortega et al., 2013); Negative binomial distribution (e.g. Castro, Cancho and Rodrigues, 2009; Ortega et al., 2014); COM-Poisson distribution (e.g. Rodrigues et al., 2009; Balakrishnan and Pal, 2012, 2013a,b,c); Compound weighted Poisson distribution (e.g. Rodrigues et al., 2011); Geometric distribution (e.g. Louzada et al., 2014);
26 Unified represantion Taking into account the relation S P (t)= S(t) i P(M = i) i=0 we can say that, both of the two classes of cure rate models (the binary and competing cause model) can be expressed as S P (t)=e[s(t) M ]=G M (S(t)), where the expectations is taken w.r.t. M and G M is the probability generating function of M (see e.g. Tsodikov et al., 2003).
27 Our semi-parametric approach In this work, we are going to study both the logistic-binary mixture and the logistic-poisson mixture models. The common hazard function h(t) of the time-to-event due to ith competing cause is approximated by a piecewise linear function. h( t ) φ N φn 1 φ 2 φ 1 φ 0 τ 0 τ1 τ 2 τ N 1 τ N t
28 Our semi-parametric approach The piecewise linear approximation (PLA), for the common hazard function h(t) of W i, is based on the following decisions: the number of lines, N, to be used; the selection of cut points τ 0 < τ 1 <...<τ N 1 < τ N for forming the line segments. Of course, we further assume that the PLA is continuous function at cut points.
29 Hence, Our semi-parametric approach h L (t)= N (c j + s j t)i [τj 1,τ j ](t), j=1 where I [τj 1,τ j ](t)=1 if and only if t [τ j 1,τ j ] and c j = φ j τ j φ j φ j 1 τ j τ j 1, s j = φ j φ j 1 τ j τ j 1, for j= 1,2,...,N, with φ j = c j + s j τ j and φ 0 = c 1 + s 1 τ 0. h( t ) φ N φn 1 φ 2 φ 1 φ 0 τ 0 τ1 τ 2 τ N 1 τ N t
30 Data and Estimation The MLE and EM algorithm is carried out for estimating the parameters of the model. We consider the scenario where the time-to-event is subject to noninformative random right censoring.
31 Data and Estimation : non-cured/ susceptibles : cured start time end of study
32 Data and Estimation Denoting with C i and T i the censoring time and lifetime of the ith individual, respectively, we then observe and δ i = I(T i C i ), i.e. δ i = Y i = min{t i,c i } { 1, if Yi is a time-to-event 0, if Y i is a censoring time,,i=1,2,...,n.
33 Data and Estimation ( y1, δ 1) = (exact lifetime,1) ( y2, δ 2) = (censoring time,0) ( y3, δ 3) = (censoring time,0) ( y4, δ 4) = (exact lifetime,1) ( y5, δ 5) = (censoring time,0) ( y6, δ 6) = (censoring time,0) ( y7, δ 7) = (exact lifetime,1) ( y8, δ 8) = (censoring time,0) ( y9, δ 9) = (censoring time,0) start time end of study
34 Data and Estimation: likelihood function From n pairs of times and censoring indicators (y 1,δ 1 ),...,(y n,δ n ), the likelihood function can be written as L=L(θ; x, y,δ) n f P (y i, x i ;θ) δ i S P (y i, x i ;θ) 1 δ i. i=1
35 Data and Estimation: likelihood function From n pairs of times and censoring indicators (y 1,δ 1 ),...,(y n,δ n ), the likelihood function can be written as L=L(θ; x, y,δ) Thus, the likelihood becomes n f P (y i, x i ;θ) δ i S P (y i, x i ;θ) 1 δ i. i=1 L(θ; x, y,δ) n (1 p 0 (x i,β)) δ i f U (y i, x i ;θ) δ i [p 0 (x i,β)+(1 p 0 (x i,β))s U (y i, x i ;θ)] 1 δ i, i=1 where S U and f U are the probability density and survival function of susceptibles, respectively.
36 Data and Estimation: likelihood function Let us now assume that I i = { 1, if the ith individual is susceptible 0, if the ith individual is cured,,i=1,2,...,n. Note that the values of I i for censored items are not observable.
37 Data and Estimation: likelihood function Let us now assume that I i = { 1, if the ith individual is susceptible 0, if the ith individual is cured,,i=1,2,...,n. Note that the values of I i for censored items are not observable. Then, the complete likelihood function can be written as L c = L c (θ; x, y,δ) (1 p 0 (x i,β)) f U (y i, x i ;θ) p 0 (x i,β) 1 I i i 1 i 1 i 0 [(1 p 0 (x i,β))s U (y i, x i ;θ)] I i i 0 (note that i 1 I i = 1).
38 Data and Estimation: EM algorithm For the E-step of the EM algorithm, we have where w (z) i = E[I i θ (z), O]= exp(x i β(z) )S U (y i, x;θ (z) ) 1+exp(x i β(z) )S U (y i, x;θ (z) ), Binary case : S U (t, x;θ)=s L (y i ;φ (z) ), Poisson case : S U (t, x;θ)= exp( θ(x,β)f L(t;φ)) exp( θ(x,β)). 1 exp( θ(x,β))
39 Numerical results A set of data of size n=300 is generated, in which the random variables W i follow a Weibull distribution. We assume that: 1. we have only one covariate with 4 possible values/groups (x {1,2,3,4}); 2. the censoring times C i follow Exponential distributions;
40 Numerical results A set of data of size n=300 is generated, in which the random variables W i follow a Weibull distribution. We assume that: 1. we have only one covariate with 4 possible values/groups (x {1,2,3,4}); 2. the censoring times C i follow Exponential distributions; We fix the cured proportions for the first and fourth group to 0.30 and 0.15, respectively; therefore, based on the equations 0.30= 1 1+exp(β 0 + β 1 1),0.15= 1 1+exp(β 0 + β 1 4) we have that β 0 = and β 1 =
41 Numerical results The cut points τ 0,τ 1,...,τ N 1 of the PLA are taken to be the sample percentiles of the uncensored data (i.e. 0,1/N,2/N,...,(N 1)/Nth percentiles); the last point τ N is the maximum of Y i. The initial values for β 0 and β 1 can be given by replacing the cured proportions with the observed censoring proportions; the initial values for φ 0,φ 1,...,φ N may be computed using a set of estimates through φ j = ĥ(τ j),j = 0,1,...,N; in our case, the Kaplan-Meier estimator was used. Moreover, we assume for each group that the censoring proportions exceed the cured proportions by c= 0.10 or c= 0.20.
42 Numerical results Table: Logistic-binary mixture model: Weibull distribution and low censored proportion (c =.10; n = 300) Cured proportions for Groups 1-4:(.3,.242,.192,.15) Censored proportions for Groups 1-4:(.3,.242,.192,.15)+c,c=.10 Sample size for Groups 1-4: (60,105,65,70),n=300 Parameters:β 0 =.5515,β 1 =.2958,α=1.5,γ=1.5(parameters of the distribution) ˆφ i,(ˆα, ˆγ) N (ˆβ 0, ˆβ 1 ) ˆp 0 l i=0 i=1 i=2 i=3 i=4 Mean 1 (.542,.310) (.301,.240,.189,.148) (.542,.311) (.301,.239,.188,.148) (.543,.311) (.301,.239,.188,.148) (.543,.311) (.301,.239,.188,.148) (.546,.296) (.304,.244,.194,.154) (1.504,1.510) RMSE 1 (.362,.148) (.050,.030,.031,.038) 2 (.362,.148) (.050,.030,.031,.038) 3 (.362,.148) (.050,.030,.031,.038) 4 (.362,.148) (.050,.030,.031,.038) (.377,.144) (.053,.031,.028,.035).007 Std 1 (.363,.148) (.050,.030,.031,.038) (.364,.148) (.050,.030,.031,.038) (.363,.148) (.050,.030,.031,.038) (.363,.148) (.050,.030,.031,.038) (.378,.145) (.053,.031,.028,.034) (.081,.070)
43 Data and Estimation
44 Numerical results Table: Logistic-binary mixture model: Weibull distribution and high censored proportion (c=.20; n=300) Cured proportions for Groups 1-4:(.3,.242,.192,.15) Censored proportions for Groups 1-4:(.3,.242,.192,.15)+c,c=.20 Sample size for Groups 1-4: (60,105,65,70),n=300 Parameters:β 0 =.5515,β 1 =.2958,α=1.5,γ=1.5(parameters of the distribution) ˆφ i,(ˆα, ˆγ) N (ˆβ 0, ˆβ 1 ) ˆp 0 l i=0 i=1 i=2 i=3 i=4 Mean 1 (.460,.343) (.314,.244,.187,.144) (.460,.348) (.313,.242,.185,.141) (.460,.349) (.313,.242,.185,.141) (.462,.349) (.312,.242,.185,.141) (.457,.336) (.316,.247,.190,.148) (1.526,1.503) RMSE 1 (.491,.183) (.073,.045,.039,.044) 2 (.492,.185) (.072,.044,.039,.044) 3 (.490,.185) (.072,.044,.039,.044) 4 (.489,.185) (.072,.044,.039,.044) (.503,.193) (.071,.040,.036,.043).010 Std 1 (.484,.177) (.072,.045,.039,.044) (.485,.178) (.071,.044,.039,.043) (.483,.178) (.071,.044,.039,.043) (.483,.178) (.071,.044,.039,.043) (.496,.190) (.069,.040,.036,.043) (.099,.092)
45 Data and Estimation
46 Numerical results Table: Logistic-Poisson mixture model: Weibull distribution and low censored proportion (c=.10; n=300) Cured proportions for Groups 1-4:(.3,.242,.192,.15) Censored proportions for Groups 1-4:(.3,.242,.192,.15)+c,c=.10 Sample size for Groups 1-4: (60,105,65,70),n=300 Parameters:β 0 =.5515,β 1 =.2958,α=1.5,γ=1.5(parameters of the distribution) ˆφ i,(ˆα, ˆγ) N (ˆβ 0, ˆβ 1 ) ˆp 0 l i=0 i=1 i=2 i=3 i=4 Mean 1 (.593,.307) (.292,.231,.182,.143) (.584,.303) (.294,.235,.185,.146) (.585,.303) (.294,.235,.185,.146) (.585,.303) (.294,.235,.185,.146) (.585,.303) (.294,.235,.185,.146) (1.477,1.543) RMSE 1 (.371,.143) (.051,.032,.029,.035) 2 (.365,.141) (.051,.031,.028,.035) 3 (.365,.140) (.051,.031,.028,.035) 4 (.365,.140) (.051,.031,.028,.035) (.365,.140) (.051,.031,.028,.035) (.005,.01) Std 1 (.369,.143) (.051,.030,.028,.035) (.365,.141) (.051,.030,.028,.034) (.364,.141) (.051,.030,.028,.034) (.364,.140) (.051,.030,.028,.034) (.365,.140) (.051,.030,.028,.035) (.067,.093)
47 Data and Estimation
48 Numerical results One more quantity of interest is the probability P(cured T > t)= p 0 (x β) p 0 (x β)+(1 p 0 (x β))s U (t;θ). Hence, in the next Table we mention the mean value (and sample standard deviation, in parentheses) of the above probability at the point t 0.95 for which P(cured T > t 0.95 )=0.95 (at each of the four possible values of our covariate; i.e. for x=1,2,3,4).
49 Numerical results Table: The sample mean of P(cured T > t 0.95 ) (n=300,c=.10) and for each possible value of covariate x. Weibull (α = 1,γ = 1.5)/Binary Weibull (α = 1.5,γ = 1.5)/Binary PLA PLA x N = 1 N = 2 N = 3 N = 4 Par* N = 1 N = 2 N = 3 N = 4 Par* 1.946(.02).945(.02).945(.02).946(.02).947(.02).971(.02).968(.01).967(.02).968(.02).947(.02) 2.947(.02).945(.02).945(.02).946(.02).947(.02).971(.03).967(.02).966(.02).967(.02).946(.02) 3.946(.02).944(.02).944(.02).945(.02).946(.02).970(.03).964(.02).964(.02).964(.02).945(.03) 4.944(.03).942(.03).942(.03).943(.03).944(.02).968(.04).961(.03).961(.03).961(.03).943(.03) Weibull (α = 1,γ = 1.5)/Poisson Weibull (α = 1.5,γ = 1.5)/Poisson 1.935(.04).933(.04).931(.04).931(.04).941(.03).958(.02).931(.04).932(.04).932(.04).931(.03) 2.933(.04).930(.04).928(.04).928(.04).940(.03).960(.02).928(.04).930(.04).931(.04).929(.03) 3.930(.04).926(.05).925(.05).925(.05).939(.03).962(.03).925(.05).927(.05).928(.05).926(.04) 4.926(.05).921(.05).920(.05).920(.06).936(.04).963(.03).921(.06).923(.06).924(.06).921(.05) *The respective estimations using the parametric approach.
50 Melanoma data The study took place from 1991 to 1995, and follow-up was conducted until 1998 (n=427); the data, taken from Ibrahim et al. (2005), present the survival times (in years) until the patient s death or the censoring time. In our application, the nodule category is the only covariate with 4 possible values x = 1,2,3,4 (lymph nodes involved in the disease) and the group sizes are n 1 = 111,n 2 = 137,n 3 = 87 and n 4 = 82, with respective censoring proportions 67.6%,61.3%,52.9% and 32.9%. The results from the PLA will be compared with those of the parametric approach.
51 Melanoma data Table: The estimation of the logistic-binary mixture model. N = 1 N = 2 N = 3 N = 4 N = 5 N = 6 Par 1 Par 2 p (.046).615(.046).617(.046).615(.046).615(.046).615(.046).632(.048).667(.040) p (.042).503(.042).504(.042).502(.042).503(.042).502(.042).506(.041).557(.030) p (.053).422(.053).423(.053).421(.052).421(.052).421(.052).379(.047).440(.033) p (.049).274(.049).275(.049).273(.049).273(.049).273(.049).267(.056).330(.046) φ φ φ φ φ φ φ l AIC BIC *Par 1 : the parametric Log-normal/logistic-binary mixture model of Balakrishnan and Pal (2013a) Par 2 : the parametric Weibull/logistic-binary mixture model of Balakrishnan and Pal (2013b)
52 Melanoma data P( M = 0 T > t) x=1 x=2 x=3 0.4 x=4 t Figure: The probability P(M = 0 T > t) = P(cured T > t), for each nodule category.
53 Similar works Larson and Dinse (1985) used a piecewise constant approximation for the baseline hazard function of the Cox s PH model (under a different competing cause scenario). In Lo et al. (1993), the baseline hazard was determined now by a piecewise continuous linear function (for the Bernoulli model and a different estimation approach). Chen and Ibrahim (2001) studied the Poisson mixture model where the common cdf of competing causes was approximated by a piecewise constant function (using different estimation approach). A piecewise constant approximation under Bayesian framework is met in Ibrahim, Chen and Sinha (2001) and Kim et al. (2007).
54 Conclusions The accuracy of PLA for estimating the regression coefficients is very close to that gained by a parametric approach; a similar remark made by Taylor (1995), through a KM-based non-parametric approach. The last line of our PLA seems to be significantly affected by the underlying censoring mechanism. In Lo et al. (1993) was found that the choice of cut points had a minimal effect on the results which is in accordance with our findings.
55 Conclusions The suggested non-parametric approach is quite flexible and the larger the N the more non-parametric the model is. We suggest to choose small to moderate values of N and trying to check the robustness of the estimates of the regression coefficients. No more lines can offer better approximations if the additional estimates of φ s lie very close to existing lines. Finally, we noted that more than 4 or 5 lines do not offer any improvement using either the maximum likelihood or MSE criterion.
56 Future work Future work includes: assuming that the number of competing causes follows a wider class of discrete distributions; under a fully non-parametric framework about the number of competing causes; a Cox proportional hazard model approach using different estimation methods; the asymptotic properties of these models.
57 References Balakrishnan, N. and Pal, S. (2012). EM algorithm-based likelihood estimation for some cure rate models. Journal of Statistical Theory and Practice, 6, Balakrishnan, N. and Pal, S. (2013a). Lognormal distribution and likelihood-based inference for flexible cure rate models based on COM-Poisson family. Computational Statistics & Data Analysis, 67, Balakrishnan, N. and Pal, S. (2013b). Expectation maximization-based likelihood inference for flexible cure rate models with Weibull distribution. Statistical Methods in Medical Research, available online. Balakrishnan, N. and Pal, S. (2013c). COM-Poisson cure rate models and associated likelihood-based inference with exponential and Weibull lifetimes, in Applied Reliability Engineering and Risk Analysis: Probabilistic Models and Statistical Inference (eds I. B. Frenkel, A. Karagrigoriou, A. Lisnianski and A. Kleyner), John Wiley & Sons, Chichester. Berkson, J. and Gage, R.P. (1952), Survival curves for cancer patients following treatment. Journal of the American Statistical Association, 47, Bentzen, S.M., Thames, H.D., Travis, E.L., Kian Ang, K., Van Der Schueren, E., Dewit, L. and Dixon, D.O. (1989). Direct estimation of latent time for radiation injury in late-responding normal tissues: gut, lung, and spinal cord. International journal of radiation biology, 55, Boag, J.M. (1949), Maximum likelihood estimates of the proportion of patients cured by cancer therapy. Journal of the Royal Statistical Society, Ser. B, 11, Cantor, A.B. and Shuster, J.J. (1992). Parametric versus non-parametric methods for estimating cure rates based on censored survival data. Statistics in Medicine, 11, Castro, M.D., Cancho, V.G. and Rodrigues, J. (2009). A bayesian long-term survival model parametrized in the cured fraction. Biometrical Journal, 51, Chen, M.H. and Ibrahim, J.G. (2001). Maximum likelihood methods for cure rate models with missing covariates. Biometrics, 57, Chen, M.H., Ibrahim, J.G. and Sinha, D. (1999). A new Bayesian model for survival data with a surviving fraction. Journal of the American Statistical Association, 94, Fang, H.B., Li, G. and Sun, J. (2005). Maximum Likelihood Estimation in a Semiparametric Logistic/Proportional-Hazards Mixture Model. Scandinavian Journal of Statistics, 32,
58 References Farewell, V.T. (1977a). A model for a binary variable with time-censored observations. Biometrika, 64, Farewell, V.T. (1977b). The combined effect of breast cancer risk factors. Cancer, 40, Farewell, V.T. (1982). The use of mixture models for the analysis of survival data with long-term survivors. Biometrics, 38, Farewell, V.T. (1986). Mixture models in survival analysis: Are they worth the risk? Canadian Journal of Statistics, 14, Ghitany, M. E. (1993). On the information matrix of exponential mixture models with long-term survivors. Biometrical Journal, 35, Ghitany, M.E. and Maller, R.A. (1992). Asymptotic results for exponential mixture models with long-term survivors. Statistics: A Journal of Theoretical and Applied Statistics, 23, Hashimoto, E.M., Ortega, E.M., Cordeiro, G.M. and Cancho, V.G. (2013). The Poisson Birnbaum-Saunders model with long-term survivors. Statistics, available online. Ibrahim, J.G., Chen, M.H. and Sinha, D. (2001). Bayesian semiparametric models for survival data with a cure fraction. Biometrics, 57, Ibrahim, J.G., Chen, M.H. and Sinha, D. (2005). Bayesian Survival Analysis. John Wiley & Sons, New York. Kim, S., Chen, M.H., Dey, D.K. and Gamerman, D. (2007). Bayesian dynamic models for survival data with a cure fraction. Lifetime Data Analysis, 13, Kuk, A.Y. and Chen, C.H. (1992). A mixture model combining logistic regression with proportional hazards regression. Biometrika, 79, Larson, M.G. and Dinse, G.E. (1985). A mixture model for the regression analysis of competing risks data. Applied statistics, 34,
59 References Laska, E.M. and Meisner, M.J. (1992). Nonparametric estimation and testing in a cure model. Biometrics, 48, Lo, Y.C., Taylor, J.M., McBride, W.H. and Withers, H.R. (1993). The effect of fractionated doses of radiation on mouse spinal cord. International Journal of Radiation Oncology, Biology, Physics, 27, Louzada, F., de Castro, M., Tomazella, V. and Gonzales, J.F. (2014). Modeling categorical covariates for lifetime data in the presence of cure fraction by Bayesian partition structures. Journal of Applied Statistics, 41, Maller, R.A. and Zhou, S. (1992). Estimating the proportion of immunes in a censored sample. Biometrika, 79, Maller, R.A. and Zhou, X. (1996). Survival Analysis with Long-term Survivors. John Wiley & Sons, New York. Mould, R.F. (1973). Statistical Models for Studying Long Term Survival Results following Treatment for Carcinoma of the Cervix. Ph.D. Thesis, University of London. Ortega, E.M., Barriga, G.D., Hashimoto, E.M., Cancho, V.G. and Cordeiro, G.M. (2014). A new class of survival regression models with cure fraction. Journal of Data Science, 12, Peng, Y. and Dear, K.B. (2000). A nonparametric mixture model for cure rate estimation. Biometrics, 56, Peng, Y., Dear, K.B. and Denham, J.W. (1998). A generalized F mixture model for cure rate estimation. Statistics in Medicine, 17, Rodrigues, J., Cancho, V.G., de Castro, M. and Louzada-Neto, F. (2009). On the unification of long-term survival models. Statistics & Probability Letters, 79,
60 References Rodrigues, J., de Castro, M., Balakrishnan, N. and Cancho, V.G. (2011). Destructive weighted Poisson cure rate models. Lifetime Data Analysis, 17, Struthers, C.A. and Farewell, V.T. (1989). A mixture model for time to AIDS data with left truncation and an uncertain origin. Biometrika, 76, Sy, J.P. and Taylor, J.M. (2000). Estimation in a Cox proportional hazards cure model. Biometrics, 56, Sy, J.P. and Taylor, J.M. (2001). Standard errors for the Cox proportional hazards cure model. Mathematical and Computer Modelling, 33, Taylor, J.M. (1995). Semi-parametric estimation in failure time mixture models. Biometrics, 51, Tsodikov, A.D., Ibrahim, J.G. and Yakovlev, A.Y. (2003). Estimating cure rates from survival data. Journal of the American Statistical Association, 98, Yakovlev, A.Y., Asselain, B., Bardou, V.J., Fourquet, A., Hoang, T., Rochefediere, A. and Tsodikov, A D. (1993). A simple stochastic model of tumor recurrence and its applications to data on premenopausal breast cancer. Biometrie et Analyse de Donnees Spatio-temporelles, 12, Yakovlev, A.Y., Cantor, A.B. and Shuster, J.J. (1994). Parametric versus non-parametric methods for estimating cure rates based on censored survival data. Statistics in Medicine, 13, Zhao, L., Feng, D., Bellile, E.L. and Taylor, J.M. (2014). Bayesian random threshold estimation in a Cox proportional hazards cure model. Statistics in Medicine, 33,
Parametric Maximum Likelihood Estimation of Cure Fraction Using Interval-Censored Data
Columbia International Publishing Journal of Advanced Computing (2013) 1: 43-58 doi:107726/jac20131004 Research Article Parametric Maximum Likelihood Estimation of Cure Fraction Using Interval-Censored
More informationAnalysis of Cure Rate Survival Data Under Proportional Odds Model
Analysis of Cure Rate Survival Data Under Proportional Odds Model Yu Gu 1,, Debajyoti Sinha 1, and Sudipto Banerjee 2, 1 Department of Statistics, Florida State University, Tallahassee, Florida 32310 5608,
More informationNegative Binomial Kumaraswamy-G Cure Rate Regression Model
Journal of Risk and Financial Management Article Negative Binomial Kumaraswamy-G Cure Rate Regression Model Amanda D Andrea 1,2, * ID, Ricardo Rocha 3, Vera Tomazella 1 ID and Francisco Louzada 2 1 Department
More informationEdinburgh Research Explorer
Edinburgh Research Explorer The zero-inflated promotion cure rate model applied to financial data on time-to-default Citation for published version: Oliveira Jr, M, Moreira, F & Louzada, F 2017, 'The zero-inflated
More informationAnalysis of competing risks data and simulation of data following predened subdistribution hazards
Analysis of competing risks data and simulation of data following predened subdistribution hazards Bernhard Haller Institut für Medizinische Statistik und Epidemiologie Technische Universität München 27.05.2013
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 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 informationStatistical Modeling and Analysis for Survival Data with a Cure Fraction
Statistical Modeling and Analysis for Survival Data with a Cure Fraction by Jianfeng Xu A thesis submitted to the Department of Mathematics and Statistics in conformity with the requirements for the degree
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 informationSemiparametric Models for Joint Analysis of Longitudinal Data and Counting Processes
Semiparametric Models for Joint Analysis of Longitudinal Data and Counting Processes by Se Hee Kim A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial
More information3003 Cure. F. P. Treasure
3003 Cure F. P. reasure November 8, 2000 Peter reasure / November 8, 2000/ Cure / 3003 1 Cure A Simple Cure Model he Concept of Cure A cure model is a survival model where a fraction of the population
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 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 informationBIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY
BIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY Ingo Langner 1, Ralf Bender 2, Rebecca Lenz-Tönjes 1, Helmut Küchenhoff 2, Maria Blettner 2 1
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 informationPromotion Time Cure Rate Model with Random Effects: an Application to a Multicentre Clinical Trial of Carcinoma
www.srl-journal.org Statistics Research Letters (SRL) Volume Issue, May 013 Promotion Time Cure Rate Model with Random Effects: an Application to a Multicentre Clinical Trial of Carcinoma Diego I. Gallardo
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 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 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 informationarxiv: v2 [stat.me] 19 Sep 2015
The zero-inflated cure rate regression model: Applications to fraud detection in bank loan portfolios Francisco Louzada 1 Institute of Mathematical Science and Computing at the University of São Paulo
More informationAssessing the effect of a partly unobserved, exogenous, binary time-dependent covariate on -APPENDIX-
Assessing the effect of a partly unobserved, exogenous, binary time-dependent covariate on survival probabilities using generalised pseudo-values Ulrike Pötschger,2, Harald Heinzl 2, Maria Grazia Valsecchi
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 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 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 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 informationGoodness-of-fit tests for the cure rate in a mixture cure model
Biometrika (217), 13, 1, pp. 1 17 Printed in Great Britain Advance Access publication on 31 July 216 Goodness-of-fit tests for the cure rate in a mixture cure model BY U.U. MÜLLER Department of Statistics,
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 informationBivariate Weibull Distributions Derived From Copula Functions In The Presence Of Cure Fraction And Censored Data
Journal of Data Science 14(2016), 295-316 Bivariate Weibull Distributions Derived From Copula Functions In The Presence Of Cure Fraction And Censored Data Em ılio A. Coelho-Barros 1*, Jorge Alberto Achcar
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 informationIntroduction to Reliability Theory (part 2)
Introduction to Reliability Theory (part 2) Frank Coolen UTOPIAE Training School II, Durham University 3 July 2018 (UTOPIAE) Introduction to Reliability Theory 1 / 21 Outline Statistical issues Software
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 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 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 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 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 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 informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationST5212: Survival Analysis
ST51: Survival Analysis 8/9: Semester II Tutorial 1. A model for lifetimes, with a bathtub-shaped hazard rate, is the exponential power distribution with survival fumction S(x) =exp{1 exp[(λx) α ]}. (a)
More informationCURE MODEL WITH CURRENT STATUS DATA
Statistica Sinica 19 (2009), 233-249 CURE MODEL WITH CURRENT STATUS DATA Shuangge Ma Yale University Abstract: Current status data arise when only random censoring time and event status at censoring are
More informationVertical modeling: analysis of competing risks data with a cure proportion
Vertical modeling: analysis of competing risks data with a cure proportion M A Nicolaie a and JMG Taylor b and C Legrand a a Institute of Statistics, Biostatistics and Actuarial Sciences, Catholic University
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 POISSON-EXPONENTIAL SURVIVAL DISTRIBUTION FOR LIFETIME DATA Vicente G. Cancho Francisco Louzada-Neto Gladys
More informationConstrained estimation for binary and survival data
Constrained estimation for binary and survival data Jeremy M. G. Taylor Yong Seok Park John D. Kalbfleisch Biostatistics, University of Michigan May, 2010 () Constrained estimation May, 2010 1 / 43 Outline
More informationSimple Step-Stress Models with a Cure Fraction
Simple Step-Stress Models with a Cure Fraction Nandini Kannan & Debasis Kundu Abstract In this article, we consider models for time-to-event data obtained from experiments in which stress levels are altered
More informationDuration Analysis. Joan Llull
Duration Analysis Joan Llull Panel Data and Duration Models Barcelona GSE joan.llull [at] movebarcelona [dot] eu Introduction Duration Analysis 2 Duration analysis Duration data: how long has an individual
More informationStatistical Methods for Alzheimer s Disease Studies
Statistical Methods for Alzheimer s Disease Studies Rebecca A. Betensky, Ph.D. Department of Biostatistics, Harvard T.H. Chan School of Public Health July 19, 2016 1/37 OUTLINE 1 Statistical collaborations
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 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 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 informationTHESIS for the degree of MASTER OF SCIENCE. Modelling and Data Analysis
PROPERTIES OF ESTIMATORS FOR RELATIVE RISKS FROM NESTED CASE-CONTROL STUDIES WITH MULTIPLE OUTCOMES (COMPETING RISKS) by NATHALIE C. STØER THESIS for the degree of MASTER OF SCIENCE Modelling and Data
More informationREGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520
REGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520 Department of Statistics North Carolina State University Presented by: Butch Tsiatis, Department of Statistics, NCSU
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 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 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 informationCensoring mechanisms
Censoring mechanisms Patrick Breheny September 3 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Fixed vs. random censoring In the previous lecture, we derived the contribution to the likelihood
More informationPower and Sample Size Calculations with the Additive Hazards Model
Journal of Data Science 10(2012), 143-155 Power and Sample Size Calculations with the Additive Hazards Model Ling Chen, Chengjie Xiong, J. Philip Miller and Feng Gao Washington University School of Medicine
More informationA new extended Birnbaum-Saunders model with cure fraction: classical and Bayesian approach
Communications for Statistical Applications and Methods 017, Vol. 4, No. 4, 397 419 https://doi.org/10.5351/csam.017.4.4.397 Print ISSN 87-7843 / Online ISSN 383-4757 A new extended Birnbaum-Saunders model
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 informationCox s proportional hazards model and Cox s partial likelihood
Cox s proportional hazards model and Cox s partial likelihood Rasmus Waagepetersen October 12, 2018 1 / 27 Non-parametric vs. parametric Suppose we want to estimate unknown function, e.g. survival function.
More informationThe nltm Package. July 24, 2006
The nltm Package July 24, 2006 Version 1.2 Date 2006-07-17 Title Non-linear Transformation Models Author Gilda Garibotti, Alexander Tsodikov Maintainer Gilda Garibotti Depends
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 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 informationAdvanced Methodology Developments in Mixture Cure Models
University of South Carolina Scholar Commons Theses and Dissertations 1-1-2013 Advanced Methodology Developments in Mixture Cure Models Chao Cai University of South Carolina Follow this and additional
More informationStat 542: Item Response Theory Modeling Using The Extended Rank Likelihood
Stat 542: Item Response Theory Modeling Using The Extended Rank Likelihood Jonathan Gruhl March 18, 2010 1 Introduction Researchers commonly apply item response theory (IRT) models to binary and ordinal
More informationLongitudinal + Reliability = Joint Modeling
Longitudinal + Reliability = Joint Modeling Carles Serrat Institute of Statistics and Mathematics Applied to Building CYTED-HAROSA International Workshop November 21-22, 2013 Barcelona Mainly from Rizopoulos,
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 informationTechnical Report - 7/87 AN APPLICATION OF COX REGRESSION MODEL TO THE ANALYSIS OF GROUPED PULMONARY TUBERCULOSIS SURVIVAL DATA
Technical Report - 7/87 AN APPLICATION OF COX REGRESSION MODEL TO THE ANALYSIS OF GROUPED PULMONARY TUBERCULOSIS SURVIVAL DATA P. VENKATESAN* K. VISWANATHAN + R. PRABHAKAR* * Tuberculosis Research Centre,
More informationSurvival Prediction Under Dependent Censoring: A Copula-based Approach
Survival Prediction Under Dependent Censoring: A Copula-based Approach Yi-Hau Chen Institute of Statistical Science, Academia Sinica 2013 AMMS, National Sun Yat-Sen University December 7 2013 Joint work
More informationLatent cure rate model under repair system and threshold effect
Universidade de São Paulo Biblioteca Digital da Produção Intelectual - BDPI Outros departamentos - ICMC/Outros Importação - 2015 2015 Latent cure rate model under repair system and threshold effect Journal
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationTime-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation
Time-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation H. Zhang, E. Cutright & T. Giras Center of Rail Safety-Critical Excellence, University of Virginia,
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 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 informationConstant Stress Partially Accelerated Life Test Design for Inverted Weibull Distribution with Type-I Censoring
Algorithms Research 013, (): 43-49 DOI: 10.593/j.algorithms.01300.0 Constant Stress Partially Accelerated Life Test Design for Mustafa Kamal *, Shazia Zarrin, Arif-Ul-Islam Department of Statistics & Operations
More informationPart III Measures of Classification Accuracy for the Prediction of Survival Times
Part III Measures of Classification Accuracy for the Prediction of Survival Times Patrick J Heagerty PhD Department of Biostatistics University of Washington 102 ISCB 2010 Session Three Outline Examples
More informationMixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals. John W. Mac McDonald & Alessandro Rosina
Mixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals John W. Mac McDonald & Alessandro Rosina Quantitative Methods in the Social Sciences Seminar -
More informationJoint Modeling of Longitudinal Item Response Data and Survival
Joint Modeling of Longitudinal Item Response Data and Survival Jean-Paul Fox University of Twente Department of Research Methodology, Measurement and Data Analysis Faculty of Behavioural Sciences Enschede,
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 informationSurvival analysis with long-term survivors and partially observed covariates
The Canadian Journal of Statistics Vol. 29, No. 2, 2001, Pages???-??? La revue canadienne de statistique Survival analysis with long-term survivors and partially observed covariates Meehyung CHO, Nathaniel
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 information[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements
[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements Aasthaa Bansal PhD Pharmaceutical Outcomes Research & Policy Program University of Washington 69 Biomarkers
More informationA SMOOTHED VERSION OF THE KAPLAN-MEIER ESTIMATOR. Agnieszka Rossa
A SMOOTHED VERSION OF THE KAPLAN-MEIER ESTIMATOR Agnieszka Rossa Dept. of Stat. Methods, University of Lódź, Poland Rewolucji 1905, 41, Lódź e-mail: agrossa@krysia.uni.lodz.pl and Ryszard Zieliński Inst.
More informationA Regression Model For Recurrent Events With Distribution Free Correlation Structure
A Regression Model For Recurrent Events With Distribution Free Correlation Structure J. Pénichoux(1), A. Latouche(2), T. Moreau(1) (1) INSERM U780 (2) Université de Versailles, EA2506 ISCB - 2009 - Prague
More informationMultivariate Survival Analysis
Multivariate Survival Analysis Previously we have assumed that either (X i, δ i ) or (X i, δ i, Z i ), i = 1,..., n, are i.i.d.. This may not always be the case. Multivariate survival data can arise in
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 informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
More informationGaussian Process Regression Model in Spatial Logistic Regression
Journal of Physics: Conference Series PAPER OPEN ACCESS Gaussian Process Regression Model in Spatial Logistic Regression To cite this article: A Sofro and A Oktaviarina 018 J. Phys.: Conf. Ser. 947 01005
More informationSAMPLE SIZE ESTIMATION FOR SURVIVAL OUTCOMES IN CLUSTER-RANDOMIZED STUDIES WITH SMALL CLUSTER SIZES BIOMETRICS (JUNE 2000)
SAMPLE SIZE ESTIMATION FOR SURVIVAL OUTCOMES IN CLUSTER-RANDOMIZED STUDIES WITH SMALL CLUSTER SIZES BIOMETRICS (JUNE 2000) AMITA K. MANATUNGA THE ROLLINS SCHOOL OF PUBLIC HEALTH OF EMORY UNIVERSITY SHANDE
More informationA general mixed model approach for spatio-temporal regression data
A general mixed model approach for spatio-temporal regression data Thomas Kneib, Ludwig Fahrmeir & Stefan Lang Department of Statistics, Ludwig-Maximilians-University Munich 1. Spatio-temporal regression
More informationAdaptive Prediction of Event Times in Clinical Trials
Adaptive Prediction of Event Times in Clinical Trials Yu Lan Southern Methodist University Advisor: Daniel F. Heitjan May 8, 2017 Yu Lan (SMU) May 8, 2017 1 / 19 Clinical Trial Prediction Event-based trials:
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 informationIn contrast, parametric techniques (fitting exponential or Weibull, for example) are more focussed, can handle general covariates, but require
Chapter 5 modelling Semi parametric We have considered parametric and nonparametric techniques for comparing survival distributions between different treatment groups. Nonparametric techniques, such as
More informationSimple techniques for comparing survival functions with interval-censored data
Simple techniques for comparing survival functions with interval-censored data Jinheum Kim, joint with Chung Mo Nam jinhkim@suwon.ac.kr Department of Applied Statistics University of Suwon Comparing survival
More informationLecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL
Lecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL The Cox PH model: λ(t Z) = λ 0 (t) exp(β Z). How do we estimate the survival probability, S z (t) = S(t Z) = P (T > t Z), for an individual with covariates
More informationPairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion
Pairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion Glenn Heller and Jing Qin Department of Epidemiology and Biostatistics Memorial
More informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationSurvival Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation
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 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 informationMultistate models and recurrent event models
and recurrent event models Patrick Breheny December 6 Patrick Breheny University of Iowa Survival Data Analysis (BIOS:7210) 1 / 22 Introduction In this final lecture, we will briefly look at two other
More information