A predictive risk probability approach for microarray data with survival as an endpoint.
Gene expression profiling has played an important role in cancer risk classification and has shown promising results. Since gene expression profiling often involves determination of a set of top rank genes for analysis, it is important to evaluate how modeling performance varies with the number of selected top ranked genes incorporated in the model. We used a colon data set collected at Moffitt Cancer Center as an example of the study, and ranked genes based on the univariate Cox proportional hazards model. A set of top ranked genes was selected for evaluation. The selection was done by choosing the top k ranked genes for k = 1 to 12,500. An analysis indicated a considerable variation of classification outcomes when the number of top ranked genes was changed. We developed a predictive risk probability approach to accommodate this variation by identifying a range number of top ranked genes. For each number of top ranked genes, the procedure classifies each patient as having high risk (score = 1) or low risk (score = 0). The categorizations are then averaged, giving a risk score between 0 and 1, thus providing a ranking for the patient's need for further treatment. This approach was applied to the colon data set and demonstrated the strength of this approach by three criteria: First, a univariate Cox proportional hazards model showed a highly statistically significant level (log-rank chi(2) statistics = 110 with p-value <10(-16)) for the predictive risk probability classification. Second, the survival tree model used the risk probability to partition patients into five risk groups showing a good separation of survival curves (log-rank chi(2) statistics = 215). In addition, utilization of the risk group status identified a small set of risk genes that may be practical for biological validation. Third, analysis of resampling the risk probability suggested the variation pattern of the log-rank chi(2) in the colon cancer data set was unlikely caused by chance.
Chen
Dung-Tsa
DT
Biostatistics Division, Moffitt Cancer Center & Research Institute, University of South Florida, Tampa, Florida 33612, USA. Dung-Tsa.Chen@moffitt.org
Schell
Michael J
MJ
Chen
James J
JJ
Fulp
William J
WJ
Eschrich
Steven
S
Yeatman
Timothy
T
D003110:Colonic Neoplasms / Q000209:etiology* / Q000235:genetics / Q000401:mortality; D023381:Endpoint Determination; D020869:Gene Expression Profiling; D006801:Humans; D020411:Oligonucleotide Array Sequence Analysis; D011336:Probability; D016016:Proportional Hazards Models; D012306:Risk
J Clin Oncol. 2005 May 20;23(15):3526-35
15908663
Cell Cycle. 2003 Sep-Oct;2(5):415-7
12963829
N Engl J Med. 2002 Jun 20;346(25):1937-47
12075054
Proc Natl Acad Sci U S A. 2003 Jul 8;100(14):8418-23
12829800
Nucleic Acids Res. 2003 Feb 15;31(4):e15
12582260
Bioinformatics. 2002 Dec;18(12):1625-32
12490447
J Clin Oncol. 2005 Mar 10;23(8):1631-5
15755970
Proc Natl Acad Sci U S A. 2001 Sep 25;98(20):11462-7
11562467
PLoS Biol. 2004 Apr;2(4):E108
15094809
Proc Natl Acad Sci U S A. 2001 Apr 24;98(9):5116-21
11309499
N Engl J Med. 2002 Dec 19;347(25):1999-2009
12490681
Bioinformatics. 2002;18 Suppl 1:S120-7
12169539