Lecture 8 Stat D. Gillen


 Rose Johnston
 8 months ago
 Views:
Transcription
1 Statistics Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 8.1
2 Example of two ways to stratify Suppose a confounder C has 3 levels on which we would like to stratify when comparing λ(t E) and λ(t Ē) where E is an indicator of exposure" Dummy variable stratification: Consider x E = { 1 E 0 Ē c 2 = { 1 level 2 0 else c 3 = { 1 level 3 0 else Now suppose we fit the following proportional hazards model: λ(t) = λ 0 (t)e β E x E +β 2 c 2 +β 3 c 3 8.2
3 Example of two ways to stratify From this, we have Level 1 of C: exposed : λ(t) = unexposed : λ(t) = HR= Level 2 of C: exposed : λ(t) = unexposed : λ(t) = HR= Level 3 of C: exposed : λ(t) = unexposed : λ(t) = HR= *Key point: Stratumtostratum differences in log λ(t) for unexposed have arbitrarily different magnitudes, but they are parallel (). 8.3
4 Example of two ways to stratify True" stratification: Alternatively we could also consider the following proportional hazards model: λ(t) = λ 0i (t)e β E x E, i = 1, 2, 3 In this case we have different baseline hazards for each stratum Note: This is the implicit model for the stratified logrank test 8.4
5 Example of two ways to stratify From the true stratification" model, we have exposed : λ(t) = Level 1 of C: unexposed : λ(t) = exposed : λ(t) = Level 2 of C: unexposed : λ(t) = HR= HR= Level 3 of C: exposed : λ(t) = unexposed : λ(t) = HR= *Key point: Stratumtostratum differences in log λ(t) for unexposed are completely arbitrary (ie. no s). 8.5
6 Example of two ways to stratify Graphically, we can think of the difference as follows: Dummy Variable λ(t) = λ 0 (t)e β E x E +β 2 c 2 +β 3 c 3 True" Stratification λ(t) = λ 0i (t)e β E x E 8.6
7 The stratified proportional hazards regression model The general stratified proportional hazards regression model is: λ i (t x i, stratum g) = λ g (t) exp(β T x i ) for strata g = 1,..., G. A semiparametric model with a nonparametric baseline hazard function within each stratum Alternatively, it can be useful because the hazard can difficult to model as a function of complicated set of adjustment covariates 8.7
8 Considerations Under dummy variable stratification, the proportional stratumtostratum hazards may not be correct The confounder (adjustment/stratification variable) may be inadequately controlled, and may still confound Proportionality can be checked by examining stratumspecific log Λ(t) plots (coming up) True stratification is a more thorough adjustment as long as observations within each level are homogenous If the confounder can be measured continuously, control for confounding might be achieved with the continuous covariate adjustment 8.8
9 Considerations If the confounder is controlled using true stratification, there is no way to estimate a summary relative risk comparing two levels of the confounder True stratification generally requires more data to obtain the same precision in coefficient estimates (biasvariance tradeoff) 8.9
10 Example: 6MP Data Example: Recall that in the 6MP leukemia data (Section 1.2 in K & M): 21 patient pairs Interest on the effect of the treatment 6MP on survival, and Matching was on remission status (remstat) and hospital (pairid) Recall also, that matched pairs are like many strata Goal: Estimate the relative risk for 6MP versus placebo, adjusting for hospital and remission status 8.10
11 Example: 6MP Data We can conduct the stratified analysis accounting for matching by including a strata() term > fit < coxph( Surv( time, irelapse ) ~ sixmp + strata(pairid), data=sixmplong ) > summary( fit ) Call: coxph(formula = Surv(time, irelapse) ~ sixmp + strata(pairid), data = sixmplong) n= 42 coef exp(coef) se(coef) z Pr(> z ) sixmp **  Signif. codes: 0 Ô***Õ Ô**Õ 0.01 Ô*Õ 0.05 Ô.Õ 0.1 Ô Õ 1 exp(coef) exp(coef) lower.95 upper.95 sixmp Rsquare= (max possible= 0.5 ) Likelihood ratio test= 11.9 on 1 df, p= Wald test = 8.26 on 1 df, p= Score (logrank) test = 10.7 on 1 df, p=
12 Example: 6MP Data Conclusion: the withinstratum estimated relative risk for 6MP versus placebo is 0.17, 95% CI [0.05,.57], indicating that the risk of relapse is about 6 times higher in the placebo than in the 6MP group, after adjusting for any effects due to hospital and remission status. 8.12
13 Example: 6MP Data For comparison, let s consider the unmatched (incorrect) analysis > fit2 < coxph( Surv( time, irelapse ) ~ sixmp, data=sixmplong ) > summary( fit2 ) Call: coxph(formula = Surv(time, irelapse) ~ sixmp, data = sixmplong) n= 42 coef exp(coef) se(coef) z Pr(> z ) sixmp ***  Signif. codes: 0 Ô***Õ Ô**Õ 0.01 Ô*Õ 0.05 Ô.Õ 0.1 Ô Õ 1 exp(coef) exp(coef) lower.95 upper.95 sixmp Rsquare= (max possible= ) Likelihood ratio test= 16.4 on 1 df, p=5.26e05 Wald test = 14.5 on 1 df, p= Score (logrank) test = 17.2 on 1 df, p=3.28e
14 Example: 6MP Data (Incorrect) Conclusion: the estimated relative risk for 6MP versus placebo is 0.208, indicating that the risk of relapse is about 4.8 times higher in the placebo than in the 6MP group. Incorrect because not at all clear what the reference population is 8.14
15 Partial Likelihood for the Stratified PH Partial Likelihood for the Stratified PH Recall that the partial likelihood for the ordinary proportional hazards model is ( ) exp(β T x (j) ) L P = i R failure times j (j) exp(β T x i ) L P compares the covariate x (j) of the failed subject to the covariates x i in the risk set R (j) at time t (j) The stratified partial likelihood combines partial likelihoods L Pg across strata g = 1,..., G: L P = G g=1 L Pg 8.15
16 Partial Likelihood for the Stratified PH Partial Likelihood for the Stratified PH Alternatively, who is in the effective risk set at time t (j)? A: The subject failing at t (j) and other subjects not yet failed, but also in the same stratum (i.e., in stratum g) Suppose t (j) is a failure time for a subject in stratum g Define R g (j) to be the set of subjects in stratum g and at risk for failure at t (j) Then the stratified partial likelihood (in the absence of ties) is ( ) exp(β T x (j) ) L P = i R failure times g exp(β T x i ) j (j) 8.16
17 Assumption Example: Bone marrow transplant data Ex: Bone marrow transplant data (Section 1.3 in K & M) Question: For leukemia patients, does highrisk versus lowrisk subclassification have any prognosticative meaning or value, after adjusting for other baseline clinical information? Data: time (in days) of diseasefree survival waittime: the waiting time to transplant (days) fab: FAB grade of disease (4 or 5, coded 0 or 1) i.mtx: MTX, a graphvshost disease prophylactic treatment agep: age of the patient aged: age of the donor g: Indicator of whether patient has acute myeloid leukemia (AML) Prior model building already performed
18 Assumption Example: Bone marrow transplant data Ex: Bone marrow transplant data (Section 1.3 in K & M) Question: Adjusting for other covariates, does the effect of imtx follow a proportional hazards model? Method: If the proportional hazards holds, the log cumulative hazards for the two MTX groups, adjusting for other covariates x i, should be roughly parallel 8.18
19 Assumption Example: Bone marrow transplant data Consider strata defined by imtx = 0 or 1 and model λ i (t x i, imtx i ) = λ imtxi (t) exp(β T x i ) where x i contains the other covariates Recall from Lecture 4: If proportional hazards holds for imtx, then log Λ imtx=1 (t) = log(φ) + log Λ imtx=0 (t) where φ = exp(β T x) So, if the baseline log cumulative hazards are not parallel this is evidence of nonproportional hazards 8.19
20 Assumption Example: Bone marrow transplant data R will estimate the withinstratum baseline hazard when a stratified analysis is fit > ## > ##### > ##### Read in BMT data > ##### > ## > bmt < read.csv( " STAT255/Data/bmt.csv" ) > ## > ##### Fit the stratified cox model and plot > ##### logcum hazard estimates > ## > fit < coxph( Surv(tnodis,inodis) ~ agep*aged + g + strata(imtx), data=bmt ) > plot( survfit(fit), fun="cloglog", + xlab="time from Transplantation (logarithmic scale)", + ylab="logcumulative Hazard Function" ) 8.20
21 Assumption Example: Bone marrow transplant data Log Cumulative Hazard Function Time from Transplantation 8.21
22 Assumption Example: Bone marrow transplant data Conclusion: Hard to tell but not strong evidence of nonproportional hazards This is a good look at gross departures, but it is far from a formal test...more later! 8.22
23 β k is the change in the loghazard function, i.e. the logrelative risk, associated with a unit difference in x k, for subjects in the same stratum and holding other covariates constant. Therefore, the contrast expressed by β k is adjusted for all other covariates in the model as well as any variables contained in the stratification variable. The key is that the effect of covariate x ik is the same in each stratum. If it is not, ie. there is an interaction between the stratification variable and the predictor of interest, stratumspecific estimates should be presented 8.23
24 With p covariates in the model, plots based on a model stratified on values of new covariate x i,p+1 can help assess the validity of the proportional hazards for x i,p+1. The stratified partial likelihood model works by only allowing subjects from the same stratum to be in the same risk set. 8.24
Chapter 4 Regression Models
23.August 2010 Chapter 4 Regression Models The target variable T denotes failure time We let x = (x (1),..., x (m) ) represent a vector of available covariates. Also called regression variables, regressors,
More informationMultivariable Fractional Polynomials
Multivariable Fractional Polynomials Axel Benner September 7, 2015 Contents 1 Introduction 1 2 Inventory of functions 1 3 Usage in R 2 3.1 Model selection........................................ 3 4 Example
More informationLecture 11. Interval Censored and. DiscreteTime Data. Statistics Survival Analysis. Presented March 3, 2016
Statistics 255  Survival Analysis Presented March 3, 2016 Motivating Dan Gillen Department of Statistics University of California, Irvine 11.1 First question: Are the data truly discrete? : Number of
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 informationExtensions of Cox Model for NonProportional Hazards Purpose
PhUSE 2013 Paper SP07 Extensions of Cox Model for NonProportional Hazards Purpose Jadwiga Borucka, PAREXEL, Warsaw, Poland ABSTRACT Cox proportional hazard model is one of the most common methods used
More informationTied survival times; estimation of survival probabilities
Tied survival times; estimation of survival probabilities Patrick Breheny November 5 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Tied survival times Introduction Breslow approximation
More informationIgnoring the matching variables in cohort studies  when is it valid, and why?
Ignoring the matching variables in cohort studies  when is it valid, and why? Arvid Sjölander Abstract In observational studies of the effect of an exposure on an outcome, the exposureoutcome association
More informationJournal of Statistical Software
JSS Journal of Statistical Software January 2011, Volume 38, Issue 7. http://www.jstatsoft.org/ mstate: An R Package for the Analysis of Competing Risks and MultiState Models Liesbeth C. de Wreede Leiden
More informationLogistic regression model for survival time analysis using timevarying coefficients
Logistic regression model for survival time analysis using timevarying coefficients Accepted in American Journal of Mathematical and Management Sciences, 2016 Kenichi SATOH ksatoh@hiroshimau.ac.jp Research
More informationMachine Learning Linear Classification. Prof. Matteo Matteucci
Machine Learning Linear Classification Prof. Matteo Matteucci Recall from the first lecture 2 X R p Regression Y R Continuous Output X R p Y {Ω 0, Ω 1,, Ω K } Classification Discrete Output X R p Y (X)
More informationOptimal Treatment Regimes for Survival Endpoints from a Classification Perspective. Anastasios (Butch) Tsiatis and Xiaofei Bai
Optimal Treatment Regimes for Survival Endpoints from a Classification Perspective Anastasios (Butch) Tsiatis and Xiaofei Bai Department of Statistics North Carolina State University 1/35 Optimal Treatment
More informationA new strategy for metaanalysis of continuous covariates in observational studies with IPD. Willi Sauerbrei & Patrick Royston
A new strategy for metaanalysis of continuous covariates in observational studies with IPD Willi Sauerbrei & Patrick Royston Overview Motivation Continuous variables functional form Fractional polynomials
More informationUnderstanding the Cox Regression Models with TimeChange Covariates
Understanding the Cox Regression Models with TimeChange Covariates Mai Zhou University of Kentucky The Cox regression model is a cornerstone of modern survival analysis and is widely used in many other
More informationStat 642, Lecture notes for 04/12/05 96
Stat 642, Lecture notes for 04/12/05 96 HosmerLemeshow Statistic The HosmerLemeshow Statistic is another measure of lack of fit. Hosmer and Lemeshow recommend partitioning the observations into 10 equal
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 informationThe Ef ciency of Simple and Countermatched Nested Casecontrol 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 Casecontrol
More informationLoglinearity for Cox s regression model. Thesis for the Degree Master of Science
Loglinearity 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 informationUnobserved Heterogeneity
Unobserved Heterogeneity Germán Rodríguez grodri@princeton.edu Spring, 21. Revised Spring 25 This unit considers survival models with a random effect representing unobserved heterogeneity of frailty, a
More informationAnalysis of competing risks data and simulation of data following predened subdistribution hazards
Analysis of competing risks data and simulation of data following predened subdistribution hazards Bernhard Haller Institut für Medizinische Statistik und Epidemiologie Technische Universität München 27.05.2013
More informationIntegrated likelihoods in survival models for highlystratified
Working Paper Series, N. 1, January 2014 Integrated likelihoods in survival models for highlystratified censored data Giuliana Cortese Department of Statistical Sciences University of Padua Italy Nicola
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 informationMüller: Goodnessoffit criteria for survival data
Müller: Goodnessoffit criteria for survival data Sonderforschungsbereich 386, Paper 382 (2004) Online unter: http://epub.ub.unimuenchen.de/ Projektpartner Goodness of fit criteria for survival data
More informationPropensity Score Methods for Causal Inference
John Pura BIOS790 October 2, 2015 Causal inference Philosophical problem, statistical solution Important in various disciplines (e.g. Koch s postulates, Bradford Hill criteria, Granger causality) Good
More informationPrentice Hall Stats: Modeling the World 2004 (Bock) Correlated to: National Advanced Placement (AP) Statistics Course Outline (Grades 912)
National Advanced Placement (AP) Statistics Course Outline (Grades 912) Following is an outline of the major topics covered by the AP Statistics Examination. The ordering here is intended to define the
More informationStatistical Methods in Clinical Trials Categorical Data
Statistical Methods in Clinical Trials Categorical Data Types of Data quantitative Continuous Blood pressure Time to event Categorical sex qualitative Discrete No of relapses Ordered Categorical Pain level
More informationStandardization methods have been used in epidemiology. Marginal Structural Models as a Tool for Standardization ORIGINAL ARTICLE
ORIGINAL ARTICLE Marginal Structural Models as a Tool for Standardization Tosiya Sato and Yutaka Matsuyama Abstract: In this article, we show the general relation between standardization methods and marginal
More informationLongitudinal Modeling with Logistic Regression
Newsom 1 Longitudinal Modeling with Logistic Regression Longitudinal designs involve repeated measurements of the same individuals over time There are two general classes of analyses that correspond to
More informationCorrelation and regression
1 Correlation and regression Yongjua Laosiritaworn Introductory on Field Epidemiology 6 July 2015, Thailand Data 2 Illustrative data (Doll, 1955) 3 Scatter plot 4 Doll, 1955 5 6 Correlation coefficient,
More informationSTAT 331. Accelerated Failure Time Models. Previously, we have focused on multiplicative intensity models, where
STAT 331 Accelerated Failure Time Models Previously, we have focused on multiplicative intensity models, where h t z) = h 0 t) g z). These can also be expressed as H t z) = H 0 t) g z) or S t z) = e Ht
More informationSurvival Analysis APTS 2016/17. Ingrid Van Keilegom ORSTAT KU Leuven. Glasgow, August 2125, 2017
Survival Analysis APTS 2016/17 Ingrid Van Keilegom ORSTAT KU Leuven Glasgow, August 2125, 2017 Basic What is Survival analysis? Survival analysis (or duration analysis) is an area of statistics that and
More informationPASS Sample Size Software. Poisson Regression
Chapter 870 Introduction Poisson regression is used when the dependent variable is a count. Following the results of Signorini (99), this procedure calculates power and sample size for testing the hypothesis
More informationMonitoring clinical trial outcomes with delayed response: incorporating pipeline data in group sequential designs. Christopher Jennison
Monitoring clinical trial outcomes with delayed response: incorporating pipeline data in group sequential designs Christopher Jennison Department of Mathematical Sciences, University of Bath http://people.bath.ac.uk/mascj
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 informationTIMEDEPENDENT COVARIATES IN THE COX PROPORTIONALHAZARDS REGRESSION MODEL
Annu. Rev. Public Health. 1999. 20:145 57 Copyright c 1999 by Annual Reviews. All rights reserved TIMEDEPENDENT COVARIATES IN THE COX PROPORTIONALHAZARDS REGRESSION MODEL Lloyd D. Fisher and D. Y. Lin
More informationIntroduction to the rstpm2 package
Introduction to the rstpm2 package Mark Clements Karolinska Institutet Abstract This vignette outlines the methods and provides some examples for linkbased survival models as implemented in the R rstpm2
More informationSTAT 5500/6500 Conditional Logistic Regression for Matched Pairs
STAT 5500/6500 Conditional Logistic Regression for Matched Pairs The data for the tutorial came from support.sas.com, The LOGISTIC Procedure: Conditional Logistic Regression for Matched Pairs Data :: SAS/STAT(R)
More informationGuideline on adjustment for baseline covariates in clinical trials
26 February 2015 EMA/CHMP/295050/2013 Committee for Medicinal Products for Human Use (CHMP) Guideline on adjustment for baseline covariates in clinical trials Draft Agreed by Biostatistics Working Party
More informationAnalytic Methods for Applied Epidemiology: Framework and Contingency Table Analysis
Analytic Methods for Applied Epidemiology: Framework and Contingency Table Analysis 2014 Maternal and Child Health Epidemiology Training PreTraining Webinar: Friday, May 16 24pm Eastern Kristin Rankin,
More informationAttributable Risk Function in the Proportional Hazards Model
UW Biostatistics Working Paper Series 5312005 Attributable Risk Function in the Proportional Hazards Model Ying Qing Chen Fred Hutchinson Cancer Research Center, yqchen@u.washington.edu Chengcheng Hu
More informationBIAS OF MAXIMUMLIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY
BIAS OF MAXIMUMLIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY Ingo Langner 1, Ralf Bender 2, Rebecca LenzTönjes 1, Helmut Küchenhoff 2, Maria Blettner 2 1
More informationE509A: Principle of Biostatistics. (Week 11(2): Introduction to nonparametric. methods ) GY Zou.
E509A: Principle of Biostatistics (Week 11(2): Introduction to nonparametric methods ) GY Zou gzou@robarts.ca Sign test for two dependent samples Ex 12.1 subj 1 2 3 4 5 6 7 8 9 10 baseline 166 135 189
More informationSurvival models and health sequences
Survival models and health sequences Walter Dempsey University of Michigan July 27, 2015 Survival Data Problem Description Survival data is commonplace in medical studies, consisting of failure time information
More informationMODELING THE SUBDISTRIBUTION OF A COMPETING RISK
Statistica Sinica 16(26), 13671385 MODELING THE SUBDISTRIBUTION OF A COMPETING RISK Liuquan Sun 1, Jingxia Liu 2, Jianguo Sun 3 and MeiJie Zhang 2 1 Chinese Academy of Sciences, 2 Medical College of
More informationST495: Survival Analysis: Maximum likelihood
ST495: Survival Analysis: Maximum likelihood Eric B. Laber Department of Statistics, North Carolina State University February 11, 2014 Everything is deception: seeking the minimum of illusion, keeping
More informationMultinomial Logistic Regression Models
Stat 544, Lecture 19 1 Multinomial Logistic Regression Models Polytomous responses. Logistic regression can be extended to handle responses that are polytomous, i.e. taking r>2 categories. (Note: The word
More informationLecture 3. Truncation, lengthbias and prevalence sampling
Lecture 3. Truncation, lengthbias 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 informationGroup Sequential Methods. for Clinical Trials. Christopher Jennison, Dept of Mathematical Sciences, University of Bath, UK
Group Sequential Methods for Clinical Trials Christopher Jennison, Dept of Mathematical Sciences, University of Bath, UK http://www.bath.ac.uk/οmascj PSI, London, 22 October 2003 1 Plan of talk 1. Why
More informationPoisson Regression: Let me count the uses!
Poisson Regression: Let me count the uses! Lisa Kaltenbach November 20, 2008 Poisson Reg may analyze 1. Count data (ex. no. of surgical site infections) 2. Binary data (ex. received vaccine(yes/no) as
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 informationJournal of Statistical Software
JSS Journal of Statistical Software January 2011, Volume 38, Issue 2. http://www.jstatsoft.org/ Analyzing Competing Risk Data Using the R timereg Package Thomas H. Scheike University of Copenhagen MeiJie
More informationAssessment of time varying long term effects of therapies and prognostic factors
Assessment of time varying long term effects of therapies and prognostic factors Dissertation by Anika Buchholz Submitted to Fakultät Statistik, Technische Universität Dortmund in Fulfillment of the Requirements
More informationEffects of multiple interventions
Chapter 28 Effects of multiple interventions James Robins, Miguel Hernan and Uwe Siebert 1. Introduction The purpose of this chapter is (i) to describe some currently available analytical methods for using
More informationQuantifying the Predictive Accuracy of Time to Event Models in the Presence of Competing Risks
Quantifying the Predictive Accuracy of Time to Event Models in the Presence of Competing Risks Rotraut Schoop, Jan Beyersmann, Martin Schumacher & Harald Binder Universität Freiburg i. Br. FDMPreprint
More information8 Nominal and Ordinal Logistic Regression
8 Nominal and Ordinal Logistic Regression 8.1 Introduction If the response variable is categorical, with more then two categories, then there are two options for generalized linear models. One relies on
More informationStatistics 262: Intermediate Biostatistics Regression & Survival Analysis
Statistics 262: Intermediate Biostatistics Regression & Survival Analysis Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Introduction This course is an applied course,
More informationTests for Two Correlated Proportions in a Matched Case Control Design
Chapter 155 Tests for Two Correlated Proportions in a Matched Case Control Design Introduction A 2byM casecontrol study investigates a risk factor relevant to the development of a disease. A population
More informationBinomial Model. Lecture 10: Introduction to Logistic Regression. Logistic Regression. Binomial Distribution. n independent trials
Lecture : Introduction to Logistic Regression Ani Manichaikul amanicha@jhsph.edu 2 May 27 Binomial Model n independent trials (e.g., coin tosses) p = probability of success on each trial (e.g., p =! =
More informationLinear Regression Models P8111
Linear Regression Models P8111 Lecture 25 Jeff Goldsmith April 26, 2016 1 of 37 Today s Lecture Logistic regression / GLMs Model framework Interpretation Estimation 2 of 37 Linear regression Course started
More informationHypothesis Testing, Power, Sample Size and Confidence Intervals (Part 2)
Hypothesis Testing, Power, Sample Size and Confidence Intervals (Part 2) B.H. Robbins Scholars Series June 23, 2010 1 / 29 Outline Ztest χ 2 test Confidence Interval Sample size and power Relative effect
More informationpensim Package Example (Version 1.2.9)
pensim Package Example (Version 1.2.9) Levi Waldron March 13, 2014 Contents 1 Introduction 1 2 Example data 2 3 Nested crossvalidation 2 3.1 Summarization and plotting..................... 3 4 Getting
More informationCorrelation and Linear Regression
Correlation and Linear Regression Correlation: Relationships between Variables So far, nearly all of our discussion of inferential statistics has focused on testing for differences between group means
More informationGeneralized linear models for binary data. A better graphical exploratory data analysis. The simple linear logistic regression model
Stat 3302 (Spring 2017) Peter F. Craigmile Simple linear logistic regression (part 1) [Dobson and Barnett, 2008, Sections 7.1 7.3] Generalized linear models for binary data Beetles doseresponse example
More informationChapter 5: Logistic RegressionI
: Logistic RegressionI Dipankar Bandyopadhyay Department of Biostatistics, Virginia Commonwealth University BIOS 625: Categorical Data & GLM [Acknowledgements to Tim Hanson and Haitao Chu] D. Bandyopadhyay
More informationFRAILTY MODELS FOR MODELLING HETEROGENEITY
FRAILTY MODELS FOR MODELLING HETEROGENEITY By ULVIYA ABDULKARIMOVA, B.Sc. A Thesis Submitted to the School of Graduate Studies in Partial Fulfillment of the Requirements for the Degree Master of Science
More informationNonlinear Regression
Nonlinear Regression James H. Steiger Department of Psychology and Human Development Vanderbilt University James H. Steiger (Vanderbilt University) Nonlinear Regression 1 / 36 Nonlinear Regression 1 Introduction
More informationSTA6938Logistic Regression Model
Dr. Ying Zhang STA6938Logistic Regression Model Topic 2Multiple Logistic Regression Model Outlines:. Model Fitting 2. Statistical Inference for Multiple Logistic Regression Model 3. Interpretation of
More informationThe Accelerated Failure Time Model Under Biased. Sampling
The Accelerated Failure Time Model Under Biased Sampling Micha Mandel and Ya akov Ritov Department of Statistics, The Hebrew University of Jerusalem, Israel July 13, 2009 Abstract Chen (2009, Biometrics)
More informationM3 Symposium: Multilevel Multivariate Survival Models For Analysis of Dyadic Social Interaction
M3 Symposium: Multilevel Multivariate Survival Models For Analysis of Dyadic Social Interaction Mike Stoolmiller: stoolmil@uoregon.edu University of Oregon 5/21/2013 Outline Example Research Questions
More informationPackage CPE. R topics documented: February 19, 2015
Package CPE February 19, 2015 Title Concordance Probability Estimates in Survival Analysis Version 1.4.4 Depends R (>= 2.10.0),survival,rms Author Qianxing Mo, Mithat Gonen and Glenn Heller Maintainer
More informationThe Type I Generalized Half Logistic Survival Model
International Journal of Theoretical and Applied Mathematics 2016; 2(2): 7478 http://www.sciencepublishinggroup.com/j/ijtam doi: 10.11648/j.ijtam.20160202.17 The Type I Generalized Half Logistic urvival
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models
Introduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationCompare Predicted Counts between Groups of Zero Truncated Poisson Regression Model based on Recycled Predictions Method
Compare Predicted Counts between Groups of Zero Truncated Poisson Regression Model based on Recycled Predictions Method Yan Wang 1, Michael Ong 2, Honghu Liu 1,2,3 1 Department of Biostatistics, UCLA School
More informationSigmaplot di Systat Software
Sigmaplot di Systat Software SigmaPlot Has Extensive Statistical Analysis Features SigmaPlot is now bundled with SigmaStat as an easytouse package for complete graphing and data analysis. The statistical
More informationDouble Bootstrap Confidence Interval Estimates with Censored and Truncated Data
Journal of Modern Applied Statistical Methods Volume 13 Issue 2 Article 22 112014 Double Bootstrap Confidence Interval Estimates with Censored and Truncated Data Jayanthi Arasan University Putra Malaysia,
More informationSensitivity study of dosefinding methods
to of dosefinding methods Sarah Zohar 1 John O Quigley 2 1. Inserm, UMRS 717,Biostatistic Department, Hôpital SaintLouis, Paris, France 2. Inserm, Université Paris VI, Paris, France. NY 2009 1 / 21 to
More information9 Generalized Linear Models
9 Generalized Linear Models The Generalized Linear Model (GLM) is a model which has been built to include a wide range of different models you already know, e.g. ANOVA and multiple linear regression models
More informationSAS/STAT 13.2 User s Guide. Introduction to Survey Sampling and Analysis Procedures
SAS/STAT 13.2 User s Guide Introduction to Survey Sampling and Analysis Procedures This document is an individual chapter from SAS/STAT 13.2 User s Guide. The correct bibliographic citation for the complete
More informationarxiv: v1 [stat.me] 15 May 2011
Working Paper Propensity Score Analysis with Matching Weights Liang Li, Ph.D. arxiv:1105.2917v1 [stat.me] 15 May 2011 Associate Staff of Biostatistics Department of Quantitative Health Sciences, Cleveland
More informationUnit 9: Inferences for Proportions and Count Data
Unit 9: Inferences for Proportions and Count Data Statistics 571: Statistical Methods Ramón V. León 12/15/2008 Unit 9  Stat 571  Ramón V. León 1 Large Sample Confidence Interval for Proportion ( pˆ p)
More informationPoisson Regression. The Training Data
The Training Data Poisson Regression Office workers at a large insurance company are randomly assigned to one of 3 computer use training programmes, and their number of calls to IT support during the following
More informationRandomization Tests for Regression Models in Clinical Trials
Randomization Tests for Regression Models in Clinical Trials A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at George Mason University By Parwen
More informationSAS/STAT 13.1 User s Guide. Introduction to Survey Sampling and Analysis Procedures
SAS/STAT 13.1 User s Guide Introduction to Survey Sampling and Analysis Procedures This document is an individual chapter from SAS/STAT 13.1 User s Guide. The correct bibliographic citation for the complete
More informationLecture 8. Poisson models for counts
Lecture 8. Poisson models for counts Jesper Rydén Department of Mathematics, Uppsala University jesper.ryden@math.uu.se Statistical Risk Analysis Spring 2014 Absolute risks The failure intensity λ(t) describes
More informationModel Estimation Example
Ronald H. Heck 1 EDEP 606: Multivariate Methods (S2013) April 7, 2013 Model Estimation Example As we have moved through the course this semester, we have encountered the concept of model estimation. Discussions
More informationConfidence Intervals for the Odds Ratio in Logistic Regression with One Binary X
Chapter 864 Confidence Intervals for the Odds Ratio in Logistic Regression with One Binary X Introduction Logistic regression expresses the relationship between a binary response variable and one or more
More informationMarginal Structural Models and Causal Inference in Epidemiology
Marginal Structural Models and Causal Inference in Epidemiology James M. Robins, 1,2 Miguel Ángel Hernán, 1 and Babette Brumback 2 In observational studies with exposures or treatments that vary over time,
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 22877843 / Online ISSN 23834757 Estimation of Conditional
More informationConfidence Intervals for the Odds Ratio in Logistic Regression with Two Binary X s
Chapter 866 Confidence Intervals for the Odds Ratio in Logistic Regression with Two Binary X s Introduction Logistic regression expresses the relationship between a binary response variable and one or
More informationGeneralized Linear Modeling  Logistic Regression
1 Generalized Linear Modeling  Logistic Regression Binary outcomes The logit and inverse logit interpreting coefficients and odds ratios Maximum likelihood estimation Problem of separation Evaluating
More informationFrailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Mela. P.
Frailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Melanie M. Wall, Bradley P. Carlin November 24, 2014 Outlines of the talk
More informationLarge n normal approximations (Central Limit Theorem). xbar ~ N[mu, sigma 2 / n] (sketch a normal with mean mu and sd = sigma / root(n)).
1. 51, 52, 53 Large n normal approximations (Central Limit Theorem). xbar ~ N[mu, sigma 2 / n] (sketch a normal with mean mu and sd = sigma / root(n)). phat ~ N[p, pq / n] (sketch a normal with mean
More informationPromotion Time Cure Rate Model with Random Effects: an Application to a Multicentre Clinical Trial of Carcinoma
www.srljournal.org Statistics Research Letters (SRL) Volume Issue, May 013 Promotion Time Cure Rate Model with Random Effects: an Application to a Multicentre Clinical Trial of Carcinoma Diego I. Gallardo
More informationAnalysing categorical data using logit models
Analysing categorical data using logit models Graeme Hutcheson, University of Manchester The lecture notes, exercises and data sets associated with this course are available for download from: www.researchtraining.net/manchester
More informationStat/F&W Ecol/Hort 572 Review Points Ané, Spring 2010
1 Linear models Y = Xβ + ɛ with ɛ N (0, σ 2 e) or Y N (Xβ, σ 2 e) where the model matrix X contains the information on predictors and β includes all coefficients (intercept, slope(s) etc.). 1. Number of
More informationTRIM Workshop. Arco van Strien Wildlife statistics Statistics Netherlands (CBS)
TRIM Workshop Arco van Strien Wildlife statistics Statistics Netherlands (CBS) What is TRIM? TRends and Indices for Monitoring data Computer program for the analysis of time series of count data with missing
More informationFractional Polynomials and Model Averaging
Fractional Polynomials and Model Averaging Paul C Lambert Center for Biostatistics and Genetic Epidemiology University of Leicester UK paul.lambert@le.ac.uk Nordic and Baltic Stata Users Group Meeting,
More informationModels for Ordinal Response Data
Models for Ordinal Response Data Robin High Department of Biostatistics Center for Public Health University of Nebraska Medical Center Omaha, Nebraska Recommendations Analyze numerical data with a statistical
More informationModel Selection Procedures
Model Selection Procedures Statistics 135 Autumn 2005 Copyright c 2005 by Mark E. Irwin Model Selection Procedures Consider a regression setting with K potential predictor variables and you wish to explore
More information2/26/2017. PSY 512: Advanced Statistics for Psychological and Behavioral Research 2
PSY 512: Advanced Statistics for Psychological and Behavioral Research 2 When and why do we use logistic regression? Binary Multinomial Theory behind logistic regression Assessing the model Assessing predictors
More informationAP Statistics Unit 6 Note Packet Linear Regression. Scatterplots and Correlation
Scatterplots and Correlation Name Hr A scatterplot shows the relationship between two quantitative variables measured on the same individuals. variable (y) measures an outcome of a study variable (x) may
More information