Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis
|
|
- Sharlene Atkinson
- 5 years ago
- Views:
Transcription
1 Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis Andrea Arfè 1 Stefano Peluso 2 Pietro Muliere 1 1 Università Commerciale Luigi Bocconi 2 Università Cattolica del Sacro Cuore SISBAYES 217 workshop, Rome, Italy Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
2 Introduction In clinical prognostic research with a time-to-event outcome, occurrence of competing risks may preclude the occurrence of another event of interest. Example: Amsterdam Cohort Study of AIDS progression in HIV-infected men [Geskus et al., 23] Syncytium Inducing (SI) AIDS At risk non-si AIDS Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
3 Subdistribution functions Subdistribution function [Kalbfleisch and Prentice, 22] F (t, d) = P (T t, C = c) where: T > is the time to event onset C {1,..., K} encodes the type of occurring event Can consider covariates: F (t, c x) Note: will assume that T {1, 2, 3,...} Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
4 Outline 1 Literature review 2 The subdistribution beta-stacy process 3 Predictive characterization via reinforced urn processes 4 Theoretical properties 5 A semiparametric competing risks regression model 6 Application: analysis of data from the Amsterdam Cohort Studies 7 Future work Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
5 Literature review I Competing risks data has received widespread attention in the classical nonparametric literature, c.f. Kalbfleisch and Prentice [22], Aalen et al. [28], Andersen et al. [212], Lawless [211], Crowder [212], and Pintilie [26]. Classical Kalbfleisch and Prentice [22] estimator of the subdistribution function: F (t, c) = t Ŝ(u 1) Âc(u) u=1 where: Ŝ(t) is the Kaplan-Meier estimator of the survival function of T ; Â c (t) is the Nelson-Aalen estimator of the cause-specific cumulative hazard function for the event C = c. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
6 Literature review II In contrast with the frequentist literature, the Bayesian nonparametric literature on competing risks is still sparse. Nonparametric bayesian inference for competing risks involves a process prior on the space of subdistribution functions. Current proposals are based on either the gamma process [Ge and Chen, 212] or the beta process [Hjort, 199, De Blasi and Hjort, 27] and its extensions, such as the beta-dirichlet process [Kim et al., 212, Chae et al., 213], for subdistribution or cause-specific hazard functions. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
7 The subdistribution beta-stacy process prior Let α t,c > for all t = 1, 2,... and c = 1,..., K and suppose that τ α t, lim τ + k d= α =. (1) t,d t=1 Subdistribution beta-stacy process A random subdistribution function F is subdistribution beta-stacy if: 1 F (, c) = with probability 1 for all c = 1,..., k; 2 for all c = 1,..., k and all t 1, t 1 F (t, c) = W t,c u=1 ( 1 where F (t, c) = F (t, c) F (t 1, c) and ) k W u,d, d=1 W t = (W t,,..., W t,k ) ind Dirichlet(α t,,..., α t,k ). Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
8 Reinforced Urn Processes A reinforced urn process (RUP) is a stochastic process (X n ) n with countable state-space S and X = x S. Each state x S is associated with a Pólya urn U(x) with balls of colors represented by the elements of the finite set E; If X n = x and a ball of color c E is extracted from U(x), this is replaced in U(x) together with another ball of the same color (reiforcement). Then, X n+1 = q(x, c) S, where q(x, c) is the law of motion. Theorem (Muliere et al. [2]) If (X n ) n is recurrent, then there exists a random transition matrix Q on S conditionally on which (X n ) n is a Markov Chain with transition matrix Q. The rows of Q are independent Dirichlet processes with base measure determined by the initial composition of the urns. RUPs provide a predictive characterization of Pólya trees and discrete-time neutral-to-the-right processes [Muliere et al., 2]. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
9 Construction via reinforced urn processes I a) Patient 1 b) Patient 2 Color 2 Color 2 (,2) (1,2) (2,2) X 3 =(3,2) (,2) (1,2) (2,2) (3,2) Color 1 Color 1 (,1) (1,1) (2,1) (3,1) (,1) (1,1) X 6 =(2,1) (3,1) Color Color X =(,) X 1 =(1,) X 2 =(2,) (3,) X 4 =(,) X 5 =(1,) (2,) (3,) time 1 2 X = (, ), X 1 = (1, ), X 2 = (2, ), X 3 = (3, 2), X 4 = (, ), X 5 = (1, ), X 6 = (2, 1) (T 1, C 1 ) = (3, 2), (T 2, C 2 ) = (2, 1) 3 time Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
10 Construction via reinforced urn processes II The parameters α t,, α t,1,..., α t,k determine the initial composition of the urn U((t 1, )). Lemma (Recurrency condition) (X n ) n is recurrent if and only if lim τ τ + t=1 α t, k d= α t,d Theorem (Predictive characterization) = (2) Assume (2) holds. The sequence ((T i, C i )) i 1 generated by (X n ) n is exchangeable. Its de Finetti measure is the law of a subdistribution beta-stacy process. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
11 Centering the process Let F be a fixed subdistribution function and ω t > for all t 1. Write F sbs(ω, F ) if F has is a subdistribution beta-stacy process with parameters ) k α t,c = ω t F (t, c), α t, = ω t (1 F (t, d) d=1 Moments of F sbs(ω, F ) E[F (t, c)] = F (t, c) Var (F (t, c)) is a decreasing function of ω t Var (F (t, c)) as ω t + Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
12 Relation with other prior processes If F (t, c) sbs(ω, F ), then K d=1 F (t, d) is a random distribution function distributed according to the discrete-time beta-stacy process of Walker and Muliere [1997]. The cumulative hazards A c (t) = F (t, c)/(1 k d=1 F (t 1, d)) are independent discrete-time beta processes [Hjort, 199]. The subdistribution beta-stacy process is therefore also related to the beta-dirichlet process of Kim et al. [212], a generalization of Hjort s beta process. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
13 Posterior computations with censored data Right-censored data: T i = min(t i, R i ), C i = C i I {T i < R i }, where R i is a random censoring time, i = 1,..., n. Censoring is ignorable if i) R 1,..., R n are independent with common distribution function H, ii) conditional on F and H, (T 1, C 1 ),..., (T n, C n ) and R 1,..., R n are independent [Heitjan and Rubin, 1991]. Theorem 1 If censoring is ignorable and F sbs(ω, F ) a priori, conditionally on (T 1, C 1 ),..., (T n, C n) the posterior distribution of F is sbs(ω, F ). 2 F (t, c) F (t, c) as max(ω t ), where F is the classical Kalbfleisch-Prentice estimator of F. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
14 Semiparametric regression for competing risks A subdistribution beta-stacy mixture regression model: (T i, C i ) ind F i (i = 1,..., n) F i ind sbs(ω(θ, x i ), F ( θ, x i )) F (t, c θ, x i ) = F (1) (c θ)f (2) (t c, θ, x i ) 1 ω t (θ, x i ) = K d=1 F (t, d θ, x i ) θ π(θ) We choose: F (1) = multinomial logistic model; F (2) = Weibull regression model; regression coefficients are assigned diffuse normal priors; Weibull scale parameters are assigned non-informative Gamma(ɛ, ɛ) distributions. Inference can be easily carried out via a Metropolis-Hastings algorithm after marginalizing away the F i from the likelihood. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
15 Analysis of Amsterdam Cohort Studies data I Objective: assessing the long-term prognosis of HIV infected men with respect to the risk of non-si AIDS onset or SI onset as a function of CCR5 genotype: WW (wild type allele on both chromosomes) or WM (mutant allele on one chromosome). A total of 324 male patients were followed-up from the date of HIV infection to the earliest among the dates of AIDSs onset, SI phenotype onset, or right censoring, i.e. death, study drop-out, or end of the study period. A total of 65 (2.1%) patients had a WM genotype, while mean age at HIV infection was about 34.6 years (s.d., 7.2 years). Overall, the 324 patients accumulated person-years of follow up (minimum - maximum follow-up:.1 years years), generating 117 cases of AIDS onset and 17 cases of SI onset. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
16 Analysis of Amsterdam Cohort Studies data II AIDS (subdistr. beta Stacy model) SI (subdistr. beta Stacy model) Cumulative probability WW WM Cumulative probability WM WW Years since HIV infection Years since HIV infection AIDS (Classical estimates) SI (Classical estimates) Cumulative probability WW WM Cumulative probability WM WW Years since HIV infection Years since HIV infection Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
17 Future developments Other reinforced urn schemes could be contemplated. For example, each extracted ball may be reinforced by a general amount m > of new similar balls, or m may be random variable depending on the color of the extracted balls, as in Muliere et al. [26] A continuous-time generalization of the subdistribution beta-stacy process could be considered. We are currently developing a characterization of such process from a predictive perspective by means of the continuous-time urn models of Muliere et al. [23] and Bulla and Muliere [27] A generalization of the reinforced urn process considered in this work could be attempted to characterize a process prior on the space of transition kernels of a Markovian multistate process. Such process could be useful for the predictive Bayesian nonparametric analysis of event-history data [Aalen et al., 28] Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
18 References I Odd Aalen, Ornulf Borgan, and Hakon Gjessing. Survival and event history analysis: a process point of view. Springer Science & Business Media, New York, 28. Per Kragh Andersen, Ornulf Borgan, Richard D Gill, and Niels Keiding. Statistical models based on counting processes. Springer Science & Business Media, New York, 212. Paolo Bulla and Pietro Muliere. Bayesian nonparametric estimation for reinforced markov renewal processes. Statistical Inference for Stochastic Processes, 1(3):283 33, 27. Minwoo Chae, Rafael Weißbach, Kwang Hyun Cho, and Yongdai Kim. A mixture of beta dirichlet processes prior for bayesian analysis of event history data. Journal of the Korean Statistical Society, 42 (3): , 213. Martin J Crowder. Multivariate survival analysis and competing risks. CRC Press, Boca Raton, Florida, 212. Pierpaolo De Blasi and Nils Lid Hjort. Bayesian survival analysis in proportional hazard models with logistic relative risk. Scandinavian Journal of Statistics, 34(1): , 27. Miaomiao Ge and Ming-Hui Chen. Bayesian inference of the fully specified subdistribution model for survival data with competing risks. Lifetime Data Analysis, 18(3): , 212. Ronald B Geskus, Frank A Miedema, Jaap Goudsmit, Peter Reiss, Hanneke Schuitemaker, and Roel A Coutinho. Prediction of residual time to aids and death based on markers and cofactors. Journal of Acquired Immune Deficiency Syndromes, 32(5): , 23. Daniel F Heitjan and Donald B Rubin. Ignorability and coarse data. The Annals of Statistics, 19: , Nils Lid Hjort. Nonparametric Bayes estimators based on beta processes in models for life history data. The Annals of Statistics, 18: , 199. John D Kalbfleisch and Ross L Prentice. The statistical analysis of failure time data. John Wiley & Sons, Hoboken, New Jersey, 2nd edition edition, 22. Yongdai Kim, Lancelot James, and Rafael Weissbach. Bayesian analysis of multistate event history data: beta-dirichlet process prior. Biometrika, 99(1):127 14, 212. Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
19 References II Jerald F Lawless. Statistical models and methods for lifetime data. John Wiley & Sons, Hoboken, New Jersey, 211. P Muliere, P Secchi, and SG Walker. Urn schemes and reinforced random walks. Stochastic Processes and their Applications, 88(1):59 78, 2. Pietro Muliere, Piercesare Secchi, and Stephen G Walker. Reinforced random processes in continuous time. Stochastic Processes and their Applications, 14(1):117 13, 23. Pietro Muliere, Anna Maria Paganoni, and Piercesare Secchi. A randomly reinforced urn. Journal of Statistical Planning and Inference, 136(6): , 26. Melania Pintilie. Competing risks: a practical perspective. John Wiley & Sons, Chichester, England, 26. Stephen Walker and Pietro Muliere. Beta-stacy processes and a generalization of the pólya-urn scheme. The Annals of Statistics, 25: , Andrea Arfè (Bocconi) Subdistribution beta-stacy 8 Feb / 19
Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis
Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis Andrea Arfè Department of Decision Sciences, Bocconi University, 20136 Milan, Italy. E-mail: andrea.arfe@phd.unibocconi.it
More informationarxiv: v1 [math.st] 29 Nov 2018
Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis Andrea Arfé Stefano Peluso Pietro Muliere Accepted manuscript, Scandinavian Journal of Statistics, 2018 arxiv:1811.12304v1
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 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 informationBayesian semiparametric analysis of short- and long- term hazard ratios with covariates
Bayesian semiparametric analysis of short- and long- term hazard ratios with covariates Department of Statistics, ITAM, Mexico 9th Conference on Bayesian nonparametrics Amsterdam, June 10-14, 2013 Contents
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 informationReconstruction of individual patient data for meta analysis via Bayesian approach
Reconstruction of individual patient data for meta analysis via Bayesian approach Yusuke Yamaguchi, Wataru Sakamoto and Shingo Shirahata Graduate School of Engineering Science, Osaka University Masashi
More informationRonald Christensen. University of New Mexico. Albuquerque, New Mexico. Wesley Johnson. University of California, Irvine. Irvine, California
Texts in Statistical Science Bayesian Ideas and Data Analysis An Introduction for Scientists and Statisticians Ronald Christensen University of New Mexico Albuquerque, New Mexico Wesley Johnson University
More informationMulti-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 informationIntroduction to BGPhazard
Introduction to BGPhazard José A. García-Bueno and Luis E. Nieto-Barajas February 11, 2016 Abstract We present BGPhazard, an R package which computes hazard rates from a Bayesian nonparametric view. This
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 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 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 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 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 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 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 informationMultistate models and recurrent event models
Multistate models Multistate models and recurrent event models Patrick Breheny December 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Multistate models In this final lecture,
More informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume II: Probability Emlyn Lloyd University oflancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester - New York - Brisbane
More informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
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 informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
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 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 informationDescription Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see
Title stata.com stcrreg postestimation Postestimation tools for stcrreg Description Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see
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 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 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 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 informationBayesian nonparametrics
Bayesian nonparametrics 1 Some preliminaries 1.1 de Finetti s theorem We will start our discussion with this foundational theorem. We will assume throughout all variables are defined on the probability
More informationSUPERPOSITION OF BETA PROCESSES
SUPERPOSITION OF BETA PROCESSES Pierpaolo De Blasi 1, Stefano Favaro 1 and Pietro Muliere 2 1 Dipartiento di Statistica e Mateatica Applicata, Università degli Studi di Torino Corso Unione Sovietica 218/bis,
More informationA class of neutral to the right priors induced by superposition of beta processes
A class of neutral to the right priors induced by superposition of beta processes Pierpaolo De Blasi Stefano Favaro Pietro Muliere No. 13 December 29 www.carloalberto.org/working_papers 29 by Pierpaolo
More informationBayesian semiparametric inference for the accelerated failure time model using hierarchical mixture modeling with N-IG priors
Bayesian semiparametric inference for the accelerated failure time model using hierarchical mixture modeling with N-IG priors Raffaele Argiento 1, Alessandra Guglielmi 2, Antonio Pievatolo 1, Fabrizio
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 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 informationSimulation-based robust IV inference for lifetime data
Simulation-based robust IV inference for lifetime data Anand Acharya 1 Lynda Khalaf 1 Marcel Voia 1 Myra Yazbeck 2 David Wensley 3 1 Department of Economics Carleton University 2 Department of Economics
More informationBayesian Nonparametric Accelerated Failure Time Models for Analyzing Heterogeneous Treatment Effects
Bayesian Nonparametric Accelerated Failure Time Models for Analyzing Heterogeneous Treatment Effects Nicholas C. Henderson Thomas A. Louis Gary Rosner Ravi Varadhan Johns Hopkins University September 28,
More informationChapter 2. Data Analysis
Chapter 2 Data Analysis 2.1. Density Estimation and Survival Analysis The most straightforward application of BNP priors for statistical inference is in density estimation problems. Consider the generic
More informationDoes Cox analysis of a randomized survival study yield a causal treatment effect?
Published in final edited form as: Lifetime Data Analysis (2015), 21(4): 579 593 DOI: 10.1007/s10985-015-9335-y Does Cox analysis of a randomized survival study yield a causal treatment effect? Odd O.
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 informationEfficiency Comparison Between Mean and Log-rank Tests for. Recurrent Event Time Data
Efficiency Comparison Between Mean and Log-rank Tests for Recurrent Event Time Data Wenbin Lu Department of Statistics, North Carolina State University, Raleigh, NC 27695 Email: lu@stat.ncsu.edu Summary.
More informationBayesian Inference of Interactions and Associations
Bayesian Inference of Interactions and Associations Jun Liu Department of Statistics Harvard University http://www.fas.harvard.edu/~junliu Based on collaborations with Yu Zhang, Jing Zhang, Yuan Yuan,
More informationCOPYRIGHTED MATERIAL CONTENTS. Preface Preface to the First Edition
Preface Preface to the First Edition xi xiii 1 Basic Probability Theory 1 1.1 Introduction 1 1.2 Sample Spaces and Events 3 1.3 The Axioms of Probability 7 1.4 Finite Sample Spaces and Combinatorics 15
More informationBayesian Analysis for Markers and Degradation
Bayesian Analysis for Markers and Degradation Mei-Ling Ting Lee 12, Maria Shubina 1, and Alan Zaslavsky 13 1 Biostatistics Department, Harvard School of Public Health, Boston, USA 2 Channing Laboratory,
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 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 informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
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 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 informationFoundations of Nonparametric Bayesian Methods
1 / 27 Foundations of Nonparametric Bayesian Methods Part II: Models on the Simplex Peter Orbanz http://mlg.eng.cam.ac.uk/porbanz/npb-tutorial.html 2 / 27 Tutorial Overview Part I: Basics Part II: Models
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 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 informationBayesian Statistics. Debdeep Pati Florida State University. April 3, 2017
Bayesian Statistics Debdeep Pati Florida State University April 3, 2017 Finite mixture model The finite mixture of normals can be equivalently expressed as y i N(µ Si ; τ 1 S i ), S i k π h δ h h=1 δ h
More informationBayesian estimation of the discrepancy with misspecified parametric models
Bayesian estimation of the discrepancy with misspecified parametric models Pierpaolo De Blasi University of Torino & Collegio Carlo Alberto Bayesian Nonparametrics workshop ICERM, 17-21 September 2012
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 informationLecture 8: The Metropolis-Hastings Algorithm
30.10.2008 What we have seen last time: Gibbs sampler Key idea: Generate a Markov chain by updating the component of (X 1,..., X p ) in turn by drawing from the full conditionals: X (t) j Two drawbacks:
More informationBivariate Degradation Modeling Based on Gamma Process
Bivariate Degradation Modeling Based on Gamma Process Jinglun Zhou Zhengqiang Pan Member IAENG and Quan Sun Abstract Many highly reliable products have two or more performance characteristics (PCs). The
More informationIndividualized Treatment Effects with Censored Data via Nonparametric Accelerated Failure Time Models
Individualized Treatment Effects with Censored Data via Nonparametric Accelerated Failure Time Models Nicholas C. Henderson Thomas A. Louis Gary Rosner Ravi Varadhan Johns Hopkins University July 31, 2018
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 informationTesting Goodness-of-Fit of a Uniform Truncation Model
Testing Goodness-of-Fit of a Uniform Truncation Model Micha Mandel and Rebecca A. Betensky Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts 02115 email: mmandel@hsph.harvard.edu
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 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 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 informationBayesian Nonparametrics
Bayesian Nonparametrics Peter Orbanz Columbia University PARAMETERS AND PATTERNS Parameters P(X θ) = Probability[data pattern] 3 2 1 0 1 2 3 5 0 5 Inference idea data = underlying pattern + independent
More informationSample-weighted semiparametric estimates of cause-specific cumulative incidence using left-/interval censored data from electronic health records
1 / 22 Sample-weighted semiparametric estimates of cause-specific cumulative incidence using left-/interval censored data from electronic health records Noorie Hyun, Hormuzd A. Katki, Barry I. Graubard
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 informationNonparametric Bayes Estimator of Survival Function for Right-Censoring and Left-Truncation Data
Nonparametric Bayes Estimator of Survival Function for Right-Censoring and Left-Truncation Data Mai Zhou and Julia Luan Department of Statistics University of Kentucky Lexington, KY 40506-0027, U.S.A.
More informationBayesian Nonparametric Regression for Diabetes Deaths
Bayesian Nonparametric Regression for Diabetes Deaths Brian M. Hartman PhD Student, 2010 Texas A&M University College Station, TX, USA David B. Dahl Assistant Professor Texas A&M University College Station,
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 informationBayesian nonparametric inference for random distributions and related functions
J. R. Statist. Soc. B(1999) 61, Part 3, pp. 485^527 Bayesian nonparametric inference for random distributions and related functions Stephen G. Walker, Imperial College of Science, Technology and Medicine,
More informationModel Selection in Bayesian Survival Analysis for a Multi-country Cluster Randomized Trial
Model Selection in Bayesian Survival Analysis for a Multi-country Cluster Randomized Trial Jin Kyung Park International Vaccine Institute Min Woo Chae Seoul National University R. Leon Ochiai International
More informationMultinomial Data. f(y θ) θ y i. where θ i is the probability that a given trial results in category i, i = 1,..., k. The parameter space is
Multinomial Data The multinomial distribution is a generalization of the binomial for the situation in which each trial results in one and only one of several categories, as opposed to just two, as in
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 informationEMPIRICAL ENVELOPE MLE AND LR TESTS. Mai Zhou University of Kentucky
EMPIRICAL ENVELOPE MLE AND LR TESTS Mai Zhou University of Kentucky Summary We study in this paper some nonparametric inference problems where the nonparametric maximum likelihood estimator (NPMLE) are
More informationNonparametric Bayesian Survival Analysis using Mixtures of Weibull Distributions
Nonparametric Bayesian Survival Analysis using Mixtures of Weibull Distributions ATHANASIOS KOTTAS University of California, Santa Cruz ABSTRACT. Bayesian nonparametric methods have been applied to survival
More informationThe Metropolis-Hastings Algorithm. June 8, 2012
The Metropolis-Hastings Algorithm June 8, 22 The Plan. Understand what a simulated distribution is 2. Understand why the Metropolis-Hastings algorithm works 3. Learn how to apply the Metropolis-Hastings
More informationBayesian Nonparametric Inference Methods for Mean Residual Life Functions
Bayesian Nonparametric Inference Methods for Mean Residual Life Functions Valerie Poynor Department of Applied Mathematics and Statistics, University of California, Santa Cruz April 28, 212 1/3 Outline
More informationCASE STUDY: Bayesian Incidence Analyses from Cross-Sectional Data with Multiple Markers of Disease Severity. Outline:
CASE STUDY: Bayesian Incidence Analyses from Cross-Sectional Data with Multiple Markers of Disease Severity Outline: 1. NIEHS Uterine Fibroid Study Design of Study Scientific Questions Difficulties 2.
More informationDiscussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs
Discussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs Michael J. Daniels and Chenguang Wang Jan. 18, 2009 First, we would like to thank Joe and Geert for a carefully
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 informationBAYESIAN ESTIMATION OF LINEAR STATISTICAL MODEL BIAS
BAYESIAN ESTIMATION OF LINEAR STATISTICAL MODEL BIAS Andrew A. Neath 1 and Joseph E. Cavanaugh 1 Department of Mathematics and Statistics, Southern Illinois University, Edwardsville, Illinois 606, USA
More informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationSome Concepts of Probability (Review) Volker Tresp Summer 2018
Some Concepts of Probability (Review) Volker Tresp Summer 2018 1 Definition There are different way to define what a probability stands for Mathematically, the most rigorous definition is based on Kolmogorov
More informationEstimation of Conditional Kendall s Tau for Bivariate Interval Censored Data
Communications for Statistical Applications and Methods 2015, Vol. 22, No. 6, 599 604 DOI: http://dx.doi.org/10.5351/csam.2015.22.6.599 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation of Conditional
More informationOutline. Clustering. Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models
Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models Collaboration with Rudolf Winter-Ebmer, Department of Economics, Johannes Kepler University
More informationPractical Bayesian Quantile Regression. Keming Yu University of Plymouth, UK
Practical Bayesian Quantile Regression Keming Yu University of Plymouth, UK (kyu@plymouth.ac.uk) A brief summary of some recent work of us (Keming Yu, Rana Moyeed and Julian Stander). Summary We develops
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 informationBayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features. Yangxin Huang
Bayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features Yangxin Huang Department of Epidemiology and Biostatistics, COPH, USF, Tampa, FL yhuang@health.usf.edu January
More informationBAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA
BAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA Intro: Course Outline and Brief Intro to Marina Vannucci Rice University, USA PASI-CIMAT 04/28-30/2010 Marina Vannucci
More informationHarvard University. Harvard University Biostatistics Working Paper Series
Harvard University Harvard University Biostatistics Working Paper Series Year 2008 Paper 94 The Highest Confidence Density Region and Its Usage for Inferences about the Survival Function with Censored
More informationModels for Multivariate Panel Count Data
Semiparametric Models for Multivariate Panel Count Data KyungMann Kim University of Wisconsin-Madison kmkim@biostat.wisc.edu 2 April 2015 Outline 1 Introduction 2 3 4 Panel Count Data Motivation Previous
More informationFahrmeir: Discrete failure time models
Fahrmeir: Discrete failure time models Sonderforschungsbereich 386, Paper 91 (1997) Online unter: http://epub.ub.uni-muenchen.de/ Projektpartner Discrete failure time models Ludwig Fahrmeir, Universitat
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 informationEvaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions
Evaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions C. Xing, R. Caspeele, L. Taerwe Ghent University, Department
More informationInference on the Univariate Frailty Model: An Approach Bayesian Reference Analysis
Inference on the Univariate Frailty Model: An Approach Bayesian Reference Analysis Vera Lucia D. Tomazella and Camila Bertini Martins Universidade Federal de São Carlos, SP-Brasil (vera@ufscar.br; cacamila-bertini@yahoo.com.br)
More informationA Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness
A Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness A. Linero and M. Daniels UF, UT-Austin SRC 2014, Galveston, TX 1 Background 2 Working model
More informationBayesian Semiparametric Methods for Joint Modeling Longitudinal and Survival Data
Bayesian Semiparametric Methods for Joint Modeling Longitudinal and Survival Data Tim Hanson Division of Biostatistics University of Minnesota and Adam Branscum University of Kentucky Departments of Biostatistics,
More informationBayesian Regularization
Bayesian Regularization Aad van der Vaart Vrije Universiteit Amsterdam International Congress of Mathematicians Hyderabad, August 2010 Contents Introduction Abstract result Gaussian process priors Co-authors
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 informationFirst Aid Kit for Survival. Hypoxia cohort. Goal. DFS=Clinical + Marker 1/21/2015. Two analyses to exemplify some concepts of survival techniques
First Aid Kit for Survival Melania Pintilie pintilie@uhnres.utoronto.ca Two analyses to exemplify some concepts of survival techniques Checking linearity Checking proportionality of hazards Predicted curves:
More information