
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Oncology, Markers, Clinical Correlates |
Departments of Epidemiology and Biostatistics [D. A. V., G. H.] and Medicine [W. K. K., H. I. S.], Memorial Sloan-Kettering Cancer Center, New York, New York 10021
| ABSTRACT |
|---|
|
|
|---|
Experimental Design: With serial PSA measurements from 254 patients with androgen-independent prostate cancer, we use a time-dependent Cox model with a nonlinear log relative risk function to quantify the strength of the association between time-dependent PSA and survival. The nonlinear log relative risk function provides a flexible Cox model and a more accurate measure of the strength of association between PSA and survival.
Results: Among these 254 patients, there were 247 deaths (median survival time = 13.0 months). Median follow-up time for those alive or censored was 59.3 months (range, 36.871.3 months). An association was observed between PSA and survival (P < 0.01, two-sided test). However, time-dependent PSA explains only 17% of the variation in survival.
Conclusions: Use of this methodology demonstrates that there remains sufficient variation in survival unaccounted for by PSA measurements in this patient cohort. Other factors, perhaps unknown, exist that determine survival outcome. Consideration of PSA alone as a surrogate can produce misleading information regarding the risk of death; its use as a surrogate for survival is not warranted when designing a clinical trial in this patient population.
| INTRODUCTION |
|---|
|
|
|---|
These studies, however, have focused only on establishing a statistical association between PSA and survival without quantifying the strength of the association. This is troubling because the establishment of a statistical association through an appreciably small P value is not sufficient to justify its use as a surrogate for survival (5) . Focusing on PSA, when survival is only marginally affected by biological mechanisms (including treatment), is misguided; instead the investigative effort should be expanded to include the identification and understanding of the effect of other biological mechanisms on survival (6) . An example of this phenomenon can be seen in a recent study of hormone therapy in the noncastrate metastatic population conducted by Eisenberger et al. (7) . Specifically, the trial, using a PSA end point, demonstrated a higher proportion of patients achieving a normalization of PSA using a combined androgen blockade versus a monotherapy approach. Despite this, no difference in survival was observed between the two cohorts (7) .
In this study, our objective is to quantify the level of association between PSA and survival. To accomplish this, we built a time-dependent Cox model with a general, nonlinear log relative risk function. The patient population used for the analysis is a cohort of castrate metastatic (androgen-independent) prostate cancer patients, with progressive disease, enrolled on consecutive therapeutic trials. We chose the castrate metastatic state because it is reasonable to assume that the survival end point will be observed and that there will be a reduced confounding effect of testosterone modulation, which is known to affect serum PSA levels. For this analysis, a post-therapy change in PSA is not specifically defined; rather, we use every PSA measurement available from all patients post-therapy to estimate the log relative risk function. Our approach is population-based. PSA values are not identified by individual subjects but rather are aggregated to assess the overall relationship in the population between PSA and survival time. Using the nonlinear log relative risk Cox model, we demonstrate that there remains sufficient variation in survival unaccounted for by PSA. This suggests that focusing on PSA alone, while important, does not provide sufficient information to predict survival. As a result, we conclude that the various post-therapy change definitions of PSA such as percentage decrease from baseline (at varying landmark time points), PSA slope, PSA doubling time, and PSA velocity (3 , 7, 8, 9, 10, 11, 12, 13) should not be used as surrogates for survival time when designing and evaluating clinical trials.
| MATERIALS AND METHODS |
|---|
|
|
|---|
![]() |
To quantify the level of association between survival and PSA, we propose a model that builds upon the above time-dependent Cox model but uses a nonlinear specification of the log relative risk function (16, 17, 18, 19)
. This general, more flexible log relative risk function does not assume a linear relationship between the log relative risk of death and PSA but rather can take many forms, including that of a linear relationship if in fact that is appropriate. This has the practical advantage of allowing the Cox model to be applicable in a wider range of data sets. The model is written as follows:
![]() |
The reason for using the time-dependent Cox model with a general, nonlinear log relative risk function is to add flexibility when we use a measure of association. Model-based measures of association are influenced by (a) the strength of the relationship between the covariate and response and (b) proper model specification. In order for the measure to primarily reflect the strength of the association between time-dependent PSA values and survival, we replaced the linear log relative risk function Ztß found in the conventional Cox model with the less restrictive general log relative risk function S(Zt). Thus, in this flexible Cox model, we do not constrain the function to be linear, and furthermore, we allow the data to determine its functional form. The method of estimating the nonlinear function S(Zt) is described in the "Appendix. "
The measure that quantifies the strength of association is the R2 statistic adapted to the Cox model (20) . This statistic, similar to the R2 statistic used to summarize the goodness of fit in linear regression models, measures the amount of variability in survival explained by time-dependent PSA values (see "Appendix "). In the Cox model, R2 is a function of the relative risk function exp[S(t)] over time and is formed by the ratio of the information provided by having the time-dependent PSA values (Zt) in the model to the information provided by the model without PSA. This measure varies between 0 and 1, with high R2 values indicating that PSA is discriminant with respect to the risk of death. Heuristically, when patients whose death is imminent have high PSA values (and conversely when patients with low and moderate PSA levels do not die in the short term), a high R2 would result. All analyses were performed with S-Plus 2000 Professional software. Code is available from the authors upon request.
| RESULTS |
|---|
|
|
|---|
|
Our objective was to quantify the strength of the association between the time-dependent PSA values and survival. In the proposed Cox model, this association is described by the general, nonlinear log relative risk function S. The estimated relative risk function S is presented in Fig. 2
. It is clear that the relationship between loge(1 + PSA) and the relative risk is not linear. The flexibility of this nonlinear function allows for this situation. We can see that the relative risk increases between 0 and 1, decreases for values of loge(1 + PSA) between 1 and 3, and then increases thereafter, more steeply between 3 and 6 and then moderately between 6 and 9, before rising again. This region of decreasing relative risk was unanticipated.
|
| DISCUSSION |
|---|
|
|
|---|
In our example, we examined the relationship between PSA, modeled as a time-dependent covariate, and survival in 254 castrate metastatic (androgen-independent) cancer patients. We observed a region of decreasing relative risk (1 < loge(1 + PSA) <3), which was unanticipated; a possible biologically plausible explanation for this observation is that there is a heterogeneous group of patients within this region that contains patients with poorly differentiated metastatic tumors with low PSA values and patients with low PSA values but without metastatic disease. The former group of patients does much worse clinically and therefore causes the risk of death, relative to the other patients, to decrease. The flexibility of this nonlinear function allowed for this observation. Had we used the conventional Cox model with a linear log relative risk function, this group of patients would not have been identified.
We found that despite a statistically significant association between PSA and survival (P < 0.01, two-sided test), only 17% of the variability in survival could be explained by the time-dependent PSA values; other factors, many of which are unknown, account for the remaining variation. Until these factors are discovered and incorporated, the results suggest that PSA is not adequate as a surrogate end point for survival in this clinical state, and using it (or any of the various definitions of a post-therapy change in PSA) as a surrogate in clinical trials could lead to problems in the design and the final interpretation of outcome of these clinical trials. Despite this finding, PSA is important as a predictor for survival. This methodology could be extended by incorporating other predictors (e.g., age, pretherapy serum lactate dehydrogenase, pretherapy hemoglobin, and other factors collected serially), in addition to PSA, into the nonparametric generalized additive Cox model to determine how much variability is unexplained by such a model.
This study is limited by a relatively small number of patients who were treated heterogeneously and by the fact that we looked only at castrate metastatic (androgen-independent) prostate cancer patients. In addition, the estimation of the nonlinear log relative risk function is sensitive to the number of knots chosen for the B-spline function (see "Appendix "). In general, an increase in the number of knots used to estimate the nonlinear log relative risk function will increase the R2 value, but at the cost of increasing the complexity of the model. The goal is to balance the level of association against model parsimony, which we believe we have done.
Currently, PSA remains the best clinical marker we have, but the value of PSA as a surrogate for survival, as suggested by the results of this study, is limited. Efforts must continue to be made to find other predictors of survival.
| Appendix 1 |
|---|
|
|
|---|
Involved in the specification of this model is choosing the number of knots that will be used to estimate the cubic B-splines (23) . These knots are break points in the data where each of the B-splines will be estimated. Our goal is to produce S as a simple, smooth function of the time-dependent PSA covariate. For this data, four knots are used to specify this relative risk covariate function, and we allow the data to determine where the knots are placed (24) .
The measure of association chosen for this model is the R2 statistic of Kent and OQuigley shown below (20)
.
![]() |
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Supported in part by CaPCURE. ![]()
2 To whom requests for reprints should be addressed, at Department of Epidemiology and Biostatistics, Memorial Sloan-Kettering Cancer Center, 1275 York Avenue, Box 44, New York, NY 10021. Phone: (212) 639-2145; Fax: (212) 717-3137; E-mail: verbeld{at}mskcc.org ![]()
3 The abbreviation used is: PSA, prostate-specific antigen. ![]()
Received 1/25/02; revised 4/12/02; accepted 5/ 1/02.
| REFERENCES |
|---|
|
|
|---|
Commentary
This article has been cited by other articles:
![]() |
S. F. Slovin, A. S. Wilton, G. Heller, and H. I. Scher Time to Detectable Metastatic Disease in Patients with Rising Prostate-Specific Antigen Values following Surgery or Radiation Therapy Clin. Cancer Res., December 15, 2005; 11(24): 8669 - 8673. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. I. Scher, M. J. Morris, W. K. Kelly, L. H. Schwartz, and G. Heller Prostate Cancer Clinical Trial End Points: "RECIST"ing a Step Backwards Clin. Cancer Res., July 15, 2005; 11(14): 5223 - 5232. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. L. Kaufman, W. Wang, J. Manola, R. S. DiPaola, Y.-J. Ko, C. Sweeney, T. L. Whiteside, J. Schlom, G. Wilding, and L. M. Weiner Phase II Randomized Study of Vaccine Treatment of Advanced Prostate Cancer (E7897): A Trial of the Eastern Cooperative Oncology Group J. Clin. Oncol., June 1, 2004; 22(11): 2122 - 2132. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. I. Scher, M. Eisenberger, A. V. D'Amico, S. Halabi, E. J. Small, M. Morris, M. W. Kattan, M. Roach, P. Kantoff, K. J. Pienta, et al. Eligibility and Outcomes Reporting Guidelines for Clinical Trials for Patients in the State of a Rising Prostate-Specific Antigen: Recommendations From the Prostate-Specific Antigen Working Group J. Clin. Oncol., February 1, 2004; 22(3): 537 - 556. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. M. Sandler and M. L. DeSilvio Surrogate End Points for Prostate Cancer: What Is Prostate-Specific Antigen Telling Us? J Natl Cancer Inst, September 17, 2003; 95(18): 1352 - 1353. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |