In response to our paper concerning the distribution of tumor clonogens at the end of treatment, Dr. Yakovlev points out that the probability-generating function, and hence the distribution, can be calculated exactly, and thus we could have reached our conclusions without using computer simulation. We thank Dr. Yakovlev for making this useful observation. The possibility of doing this was also pointed out to us by Elliot Landaw, a colleague at UCLA, soon after our paper appeared. In our paper, we derived a formula for the probability of tumor cure after fractionated treatment [our Eq. (8)] that approximated the simulation results quite well and represented a clear improvement over the usual Poisson formula [our Eqs. (1) and (2)] for cure probability among tumors in which proliferation occurs during treatment. A slightly different formula for tumor cure probability can be obtained from the probability-generating function derived by Dr. Yakovlev by substituting s = 0 into his Eqs. (1), (2), (5), and (6). The only difference between our formula, which could also be obtained using generating functions (1), and Dr. Yakovlev's formula is that we used a stochastic term for the cell-killing component but a deterministic term for the growth component, whereas Dr. Yakovlev used stochastic elements for both of these components. The good agreement between our formula and the simulation results is a reflection of the stronger relative influence of cell killing than of cell proliferation for the treatment schedules we considered. Our model and Dr. Yakovlev's model are identical when there is no cell division during treatment (i.e., a = 0), in which case both formulae reduce to a binomial form. We noted in our paper [(2), p. 277 and Appendix B], as Dr. Yakovlev points out again, that the Poisson model provides an excellent approximation in this case. Thus we too were not at all surprised by the good agreement between the predictions of Poisson statistics and the numerical results in the assumed absence of tumor cell proliferation during treatment. To assist readers who might be interested in investigating the differences in tumor cure between various fractionation schemes (3) and are not familiar with generating functions, we gi e below the explicit formula for calculating the probability of tumor cure for fractionation regimens in which e ch dose fraction might be of a different size and in which the time between dose fractions could vary. The formula g ven below applies only to the situation where there is a rel tively short interval between fractions, such that the probability of a cell undergoing two successive divisions in that interval is negligibly small. We are currently investigating the more general situation where the interval between fractions could be much larger, for example, in a splitcourse radiotherapy regimen. The formula given below is derived from the generating functions given by Dr. Yakovlev. Assume there are n fractions and let p, be the probability of cell survival associated with the ith dose and ai be the probability of cell division between the ith and (i + 1)st dose. Let