Nonparametric multivariate inference on shift parameters.

RATIONALE AND OBJECTIVES Consider a study evaluating the prognostic value of prostate-specific antigen (PSA), in the presence of Gleason score, in differentiating between early and advanced prostate cancer. This data set features subjects divided into two groups (early versus advanced cancer), with one manifest variable (PSA), one covariate (Gleason score), and one stratification variable (hospital, taking three values). We present a nonparametric method for estimating a shift in median PSA score between the two groups, after adjusting for differences in Gleason score and stratifying on hospital. This method can also be extended to cases involving multivariate manifest variable. MATERIALS AND METHODS Our method uses estimating equations derived from an existing rank-based estimator of the area under the receiver operating characteristic curve (AUC). This existing AUC estimator is adjusted for stratification and for covariates. We use the adjusted AUC estimator to construct a family of tests by shifting manifest variables in one of the treatment groups by an arbitrary constant. The null hypothesis for these tests is that the AUC is half. We report the set of shift values for which the null hypothesis is not rejected as the resulting confidence region. RESULTS Simulated data show performance consistent with the distributional approximations used by the proposed methodology. This methodology is applied to two examples. In the first example, the mean difference in PSA levels between advanced and nonadvanced prostate cancer patients is estimated, controlling for Gleason score. In the second example, to assess the degree to which age and baseline tumor size are prognostic factors for breast cancer recurrence, differences in age and tumor size between subjects with recurrent and nonrecurrent breast cancer, stratified on Tamoxifen treatment and adjusting for tumor grade, are estimated. CONCLUSIONS The proposed methodology provides a nonparametric method for a statistic measuring adjusted AUC to be used to estimate shift between two manifest variables.

[1]  K. Zou,et al.  Statistical validation based on parametric receiver operating characteristic analysis of continuous classification data. , 2003, Academic radiology.

[2]  Jana Jurečková,et al.  Asymptotic Linearity of a Rank Statistic in Regression Parameter , 1969 .

[3]  Subha Madhavan,et al.  G-DOC: a systems medicine platform for personalized oncology. , 2011, Neoplasia.

[4]  D. Bamber.  The area above the ordinal dominance graph and the area below the receiver operating characteristic graph , 1975 .

[5]  A. Renshaw,et al.  Biochemical Outcome after radical prostatectomy, external beam Radiation Therapy, or interstitial Radiation therapy for clinically localized prostate cancer , 1998 .

[6]  D. Mehrotra,et al.  Rank-Based Analyses of Stratified Experiments: Alternatives to the van Elteren Test , 2010 .

[7]  Mirko Schiemann,et al.  Vascular guide wire navigation with a magnetic guidance system: experimental results in a phantom. , 2004, Radiology.

[8]  J. Hanley,et al.  The meaning and use of the area under a receiver operating characteristic (ROC) curve. , 1982, Radiology.

[9]  M. J. van de Vijver,et al.  Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. , 2006, Journal of the National Cancer Institute.

[10]  H. B. Mann,et al.  On a Test of Whether one of Two Random Variables is Stochastically Larger than the Other , 1947 .

[11]  Robert V. Hogg,et al.  Discriminant Analysis Based on Ranks , 1978 .

[12]  Kyle J Myers,et al.  The importance of ROC data. , 2011, Academic radiology.

[13]  G. Koch,et al.  Stratified Multivariate Mann–Whitney Estimators for the Comparison of Two Treatments with Randomization Based Covariance Adjustment , 2011 .

[14]  C. Metz.  Basic principles of ROC analysis. , 1978, Seminars in nuclear medicine.

[15]  K. Zou,et al.  Comparison of adjustment methods for stratified two‐sample tests in the context of ROC analysis , 2012, Biometrical journal. Biometrische Zeitschrift.

[16]  A James O'Malley,et al.  Bayesian multivariate hierarchical transformation models for ROC analysis , 2006, Statistics in medicine.

[17]  J R Fielding,et al.  Bayesian regression methodology for estimating a receiver operating characteristic curve with two radiologic applications: prostate biopsy and spiral CT of ureteral stones. , 2001, Academic radiology.

[18]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[19]  C. Metz,et al.  Maximum likelihood estimation of receiver operating characteristic (ROC) curves from continuously-distributed data. , 1998, Statistics in medicine.