Survival models and health sequences
|
|
- Amice Small
- 6 years ago
- Views:
Transcription
1 Survival models and health sequences Walter Dempsey University of Michigan July 27, 2015
2 Survival Data Problem Description Survival data is commonplace in medical studies, consisting of failure time information for each patient i, T i. Many studies require patient monitoring, generating a series of measurements of a health process, Y i. Acompleteuncensoredobservationforonepatientconsistsofthefollowingtriple: (Y, T, t) Y health process (ex. CD4 cell count, Prothrombin ndex, Quality of Life ndex) T failure time as measured from recruitment t appointment schedule Joint modeling of the repeated measurements and survival data is necessary in order to gauge the e ect of covariates on the response. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
3 Survival Process Problem Description Survival Process AstochasticprocessY i (t) :R!Rwhere R is the state space. Typically Y i (t) is state of health or quality of life of patient i at time t, typically measured since recruitment ex. Simple Survival Process, state space is R = {0, 1} is su cient Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
4 Survival Process Problem Description Survival Process AstochasticprocessY i (t) :R!Rwhere R is the state space. Typically Y i (t) is state of health or quality of life of patient i at time t, typically measured since recruitment ex. Simple Survival Process, state space is R = {0, 1} is su cient Flatlining An absorbing state [ 2Rsuch that Y (t) =[ ) Y (t 0 )=[ for all t 0 > t. The survival time, T, is the time to failure: T i =sup{t : Y i (t) 6= [} Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
5 Appointment Schedule Problem Description The appointment schedule is a random subset t [0, T ), which is informative for survival regardless of health trajectory Sequential Conditional ndependence Consider patient with k appointments t (k) =(t 0 <...<t k 1 ). The sequence Y [t (k) ] may a ect the scheduled appointment date t k.thesequential conditional independence assumption states t k? Y (T, t (k), Y [t (k) ]). (1) (Example) Half of patients had appointment within the last three weeks of life! unclear if assumption is violated by patient-initiated appointments (71 within 10 days prior to death). Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
6 Problem Description Motivating Example : Prednisone Case Study From 1962 to 1969, 488 patients with histologically verified liver cirrhosis at hospitals in Copenhagen were randomly assigned to the two treatment arms. 251 and 237 patients in the prednisone and placebo treatment groups respectively 292 uncensored records ) 40% patients have censored records (45% of total measurements correspond to censored patient records) The purpose of the trial was to ascertain whether prednisone prolonged survival for patients with cirrhosis. Prothrombin ndex Time Since Recruitment Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
7 Problem Description Motivating Example : Prednisone Case Study Table 1 shows average Y -values by T and t. Thecellaverages8 266non-independent, highly variable measurements. Table 1 : Average prothrombin levels indexed by T and t Survival Time t after recruitment (yrs) time (T ) Suggests reverse-time is a more e ective way of organizing the data to display main trends in the mean response. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
8 Revival Process Section 2 Revival Process Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
9 The Revival Process Revival Process Revival Model Assuming the survival time is finite with probability one, wecandefinetherevival Process Z i (s) =Y i (T i s) The revival process is the re-aligned health process. This is well-defined conditional on the survival time and provides an invertible mapping from the survival process to the joint survival time and revival process. Y $ (Z, T ) We consider this joint process in lieu of the survival process as we hope it provides e ective alignment of patient records for comparison and signal extraction. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
10 The Revival Process Revival Process Revival Model Table 2 confirms there is considerable excess variation associated with rows (116.8) and with the reverse-time factor (77.8) but not so much with columns. Table 2 : ANOVA decomposition for Table 1 Source U/V kp UY k 2 kp VY k 2 d.f. M.S. Diagonal (R + C + D)/(R + C) Column (R + C + D)/(R + D) Row (R + C + D)/(C + D) Residual RC/(R + C + D) ) Time reversal yields e ective alignment of patient records for comparison and signal extraction. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
11 Covariates Revival Process Revival Model Definition (Temporal and time-evolving variable) A temporal variable x is a function defined for every t 0. t is a covariate if is a function on the units (entire function is determined at baseline) Typically implies x constant, but there are exceptions (eg. patient s age) A time-evolving variable x 0 is a temporal variable that is not a function on the units (eg. marital status, quality of life, and air quality) Every time-evolving variable is necessarily part of the response process The joint distribution of time-evolving variables and survival time may be used to predict survival time beyond t whose status is known at times t prior to t. Probabilistic prediction is not possible without the requisite mathematical structure of -fields. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
12 Revival Process Revival Model Treatment e ect: definition and estimation For patient i, leta i (t) be the treatment arm scheduled for patient i at time t. Null level required for times at and before recruitment, t apple 0. Entire temporal trajectory, a i (t), for t 0isspecifiedatbaselineanddeterminedby randomization (ie. it is a time-dependent covariate). Let ā i (s) =a i (T s) bethetreatmentarmexpressedinrevivaltime. Z? T ā because T is a function of ā Lack of interference () :thetreatmentassignedtooneindividualhasnoe ecton the response distribution for other individuals. Lack of interference () :thetreatmentprotocolatonetimehasnoe ectonthe response distribution at other times. Z[s]? ā ā[s] Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
13 Revival Process Revival Model Treatment e ect: definition and estimation Consider two patients, one in each treatment arm, a i (t) =ā i (T i t)=1, a j (t) =ā j (T j t)=0 such that x i = x j. The conventional treatment definitions are 10(t) =E(Y i (t) E(Y j (t)) or 0 10(t) =E(Y i (t) E(Y j (t) T i, T j > t) Supposing the dependence on T is additive, the di erence of means at revival time s E(Z i (s) T ) E(Z j (s) T )= 10(s)+ (T i ) (T j ) contains both a treatment e ect and an e ect due to the di erence in survival times. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
14 Conditional Distribution Survival Prediction The joint density of (T, t (k), Y [t (k) ]at(t, t (k), y) is a product of three factors f (t) Y j<k p(t j, y j H j, T = t) = f (t) Y j<k p(y j H j, T = t) Y j<k p(t j H j, T = t) = f (t) g k (y; t t (k) t) Y j<k p(t j H j, T = t), (2) where f = F 0 is the survival density, and H j is the observed history (t (j), Y [t (j) ]) at time t j 1. We assume the appointment schedule is uninformative for prediction in the sense that for t k 1 < t k < t. p(t k H k, T = t) =p(t k H k, T = 1) (3) Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
15 Parameter Estimation Likelihood Factorization : Survival distribution specification Re-alignment by the survival time requires the survival time to be finite with probability one. pr(t < 1) =1 Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
16 Parameter Estimation Likelihood Factorization : Survival distribution specification Re-alignment by the survival time requires the survival time to be finite with probability one. pr(t < 1) =1 Harmonic Process The harmonic process is an exchangeable survival process defined by two non-negative parameters (, ). T i Exp ( ( (1 + ) ( ))), where ( ) is the derivative of the log-gamma function. Given unique survival times T 1 <...<T k the conditional hazard is the product of a continuous and discrete component. The continuous component is H(t) = X i:t i applet T i T i 1 R ] (T i 1 )+ + t T j R ] (T j )+, where R ] (t) =#{i : T i < t}. Thediscretecomponentis Y j:t j applet r j + r j + d j + Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
17 ncomplete Records Parameter Estimation On the assumption that censoring is uninformative, the joint density is Z f (t)p(t c t)g(y; t t t)dt t c where t c = t \ [0, c] andy [t c]=y. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
18 ncomplete Records Parameter Estimation On the assumption that censoring is uninformative, the joint density is Z f (t)p(t c t)g(y; t t t)dt where t c = t \ [0, c] andy [t c]=y. mputation t c mpute survival times T 0 using the conditional survival distribution f T Y [t c], t c, T > c; ˆu, ˆ / f (T ; ˆ )g(y [t c]; T t c; ˆu)1[T > c] The log-likelihood component associated with the imputed, uncensored record is given by log g(y; T 0 t c; )+logf (T 0 ; ) so parameter estimation after imputation is straightforward. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
19 A worked example : cirrhosis study Section 5 A worked example : cirrhosis study Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
20 A worked example : cirrhosis study Motivating Example: Continued Figure 1 : Smoothed Mean Curve for Prothrombin ndex Aligned by Time until Failure Control Prednisone Prothrombin ndex Time Until Failure Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
21 A worked example : cirrhosis study Motivating Example: Continued We suggest the following model based on Figure 1, E(Z i (s) T )= + āi (s) + 0T i + 1s + s log(s + ) cov(z i (s), Z j (s 0 ) T )= 1 2 ijk 1(s, s 0 ; )+ 2 2 ij ij ss 0 s s K 1(s, s 0 0 ; )=exp Table 3 : Coe cients for revival model Censored records Uncensored records Covariate Coef. S.E. Ratio Coef. S.E. Ratio T? Null Treatment Control Prednizone Survival (T ) Revival (s) log(s + ) Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
22 A worked example : cirrhosis study Motivating Example: Continued Table 4 : Variance components Censored records Uncensored records Coe cient S.E. Coe cient S.E. AR Patient White Noise Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
23 A worked example : cirrhosis study E ect of prothrombin on prognosis Over a period of 5 years and one month following recruitment, patient u had eight appointments with prothrombin values as follows: t u (days) Y u[t u] This is the record for patient 402 who was assigned to prednisone and was subsequently censored at 2661 days. The log density at y u is a quadratic form depending on t only through µ. h(t, y u)=const (y u µ) 0 1 (y u µ)/2 Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
24 A worked example : cirrhosis study E ect of prothrombin on prognosis This estimated factor is shown in Fig. 2a for three versions of the record in which the final prothrombin value is 59, 69 or 79. Figure 2 : Three versions of the record for patient 402: log modification factors for the predictive survival density (left panel) and hazard functions (right panel). Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
25 Summary A worked example : cirrhosis study The health sequence is regarded as a random process in its own right, not as a time-dependent covariate governing survival. To a substantial extent, the model for survival time is decoupled from the revival model for the behaviour of the health sequence in reverse time. Realignment implies that value Y i (0) at recruitment must not be treated as a covariate, but as an integral part of the response sequence. f they were available, values prior to recruitment could also be used. The definition of a treatment e ect is not the usual one because the natural way to compare the records for two individuals is not at a fixed time following recruitment, but at a fixed revival time. The treatment value need not be constant in revival time. The predictive value of a partial health sequence for subsequent survival emerges naturally from the joint survival-revival distribution. n particular, the conditional hazard given the finite sequence of earlier values is typically not constant during the subsequent inter-appointment period. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
26 Summary A worked example : cirrhosis study Records cannot be aligned until the patient dies, which means that the revival process is not observable component-wise until T is known. As a result, the likelihood analysis for incomplete records is technically more complicated. The omission of incomplete records from the revival likelihood does not lead to bias in estimation, but it does lead to ine ciency, which could be substantial if the majority of records are incomplete. The principal assumption, that appointment dates be uninformative for subsequent survival, does not a ect likelihood calculations, but it does a ect prognosis calculations for individual patients. For that reason, it is advisable to label all appointments as scheduled or unscheduled. Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
27 A worked example : cirrhosis study Thank you Walter Dempsey (University of Michigan) Survival models and health sequences July 27, / 24
Survival models and health sequences
Survival models and health sequences Peter McCullagh July 3, 2014 Abstract Survival studies often generate not only a survival time for each patient but also a sequence of health measurements at annual
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 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 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 informationSTAT 526 Spring Final Exam. Thursday May 5, 2011
STAT 526 Spring 2011 Final Exam Thursday May 5, 2011 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
More informationSimple techniques for comparing survival functions with interval-censored data
Simple techniques for comparing survival functions with interval-censored data Jinheum Kim, joint with Chung Mo Nam jinhkim@suwon.ac.kr Department of Applied Statistics University of Suwon Comparing survival
More informationAnalysis of Time-to-Event Data: Chapter 6 - Regression diagnostics
Analysis of Time-to-Event Data: Chapter 6 - Regression diagnostics Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/25 Residuals for the
More informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationSample size re-estimation in clinical trials. Dealing with those unknowns. Chris Jennison. University of Kyoto, January 2018
Sample Size Re-estimation in Clinical Trials: Dealing with those unknowns Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj University of Kyoto,
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 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 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 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 informationBIOS 2083: Linear Models
BIOS 2083: Linear Models Abdus S Wahed September 2, 2009 Chapter 0 2 Chapter 1 Introduction to linear models 1.1 Linear Models: Definition and Examples Example 1.1.1. Estimating the mean of a N(μ, σ 2
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 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 informationSTAT 525 Fall Final exam. Tuesday December 14, 2010
STAT 525 Fall 2010 Final exam Tuesday December 14, 2010 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
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 informationBIOS 2083 Linear Models c Abdus S. Wahed
Chapter 5 206 Chapter 6 General Linear Model: Statistical Inference 6.1 Introduction So far we have discussed formulation of linear models (Chapter 1), estimability of parameters in a linear model (Chapter
More informationLecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL
Lecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL The Cox PH model: λ(t Z) = λ 0 (t) exp(β Z). How do we estimate the survival probability, S z (t) = S(t Z) = P (T > t Z), for an individual with covariates
More informationThe Ef ciency of Simple and Countermatched Nested Case-control Sampling
Published by Blackwell Publishers Ltd, 108 Cowley Road, Oxford OX4 1JF, UK and 350 Main Street, Malden, MA 02148, USA Vol 26: 493±509, 1999 The Ef ciency of Simple and Countermatched Nested Case-control
More information1.4 Properties of the autocovariance for stationary time-series
1.4 Properties of the autocovariance for stationary time-series In general, for a stationary time-series, (i) The variance is given by (0) = E((X t µ) 2 ) 0. (ii) (h) apple (0) for all h 2 Z. ThisfollowsbyCauchy-Schwarzas
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
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 9. Statistics Survival Analysis. Presented February 23, Dan Gillen Department of Statistics University of California, Irvine
Statistics 255 - Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 9.1 Survival analysis involves subjects moving through time Hazard may
More informationWritten Exam (2 hours)
M. Müller Applied Analysis of Variance and Experimental Design Summer 2015 Written Exam (2 hours) General remarks: Open book exam. Switch off your mobile phone! Do not stay too long on a part where you
More informationβ j = coefficient of x j in the model; β = ( β1, β2,
Regression Modeling of Survival Time Data Why regression models? Groups similar except for the treatment under study use the nonparametric methods discussed earlier. Groups differ in variables (covariates)
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 informationDynamic Prediction of Disease Progression Using Longitudinal Biomarker Data
Dynamic Prediction of Disease Progression Using Longitudinal Biomarker Data Xuelin Huang Department of Biostatistics M. D. Anderson Cancer Center The University of Texas Joint Work with Jing Ning, Sangbum
More informationThe coxvc_1-1-1 package
Appendix A The coxvc_1-1-1 package A.1 Introduction The coxvc_1-1-1 package is a set of functions for survival analysis that run under R2.1.1 [81]. This package contains a set of routines to fit Cox models
More informationLocal Likelihood Bayesian Cluster Modeling for small area health data. Andrew Lawson Arnold School of Public Health University of South Carolina
Local Likelihood Bayesian Cluster Modeling for small area health data Andrew Lawson Arnold School of Public Health University of South Carolina Local Likelihood Bayesian Cluster Modelling for Small Area
More informationExtensions of Cox Model for Non-Proportional Hazards Purpose
PhUSE Annual Conference 2013 Paper SP07 Extensions of Cox Model for Non-Proportional Hazards Purpose Author: Jadwiga Borucka PAREXEL, Warsaw, Poland Brussels 13 th - 16 th October 2013 Presentation Plan
More informationST 732, Midterm Solutions Spring 2019
ST 732, Midterm Solutions Spring 2019 Please sign the following pledge certifying that the work on this test is your own: I have neither given nor received aid on this test. Signature: Printed Name: There
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 informationCausal Hazard Ratio Estimation By Instrumental Variables or Principal Stratification. Todd MacKenzie, PhD
Causal Hazard Ratio Estimation By Instrumental Variables or Principal Stratification Todd MacKenzie, PhD Collaborators A. James O Malley Tor Tosteson Therese Stukel 2 Overview 1. Instrumental variable
More informationGeneral Linear Model: Statistical Inference
Chapter 6 General Linear Model: Statistical Inference 6.1 Introduction So far we have discussed formulation of linear models (Chapter 1), estimability of parameters in a linear model (Chapter 4), least
More informationAn application of the GAM-PCA-VAR model to respiratory disease and air pollution data
An application of the GAM-PCA-VAR model to respiratory disease and air pollution data Márton Ispány 1 Faculty of Informatics, University of Debrecen Hungary Joint work with Juliana Bottoni de Souza, Valdério
More informationEstimating the Mean Response of Treatment Duration Regimes in an Observational Study. Anastasios A. Tsiatis.
Estimating the Mean Response of Treatment Duration Regimes in an Observational Study Anastasios A. Tsiatis http://www.stat.ncsu.edu/ tsiatis/ Introduction to Dynamic Treatment Regimes 1 Outline Description
More informationBIOSTATISTICAL METHODS
BIOSTATISTICAL METHODS FOR TRANSLATIONAL & CLINICAL RESEARCH Cross-over Designs #: DESIGNING CLINICAL RESEARCH The subtraction of measurements from the same subject will mostly cancel or minimize effects
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 informationAnalyzing Pilot Studies with Missing Observations
Analyzing Pilot Studies with Missing Observations Monnie McGee mmcgee@smu.edu. Department of Statistical Science Southern Methodist University, Dallas, Texas Co-authored with N. Bergasa (SUNY Downstate
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 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 informationI i=1 1 I(J 1) j=1 (Y ij Ȳi ) 2. j=1 (Y j Ȳ )2 ] = 2n( is the two-sample t-test statistic.
Serik Sagitov, Chalmers and GU, February, 08 Solutions chapter Matlab commands: x = data matrix boxplot(x) anova(x) anova(x) Problem.3 Consider one-way ANOVA test statistic For I = and = n, put F = MS
More informationleebounds: Lee s (2009) treatment effects bounds for non-random sample selection for Stata
leebounds: Lee s (2009) treatment effects bounds for non-random sample selection for Stata Harald Tauchmann (RWI & CINCH) Rheinisch-Westfälisches Institut für Wirtschaftsforschung (RWI) & CINCH Health
More informationRank preserving Structural Nested Distribution Model (RPSNDM) for Continuous
Rank preserving Structural Nested Distribution Model (RPSNDM) for Continuous Y : X M Y a=0 = Y a a m = Y a cum (a) : Y a = Y a=0 + cum (a) an unknown parameter. = 0, Y a = Y a=0 = Y for all subjects Rank
More informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
More informationCHL 5225 H Crossover Trials. CHL 5225 H Crossover Trials
CHL 55 H Crossover Trials The Two-sequence, Two-Treatment, Two-period Crossover Trial Definition A trial in which patients are randomly allocated to one of two sequences of treatments (either 1 then, or
More informationTime Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley
Time Series Models and Inference James L. Powell Department of Economics University of California, Berkeley Overview In contrast to the classical linear regression model, in which the components of the
More informationCensoring mechanisms
Censoring mechanisms Patrick Breheny September 3 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Fixed vs. random censoring In the previous lecture, we derived the contribution to the likelihood
More 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 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 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 informationA Generalized Global Rank Test for Multiple, Possibly Censored, Outcomes
A Generalized Global Rank Test for Multiple, Possibly Censored, Outcomes Ritesh Ramchandani Harvard School of Public Health August 5, 2014 Ritesh Ramchandani (HSPH) Global Rank Test for Multiple Outcomes
More informationGroup Sequential Tests for Delayed Responses. Christopher Jennison. Lisa Hampson. Workshop on Special Topics on Sequential Methodology
Group Sequential Tests for Delayed Responses Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Lisa Hampson Department of Mathematics and Statistics,
More informationChapter 7: Hypothesis testing
Chapter 7: Hypothesis testing Hypothesis testing is typically done based on the cumulative hazard function. Here we ll use the Nelson-Aalen estimate of the cumulative hazard. The survival function is used
More informationA Framework for Modeling Positive Class Expansion with Single Snapshot
A Framework for Modeling Positive Class Expansion with Single Snapshot Yang Yu and Zhi-Hua Zhou LAMDA Group National Key Laboratory for Novel Software Technology Nanjing University, China Motivating task
More informationCausal Inference. Miguel A. Hernán, James M. Robins. May 19, 2017
Causal Inference Miguel A. Hernán, James M. Robins May 19, 2017 ii Causal Inference Part III Causal inference from complex longitudinal data Chapter 19 TIME-VARYING TREATMENTS So far this book has dealt
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 informationStock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis
Stock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis Michael P. Babington and Javier Cano-Urbina August 31, 2018 Abstract Duration data obtained from a given stock of individuals
More informationA. Motivation To motivate the analysis of variance framework, we consider the following example.
9.07 ntroduction to Statistics for Brain and Cognitive Sciences Emery N. Brown Lecture 14: Analysis of Variance. Objectives Understand analysis of variance as a special case of the linear model. Understand
More informationResiduals and model diagnostics
Residuals and model diagnostics Patrick Breheny November 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/42 Introduction Residuals Many assumptions go into regression models, and the Cox proportional
More informationANOVA TESTING 4STEPS. 1. State the hypothesis. : H 0 : µ 1 =
Introduction to Statistics in Psychology PSY 201 Professor Greg Francis Lecture 35 ANalysis Of VAriance Ignoring (some) variability TESTING 4STEPS 1. State the hypothesis. : H 0 : µ 1 = µ 2 =... = µ K,
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 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 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 informationPANEL DATA RANDOM AND FIXED EFFECTS MODEL. Professor Menelaos Karanasos. December Panel Data (Institute) PANEL DATA December / 1
PANEL DATA RANDOM AND FIXED EFFECTS MODEL Professor Menelaos Karanasos December 2011 PANEL DATA Notation y it is the value of the dependent variable for cross-section unit i at time t where i = 1,...,
More informationEconometrics Midterm Examination Answers
Econometrics Midterm Examination Answers March 4, 204. Question (35 points) Answer the following short questions. (i) De ne what is an unbiased estimator. Show that X is an unbiased estimator for E(X i
More informationQuantitative Techniques - Lecture 8: Estimation
Quantitative Techniques - Lecture 8: Estimation Key words: Estimation, hypothesis testing, bias, e ciency, least squares Hypothesis testing when the population variance is not known roperties of estimates
More information4. Comparison of Two (K) Samples
4. Comparison of Two (K) Samples K=2 Problem: compare the survival distributions between two groups. E: comparing treatments on patients with a particular disease. Z: Treatment indicator, i.e. Z = 1 for
More informationAkaike criterion: Kullback-Leibler discrepancy
Model choice. Akaike s criterion Akaike criterion: Kullback-Leibler discrepancy Given a family of probability densities {f ( ; ), 2 }, Kullback-Leibler s index of f ( ; ) relativetof ( ; ) is Z ( ) =E
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 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 informationSSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that need work
SSUI: Presentation Hints 1 Comparing Marginal and Random Eects (Frailty) Models Terry M. Therneau Mayo Clinic April 1998 SSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that
More informationLecture 6: Single-classification multivariate ANOVA (k-group( MANOVA)
Lecture 6: Single-classification multivariate ANOVA (k-group( MANOVA) Rationale and MANOVA test statistics underlying principles MANOVA assumptions Univariate ANOVA Planned and unplanned Multivariate ANOVA
More informationPublished online: 10 Apr 2012.
This article was downloaded by: Columbia University] On: 23 March 215, At: 12:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationVIII. ANCOVA. A. Introduction
VIII. ANCOVA A. Introduction In most experiments and observational studies, additional information on each experimental unit is available, information besides the factors under direct control or of interest.
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 informationDESAIN EKSPERIMEN BLOCKING FACTORS. Semester Genap 2017/2018 Jurusan Teknik Industri Universitas Brawijaya
DESAIN EKSPERIMEN BLOCKING FACTORS Semester Genap Jurusan Teknik Industri Universitas Brawijaya Outline The Randomized Complete Block Design The Latin Square Design The Graeco-Latin Square Design Balanced
More information18 Bivariate normal distribution I
8 Bivariate normal distribution I 8 Example Imagine firing arrows at a target Hopefully they will fall close to the target centre As we fire more arrows we find a high density near the centre and fewer
More informationWe begin by thinking about population relationships.
Conditional Expectation Function (CEF) We begin by thinking about population relationships. CEF Decomposition Theorem: Given some outcome Y i and some covariates X i there is always a decomposition where
More informationOptimization. The value x is called a maximizer of f and is written argmax X f. g(λx + (1 λ)y) < λg(x) + (1 λ)g(y) 0 < λ < 1; x, y X.
Optimization Background: Problem: given a function f(x) defined on X, find x such that f(x ) f(x) for all x X. The value x is called a maximizer of f and is written argmax X f. In general, argmax X f may
More informationMIXED MODELS THE GENERAL MIXED MODEL
MIXED MODELS This chapter introduces best linear unbiased prediction (BLUP), a general method for predicting random effects, while Chapter 27 is concerned with the estimation of variances by restricted
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 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 informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
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 informationGeneralizing the MCPMod methodology beyond normal, independent data
Generalizing the MCPMod methodology beyond normal, independent data José Pinheiro Joint work with Frank Bretz and Björn Bornkamp Novartis AG Trends and Innovations in Clinical Trial Statistics Conference
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 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 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 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 informationWhat is Experimental Design?
One Factor ANOVA What is Experimental Design? A designed experiment is a test in which purposeful changes are made to the input variables (x) so that we may observe and identify the reasons for change
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 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 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 informationPh.D. Preliminary Examination Statistics June 2, 2014
Ph.D. Preliminary Examination Statistics June, 04 NOTES:. The exam is worth 00 points.. Partial credit may be given for partial answers if possible.. There are 5 pages in this exam paper. I have neither
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 informationDepartment of Statistical Science FIRST YEAR EXAM - SPRING 2017
Department of Statistical Science Duke University FIRST YEAR EXAM - SPRING 017 Monday May 8th 017, 9:00 AM 1:00 PM NOTES: PLEASE READ CAREFULLY BEFORE BEGINNING EXAM! 1. Do not write solutions on the exam;
More informationStatistics and Probability Letters. Using randomization tests to preserve type I error with response adaptive and covariate adaptive randomization
Statistics and Probability Letters ( ) Contents lists available at ScienceDirect Statistics and Probability Letters journal homepage: wwwelseviercom/locate/stapro Using randomization tests to preserve
More information