An improved approximation to the exact solution of the two-stage clonal growth model of cancer

Multistage clonal growth models are of interest for cancer risk assessment because they can explicitly incorporate data on cell replication. Both approximate and exact formulations of the two stage growth model have been described. The exact solution considers the conditional probability of tumors arising in previously tumor-free animals; the approximate solution estimates total probability of tumor formation. The exact solution is much more computationally intensive when time-dependent cell growth parameters are included. The approximate solution deviates from the exact solution at high incidences and probabilities of tumor. This report describes a computationally tractable, 'improved approximation' to the exact solution. Our improved approximation includes a correction term to adjust the unconditional expectation of intermediate cells based on the time history of formation of intermediate cells by mutation of normal cells (recruitment) or by cell division in the intermediate cell population (expansion). The improved approximation provided a much better match to the exact solution than the approximate solution for a wide range of parameter values. The correction term also appears to provide insight into the biological factors that contribute to the variance of the expectation for the number of intermediate cells over time.