Longitudinal breast density as a marker of breast cancer risk
|
|
- Magdalen Martin
- 5 years ago
- Views:
Transcription
1 Longitudinal breast density as a marker of breast cancer risk C. Armero (1), M. Rué (2), A. Forte (1), C. Forné (2), H. Perpiñán (1), M. Baré (3), and G. Gómez (4) (1) BIOstatnet and Universitat de València, (2) BIOstatnet and Universitat de Lleida, (3) Epidemiology and Assessment Unit UDIAT-Diagnostic Centre, Corporació Sanitària Parc Taulí, Sabadell, (4) BIOstatnet and Universitat Politècnica de Catalunya Joint Modeling and Beyond Universiteit Hasselt, April 14-15, 2016
2 Outline Background Scientific objective and study population Statistical methods Results Estimation Prediction Conclusions 2 / 26
3 Sección 1 Background 3 / 26
4 Background Breast density (BD) is a characteristic of the breast tissue that is reflected in mammograms. For most women, breast density decreases with age. The ordinal BI-RADS scale (American College of Radiology, 2013) is a commonly tool used to measure breast density. It categorizes breast density in four groups a (1): almost entirely fatty (low density), b (2): scattered fibroglandular densities (medium density), c (3): heterogeneously dense (high density), and d (4): extremely dense (very high density) Several studies have shown that high breast density is associated with an increased breast cancer (BC) risk. Some research (Huo et al., 2014) concluded that longitudinal measurement of mammographic density might be used to personalize breast cancer screening (Biomarker). 4 / 26
5 Sección 2 Scientific objective and material 5 / 26
6 Objective and material Objective: To study the association structure and intensity between mammographic longitudinal measures of breast density and breast cancer risk. Study population: 13,685 women attending a population-based screening program in the BC early-detection program in the Vallès Occidental East area in Catalonia (Spain), between October 1995 and June Women were followed for vital status or possible diagnosis of BC until December, 2013 (Baré et al., 2003, 2006, 2008). At study entry years old. No history of BC. Invited for a biennial mammographic exam. 6 / 26
7 Study variables and scales: Breast density and breast cancer data Data from the 13,685 women attending our population-based screening program. 81,621 mammographic measures of breast density with the subsequent identification and age of the woman. Survival time: Time from the initial age at entrance to the screening program until BC-diagnosis (Administrative censoring). Time scale: Age (50-69 years). Time zero: 50 years. Women who entered the programme after this age produced left truncated data. Baseline covariates: First-degree relative with BC: Yes-No. Prior breast procedures (breast biopsy, fine needle aspiration, cyst aspiration, breast reconstruction, lumpectomy and surgical treatment): Yes-No. 7 / 26
8 Data description women and mammograms Count 1000 BI RADS Number of screenings per women Years Prior breast procedures: 7.1 % (no diagnosis BC) and 13.2 % (diagnosis BC) First-degree relative with BC: 5.2 % (no diagnosis BC) and 9.8 % (diagnosis BC) 431 women diagnosed with BC (3.15 %) 8 / 26
9 Sección 3 Statistical methods 9 / 26
10 Statistical methodology Elements: Longitudinal data: Breast density data Survival data: Time from the initial age of the screening program to BC detection. Left truncation and administrative censoring. Objective: We want to construct a survival model to assess the risk of a diagnosis of breast cancer by means of baseline information and breast density temporal information. Put a joint model in your problem 10 / 26
11 Objective and material and methods Joint models Different approaches: conditional models, shared-parameter models, random-effects models, and joint latent class models. Shared-parameter models: 1 Longitudinal model for the longitudinal data y. 2 Survival model for the survival data s. 3 Connection between both models by means of a common vector of random-effects, b, which endows them with conditional independence. Frequentist f (y, s, b) = f (y, s b) f (b) Bayesian = f (y b) f (s b) f (b) f (y, s, b, θ, φ) = f (y, s b, θ, φ) f (b, θ, φ) = f (y b, θ, φ) f (s b, θ, φ) f (b θ, φ) π(θ, φ) = f (y b, θ) f (s b, θ) f (b φ) π(θ, φ) where θ and φ are the vector of parameters and hyperparameters, respectively. 11 / 26
12 The joint model, I Longitudinal submodel: A cumulative logit model for the longitudinal ordinal measures of breast density based on the idea of a continuous latent variable (Lunn et al., 2001; Luo et al., 2014). y ij breast density in the BI-RADS scale, y ij {1, 2, 3, 4}, of woman i, i = 1,..., n, at time t ij (age 50 + t ij ), j = 1,..., n i. We assume an underlying continuous latent variable yij that determines the breast density BI-RADS category of woman i at time t ij. Relationship between y ij and y ij y ij = k y ij (γ k 1, γ k ], y ij Logistic(m ij, s = 1) logit P(y ij > k) = m ij γ k, m it = β 0 + b i0 + (β 1 + b i1 )t, with = γ 0 < γ 1 < γ 2 = 0 < γ 3 < γ 4 = unknown cutpoints, β 0, β 1 unknown parameters and (b 0i, b 1i subject specific random effects). 12 / 26
13 The joint model, II Survival submodel: A left-truncated Cox proportional hazard model for the time-to-bc detection that incorporates information from the longitudinal process. T i, i = 1,..., n observed BC detection time for the i-th woman, T i = min(ti, C i), Ti is the true failure time and C i the right-censoring time. The event indicator, δ i = I(Ti C i). Hazard function of Ti in terms of the left truncated Cox proportional hazard model h i (t x i, θ is, t i > a i ) = h 0 (t λ, η 0 ) exp{η 1 Famhist i +η 2 Brstproc i +α m it }, t > a i h 0 (t λ, η 0 ) = λt λ 1 e η 0 : baseline risk function of a Weibull distribution; x i : baseline covariates, First-degree relative with BC (Famhist) and prior breast procedures (Brstproc) with regression coefficients η 1 and η 2 ; α: effect of the individual trajectory of breast density of a woman over their BC risk in terms of the latent BD mean; a i : age over 50 at which woman i enters the screening program thus providing a left truncated time (Uzunogullari et al., 1992); θ is = (λ, η 0, η 1, η 2, α, θ il ) T, with θ il = (β 0, β 1, b 0i, b 1i ) T. 13 / 26
14 The joint model, III 1 Longitudinal submodel: A cumulative logit model for the longitudinal ordinal measures of breast density based on the idea of a continuous latent variable (Lunn et al., 2001; Luo et al., 2014), and y ij = k y ij (γ k 1, γ k ], logit P(y ij > k) = m ij γ k, m it = β 0 + b i0 + (β 1 + b i1)t, 2 Survival submodel: A left-truncated Weibull proportional hazard model for the time-to-bc detection, that incorporates information from the longitudinal process. h i(t x i, θ is, ti > a i) = h 0(t λ, η 0) exp{η 1Famhist i + η 2Brstproc i + α m i(t)} = h 0(t λ, η 0) exp{η 1Famhist i + η 2Brstproc i + α(β 0 + b i0 + (β 1 + b i1)t)}, t 3 Both processes are connected through a shared vector of random effects, which, in the presence of covariates and parameters, endows both processes with conditional independence (Rizopoulos, 2012). 14 / 26
15 Bayes Inference Bayesian Inference: Complete specification of the joint model needs the elicitation of a prior distribution for the subsequent parameters and hyperparameters. Prior independence as a default specification and wide proper prior distributions with the aim of giving all inferential prominence to the data. Normal distributions for the regression coefficients of the longitudinal and survival submodels and the association coefficient. Truncated Normal distributions for the cutpoints of the latent scale, = γ 0 < γ 1 < γ 2 = 0 < γ 3 < γ 4 = selected to provide the same prior probability to each response category. Gamma distribution, Ga(1, 1), for the parameter of the baseline hazard function because it mimics a constant baseline hazard function (Guo et al., ) Uniform hyperdistribution distribution Un(0, 10) for the standard deviations σ 0 and σ 1 of the normal random effects b i0 N(0, σ 0) and b i1 N(0, σ 1). 15 / 26
16 Sección 4 Results 16 / 26
17 Posterior distribution Longitudinal submodel: A cumulative logit model for the longitudinal ordinal measures of breast density based on the idea of a continuous latent variable (Lunn et al., 2001; Luo et al., 2014). y ij = k y ij (γ k 1, γ k ], logit P(y ij > D k ) = m ij γ k, m it = β 0 + b i0 + (β 1 + b i1 )t, Survival submodel: A left-truncated Weibull proportional hazard model for the time-to-bc detection, that incorporates information from the longitudinal process. h i (t x i, θ is, t i > a i ) = h 0 (t λ, η 0 ) exp{η 1 Famhist i + η 2 Brstproc i + α m it } = h 0 (t λ, η 0 ) exp{η 1 Famhist i + η 2 Brstproc i + α(β 0 + b i0 + (β 1 + b i1 )t)}, t > a i Posterior distribution: Approximated by Markov Chain Monte Carlo methods. Mean Sd 2.5 % Median 97.5 % P( > 0 D) β β σ σ γ γ η η η α λ / 26
18 Posterior distribution for the probabilities associated to breast density categories P(y it = a, b, c, d θ, b i, φ) depends on θ and b i. Posterior distribution for the probability associated to each breast density categories for a generic woman i of the same population at time t, π(p(y it = a, b, c, d θ, b i, φ) Data) approximated from f (b i φ) and π(θ, φ Data). Posterior mean 95 % credible interval of π(p(y it = a, b, c, d θ, b i, φ) Data) 0.6 Probability a b c d Age (years) 18 / 26
19 Posterior predictive distribution for the latent breast density Conditional latent breast density (y it θ, b i, φ) Logistic(m it, s = 1). Posterior predictive distribution of the latent breast density of a generic woman i at t years over 50, (y it Data), (50, 55, 60, 65 and 70 in the graphic) Years y* Age (years) Probability a b c d 19 / 26
20 Posterior distribution for the hazard ratios, I Posterior distribution for the hazard ratio (HR): HR of BC diagnosis in women with prior breast procedures versus women without prior breast procedures. HR of BC diagnosis in women with family history of BC versus women without family history of BC. HR of BC diagnosis in women with family history of BC and prior breast procedures versus women without family history of BC and without prior breast procedures var Brstproc Famhist Famhist and Brstproc Hazard Ratio 20 / 26
21 Posterior distribution for the hazard ratios, II Posterior distribution of the HR of a BC diagnosis for a generic woman i with breast density in category k at time t versus this generic woman with breast density in category k at t is NOT POSSIBLE because h i (t x i, θ is, t i h i (t x i, θ is, t i > a i ) = h 0 (t λ, η 0 ) exp{η 1 Famhist i + η 2 Brstproc i + α m it }, YES > a i ) = h 0 (t λ, η 0 ) exp{η 1 Famhist i + η 2 Brstproc i + α y it }, NO We have constructed an ad hoc procedure to represent each value of the breast density at time t in terms of a particular value of the true latent breast density at t ordinal y it latent yit true latent m it (yit = m it + ɛ it ) t = 0 BI-RADS b c d a (1.450, 1.904) (1.959, 3.207) (2.437, 4.684) b (1.351, 1.685) (1.681, 2.460) c (1.244, 1.460) t = 10 BI-RADS b c d a (1.464, 1.937) (2.054, 3.483) (2.569, 5.132) b (1.403, 1.798) (1.754, 2.649) c (1.251, 1.473) 21 / 26
22 Dynamic estimation and prediction, I Bayesian reasoning Estimation of the conditional survival probability that a woman i with longitudinal information Y i at t i,ni and free of a diagnosis of breast cancer at t i,ni is also free at t, P(T i > t θ, b, φ, Y i, x i, T i > t i,ni ), by means of its posterior distribution π(p(t i > t θ, b, φ, Y i, x i, T i > t i,ni ) Data), t > t i,ni that we approach from the posterior distribution π(θ, b, φ Y i, T i > t i,ni, Data). Prediction of the longitudinal breast cancer measurement of a given woman in a scheduled time t ni +1, as the predictive distribution P(y i,tni +1 = a, b, c, d Y i, x i, Data), k = 1, 2, 3, 4 computed from The conditional probability P(y i,tni +1 = a, b, c, d θ, b, φ, Y i, x i, T i > t i,ni ), and The posterior distribution π(θ, b, φ Y i, T i > t i,ni, Data) 22 / 26
23 Dynamic prediction and prediction, II Both posterior distributions π(p(t i > t θ, b, φ, Y i, x i, T i > t i,ni ) Data), t > t i,ni P(y i,tni +1 = a, b, c, d Y i, x i, Data), k = 1, 2, 3, 4 based on π(θ, b, φ Y i, T i > t i,ni, Data) and could be applied to women in the study as well as women not in the study. But we only have an approximate sample from π(θ, b, φ Data). Consequently, if woman i is In the study and Data includes Y i, x i, T i > t i,ni then π(θ, b, φ Data) = π(θ, b, φ Y i, T i > t i,ni, Data) In the study but with new longitudinal and survival information not included in Data then Y i, x i, T i > t i,ni then π(θ, b, φ Data) π(θ, b, φ Y i, T i > t i,ni, Data) Update the posterior distribution with the new information. Sequential MCMC methods in progress. Not in the study with sequential interest we need to update π(θ, b, φ, b i Y i, T i > t i,ni, Data) Update the posterior distribution with the new information. Sequential MCMC methods in progress. 23 / 26
24 Dynamic estimation and prediction, III Posterior distribution: π(p(t i > t θ, b, φ, Y i, x i, T i > t i,ni ) D), t > t i,ni f (y i,tni +1 Y i, x i, D) Posterior mean and 95 % credible interval of the survival probability distribution and posterior mean of the predictive mammogram breast density of woman id17540 (No family history of BC and no prior breast procedures) with a given follow-up which includes all her historical longitudinal mammogram and covariates. 24 / 26
25 Sección 5 Conclusions 25 / 26
26 Conclusions Joint models of longitudinal and survival data are a suitable tool for analyzing the relationship between mammographic breast density and breast cancer risk. Our joint model is a good starting point for learning about the problem but we need to introduce more flexibility and work more on the problem: Non-linear trajectories in the longitudinal mammogram BD modeling. Leave the Weibull baseline risk function and explore other proposals which allow for multimodal or heavy tailed patterns. Assess the performance of the biomarker: Sensitivity, specificity and time-dependent receiver operating characteristic (ROC) curves. Sensitivity of the prior distributions, in particular the hyperpriors for the random effects. Improve our knowledge about latent variables in ordinal models. 26 / 26
CASE 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 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 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 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 informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
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 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 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 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 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 informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
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 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 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 informationDynamic Scheduling of the Upcoming Exam in Cancer Screening
Dynamic Scheduling of the Upcoming Exam in Cancer Screening Dongfeng 1 and Karen Kafadar 2 1 Department of Bioinformatics and Biostatistics University of Louisville 2 Department of Statistics University
More informationBayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information
Statistics Preprints Statistics 8-2010 Bayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information Jiqiang Guo Iowa State University, jqguo@iastate.edu
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 informationNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Approximate Bayesian Inference for nonhomogeneous Poisson processes with application to survival analysis by Rupali Akerkar, Sara Martino and Håvard Rue PREPRINT
More informationBayesian (conditionally) conjugate inference for discrete data models. Jon Forster (University of Southampton)
Bayesian (conditionally) conjugate inference for discrete data models Jon Forster (University of Southampton) with Mark Grigsby (Procter and Gamble?) Emily Webb (Institute of Cancer Research) Table 1:
More informationBayesian Inference in GLMs. Frequentists typically base inferences on MLEs, asymptotic confidence
Bayesian Inference in GLMs Frequentists typically base inferences on MLEs, asymptotic confidence limits, and log-likelihood ratio tests Bayesians base inferences on the posterior distribution of the unknowns
More informationNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Investigating posterior contour probabilities using INLA: A case study on recurrence of bladder tumours by Rupali Akerkar PREPRINT STATISTICS NO. 4/2012 NORWEGIAN
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 informationBayesian Learning. HT2015: SC4 Statistical Data Mining and Machine Learning. Maximum Likelihood Principle. The Bayesian Learning Framework
HT5: SC4 Statistical Data Mining and Machine Learning Dino Sejdinovic Department of Statistics Oxford http://www.stats.ox.ac.uk/~sejdinov/sdmml.html Maximum Likelihood Principle A generative model for
More informationVCMC: Variational Consensus Monte Carlo
VCMC: Variational Consensus Monte Carlo Maxim Rabinovich, Elaine Angelino, Michael I. Jordan Berkeley Vision and Learning Center September 22, 2015 probabilistic models! sky fog bridge water grass object
More informationIncorporating unobserved heterogeneity in Weibull survival models: A Bayesian approach
Incorporating unobserved heterogeneity in Weibull survival models: A Bayesian approach Catalina A. Vallejos 1 Mark F.J. Steel 2 1 MRC Biostatistics Unit, EMBL-European Bioinformatics Institute. 2 Dept.
More informationPhysician Performance Assessment / Spatial Inference of Pollutant Concentrations
Physician Performance Assessment / Spatial Inference of Pollutant Concentrations Dawn Woodard Operations Research & Information Engineering Cornell University Johns Hopkins Dept. of Biostatistics, April
More informationSTA 216, GLM, Lecture 16. October 29, 2007
STA 216, GLM, Lecture 16 October 29, 2007 Efficient Posterior Computation in Factor Models Underlying Normal Models Generalized Latent Trait Models Formulation Genetic Epidemiology Illustration Structural
More informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
More informationMultivariate spatial modeling
Multivariate spatial modeling Point-referenced spatial data often come as multivariate measurements at each location Chapter 7: Multivariate Spatial Modeling p. 1/21 Multivariate spatial modeling Point-referenced
More informationSTAT 499/962 Topics in Statistics Bayesian Inference and Decision Theory Jan 2018, Handout 01
STAT 499/962 Topics in Statistics Bayesian Inference and Decision Theory Jan 2018, Handout 01 Nasser Sadeghkhani a.sadeghkhani@queensu.ca There are two main schools to statistical inference: 1-frequentist
More informationMultilevel Statistical Models: 3 rd edition, 2003 Contents
Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction
More 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 informationBayesian Inference. Chapter 4: Regression and Hierarchical Models
Bayesian Inference Chapter 4: Regression and Hierarchical Models Conchi Ausín and Mike Wiper Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative
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 informationA Bayesian Nonparametric Model for Predicting Disease Status Using Longitudinal Profiles
A Bayesian Nonparametric Model for Predicting Disease Status Using Longitudinal Profiles Jeremy Gaskins Department of Bioinformatics & Biostatistics University of Louisville Joint work with Claudio Fuentes
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 informationAn Algorithm for Bayesian Variable Selection in High-dimensional Generalized Linear Models
Proceedings 59th ISI World Statistics Congress, 25-30 August 2013, Hong Kong (Session CPS023) p.3938 An Algorithm for Bayesian Variable Selection in High-dimensional Generalized Linear Models Vitara Pungpapong
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 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 informationReview. Timothy Hanson. Department of Statistics, University of South Carolina. Stat 770: Categorical Data Analysis
Review Timothy Hanson Department of Statistics, University of South Carolina Stat 770: Categorical Data Analysis 1 / 22 Chapter 1: background Nominal, ordinal, interval data. Distributions: Poisson, binomial,
More informationSTAC51: Categorical data Analysis
STAC51: Categorical data Analysis Mahinda Samarakoon January 26, 2016 Mahinda Samarakoon STAC51: Categorical data Analysis 1 / 32 Table of contents Contingency Tables 1 Contingency Tables Mahinda Samarakoon
More informationLatent Variable Models for Binary Data. Suppose that for a given vector of explanatory variables x, the latent
Latent Variable Models for Binary Data Suppose that for a given vector of explanatory variables x, the latent variable, U, has a continuous cumulative distribution function F (u; x) and that the binary
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 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 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 informationA short introduction to INLA and R-INLA
A short introduction to INLA and R-INLA Integrated Nested Laplace Approximation Thomas Opitz, BioSP, INRA Avignon Workshop: Theory and practice of INLA and SPDE November 7, 2018 2/21 Plan for this talk
More informationThe Gaussian process latent variable model with Cox regression
The Gaussian process latent variable model with Cox regression James Barrett Institute for Mathematical and Molecular Biomedicine (IMMB), King s College London James Barrett (IMMB, KCL) GP Latent Variable
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 informationMonte Carlo in Bayesian Statistics
Monte Carlo in Bayesian Statistics Matthew Thomas SAMBa - University of Bath m.l.thomas@bath.ac.uk December 4, 2014 Matthew Thomas (SAMBa) Monte Carlo in Bayesian Statistics December 4, 2014 1 / 16 Overview
More informationMarkov Chain Monte Carlo in Practice
Markov Chain Monte Carlo in Practice Edited by W.R. Gilks Medical Research Council Biostatistics Unit Cambridge UK S. Richardson French National Institute for Health and Medical Research Vilejuif France
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 informationUsing Mixture Latent Markov Models for Analyzing Change in Longitudinal Data with the New Latent GOLD 5.0 GUI
Using Mixture Latent Markov Models for Analyzing Change in Longitudinal Data with the New Latent GOLD 5.0 GUI Jay Magidson, Ph.D. President, Statistical Innovations Inc. Belmont, MA., U.S. statisticalinnovations.com
More informationBayesian Modeling of Accelerated Life Tests with Random Effects
Bayesian Modeling of Accelerated Life Tests with Random Effects Ramón V. León Avery J. Ashby Jayanth Thyagarajan Joint Statistical Meeting August, 00 Toronto, Canada Abstract We show how to use Bayesian
More informationEstimation of Optimally-Combined-Biomarker Accuracy in the Absence of a Gold-Standard Reference Test
Estimation of Optimally-Combined-Biomarker Accuracy in the Absence of a Gold-Standard Reference Test L. García Barrado 1 E. Coart 2 T. Burzykowski 1,2 1 Interuniversity Institute for Biostatistics and
More informationJoint longitudinal and time-to-event models via Stan
Joint longitudinal and time-to-event models via Stan Sam Brilleman 1,2, Michael J. Crowther 3, Margarita Moreno-Betancur 2,4,5, Jacqueline Buros Novik 6, Rory Wolfe 1,2 StanCon 2018 Pacific Grove, California,
More informationIntroduction to Bayesian Methods. Introduction to Bayesian Methods p.1/??
to Bayesian Methods Introduction to Bayesian Methods p.1/?? We develop the Bayesian paradigm for parametric inference. To this end, suppose we conduct (or wish to design) a study, in which the parameter
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 informationBayesian Hierarchical Models
Bayesian Hierarchical Models Gavin Shaddick, Millie Green, Matthew Thomas University of Bath 6 th - 9 th December 2016 1/ 34 APPLICATIONS OF BAYESIAN HIERARCHICAL MODELS 2/ 34 OUTLINE Spatial epidemiology
More informationBayesian Analysis for Step-Stress Accelerated Life Testing using Weibull Proportional Hazard Model
Noname manuscript No. (will be inserted by the editor) Bayesian Analysis for Step-Stress Accelerated Life Testing using Weibull Proportional Hazard Model Naijun Sha Rong Pan Received: date / Accepted:
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 informationSpatial Statistics with Image Analysis. Outline. A Statistical Approach. Johan Lindström 1. Lund October 6, 2016
Spatial Statistics Spatial Examples More Spatial Statistics with Image Analysis Johan Lindström 1 1 Mathematical Statistics Centre for Mathematical Sciences Lund University Lund October 6, 2016 Johan Lindström
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationHybrid Censoring; An Introduction 2
Hybrid Censoring; An Introduction 2 Debasis Kundu Department of Mathematics & Statistics Indian Institute of Technology Kanpur 23-rd November, 2010 2 This is a joint work with N. Balakrishnan Debasis Kundu
More informationStat 5421 Lecture Notes Proper Conjugate Priors for Exponential Families Charles J. Geyer March 28, 2016
Stat 5421 Lecture Notes Proper Conjugate Priors for Exponential Families Charles J. Geyer March 28, 2016 1 Theory This section explains the theory of conjugate priors for exponential families of distributions,
More informationApproximate Bayesian computation for spatial extremes via open-faced sandwich adjustment
Approximate Bayesian computation for spatial extremes via open-faced sandwich adjustment Ben Shaby SAMSI August 3, 2010 Ben Shaby (SAMSI) OFS adjustment August 3, 2010 1 / 29 Outline 1 Introduction 2 Spatial
More informationReview: Probabilistic Matrix Factorization. Probabilistic Matrix Factorization (PMF)
Case Study 4: Collaborative Filtering Review: Probabilistic Matrix Factorization Machine Learning for Big Data CSE547/STAT548, University of Washington Emily Fox February 2 th, 214 Emily Fox 214 1 Probabilistic
More informationNearest Neighbor Gaussian Processes for Large Spatial Data
Nearest Neighbor Gaussian Processes for Large Spatial Data Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns
More informationFast Likelihood-Free Inference via Bayesian Optimization
Fast Likelihood-Free Inference via Bayesian Optimization Michael Gutmann https://sites.google.com/site/michaelgutmann University of Helsinki Aalto University Helsinki Institute for Information Technology
More informationSeparate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages
Separate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages Xu GUO and Bradley P. CARLIN Many clinical trials and other medical and reliability studies generate both
More informationBuilding a Prognostic Biomarker
Building a Prognostic Biomarker Noah Simon and Richard Simon July 2016 1 / 44 Prognostic Biomarker for a Continuous Measure On each of n patients measure y i - single continuous outcome (eg. blood pressure,
More informationBayesian philosophy Bayesian computation Bayesian software. Bayesian Statistics. Petter Mostad. Chalmers. April 6, 2017
Chalmers April 6, 2017 Bayesian philosophy Bayesian philosophy Bayesian statistics versus classical statistics: War or co-existence? Classical statistics: Models have variables and parameters; these are
More informationContinuous Time Survival in Latent Variable Models
Continuous Time Survival in Latent Variable Models Tihomir Asparouhov 1, Katherine Masyn 2, Bengt Muthen 3 Muthen & Muthen 1 University of California, Davis 2 University of California, Los Angeles 3 Abstract
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 informationHastings-within-Gibbs Algorithm: Introduction and Application on Hierarchical Model
UNIVERSITY OF TEXAS AT SAN ANTONIO Hastings-within-Gibbs Algorithm: Introduction and Application on Hierarchical Model Liang Jing April 2010 1 1 ABSTRACT In this paper, common MCMC algorithms are introduced
More informationCIMAT Taller de Modelos de Capture y Recaptura Known Fate Survival Analysis
CIMAT Taller de Modelos de Capture y Recaptura 2010 Known Fate urvival Analysis B D BALANCE MODEL implest population model N = λ t+ 1 N t Deeper understanding of dynamics can be gained by identifying variation
More informationBayesian Inference. Chapter 4: Regression and Hierarchical Models
Bayesian Inference Chapter 4: Regression and Hierarchical Models Conchi Ausín and Mike Wiper Department of Statistics Universidad Carlos III de Madrid Advanced Statistics and Data Mining Summer School
More informationLearning Bayesian Networks for Biomedical Data
Learning Bayesian Networks for Biomedical Data Faming Liang (Texas A&M University ) Liang, F. and Zhang, J. (2009) Learning Bayesian Networks for Discrete Data. Computational Statistics and Data Analysis,
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee and Andrew O. Finley 2 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department
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 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 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 informationOptimal rules for timing intercourse to achieve pregnancy
Optimal rules for timing intercourse to achieve pregnancy Bruno Scarpa and David Dunson Dipartimento di Statistica ed Economia Applicate Università di Pavia Biostatistics Branch, National Institute of
More informationIntegrated Non-Factorized Variational Inference
Integrated Non-Factorized Variational Inference Shaobo Han, Xuejun Liao and Lawrence Carin Duke University February 27, 2014 S. Han et al. Integrated Non-Factorized Variational Inference February 27, 2014
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationThe Weibull Distribution
The Weibull Distribution Patrick Breheny October 10 Patrick Breheny University of Iowa Survival Data Analysis (BIOS 7210) 1 / 19 Introduction Today we will introduce an important generalization of the
More informationFoundations of Statistical Inference
Foundations of Statistical Inference Julien Berestycki Department of Statistics University of Oxford MT 2015 Julien Berestycki (University of Oxford) SB2a MT 2015 1 / 16 Lecture 16 : Bayesian analysis
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 informationDensity Estimation. Seungjin Choi
Density Estimation Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr http://mlg.postech.ac.kr/
More informationFitting the Bartlett-Lewis rainfall model using Approximate Bayesian Computation
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Fitting the Bartlett-Lewis rainfall model using Approximate Bayesian
More informationOn Bayesian Computation
On Bayesian Computation Michael I. Jordan with Elaine Angelino, Maxim Rabinovich, Martin Wainwright and Yun Yang Previous Work: Information Constraints on Inference Minimize the minimax risk under constraints
More informationEstimating the marginal likelihood with Integrated nested Laplace approximation (INLA)
Estimating the marginal likelihood with Integrated nested Laplace approximation (INLA) arxiv:1611.01450v1 [stat.co] 4 Nov 2016 Aliaksandr Hubin Department of Mathematics, University of Oslo and Geir Storvik
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 informationSpatial discrete hazards using Hierarchical Bayesian Modeling
Spatial discrete hazards using Hierarchical Bayesian Modeling Mathias Graf ETH Zurich, Institute for Structural Engineering, Group Risk & Safety 1 Papers -Maes, M.A., Dann M., Sarkar S., and Midtgaard,
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationBayesian Estimation for the Generalized Logistic Distribution Type-II Censored Accelerated Life Testing
Int. J. Contemp. Math. Sciences, Vol. 8, 2013, no. 20, 969-986 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2013.39111 Bayesian Estimation for the Generalized Logistic Distribution Type-II
More informationA Hierarchical Mixture Dynamic Model of School Performance in the Brazilian Mathematical Olympiads for Public Schools (OBMEP)
A of School Performance in the ian Mathematical Olympiads for Public Schools (OBMEP) Alexandra M. Schmidt IM - UFRJ, Homepage: www.dme.ufrj.br/ alex joint work with Caroline P. de Moraes and Helio S. Migon
More informationStable Limit Laws for Marginal Probabilities from MCMC Streams: Acceleration of Convergence
Stable Limit Laws for Marginal Probabilities from MCMC Streams: Acceleration of Convergence Robert L. Wolpert Institute of Statistics and Decision Sciences Duke University, Durham NC 778-5 - Revised April,
More informationPetr Volf. Model for Difference of Two Series of Poisson-like Count Data
Petr Volf Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodárenskou věží 4, 182 8 Praha 8 e-mail: volf@utia.cas.cz Model for Difference of Two Series of Poisson-like
More informationA Study of Joinpoint Models for Longitudinal Data
UNLV Theses, Dissertations, Professional Papers, and Capstones 8-1-2014 A Study of Joinpoint Models for Longitudinal Data Libo Zhou University of Nevada, Las Vegas, zhoul@unlv.nevada.edu Follow this and
More informationarxiv: v1 [stat.me] 22 Mar 2015
Personalized Screening Intervals for Biomarkers using Joint Models for Longitudinal and Survival Data arxiv:1503.06448v1 [stat.me] 22 Mar 2015 Dimitris Rizopoulos 1,, Jeremy M.G. Taylor 2, Joost van Rosmalen
More information