Strategies for power calculations in predictive biomarker studies in survival data

Biomarkers and genomic signatures represent potentially predictive tools for precision medicine. Validation of predictive biomarkers in prospective or retrospective studies requires statistical justification of power and sample size. However, the design of these studies is complex and the statistical methods and associated software are limited, especially in survival data. Herein, we address common statistical design issues relevant to these two types of studies and provide guidance and a general template for analysis.

A statistical interaction effect in the Cox proportional hazards model is used to describe predictive biomarkers. The analytic form by Peterson et al. and Lachin is utilized to calculate the statistical power for both prospective and retrospective studies.

We demonstrate that the common mistake of using only Hazard Ratio's Ratio (HRR) or two hazard ratios (HRs) can mislead power calculations. We establish that the appropriate parameter settings for prospective studies require median survival time (MST) in 4 subgroups (treatment and control in positive biomarker, treatment and control in negative biomarker). For the retrospective study which has fixed survival time and censored status, we develop a strategy to harmonize the hypothesized parameters and the study cohort. Moreover, we provide an easily-adapted R software application to generate a template of statistical plan for predictive biomarker validation so investigators can easily incorporate into their study proposals.

Our study provides guidance and software to help biostatisticians and clinicians design sound clinical studies for testing predictive biomarkers.