Predicting Survival Outcomes Based on Compound Covariate Method under Cox Proportional Hazard Models with Microarrays

Similar documents
Survival Prediction Based on Compound Covariate under Cox Proportional Hazard Models

Chapter 3 Supplemental Text Material

Chapter 7. Bounds for weighted sums of Random Variables

β (cf Khan, 2006). In this model, p independent

Introduction to mathematical Statistics

Regression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers

Random variables and sampling theory

Density estimation II

Inference on Parameters of a Shifted and Rescaled Wiener Process Based on Nonidentically Distributed Observations

An Extended Mixture Inverse Gaussian Distribution

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Differential Entropy 吳家麟教授

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Bayesian Classification. CS690L Data Mining: Classification(2) Bayesian Theorem: Basics. Bayesian Theorem. Training dataset. Naïve Bayes Classifier

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

Z = = = = X np n. n n. npq. npq pq

MATH2999 Directed Studies in Mathematics Matrix Theory and Its Applications

Section 7.2 Two-way ANOVA with random effect(s)

Stats & Summary

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

Multivariate Transformation of Variables and Maximum Likelihood Estimation

Chapter 2 Intro to Math Techniques for Quantum Mechanics

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

Chapter Two. An Introduction to Regression ( )

Bayes Estimator for Exponential Distribution with Extension of Jeffery Prior Information

Linear Regression with One Regressor

Almost Unbiased Estimation of the Poisson Regression Model

LINEAR REGRESSION ANALYSIS

Bayes (Naïve or not) Classifiers: Generative Approach

Objectives of Multiple Regression

MTH 146 Class 7 Notes

VOL. 3, NO. 11, November 2013 ISSN ARPN Journal of Science and Technology All rights reserved.

Simple Linear Regression

Special Instructions / Useful Data

ENGI 3423 Simple Linear Regression Page 12-01

CLASS NOTES. for. PBAF 528: Quantitative Methods II SPRING Instructor: Jean Swanson. Daniel J. Evans School of Public Affairs

Lecture 8: Linear Regression

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

The Probability Distribution Model of Wind Speed over East Malaysia

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

Midterm Exam 1, section 1 (Solution) Thursday, February hour, 15 minutes

3. Basic Concepts: Consequences and Properties

Midterm Exam 1, section 2 (Solution) Thursday, February hour, 15 minutes

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Chapter 4: Distributions

Kernel-based Methods and Support Vector Machines

On Several Inequalities Deduced Using a Power Series Approach

CS 2750 Machine Learning. Lecture 8. Linear regression. CS 2750 Machine Learning. Linear regression. is a linear combination of input components x

BAYESIAN INFERENCES FOR TWO PARAMETER WEIBULL DISTRIBUTION

Simple Linear Regression - Scalar Form

Design of Bayesian MDS Sampling Plan Based on the Process Capability Index

CHAPTER VI Statistical Analysis of Experimental Data

Faculty Research Interest Seminar Department of Biostatistics, GSPH University of Pittsburgh. Gong Tang Feb. 18, 2005

Chapter Gauss-Seidel Method

Statistics: Unlocking the Power of Data Lock 5

Generative classification models

The number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter

Lower Bounds of the Kirchhoff and Degree Kirchhoff Indices

COV. Violation of constant variance of ε i s but they are still independent. The error term (ε) is said to be heteroscedastic.

Solutions Manual for Polymer Science and Technology Third Edition

ESTIMATION OF PARAMETERS OF THE MARSHALL-OLKIN EXTENDED LOG-LOGISTIC DISTRIBUTION FROM PROGRESSIVELY CENSORED SAMPLES

Reaction Time VS. Drug Percentage Subject Amount of Drug Times % Reaction Time in Seconds 1 Mary John Carl Sara William 5 4

Accelerated Life Test Sampling Plans under Progressive Type II Interval Censoring with Random Removals

STRONG CONSISTENCY FOR SIMPLE LINEAR EV MODEL WITH v/ -MIXING

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]

Testing for multivariate normality of disturbances in the multivariate linear regression model Yan Su a, Shao-Yue Kang b

1 Solution to Problem 6.40

Study of Correlation using Bayes Approach under bivariate Distributions

A Technique for Constructing Odd-order Magic Squares Using Basic Latin Squares

1. The weight of six Golden Retrievers is 66, 61, 70, 67, 92 and 66 pounds. The weight of six Labrador Retrievers is 54, 60, 72, 78, 84 and 67.

Bayesian Inferences for Two Parameter Weibull Distribution Kipkoech W. Cheruiyot 1, Abel Ouko 2, Emily Kirimi 3

Correlation and Simple Linear Regression

Dimensionality Reduction and Learning

Asymptotic Statistical Analysis on Special Manifolds (Stiefel and Grassmann manifolds) Yasuko Chikuse Faculty of Engineering Kagawa University Japan

Previous lecture. Lecture 8. Learning outcomes of this lecture. Today. Statistical test and Scales of measurement. Correlation

Cooper and McGillem Chapter 4: Moments Linear Regression

9.1 Introduction to the probit and logit models

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

Likelihood Ratio, Wald, and Lagrange Multiplier (Score) Tests. Soccer Goals in European Premier Leagues

Empirical likelihood ratio tests with power one

STK3100 and STK4100 Autumn 2017

Simulation Output Analysis

MS exam problems Fall 2012

Modeling uncertainty using probabilities

Relations to Other Statistical Methods Statistical Data Analysis with Positive Definite Kernels

PARAMETER ESTIMATION OF GEOGRAPHICALLY WEIGHTED MULTIVARIATE t REGRESSION MODEL

EVALUATING COMPARISON BETWEEN CONSISTENCY IMPROVING METHOD AND RESURVEY IN AHP

TESTS FOR COVARIANCE MATRICES WITH HIGH-DIMENSIONAL DATA. Saowapha Chaipitak

CHAPTER 2. = y ˆ β x (.1022) So we can write

22 Nonparametric Methods.

ICS141: Discrete Mathematics for Computer Science I

under the curve in the first quadrant.

Differential Method of Thin Layer for Retaining Wall Active Earth Pressure and Its Distribution under Seismic Condition Li-Min XU, Yong SUN

Web Caching and Content Distribution: A View From the Interior

Some Notes on the Probability Space of Statistical Surveys

Reliability-based Structural Design

Evaluation of the Singular Stress Field. of the Crack Tip of the Infinite Elastic Plate

STA302/1001-Fall 2008 Midterm Test October 21, 2008

Transcription:

Predctg Survvl Outcomes Bsed o Compoud Covrte Method uder Cox Proportol Hzrd Models wth Mcrorrys PLoS ONE 7(10). do:10.1371/ourl.poe.0047627. http://dx.plos.org/10.1371/ourl.poe.0047627 Tkesh Emur Grdute Isttute of Sttstcs, Ntol Cetrl Uversty, Tw Jot work wth Dr. Che, Y-Hu d Dr. Che, Hsu-Yu, Isttute of Sttstcl Scece, Acdem Sc

Outles 1. Itro: Survvl lyss wth mcrorrys 2. Exstg methods (Rdge regresso & Lsso) 3. A method kow s compoud covrte predcto (ths method s ot studed suffcet detls lterture) 4. Refemet of compoud covrte method (Proposed method) 5. Comprso wth exstg methods (by rel dt)

