A copula approach for detecting prognostic genes associated with survival outcome in microarray studies

A challenging and crucial issue in clinical studies in cancer involving gene microarray experiments is the discovery, among a large number of genes, of a relatively small panel of genes whose elements are associated with a relevant clinical outcome variable such as time-to-death or time-to-recurrence of disease. A semiparametric approach, using dependence functions known as copulas, is considered to quantify and estimate the pairwise association between the outcome and each gene expression. These time-to-event type endpoints are typically subject to censoring as not all events are realized at the time of the analysis. Furthermore, given that the total number of genes is typically large, it is imperative to control a relevant error rate in any gene discovery procedure. The proposed method addresses the two aforementioned issues by direct incorporation of the censoring mechanism and by appropriate statistical adjustment for multiplicity. The performance of the proposed method is studied through simulation and illustrated with an application using a case study in lung cancer.