A note on the decomposition of number of life years lost according to causes of death
|
|
- Damon Jeffery Gilmore
- 5 years ago
- Views:
Transcription
1 A note on the decomposition of number of life years lost according to causes of death Per Kragh Andersen Research Report 12/2 Department of Biostatistics University of Copenhagen
2 A note on the decomposition of number of life years lost according to causes of death Per Kragh Andersen, Department of Biostatistics, University of Copenhagen, Ø. Farimagsgade 5 PB299, DK-114 Copenhagen K, Denmark. pka@biostat.ku.dk Abstract The standard competing risks model is studied and we show that the cause j cumulative incidence function integrated from to τ has a natural interpretation as the expected number of life years lost due to cause j before time τ. This is analogous to the τ-restricted mean life time which is the survival function integrated from to τ. The large sample properties of a non-parametric estimator are outlined, and the method is exemplified using a standard data set on survival with malignant melanoma. It is discussed how the number of years lost may be related to subject-specific explanatory variables in a regression model based on pseudo-observartions. The method is contrasted to cause-specific measures of life years lost used in demography. Key words: demography; life expectancy; Survival analysis. 1
3 Introduction In demography the expected life time is a useful summary of the age-specific mortality experience (e.g., Preston et al., 21). If α(t) is the age-specific mortality rate (hazard function) then t S(t) = exp( α(u)du) is the corresponding survival function, and the life expectancy at birth is e(, ) = S(t)dt. In practice, is represented by the maximally attainable life length, say ω, and age-specific mortality rates between ages and ω are needed to estimate the expected life time. When death can be caused by a number, k, of exclusive and exhaustive causes j = 1,...,k, the survival function is k t S(t) = exp( α j (u)du) j=1 where α j ( ) is the cause-j-specific hazard function. Thereby, the life expectancy at birth can be written as e(, ) = S 1 (t) S k (t)dt with S j (t) = exp( t α j(u)du). So-called cause-j-deleted life-tables are 2
4 then defined via the survival function S(t)/S j (t) = exp( l j t α l (u)du) and under independence of causes of death, this is the survival function in a population where cause j is no longer operating and where (by independence ) the rate of death from cause l j is still given by the cause-specific hazard α l (t). Following those lines, life expectancy without cause j e j (, ) = S(t)/S j (t)dt and number of years lost due to cause j e j (, ) e(, ) can be defined (e.g., Beltrán-Sánchez et al., 28). In biostatistics, the assumption of independent competing risks has to a large extent been abandoned by now because of the recognition of the fact that data from this world (where all k causes of death are operating) provide little information on how the mortality experience would be in a hypothetical world where cause j has been eliminated (e.g., Andersen and Keiding, 212). It is therefore of interest to develop alternative measures of the impact on life expectancy of given causes of death which do not rely on such extrapolations. This is what we will do in this note and, as we shall see, a very simple approach to the question can be taken. 3
5 Restricted mean life time and number of life years lost In biostatistics the expected life time can rarely be estimated reliably because of right-censoring and as an alternative summary of the mortality experience, the τ-restricted mean life time is some times used (e.g., Irwin, 1949; Karrison, 1987; Andersen et al., 24). This is defined as e(,τ) = S(t)dt, (1) thus for τ = ω =, this is simply the expected life time at birth. The restricted mean life time is the expectation E(T τ) of the random variable T τ where T is the (unrestricted) life time. Therefore, the difference L(,τ) = τ E(T τ) (2) is the expected number of years lost before time τ. In a competing risks situation, the cause j cumulative incidence is the probability F j (t) = P(T t,d = j) = t S(u )α j (u)du of failing from cause j before time t (where the random variable D denotes the cause of failure with possible values j = 1,...,k). For all values of t we have k S(t) + F j (t) = 1 j=1 4
6 and, therefore, the expected number of years lost before time τ can be written as a sum L(,τ) = τ k S(t)dt = F j (t)dt j=1 with a term for each of the failure causes. We will now show that each term L j (,τ) = F j (t)dt has a quite simple interpretation. The competing risks model can conveniently be represented by a multistate model for a stochastic process, X(t), with states : alive (the initial state) and j: dead by cause j, j = 1,...,k (the final, absorbing states) (e.g., Andersen et al., 1993, Section IV.4). The j transition intensity is the cause-specific hazard α j (t) and the state occupation probabilities are P(X(t) = ) = S(t) and P(X(t) = j) = F j (t),j = 1,...,k. The failure time is the (proper) random variable, T = inf t {X(t) }, denoting the time of exit from the initial state whereas the times T (j) = inf t {X(t) = j} of entry into the final state j = 1,...,k are all improper random variables because P(T (j) = ) = 1 F j ( ) >. However, for any τ, the random variable T (j) τ is proper with an expected value given by E(T (j) τ) = ts(t)α j (t)dt + τ(1 F j (τ)) 5
7 which by partial integration equals τ F j(t)dt. This leads to the result L j (,τ) = F j (t)dt = E(τ T (j) τ). (3) Thereby, the total number of years lost before time τ can be decomposed into a sum over failure causes j of terms which can be interpreted as the number of years lost due to that cause. Restricting attention to subjects still alive at, say t, is simply handled by considering the conditional survival function S(t t ) = S(t)/S(t ) and the conditional cumulative incidences F j (t t ) = F j (t)/s(t ) given survival till t. Thereby, for those still alive at time t, the expected life time before τ > t e(t,τ) = t S(t t )dt and the expected number of years lost before τ due to cause j L j (t,τ) = t F j (t t )dt can be defined. Graphically, e(,τ) and L j (,τ),j = 1,...,k can be represented, as follows. Suppose for simplicity that k = 2 causes are considered. The cumulative incidences may be presented on a stacked plot, i.e., F 1 (t) and F 1 (t) + F 2 (t) are plotted against t and on such a plot, L 1 (,τ) is the area under the lower curve between and τ, and L 2 (,τ) the area between the two curves between and τ. Finally, e(,τ) is the area between the upper 6
8 curve and 1 between and τ (see Figure 1 for an example). Non-parametric estimation of the survival and cumulative incidence functions based on right-censored (and possibly, left-truncated) survival data is straightforward using the Kaplan-Meier and Aalen-Johansen estimators. Thereby, also the quantities e(t,τ) and L j (t,τ) are easily estimable. For the former, asymptotic properties as the number, n of subjects increases were presented by Andersen et al. (1993, Example IV.3.8). For the latter, asymptotic properties are outlined in the Appendix. In the case of no censoring, i.e. all failures before time τ are observed, both the Kaplan-Meier and the Aalen-Johansen estimators are simple relative frequencies, thus in this case the estimators of (1) and (3) are simply ê(,τ) = 1 (T i τ) n i and L j (,τ) = 1 I(D i = j)(τ T i τ). (4) n i Here, T i is the time of failure for subject i = 1,...,n and D i is the cause of failure indicator for subject i = 1,...,n with possible values j = 1,...,k. If interest focuses on how the parameters e(,τ) and L j (,τ) relate to individual explanatory variables (such as gender, ethnic groups, occupational groups etc.) then regression analysis is straightforward in the case of no censoring. In the more realistic case where some observations are rightcensored (before time τ), one possible approach would be to use so-called 7
9 pseudo-observations (e.g., Andersen et al., 23, 25). If θ denotes one of these parameters and θ the estimator based on the full sample then the ith pseudo-observation for θ is θ i = n θ (n 1) θ ( i),i = 1,...,n. (5) Here, θ ( i) is the estimator applied to the sample (of size n 1) obtained by taking subject i out of the full sample. In this way the association between the parameter and explanatory variables, Z i can be ascertained via generalized estimating equations with θ i as the outcome variable (Andersen et al., 23). The method is implemented as the function pseudoyl in the R-package pseudo available at CRAN. Example: Survival with malignant melanoma Andersen et al. (1993, Example I.3.1) presented a study of survival with skin cancer (malignant melanoma). Here, 25 patients had radical surgery at Odense University Hospital, Denmark between 1962 and At the closing date of the study, 57 patients had died from the disease (cause j = 1), 14 had died from other causes (cause j = 2) and 134 were still alive. Figure 1 shows the stacked plot of the cumulative incidences of dying from the two causes in question. For τ = 1 years the restricted mean life time is ê(, 1) = 7.51 (SD=.245) while L 1 (, 1) = 2.5 (SD=.232) years were lost due to deaths from the disease and L 2 (, 1) =.485 (SD=.133) years were lost due to deaths from other causes. The corresponding areas on the stacked plot can 8
10 be seen in Figure 1. Figure 1 around here. In regression models using pseudo-observations it was found that men on average lost 1.21 (SD=.488) more years than women because of the disease and.211 (SD=.284) more years than women due to other causes. The estimated intercepts in these models (estimated numbers of years lost for women) were for deaths from the disease and.44 for deaths from other causes. Since 79/25=38.5% of the patients were men, these results are well in line with the overall numbers of years lost (since and ). The effect of tumor thickness on the number of years lost due to the disease was.361 years (SD=.82) for each mm and the effect of age on the number of years lost due to other causes was.269 years (SD=.82) for every 1 years of age. Discussion. Accounting for standard mortality In demography, number of life years lost is usually related to life expectancy based on standard mortality tables (e.g., Vaupel and Canudas-Romo, 23, Aragón et al., 28, Shkolnikov et al., 211). Thus, for uncensored data the expected number of life years lost is defined as n L(,ω) = D i e (T i,ω)/ i=1 i D i where D i = n and e (T i,ω) is the mean residual life time at the time (age) T i of death, calculated based on standard mortality rates. For cause 9
11 j = 1,...,k the number of life years lost due to that cause is similarly defined as n L j (,ω) = I(D i = j)e (T i,ω)/ i=1 i I(D i = j). Thus, the mean residual life times at time (age) of death from cause j are averaged over those who died from that cause. In the parlance of Gardner and Sanborn (199) this is a measure of Years of potential life lost. Replacing the mean residual life time in these expressions by the difference max(τ T i, ), a measure of Premature (to τ) years of potential life lost is obtained. This is seen to resemble (4) (for uncensored data), however, we divide by n when the cause-specific measure is calculated while, in demography, division by the number i I(D i = j) of deaths from cause j (before τ) is done. As we have seen, our measure of the expected number of life years lost due to cause j before τ has an interpretation as the parameter E(τ T (j) τ) from the population where all k causes of death are operating and it does not involve standard mortality tables. In biostatistics a framework where it is feasible to account for standard mortality is that of relative survival or excess hazard models. Here, the mortality rate, α(t) for a subject is written as a sum α(t) = α (t) + α (t) of a standard mortality rate, α (t) (assumed known) and an excess mortality, α (t). Often, this model is used in the absence of cause of death information, however, formally this is the hazard rate in a competing risks model with two causes of failure in which case α (t) is interpreted as the mortality rate 1
12 from the disease while α (t) is the mortality rate from other causes. Further, if cause of death information is available and if cause-specific standard mortality rates, α j(t),j = 1,...,k are given with α (t) = j α j(t) then an excess mortality competing risks model for the cause-specific hazard α j (t) = α j (j) + α j(t) can be studied. This can be seen as a competing risks model with k + 1 cause-specific hazards: α (t) for dying from disease-unrelated causes and α j (t) = α j (t) α (t) for dying from cause j,j = 1,...,k due to the disease. For this competing risks model we have the relation S(t) + j F j (t) + F (t) = 1 (with obvious notation) and the expected total number of life years lost before time τ L(,τ) = τ S(t)dt can be decomposed into L(,τ) = j F j (t)dt + F (t)dt. Here, the interpretation is that F (t)dt is the unavoidable expected number of life years lost due to the standard mortality while F j(t)dt represents 11
13 the number of years lost from cause j due to the disease. Note that, in L(,τ) the quantity L (,τ) = F (t)dt = t S(u)α (u)dudt appeared. This represents the standard number of life years lost when the excess mortality rates α j (t),j = 1,...,k are also operating and replaces the quantity t S (u)α (u)dudt = τ S (t)dt used in demography. The latter would be the number of life years lost before τ in a population where α (t) is the only hazard rate operating, a quantity which does not play a role here since all our computations relate to this world where all causes of death are present. Acknowledgements The research was supported by the Danish Natural Science Research Council [grant number Point process modelling and statistical inference ]. Comments from Niels Keiding, Department of Biostatistics, University of Copenhagen are greatly appreciated. We thank Maja Pohar Perme, Institute of Biostatistics and Medical Informatics, University of Ljubljana for updating the pseudo package. 12
14 References Andersen, PK, Borgan, Ø, Gill, RD, Keiding, N (1993) Statistical Models Based on Counting Processes. Springer-Verlag, New York. Andersen, PK, Hansen, MG, Klein, JP (24). Regression analysis of restricted mean life time using pseudo-observations. Lifetime Data Analysis 1, Andersen, PK, Keiding, N. (212). Interpretability and importance of functionals in competing risks and multi-state models. Statistics in Medicine (in press). Andersen, PK, Klein, JP, Rosthøj, S (23). Generalized linear models for correlated pseudo-observations, with applications to multi-state models. Biometrika 9, Aragón, TJ, Lichtensztajn, DY, Katcher, BS, Reiter, R, Katz, MH (28). Calculating expected years of life lost for assessing local ethnic disparities in causes of premature death. BMC Public Health 8, doi:1.1186/ Beltrán-Sánchez, H, Preston, SH, Canudas-Romo, V (28). An integrated approach to cause-of-death analysis: cause-deleted life tables and decompositions of life expectancy. Demographic Research 19,
15 Gardner, JW, Sanborn, JS (199). Years of potential life lost (YPLL) What does it measure? Epidemiology 1, Irwin, JO. (1949). The standard error of an estimate of expectation of life, with special reference to expectation of tumourless life in experiments with mice. Journal of Hygiene 47, Karrison, T. (1987). Restricted mean life with adjustment for covariates. Journal of the American Statistical Association 82, Preston, SH, Heuveline, P, Guillot, M (21). Demography. Measuring and Modeling Population Processes. Blackwell, Oxford. Shkolnikov, VM, Andreev, EM, Zhang, Z, Oeppen, J, Vaupel, JW. (211). Losses of expected lifetime in the United States and other developed countries: Methods and empirical analyses. Demography 48, Vaupel, JW, Canudas-Romo, V (23). Decomposing change in life expectancy: A bouquet of formulas in honor of Nathan Keifitz s 9th birthday. Demography 4,
16 Appendix: Asymptotics ( F j (t) F j (t))dt = = ( t ( t t ( (Ŝ(u) S(u))dÂj(u) + S(u)(dÂj(u) da j (u)) ) dt ) Ŝ(u) ( S(u) 1)S(u)dÂj(u) + S(u)(dÂj(u) da j (u)) dt u J(x) Y (x) dm(x)s(u)dâj(u) + S(u) J(u) ) Y (u) dm j(u) dt (Duhamel s equation and consistency of Ŝ and Âj) = ( t t (changing order of integration) u df j (x) J(u) ) J(u) dm(u) + S(u) Y (u) Y (u) dm j(u) dt = ( t (S(u) F j (t) + F j (u)) J(u) Y (u) dm j(u) (F j (t) F j (u)) J(u) Y (u) dm j(u) ) dt. In large samples, a process of the form n t f t(u)j(u)dm(u)/y (u) with a deterministic f t ( ) behaves like a mean zero Gaussian process U(t) with a covariance σ(s,t) which may be estimated by σ(s,t) = s t f s (u)f t (u) dn(u). Y 2 (u) Therefore, we have in large samples n ( F j (t) F j (t))dt U j (t)dt + U j (t)dt 15
17 where U j and U j are independent. It follows that L j (,τ) is asymptotically normal and that the variance may be estimated by ( σ j (s,t) + σ j (s,t))dsdt with f j t (u) = S(u) F j (t) + F j (u) and f j t (u) = F j (t) + F j (u). Thus, σ j (s,t) = s t (Ŝ(u) F j (s) + F j (u))(ŝ(u) F j (t) + F j (u)) dn Y (u) 2 j (u) and σ j (s,t) = s t ( F j (u) F j (s))( F j (u) F j (t)) d N Y (u) 2 r (u), r j where, N r ( ) counts failures from cause r. 16
18 Cumulative incidences YEARS Figure 1: Malignant melanoma study: Stacked plot of cumulative incidences of dying from the disease (dashed curve) and from other causes (difference between solid and dashed curves). 17
19 Research Reports available from Department of Biostatistics Department of Biostatistics University of Copenhagen Øster Farimagsgade 5 P.O. Box Copenhagen K Denmark 1/1 Andersen, P.K. & Skrondal, A. Biological interaction from a statistical point of view. 1/2 Rosthøj, S., Keiding, N., Schmiegelow, K. Application of History-Adjusted Marginal Structural Models to Maintenance Therapy of Children with Acute Lymphoblastic Leukaemia. 1/3 Parner, E.T. & Andersen, P.K. Regression analysis of censored data using pseudoobservations. 1/4 Nielsen, T. & Kreiner, S. Course Evaluation and Development: What can Learning Styles Contribute? 1/5 Lange, T. & Hansen, J.V. Direct and Indirect Effects in a Survival Context. 1/6 Andersen, P.K. & Keiding, N. Interpretability and importance of functionals in competing risks and multi-state models. 1/7 Gerds, T.A., Kattan, M.W., Schumacher, M. & Yu, C. Estimating a time-dependent concordance index for survival prediction models with covariate dependent censoring. 1/8 Mogensen, U.B., Ishwaran, H. & Gerds, T.A. Evaluating random forests for survival analysis using prediction error curves. 11/1 Kreiner, S. Is the foundation under PISA solid? A critical look at the scaling model underlying international comparisons of student attainment. 11/2 Andersen, P.K. Competing risks in epidemiology: Possibilities and pitfalls. 11/3 Holst, K.K. Model diagnostics based on cumulative residuals: The R-package gof. 11/4 Holst, K.K. & Budtz-Jørgensen, E. Linear latent variable models: The lava-package. 11/5 Holst, K.K., Budtz-Jørgensen, E. & Knudsen, G.M. A latent variable model with mixed binary and continuous response variables. 11/6 Cortese, G., Gerds, T.A., &Andersen, P.K. Comparison of prediction models for competing risks with time-dependent covariates.
20 11/7 Scheike, T.H. & Sun, Y. On Cross-Odds Ratio for Multivariate Competing Risks Data. 11/8 Gerds, T.A., Scheike T.H. & Andersen P.K. Absolute risk regression for competing risks: interpretation, link functions and prediction. 11/9 Scheike, T.H., Maiers, M.J., Rocha, V. & Zhang, M. Competing risks with missing covariates: Effect of haplotypematch on hematopoietic cell transplant patients. 12/1 Keiding, N., Hansen, O.K.H., Sørensen, D.N. & Slama, R. The current duration approach to estimating time to pregnancy. 12/2 Andersen, P.K. A note on the decomposition of number of life years lost according to causes of death
Multi-state models: prediction
Department of Medical Statistics and Bioinformatics Leiden University Medical Center Course on advanced survival analysis, Copenhagen Outline Prediction Theory Aalen-Johansen Computational aspects Applications
More informationMulti-state Models: An Overview
Multi-state Models: An Overview Andrew Titman Lancaster University 14 April 2016 Overview Introduction to multi-state modelling Examples of applications Continuously observed processes Intermittently observed
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 informationDAGStat Event History Analysis.
DAGStat 2016 Event History Analysis Robin.Henderson@ncl.ac.uk 1 / 75 Schedule 9.00 Introduction 10.30 Break 11.00 Regression Models, Frailty and Multivariate Survival 12.30 Lunch 13.30 Time-Variation and
More informationDefinitions and examples Simple estimation and testing Regression models Goodness of fit for the Cox model. Recap of Part 1. Per Kragh Andersen
Recap of Part 1 Per Kragh Andersen Section of Biostatistics, University of Copenhagen DSBS Course Survival Analysis in Clinical Trials January 2018 1 / 65 Overview Definitions and examples Simple estimation
More informationPhD course in Advanced survival analysis. One-sample tests. Properties. Idea: (ABGK, sect. V.1.1) Counting process N(t)
PhD course in Advanced survival analysis. (ABGK, sect. V.1.1) One-sample tests. Counting process N(t) Non-parametric hypothesis tests. Parametric models. Intensity process λ(t) = α(t)y (t) satisfying Aalen
More informationUnderstanding product integration. A talk about teaching survival analysis.
Understanding product integration. A talk about teaching survival analysis. Jan Beyersmann, Arthur Allignol, Martin Schumacher. Freiburg, Germany DFG Research Unit FOR 534 jan@fdm.uni-freiburg.de It is
More informationMultistate models in survival and event history analysis
Multistate models in survival and event history analysis Dorota M. Dabrowska UCLA November 8, 2011 Research supported by the grant R01 AI067943 from NIAID. The content is solely the responsibility of the
More informationPricing and Risk Analysis of a Long-Term Care Insurance Contract in a non-markov Multi-State Model
Pricing and Risk Analysis of a Long-Term Care Insurance Contract in a non-markov Multi-State Model Quentin Guibert Univ Lyon, Université Claude Bernard Lyon 1, ISFA, Laboratoire SAF EA2429, F-69366, Lyon,
More informationSurvival Analysis I (CHL5209H)
Survival Analysis Dalla Lana School of Public Health University of Toronto olli.saarela@utoronto.ca January 7, 2015 31-1 Literature Clayton D & Hills M (1993): Statistical Models in Epidemiology. Not really
More informationUNIVERSITY OF CALIFORNIA, 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 informationChapter 4 Fall Notations: t 1 < t 2 < < t D, D unique death times. d j = # deaths at t j = n. Y j = # at risk /alive at t j = n
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 4 Fall 2012 4.2 Estimators of the survival and cumulative hazard functions for RC data Suppose X is a continuous random failure time with
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 informationA GENERALIZED ADDITIVE REGRESSION MODEL FOR SURVIVAL TIMES 1. By Thomas H. Scheike University of Copenhagen
The Annals of Statistics 21, Vol. 29, No. 5, 1344 136 A GENERALIZED ADDITIVE REGRESSION MODEL FOR SURVIVAL TIMES 1 By Thomas H. Scheike University of Copenhagen We present a non-parametric survival model
More informationEmpirical Processes & Survival Analysis. The Functional Delta Method
STAT/BMI 741 University of Wisconsin-Madison Empirical Processes & Survival Analysis Lecture 3 The Functional Delta Method Lu Mao lmao@biostat.wisc.edu 3-1 Objectives By the end of this lecture, you will
More informationFaculty of Health Sciences. Regression models. Counts, Poisson regression, Lene Theil Skovgaard. Dept. of Biostatistics
Faculty of Health Sciences Regression models Counts, Poisson regression, 27-5-2013 Lene Theil Skovgaard Dept. of Biostatistics 1 / 36 Count outcome PKA & LTS, Sect. 7.2 Poisson regression The Binomial
More informationOn the Breslow estimator
Lifetime Data Anal (27) 13:471 48 DOI 1.17/s1985-7-948-y On the Breslow estimator D. Y. Lin Received: 5 April 27 / Accepted: 16 July 27 / Published online: 2 September 27 Springer Science+Business Media,
More informationRobust estimates of state occupancy and transition probabilities for Non-Markov multi-state models
Robust estimates of state occupancy and transition probabilities for Non-Markov multi-state models 26 March 2014 Overview Continuously observed data Three-state illness-death General robust estimator Interval
More informationConcordance for prognostic models with competing risks
Concordance for prognostic models with competing risks Marcel Wolbers Michael T. Koller Jacquline C. M. Witteman Thomas A. Gerds Research Report 13/3 Department of Biostatistics University of Copenhagen
More informationPh.D. course: Regression models. Introduction. 19 April 2012
Ph.D. course: Regression models Introduction PKA & LTS Sect. 1.1, 1.2, 1.4 19 April 2012 www.biostat.ku.dk/~pka/regrmodels12 Per Kragh Andersen 1 Regression models The distribution of one outcome variable
More informationLog-linearity for Cox s regression model. Thesis for the Degree Master of Science
Log-linearity for Cox s regression model Thesis for the Degree Master of Science Zaki Amini Master s Thesis, Spring 2015 i Abstract Cox s regression model is one of the most applied methods in medical
More informationPh.D. course: Regression models. Regression models. Explanatory variables. Example 1.1: Body mass index and vitamin D status
Ph.D. course: Regression models Introduction PKA & LTS Sect. 1.1, 1.2, 1.4 25 April 2013 www.biostat.ku.dk/~pka/regrmodels13 Per Kragh Andersen Regression models The distribution of one outcome variable
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 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 informationADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES. Cox s regression analysis Time dependent explanatory variables
ADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES Cox s regression analysis Time dependent explanatory variables Henrik Ravn Bandim Health Project, Statens Serum Institut 4 November 2011 1 / 53
More informationDepartment of Biostatistics University of Copenhagen
Cox regression with missing covariate data using a modified partial likelihood method Torben Martinussen Klaus K. Holst Thomas Scheike Research Report 14/7 Department of Biostatistics University of Copenhagen
More informationStatistics 262: Intermediate Biostatistics Non-parametric Survival Analysis
Statistics 262: Intermediate Biostatistics Non-parametric Survival Analysis Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Overview of today s class Kaplan-Meier Curve
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationA multi-state model for the prognosis of non-mild acute pancreatitis
A multi-state model for the prognosis of non-mild acute pancreatitis Lore Zumeta Olaskoaga 1, Felix Zubia Olaskoaga 2, Guadalupe Gómez Melis 1 1 Universitat Politècnica de Catalunya 2 Intensive Care Unit,
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 informationModeling Prediction of the Nosocomial Pneumonia with a Multistate model
Modeling Prediction of the Nosocomial Pneumonia with a Multistate model M.Nguile Makao 1 PHD student Director: J.F. Timsit 2 Co-Directors: B Liquet 3 & J.F. Coeurjolly 4 1 Team 11 Inserm U823-Joseph Fourier
More informationA multistate additive relative survival semi-markov model
46 e Journées de Statistique de la SFdS A multistate additive semi-markov Florence Gillaizeau,, & Etienne Dantan & Magali Giral, & Yohann Foucher, E-mail: florence.gillaizeau@univ-nantes.fr EA 4275 - SPHERE
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 informationMultistate Modelling Vertical Transmission and Determination of R 0 Using Transition Intensities
Applied Mathematical Sciences, Vol. 9, 2015, no. 79, 3941-3956 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.52130 Multistate Modelling Vertical Transmission and Determination of R 0
More informationProduct-limit estimators of the gap time distribution of a renewal process under different sampling patterns
Product-limit estimators of the gap time distribution of a renewal process under different sampling patterns arxiv:13.182v1 [stat.ap] 28 Feb 21 Richard D. Gill Department of Mathematics University of Leiden
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 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 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 informationPart IV Extensions: Competing Risks Endpoints and Non-Parametric AUC(t) Estimation
Part IV Extensions: Competing Risks Endpoints and Non-Parametric AUC(t) Estimation Patrick J. Heagerty PhD Department of Biostatistics University of Washington 166 ISCB 2010 Session Four Outline Examples
More informationChapter 7 Fall Chapter 7 Hypothesis testing Hypotheses of interest: (A) 1-sample
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 7 Fall 2012 Chapter 7 Hypothesis testing Hypotheses of interest: (A) 1-sample H 0 : S(t) = S 0 (t), where S 0 ( ) is known survival function,
More informationEstimating transition probabilities for the illness-death model The Aalen-Johansen estimator under violation of the Markov assumption Torunn Heggland
Estimating transition probabilities for the illness-death model The Aalen-Johansen estimator under violation of the Markov assumption Torunn Heggland Master s Thesis for the degree Modelling and Data Analysis
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 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 informationAccelerated Failure Time Models: A Review
International Journal of Performability Engineering, Vol. 10, No. 01, 2014, pp.23-29. RAMS Consultants Printed in India Accelerated Failure Time Models: A Review JEAN-FRANÇOIS DUPUY * IRMAR/INSA of Rennes,
More informationJournal of Statistical Software
JSS Journal of Statistical Software January 2011, Volume 38, Issue 2. http://www.jstatsoft.org/ Analyzing Competing Risk Data Using the R timereg Package Thomas H. Scheike University of Copenhagen Mei-Jie
More informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
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 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 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 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 informationRene Tabanera y Palacios 4. Danish Epidemiology Science Center. Novo Nordisk A/S Gentofte. September 1, 1995
Estimation of variance in Cox's regression model with gamma frailties. Per Kragh Andersen 2 John P. Klein 3 Kim M. Knudsen 2 Rene Tabanera y Palacios 4 Department of Biostatistics, University of Copenhagen,
More informationST745: Survival Analysis: Nonparametric methods
ST745: Survival Analysis: Nonparametric methods Eric B. Laber Department of Statistics, North Carolina State University February 5, 2015 The KM estimator is used ubiquitously in medical studies to estimate
More informationMultivariate Survival Data With Censoring.
1 Multivariate Survival Data With Censoring. Shulamith Gross and Catherine Huber-Carol Baruch College of the City University of New York, Dept of Statistics and CIS, Box 11-220, 1 Baruch way, 10010 NY.
More informationNonparametric two-sample tests of longitudinal data in the presence of a terminal event
Nonparametric two-sample tests of longitudinal data in the presence of a terminal event Jinheum Kim 1, Yang-Jin Kim, 2 & Chung Mo Nam 3 1 Department of Applied Statistics, University of Suwon, 2 Department
More informationEstimation for Modified Data
Definition. Estimation for Modified Data 1. Empirical distribution for complete individual data (section 11.) An observation X is truncated from below ( left truncated) at d if when it is at or below d
More informationPhilosophy and Features of the mstate package
Introduction Mathematical theory Practice Discussion Philosophy and Features of the mstate package Liesbeth de Wreede, Hein Putter Department of Medical Statistics and Bioinformatics Leiden University
More informationGoodness-Of-Fit for Cox s Regression Model. Extensions of Cox s Regression Model. Survival Analysis Fall 2004, Copenhagen
Outline Cox s proportional hazards model. Goodness-of-fit tools More flexible models R-package timereg Forthcoming book, Martinussen and Scheike. 2/38 University of Copenhagen http://www.biostat.ku.dk
More informationSurvival Analysis for Case-Cohort Studies
Survival Analysis for ase-ohort Studies Petr Klášterecký Dept. of Probability and Mathematical Statistics, Faculty of Mathematics and Physics, harles University, Prague, zech Republic e-mail: petr.klasterecky@matfyz.cz
More informationUSING MARTINGALE RESIDUALS TO ASSESS GOODNESS-OF-FIT FOR SAMPLED RISK SET DATA
USING MARTINGALE RESIDUALS TO ASSESS GOODNESS-OF-FIT FOR SAMPLED RISK SET DATA Ørnulf Borgan Bryan Langholz Abstract Standard use of Cox s regression model and other relative risk regression models for
More informationLecture 7 Time-dependent Covariates in Cox Regression
Lecture 7 Time-dependent Covariates in Cox Regression So far, we ve been considering the following Cox PH model: λ(t Z) = λ 0 (t) exp(β Z) = λ 0 (t) exp( β j Z j ) where β j is the parameter for the the
More informationRegression analysis of interval censored competing risk data using a pseudo-value approach
Communications for Statistical Applications and Methods 2016, Vol. 23, No. 6, 555 562 http://dx.doi.org/10.5351/csam.2016.23.6.555 Print ISSN 2287-7843 / Online ISSN 2383-4757 Regression analysis of interval
More informationMaximum likelihood estimation for Cox s regression model under nested case-control sampling
Biostatistics (2004), 5, 2,pp. 193 206 Printed in Great Britain Maximum likelihood estimation for Cox s regression model under nested case-control sampling THOMAS H. SCHEIKE Department of Biostatistics,
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 informationA Bayesian Nonparametric Approach to Causal Inference for Semi-competing risks
A Bayesian Nonparametric Approach to Causal Inference for Semi-competing risks Y. Xu, D. Scharfstein, P. Mueller, M. Daniels Johns Hopkins, Johns Hopkins, UT-Austin, UF JSM 2018, Vancouver 1 What are semi-competing
More informationfor Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park
Threshold Regression for Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park MLTLEE@UMD.EDU Outline The proportional hazards (PH) model is widely used in analyzing time-to-event data.
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
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 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 Regression Models
Survival Regression Models David M. Rocke May 18, 2017 David M. Rocke Survival Regression Models May 18, 2017 1 / 32 Background on the Proportional Hazards Model The exponential distribution has constant
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 informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
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 Time-to-Event Data: Chapter 2 - Nonparametric estimation of functions of survival time
Analysis of Time-to-Event Data: Chapter 2 - Nonparametric estimation of functions of survival time Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term
More informationNonparametric Model Construction
Nonparametric Model Construction Chapters 4 and 12 Stat 477 - Loss Models Chapters 4 and 12 (Stat 477) Nonparametric Model Construction Brian Hartman - BYU 1 / 28 Types of data Types of data For non-life
More informationStatistics in medicine
Statistics in medicine Lecture 4: and multivariable regression Fatma Shebl, MD, MS, MPH, PhD Assistant Professor Chronic Disease Epidemiology Department Yale School of Public Health Fatma.shebl@yale.edu
More informationPackage threg. August 10, 2015
Package threg August 10, 2015 Title Threshold Regression Version 1.0.3 Date 2015-08-10 Author Tao Xiao Maintainer Tao Xiao Depends R (>= 2.10), survival, Formula Fit a threshold regression
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 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 informationYou know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?
You know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?) I m not goin stop (What?) I m goin work harder (What?) Sir David
More informationTMA 4275 Lifetime Analysis June 2004 Solution
TMA 4275 Lifetime Analysis June 2004 Solution Problem 1 a) Observation of the outcome is censored, if the time of the outcome is not known exactly and only the last time when it was observed being intact,
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 informationDISCRETE PROBABILITY DISTRIBUTIONS
DISCRETE PROBABILITY DISTRIBUTIONS REVIEW OF KEY CONCEPTS SECTION 41 Random Variable A random variable X is a numerically valued quantity that takes on specific values with different probabilities The
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 informationA multi-state model for the prognosis of non-mild acute pancreatitis
A multi-state model for the prognosis of non-mild acute pancreatitis Lore Zumeta Olaskoaga 1, Felix Zubia Olaskoaga 2, Marta Bofill Roig 1, Guadalupe Gómez Melis 1 1 Universitat Politècnica de Catalunya
More informationSmoothing the Nelson-Aalen Estimtor Biostat 277 presentation Chi-hong Tseng
Smoothing the Nelson-Aalen Estimtor Biostat 277 presentation Chi-hong seng Reference: 1. Andersen, Borgan, Gill, and Keiding (1993). Statistical Model Based on Counting Processes, Springer-Verlag, p.229-255
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 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 information2011/04 LEUKAEMIA IN WALES Welsh Cancer Intelligence and Surveillance Unit
2011/04 LEUKAEMIA IN WALES 1994-2008 Welsh Cancer Intelligence and Surveillance Unit Table of Contents 1 Definitions and Statistical Methods... 2 2 Results 7 2.1 Leukaemia....... 7 2.2 Acute Lymphoblastic
More informationChapter 2 - Survival Models
2-1 Chapter 2 - Survival Models Section 2.2 - Future Lifetime Random Variable and the Survival Function Let T x = ( Future lifelength beyond age x of an individual who has survived to age x [measured in
More information1 The problem of survival analysis
1 The problem of survival analysis Survival analysis concerns analyzing the time to the occurrence of an event. For instance, we have a dataset in which the times are 1, 5, 9, 20, and 22. Perhaps those
More informationFrailty Modeling for clustered survival data: a simulation study
Frailty Modeling for clustered survival data: a simulation study IAA Oslo 2015 Souad ROMDHANE LaREMFiQ - IHEC University of Sousse (Tunisia) souad_romdhane@yahoo.fr Lotfi BELKACEM LaREMFiQ - IHEC University
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 informationDirect likelihood inference on the cause-specific cumulative incidence function: a flexible parametric regression modelling approach
Direct likelihood inference on the cause-specific cumulative incidence function: a flexible parametric regression modelling approach Sarwar I Mozumder 1, Mark J Rutherford 1, Paul C Lambert 1,2 1 Biostatistics
More informationSupport Vector Hazard Regression (SVHR) for Predicting Survival Outcomes. Donglin Zeng, Department of Biostatistics, University of North Carolina
Support Vector Hazard Regression (SVHR) for Predicting Survival Outcomes Introduction Method Theoretical Results Simulation Studies Application Conclusions Introduction Introduction For survival data,
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 informationProduct-limit estimators of the survival function with left or right censored data
Product-limit estimators of the survival function with left or right censored data 1 CREST-ENSAI Campus de Ker-Lann Rue Blaise Pascal - BP 37203 35172 Bruz cedex, France (e-mail: patilea@ensai.fr) 2 Institut
More informationOn graphical tests for proportionality of hazards in two samples
On graphical tests for proportionality of hazards in two samples Technical Report No. ASU/2014/5 Dated: 19 June 2014 Shyamsundar Sahoo, Haldia Government College, Haldia and Debasis Sengupta, Indian Statistical
More informationLecture 2: Martingale theory for univariate survival analysis
Lecture 2: Martingale theory for univariate survival analysis In this lecture T is assumed to be a continuous failure time. A core question in this lecture is how to develop asymptotic properties when
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More information