Power calculations for survival analyses via Monte Carlo estimation.

Abstract

Power calculations can be a useful step in the design of epidemiologic studies. For occupational and environmental cohort studies, however, the calculation of statistical power has been difficult because researchers are often interested in situations where exposure assignment is time-dependent, and in research questions that pertain to cumulative exposure-mortality trends evaluated with statistical methods for survival analysis. These conditions are not easily accommodated by available software or published formulas for power calculation. Monte Carlo methods can be used to estimate statistical power for survival analyses. Simple computer programs are presented to illustrate this approach. We show that, for the simple case of a randomized clinical trial involving a dichotomous exposure, the results of power calculations derived via this Monte Carlo approach conform to values derived using a previously published formula. We then illustrate how the Monte Carlo approach may be extended to obtain estimates of statistical power for analyses of cumulative exposure-mortality trends under conditions more typical of occupational cohort studies. The Monte Carlo approach provides a way to perform power calculations for a wide range of study conditions. The approach illustrated in this study should simplify the task of calculating power for survival analyses, particularly in epidemiologic research on occupational cohorts.