Survvl predcto Clcl chrcterstcs (Age / Stge / Tumor type, etc. 5 yer survvl probblty Clssfcto ( Hgh-rsk / Low-rsk ) (crter for chemotherpy ) Brest ccer ptet Cox proportol hzrdmodel: h( t x ) 0( )exp( β h t x x (Age, Stge, Tumor type, etc.) β βˆ : Prtl lkelhood estmtes from dt ), Clssfc to : βˆ x c (Low - rsk) ; βˆ x c (Hgh - rsk)

Survvl predcto wth mcrorrys Mcrorrys (v t Veer et l., 2002 Nture) Brest ccer ptet x ( x 1,..., x550 1,2,...,98 ), Mcrorrys s useful for predctg brest ccer ptets (Jese et l., 2002; v t Veer et l., 2002; Vver et l., 2002; Zho et l., 2011) Cox proportol hzrd model: h( t x ) h0 ( t)exp( βx ), Dffcult to get β βˆ due to hgh - dmesolty ( p )

Avlble methods wth mcrorry Lsso (Cox-regresso wth L_1 plty) Tbshr (1997), Gu & L (2005), Segl (2006) Rdge regresso (Cox-regresso wth L_2 pelty) Verve & Howelge(1994), Zho et l. (2011) Uvrte selecto v Cox-regresso Jesse et l. (2002), Che et l. (2007) Cluster lyss v t Veer et l., 2002; medcl studes Others (PC, supervsed PC, prtl lese squre, etc.) Amog my methods, rdge regresso hs the overll-best predcto power (Bovelstd et l., 2007; v Weerge e l., 2009; Bovelstd d Borg, 2011)

Two obectves of our study: 1. Study compoud covrte predcto (Tukey 1993) I survvl dt, compoud covrte s emprclly used: (Beer et l., 2002; Che et l., 2007; Rdmcher et l, 2002; Mtsu, 2006) * But, less studed the sttstcl lterture * So, ts comprtve performce s ukow 2. Propose to refe compoud covrte predcto v Shrkge techque

Survvl dt : {( t,, x ); 1,..., } t x ( x 1 1 f deth,,..., x p ), p Set up : ether tme to deth or tme to cesorg Exmple: 1 f cesorg Brest ccer dt (v Houwlge et l. 2006) =295, p=4919, Cesored proporto = 73% Lug ccer dt (Che et l., 2007) =125, p=672, Cesored proporto = 70% Dt lyss (lter) t t ( 1) ( 0 ) Observto Perod Cesored

Prtl lkelhood Cox regresso wth p> Cox proportol hzrd model: h( t x If p, the mxmum s ot uque Pelzed prtl lkelhood (well-kow methods) Eve f ) h L 1 log L ( t)exp( βx 1 lr exp( βx ) exp( βx ) h ( t)exp( x 0 0 1 1 p p 1 p (Lsso) p 1, where, the mxmum s uque. { l : t ( 0 s determed by cross-vldto, Verve & Houwelge, 1993) l ) log L 1 (Cox - Rdge x R ( / 2) p 1 2 regresso ) l t } )

Compoud covrte predcto Uvrte Cox regresso Pr( t t t dt t t, x ) / dt h0 ( t)exp( x ) for 1,..., p A collecto of p uvrte lkelhood estmtors βˆ(0) ( ˆ,..., ˆ ) 1 p where ˆ rg mx L 0, ( ), d where L 0, ( ) 1 lr exp( x ) exp( x l ) d R { l : t l t } Compoud covrte predcto βˆ (0) x c (Low - rsk) ; βˆ (0) x c (Hgh - rsk)

Refg Compoud covrte predcto Compoud covrte predcto uses mrgl (uvrte) lkelhood oly: 0 0 βˆ (0) ( ˆ 1,..., ˆ p ) where βˆ(0) rg mx L rg mx L, ( ) We try to ehce predcto power by corportg multvrte lkelhood formto exp( βx 1 L 1 exp( βx lr l ) Ide: Mxture of Uvrte d multvrte lkelhood l log L 1 where [0,1] s prespecfed â s determedby cross vldto (Verve& Houwelge,1993) β ) (1 )log L 0 β p 1

CompoudShrkge estmtor : where l log L 1 (1 )log L βˆ( ) rgmxl 0,

Theoretcl results Asymptotc ormlty ( βˆ( ˆ) β0) N( 0, Σ( β0)) Plug- vrce estmtor Σ ˆ (ˆ( β ˆ)) A Σ h V U 1 A { V / } A 1 h 2 2 1 ( β){ d CV ( ) / d } h /, where CV ( ) Estmtg fucto of, V observed Fsher formto I Score fucto *Resoble performce eve whe p >. U p

Comprso wth rel dt Dt: Lug ccer dt (Che et l., 2007 NEJM) β ˆ =63, p=97 trg dt compoud covrte compoud shrkge Rdge regresso Lsso Predct Low-rsk =62, p=97 test dt { x 1,...,62 } Hgh-rsk Predcto : wherec s βˆ x c (Low - rsk) ; the med of {ˆ βx βˆ x, 1,...,} c (Hgh - rsk),

0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 Survvl curves for Hgh vs. Low rsk groups for =62 testg dt; p-vlue for testg the equlty of two groups Compoud covrte Compoud shrkge p-vlue = 0.076 p-vlue = 0.179 0 10 20 30 40 0 10 20 30 40 Rdge regresso Lsso p-vlue = 0.923 p-vlue = 0.607 0 10 20 30 40 0 10 20 30 40 50

0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 Survvl curves for Hgh, Medum, Low rsk groups for =62 testg dt; p-vlue for testg the equlty of two groups Compoud covrte Compoud shrkge p-vlue = 0.547 p-vlue = 0.089 0 10 20 30 40 0 10 20 30 40 Rdge regresso Lsso p-vlue = 0.742 p-vlue = 0.547 0 10 20 30 40 0 10 20 30 40

Summry of dt lyss The compoud covrte method s best terms of the bry (good/poor) clssfcto of ptets survvl prospect. O the other hd, the three survvl curves re bestseprted by the proposed (compoud shrkge) method Overll rkg of ptets rsk my be best predcted by the proposed method Thk you for your tteto