Breast cancer tumor growth estimated through mammography screening data
 Harald WeedonFekjær^{1, 2}Email author,
 Bo H Lindqvist^{3},
 Lars J Vatten^{4},
 Odd O Aalen^{2} and
 Steinar Tretli^{1, 4}
DOI: 10.1186/bcr2092
© WeedonFekjær et al.; licensee BioMed Central Ltd. 2008
Received: 27 August 2007
Accepted: 08 May 2008
Published: 08 May 2008
Abstract
Introduction
Knowledge of tumor growth is important in the planning and evaluation of screening programs, clinical trials, and epidemiological studies. Studies of tumor growth rates in humans are usually based on small and selected samples. In the present study based on the Norwegian Breast Cancer Screening Program, tumor growth was estimated from a large population using a new estimating procedure/model.
Methods
A likelihoodbased estimating procedure was used, where both tumor growth and the screen test sensitivity were modeled as continuously increasing functions of tumor size. The method was applied to cancer incidence and tumor measurement data from 395,188 women aged 50 to 69 years.
Results
Tumor growth varied considerably between subjects, with 5% of tumors taking less than 1.2 months to grow from 10 mm to 20 mm in diameter, and another 5% taking more than 6.3 years. The mean time a tumor needed to grow from 10 mm to 20 mm in diameter was estimated as 1.7 years, increasing with age. The screen test sensitivity was estimated to increase sharply with tumor size, rising from 26% at 5 mm to 91% at 10 mm. Compared with previously used Markov models for tumor progression, the applied model gave considerably higher model fit (85% increased predictive power) and provided estimates directly linked to tumor size.
Conclusion
Screening data with tumor measurements can provide populationbased estimates of tumor growth and screen test sensitivity directly linked to tumor size. There is a large variation in breast cancer tumor growth, with faster growth among younger women.
Introduction
Mammography screening is now an established part of the health service in developed countries. There is, however, still an ongoing discussion related to optimizing mammography screening, including determining optimal time intervals between screenings and which age groups to invite. For these decisions, adequate estimates of breast cancer tumor growth and screening test sensitivity (STS) are crucial. In addition, better knowledge of tumor growth will benefit the evaluation of screening programs [1], as well as the interpretation of clinical trials and epidemiological studies. There are some observational studies of patients that were initially overlooked at earlier mammograms [2–4] or were refused treatment [2, 3], but these studies are small and are probably influenced by length of time bias, since slowgrowing tumors spend relatively longer times in preclinical stages visible on mammograms. To our knowledge, no largescale populationbased clinical observational studies of untreated cancers have therefore been performed as cancers are usually treated in populations with good cancer surveillance.
Tumor growth can also be indirectly observed as tumor progression, estimated from variations in cancer incidence in screening trials or programs. These studies [1] are usually analyzed using Markov models [5, 6], where the mean time for a breast cancer tumor to growth from screeningdetectable size to clinical detection without screening – the socalled mean sojourn time – and the STS are estimated. The Markov model, however, has no separate variable for individual variation, and the estimated variables are highly correlated with contributions from both the underlying biological processes and the given screening method. The estimated parameters therefore have no explicit relation to the biological process of tumor growth, and are often difficult to compare between different countries, as the STS is defined as 'the proportion of cancers detected at screening among screening detectable cancers', using the evaluated procedural as its own reference.
Tumor growth can be estimated by comparing tumor sizes from clinicaldetected and screeningdetected cases, but the applied statistical models only partly utilize these data. Chen and colleagues [7] used tumor size in a classical Markov model, and van Oortmarssen and colleagues [8] used tumor size in a simulation approach – but both studies only categorized tumor size into two or three groups. On the contrary, some clinical observation studies fully utilize tumor size measurements with tumor growth modeled as a continuous function of tumor size [2, 9], but these studies of nontreated or overlooked cancers are small and the results may not be valid due to either selection bias or length of time bias.
The aim of the present study was to utilize modern computer power on data from a populationbased screening program, with precise standardized measurements of tumor size, to reliably estimate tumor growth and STS.
Materials and methods
Setting: data
In 1995 the Norwegian Government initiated an organized populationbased service screening program [10], in which mammography results and interval cancer cases are carefully registered by the Cancer Registry of Norway. The Norwegian Breast Cancer Screening Program (NBCSP) originally included four counties. Other counties were subsequently included, and by 2004 the screening program achieved nationwide coverage. All women between 50 and 69 years of age receive a written invitation biannually, and twoview mammograms from participating women are independently evaluated by two readers.
To make the results comparable with estimates provided in previous studies [5, 12–15], all cases of ductal carcinoma in situ (DCIS) – a noninvasive form of breast tumor – were included. In addition, estimates were also deduced excluding DCIS cases to check the potential effect of DCIS cases. Several tumors detected at the same time in one woman were counted as one case, with size measurements given for the largest tumor. Only new primary breast cancers were included in this study.
In the NBCSP, tumor measurements are performed on pathological sections after surgery, and tumors are measured diagonally between the outer edges. All measurements were performed in a standardized manner according to specifications given in a national quality assurance manual. Tumor size measurements were available for 92% of the cancers detected at screening. There were several reasons for missing tumor measurements: some tumors were torn up at the surgical operation before tumor measurements were taken, others were difficult to measure on pathological crosssections, and some tumors had grown into the outer skin. In addition, a substantial part had received tumorreducing treatment before the pathological material was removed. Tumors of unknown size are therefore probably somewhat different from tumors with an observed tumor size. Patients who received tumorreducing treatment will typically have larger tumors, which in practice could have biased our estimates – leading to higher growth rates. Sensitivity analyses related to possible bias in tumor sizes were therefore performed.
Tumor size measurements of clinical breast cancers that emerge without screening are needed for the tumor growth model suggested in the present article. Since women who do not attend screening represent a selected group, possibly with different alertness to early symptoms, tumor size measurements made before the start of the official screening program were used. The Cancer Registry of Norway did not receive reliable information on tumor size prior to the official screening program. At Haukeland University Hospital (covering Bergen, Norway's second largest city), however, a good database for tumor measurements of clinical invasive breast cancers exists [16]. We were able to use these data, where 503 women aged 50 to 69 years were diagnosed with breast cancer between 1985 and 1994. Among these cases, 433 women (86%) had registered tumor measurements in millimeters. A comparison of tumor measurements found at screening and in the Haukeland University Hospital database of clinically detected cases is shown in Figure 2.
Growth model specification
Several studies have examined growth curves, both in general and for human breast cancer tumors in particular. The conclusion has often been that the growth curves can be described by either a logistic function [17] or a Gompertz function [9, 18]. For the range of tumor sizes that are relevant for screening, there are only minor differences between logistic and Gompertz growth given probable parameters. Spratt and colleagues used a variant of the general logistic growth curve with a maximum tumor volume of 40 cell doublings, equaling a ball of 128 mm in diameter, after testing several models on a clinical dataset that mostly consisted of overlooked tumors [9, 19]. To make the comparison with Spratt and colleagues' observations [9, 19], we used the same variant of the lognormal logistic growth model in the present study. This implies an almost exponentional growth for the smallest tumors, with decelerating growth as the tumors approaches the maximum of 128 mm in diameter (see examples in Figure 3). In addition to the chosen model, model fits for several alternative choices of maximum tumor volume were evaluated, with moderate effects on the estimated values.
where κ_{ i }is a log normally distributed growth rate with mean α_{1} and variance α_{2}, V_{max} is the maximum tumor volume (set to a tumor of 128 mm in diameter), and V_{cell} is the volume of one cell. (As all calculations in the present paper use a relative cancer time, the choice of V_{cell} does not affect the given estimates.)
Overall, this can be seen as a mixed effects model with individual logistic growth curves and a lognormally distributed random effect.
As tumor measurements in the NBCSP are the maximum diameter, the real tumor volume will in practice be smaller. The most important part of the model, however, is the modeled growth curve, and sensitivity analyses show little effects of a general reduction in modeled tumor volume as a function of tumor diameter.
Screening test sensitivity model specification
where β_{1} defines how fast tumor sensitivity increases by tumor size and β_{2} relates STS to tumor size, with β_{2} = 0 equaling S(0) = 0.5 (places the sensitivity curve in relation to tumor size).
Parameter estimation
Since mammography screening detects a higher proportion of the larger prevalent tumors compared with the smaller prevalent tumors, the pool of undiagnosed tumors is expected to have a clear overrepresentation of small tumors shortly after screening. One would suspect this could lead to relatively small tumors detected shortly after screening, followed by gradually increasing tumor sizes with the time since last screening. This trend is severely damped, however, as each tumor before detection must reach a certain individual size to produce sufficient symptoms to alarm the woman. In practice, the relationship between tumor size and clinical detection results in only a vague trend in interval cancer tumor sizes by time since screening (correlation = 0.01 in the NBCSP), whereas the number of interval cancers increases sharply. We have therefore chosen to disregard the size distribution of interval cancers, and build our estimation procedure on the observed frequency of interval cancers by the time since screening, the number of cases found at screening, the tumor size distribution of screening cancers, the assumed background incidence, and the size distribution of clinical tumors without screening (based on historical data).
Combining these data with our model, the model parameters can be estimated by maximum likelihood calculation. As the full likelihood includes several integrals, the actual maximum likelihood calculations are performed discretely, grouping the data into sufficiently small time and tumor size intervals.
where i is an indicator for the size group, sn is the number of screened women, sc_{ i }is the number of screening cases in size group i, and sp_{ i }is the probability of a woman having a tumor in size group i at screening, given the parameter set {α_{1}, α_{2}, β_{1}, β_{2}}.
where ic_{ j }is the observed number of cancers j months after screening and ie_{ j }is the expected number of cancers j months after screening, given the parameter set {α_{1}, α_{2}, β_{1}, β_{2}}.
The probability of finding a cancer in a given size group at screening (sp_{ i }) and the expected number of interval cases (ie_{ j }), given a set of known parameters, (α_{1}, α_{2}, β_{1}, β_{2}), are therefore needed for the estimation of model parameters. There is no available knowledge regarding the number of tumors initiated at different ages that have the potential of becoming screening or clinically detected cancers later on. The expected number of cases given a known tumor growth rate cannot therefore be deduced directly. It is possible, however, to calculate the expected number of cases at screening using back calculations from the expected number of clinical cancers seen without screening. This idea is not unlike the theory behind Markov models of cancer screening [12], utilizing known quantities regarding the expected number of future cancers to calculate the expected number of cases at an earlier stage.
where p_{ g }is the relative proportion of breast cancers of size g without screening.
where PYR_{ j }is the number of person years in interval j and fs_{j,g}is the probability that a clinical cancer in size group g would have been found if screened j months earlier.
In practice, both P(tumor of size g was of size i f months earlier α1, α_{2}) in equation (8) and P(tumor of size gs was of size g j months earlier α1, α_{2}) in equation (10) can be calculated in the following three stages. First, by rearranging the growth formula equation (1), expressing earlier tumor size as a function of present tumor size and tumor growth rate (κ_{ i }). Then calculating upper and lower limits for tumor growth (κ_{ i }), constituting the requested probability. Finally, calculating the probability for tumor growth within the given limits using the lognormal distribution and assumed growth parameters {α_{1}, α_{2}}.
Combining these formulas, maximum likelihood estimates of the observed dataset can be deduced by numerically maximizing the loglikelihood.
Modeling choices: specifications
While the number of cancers in the interval between screenings can be observed directly, the expected number without screening has to be estimated. As the NBCSP offers screening to all women in a defined population, no parallel control group is available to carry out this estimation. In addition, commitment to screening can, and probably does, vary with individual risk factors, so those who do not attend are not a suitable control group either.
The background incidence was therefore calculated from historical data combined with an estimated time trend. In practice, data from 1990 to 1994 were used with time trend estimates from an ageperiod cohort model with additional screening parameters [21]. Incidence rates vary among age groups and counties, and the estimate was therefore weighted by the number of personyears in each combination of age group and county. Further, it may be a problem that the sharply increased use of hormone replacement therapy (HRT) in the 1990s [22] has influenced the historical time trends in breast cancer incidence. HRT is known to increase breast cancer risk [23], and Bakken and colleagues [22] found a relative risk of breast cancer of 2.1 for current versus never users in Norway. Combining sales figures with risk estimates, Bakken and colleagues estimated the proportion of breast cancer cases that could be attributed to HRT use as 27% among Norwegian women 45 to 64 years of age. HRT use increased sharply from the period that was used to calculate the expected incidence without screening (1990 to 1994) to our estimating period (1996 to 2002). Therefore 21% was added to the estimated background incidence (when otherwise not noted), on the basis of information regarding increased breast cancer risk and HRT sale figures found in Bakken and colleagues [22]. With this correction, the expected incidence without screening was estimated as 190 cases/100,000 personyears for women 50 to 59 years of age, and as 219 cases/100,000 personyears for women 60 to 69 years of age.
When calculating the expected number of cases at screening, we cannot include an infinite number of future time intervals. We therefore limited the growth rates to realistic levels given the women's current age, and reweighted the distribution. Experiences with different limits show that the choice of growth limit had little effect on the estimated values.
Statistical calculations
All calculations, simulations, and plots were performed using the R statistical package [24]. Data were transformed from the Norwegian breast cancer screening database and were summarized using a combination of SQL commands and the statistical package SPLUS (Insightful, Seattle, USA). To double check the implemented R functions, new datasets were simulated and the results compared with the expected number of cases.
Maximum likelihood estimates were found by optimization over all four parameters simultaneously, using the optimize function found in the R package [24]. For these calculations, time intervals of 1 month were used. Tumor sizes were categorized to 1 mm, 2 mm, 5 mm, 10 mm, 15 mm, ..., 100+ mm, as the background data revealed that many pathologists approximated tumor sizes to the nearest 5 mm, 10 mm, 15 mm, ..., 100 mm (data not shown). To look at possible age differences, estimates were calculated separately for women aged 50 to 59 years and women aged 60 to 69 years, in addition to all age groups combined. Calculations were very computer intensive, with a huge number of probability calculations needed to calculate the expected number of cases for a given parameter set.
The main estimates are presented with (pointwise) confidence intervals showing their (random) uncertainty. Robust 95% confidence intervals were calculated by 1,000 smoothed biascorrected parametric bootstrap replications [25], resampling all of the observed data except the assumed breast cancer incidence without screening. Simulations were used to deduce the overall STS and the mean sojourn time. As a validation of the model fit, observed values versus expected values were plotted. In addition, the traditional Markov model of breast cancer screening [5, 12, 26] was compared with the new method using onefifth holdout crossvalidation, measuring the weighted mean square differences. For evaluation of crossvalidation results, P values calculated from 50 parametric bootstrap replications were used.
Results
Parameter estimates
For all age groups combined, model parameters were estimated as {α_{1}, α_{2}, β_{1}, β_{2}} = {1.07, 1.31, 1.47, 6.51}, while the two age groups 50 to 59 years and 60 to 69 years gave estimates of {1.38, 1.36, 1.50, 6.33} and {0.70, 1.18, 1.46, 6.65}, respectively. While parameters are hard to interpret and compare, several relevant quantities can be deduced once parameters are estimated.
Estimated tumor growth
Summary of results with 95% biascorrected bootstrap confidence intervals
With supposed higher background incidence due to increased hormone therapy use (+21%)  Combined estimate (50 to 69 years) with nonadjusted background incidence  

Women aged 50 to 59 years  Women aged 60 to 69 years  Combined estimate (50 to 69 years)  
Time taken from 10 mm to 20 mm (years)  
Mean  1.4 (1.1, 1.5)  2.1 (1.8, 2.4)  1.7 (1.5, 1.8)  2.0 (1.8, 2.2) 
Standard deviation  1.9 (1.7, 2.2)  2.4 (2.1, 2.7)  2.2 (2.0, 2.4)  2.7 (2.5, 2.9) 
Volume doubling time at 15 mm (days)  
25th percentile  29 (19, 36)  65 (47, 79)  41 (32, 48)  29 (21, 35) 
50th percentile  73 (56, 86)  143 (116, 165)  99 (84, 111)  94 (77, 107) 
75th percentile  180 (148, 205)  308 (253, 352)  234 (204, 259)  287 (243, 322) 
Screening test sensitivity  
5 mm  29 (21, 36)  24 (18, 30)  26 (22, 31)  26 (22, 31) 
10 mm  92 (88, 99)  91 (86, 98)  91 (88, 96)  91 (87, 96) 
Indicators of potential screening efficacy  
Mean sojourn time (years)  2.3 (2.0, 2.6)  3.5 (3.1, 3.9)  2.8 (2.6, 3.1)  3.4 (3.1, 3.6) 
Proportion of tumors visible on screening  0.95 (0.94, 0.96)  0.95 (0.95, 0.96)  0.95 (0.95, 0.96)  0.95 (0.95, 0.96) 
Estimates of tumor growth rates compared with Spratt and colleagues' [9] estimates based on overlooked and nontreated cancers
Percentile  Growth parameter (κ_{ i }in equation (1))  Time (years) tumor takes to grow from 10 mm to 20 mm  

Estimate  Spratt et al.  Estimate  Spratt et al.  
1st  0.2  0.2  10.9  12.9 
5th  0.4  0.6  6.3  4.6 
25th  1.3  1.7  2.0  1.5 
50th  3.0  3.2  0.9  0.8 
75th  7.2  5.2  0.4  0.5 
95th  25.4  11.8  0.1  0.2 
99th  61.9  33.7  0.0  0.1 
Estimated screening test sensitivity
The mammography STS was estimated to increase sharply from around 2 mm to 12 mm, with the STS reaching 26% at 5 mm and 91% at 10 mm (Figure 4b). There was no significant difference in the estimated STS between the two age groups (P = 0.83 for the STS at 5 mm).
Overall screening test sensitivity and mean sojourn time
Using simulations to combine the STS and the given distribution of clinical tumors, we found that nearly all cancers were likely to be visible at screening before reaching clinical detection (Table 1). Defining the mean sojourn time as the time tumors are visible at screening before clinical detection, these cancers have a mean sojourn time of 3.0 years – resulting in an overall mean sojourn time of 2.9 years for all cases. In older women the mean sojourn time was estimated to be significantly longer. There were large variations in the sojourn time, and the standard deviation was estimated as 5.0 years, indicating that the Markov model (which equals the mean sojourn time and standard deviation) does not allow for enough individual variation in growth rates.
Model fit
A more exponential tumor growth curve modeled through a higher maximum tumor volume weakened the overall model fit (data not shown), supporting the assumption that the doubling time of the tumor volume may increase with increasing tumor size (as assumed by the logistic model). To explore possible biases due to missing tumor measurements at screening, we applied several different assumptions regarding the true tumor diameter of the unknown tumors, revealing very stable parameter estimates (data not shown).
Discussion
The present study introduces a new way of modeling cancer growth and STS, based on data from a large screening program. Tumor growth was estimated to vary greatly between individual tumors, with tumors taking a mean time of 1.7 years to grow from 10 mm to 20 mm in diameter. The STS was estimated to increase rapidly with tumor size, from 26% at 5 mm to 91% at 10 mm.
Applied to the NBCSP data, the new model gives a very good model fit, and a significantly better predictive power than the previously used Markov model [26]. Certain aspects of the model need further investigation, however, and some have argued that cancer growth either follows exponential [27] or Gompertz [18] growth functions, and not the assumed logistic growth curve [19]. The practical difference between the logistic and Gompertz curve is relatively small, but an exponential growth curve could probably alter the results significantly. Mathematically, a logistic function with very large maximum tumor volume almost equals the exponential curve. Several alternative levels of maximum tumor volume were therefore tested, giving weaker model fit as the maximum tumor volumes increased, thereby strengthening our assumption of a bounded growth function (rather than an exponential growth function).
Another possible objection to the model is that the STS is assumed to always increase towards 100% as the tumor size increases, while some cancers probably never become visible on mammograms [28]. To test this alternative hypothesis, a threeparameter STS function with a parameter for maximum STS was tested, giving no indication of a lower maximum STS level. To limit the complexity of the estimation procedure and the presentation of the new model, only data from the first screening round were used in this study. Data from subsequent screening rounds were still available, and while the model predicted a 71% decline in detected cancers from first screening to second screening, the observed decline was only 46%. This is a considerable predicted–observed difference, and the NBCSP generally has shown a surprisingly high cancer rate at the second screening. In addition to possible problems with the model itself, this can be an effect of changes in HRT use in the study period (increasing the general breast cancer risk), of increased sensitivity in the second round due to use of earlier mammograms, of better training of staff with time, or of an overrepresentation of communities with high cancer risk in the second screening round.
Even with a highquality cancer registry, problems with the applied data may cause more bias to the estimated values than the applied model assumptions. Studying the fit of the new model (Figure 5), there are some signs of discrepancy in the last half of the interval following screening, with too many observed cases. This may be an effect of unregistered opportunistic screening, since opportunistic screening has been available at many private institutions, and cancers detected outside the NBCSP have in practice been registered as interval cancers. Unfortunately, no detailed information is available on the extent of opportunistic screening in the different age groups, and there is no precise information on whether interval cancers have been detected by opportunistic screening or clinical symptoms. Preliminary studies by the Norwegian Cancer Registry indicate that approximately 10% of the NBCSP's invited women are screened outside the program each year. This percentage may, however, be lower since the level of opportunistic screening may be higher among nonattendees of the public screening. Preliminary attempts to estimate the level of opportunistic screening, and to correct the estimated STS and growth rates, indicate little bias in the estimated mean cancer growth and the STS, while the variation in cancer growth rates decreased substantially.
Another problem can be the assumed background incidence without screening, as the estimates changed somewhat (Table 2) when removing the correction for a probable higher background incidence due to increased use of HRT [22]. The correction probably improves the estimates, but there is uncertainty. Based on typical user patterns, it is possible that HRT use could have been higher than assumed among woman 50 to 59 years of age, and somewhat lower among the 60 to 69 years age group. The correction may therefore be too small for the younger age group and too strong for the older age group. In addition, HRT use fluctuated during the study period, and may have influenced the cancer incidence, the STS and the tumor growth in different ways. Most importantly, HRT use is known to reduce the STS [29–31], at least partially due to increased breast density. Since HRT use has been quite common in Norway, the STS may have been even higher with moderate levels of HRT use. Tumor growth estimates may also be affected by HRT use, and both the STS and tumor growth estimates should be viewed in the light of the relatively high HRT use during the study period.
Overall, sensitivity investigations indicate that the new model is probably less vulnerable to several potential biases than the Markov model [5, 12], possibly as a result of more utilized data. The model is substantially different from the Markov model, rendering direct comparisons difficult, but the slightly different overall screening efficiency indicators confirm the estimated mean sojourn time and STS from other studies [32, 33], with a shorter mean sojourn time and a higher STS than found in a recent Norwegian study [26].
More importantly, the new model estimates tumor growth directly connected to tumor measurements, similar to the earlier nonpopulationbased studies of overlooked cancers [4, 9], but using a much larger populationbased material. The results confirm previously reported growth rates (Table 2), large variances in tumor growth, and a probable bounded growth function, suggesting less selection bias in studies of (earlier) overlooked cancers than previously assumed [9].
Earlier studies have shown decreasing tumor progression and higher STS with increasing age [32, 33]. The present study confirms the previously reported decrease in tumor growth with age, but we found no trend in STS associated with age. This is surprising, but very few new breast cancers were diagnosed in the first months after screening among women 50 to 59 years of age (Figure 1), indicating a surprisingly high STS for the younger age group. An investigation of which aspects of the data influence the various parameters revealed that differences in tumor size between screening and clinically detected tumors are vital for STS estimates. In the Norwegian screening program there is little difference in screening and clinical tumor sizes between the two age groups, a fact that indicates small differences in STS by age. This could of course be an artifact due to the modeling, but could also be an effect of very different recall rates in the two age groups [34]. Indeed, the issue clearly motivates further examinations of the STS among younger women.
Compared with studies of overlooked cancers and with studies of women who refused treatment, the populationbased approach greatly increases the number of observed cases and applies data that are probably less biased. Generally, this model combines many of the advantages of the large populationbased Markov methods [5, 12], with more specific tumor growth estimates found in clinical studies of overlooked cancers. This makes the model suited for both optimizing screening designs and cost–benefit analyses.
By combining the present tumor growth and STS estimates with death rates, different screening designs may be evaluated even more efficiently than seen in earlier studies. Plevritis presented an advanced simulation approach with continuous tumor growth [35]. The approach was based on similar tumor size backcalculation techniques to those used in our study, but the lack of estimates probably limited the practical impact of that study. In practice, Markov models have often been used to evaluate screening designs [36], but without a direct link to tumor size it is difficult to separate and compare the mean sojourn time and the STS between screening programs [26]. In recent years, more advanced simulation models have been suggested – as seen in the US National Cancer Institute Cancer Intervention and Surveillance Modeling Network [1, 37] – further emphasizing the need for precise tumor growth estimates.
Estimates with and without ductal carcinoma in situ (DCIS)
Data used  Time taken (years) to grow from 10 mm to 20 mm  Volume doubling time (days) at 15 mm  Screening test sensitivity  Indicators of potential screening efficacy  

Mean  Standard deviation  25th percentile  50th percentile  75th percentile  5 mm  10 mm  Mean sojourn time (years)  Proportion of tumors visible on screening  
With ductal carcinoma in situ  1.7  2.2  41  99  234  26  91  2.9  0.95 
Without ductal carcinoma in situ  1.5  1.9  44  96  209  23  91  2.5  0.95 
Conclusion
To summarize, tumor growth and STS estimates can be directly linked to tumor size in a full population study, resulting in very useful growth estimates directly connected to a biologically relevant measure. Tumor growth seems to vary greatly between tumors, with higher growth rates among younger women. Most tumors become visible at screening when they reach a diameter of 5 mm to 10 mm.
Abbreviations
 DCIS:

ductal carcinoma in situ
 HRT:

hormone replacement therapy
 NBCSP:

Norwegian Breast Cancer Screening Program
 STS:

screening test sensitivity.
Declarations
Acknowledgements
The authors would like to thank HansOlav Adami, Astri Syse and Tom K Grimsrud for valuable comments on the article, Solveig Hofvind and Wenche Melbye for sharing their excellent knowledge about the NBCSP and the corresponding database, and Bjørn Ove Mæhle for allowing us to use the Haukeland dataset on clinical tumor sizes. The project has received financial support from the Norwegian Cancer Society (Grant E03031/001), and is connected to the BMMS research group at University of Oslo, Norway.
Authors’ Affiliations
References
 Berry DA, Cronin KA, Plevritis SK, Fryback DG, Clarke L, Zelen M, Mandelblatt JS, Yakovlev AY, Habbema JD, Feuer EJ: Effect of screening and adjuvant therapy on mortality from breast cancer. N Engl J Med. 2005, 353: 17841792. 10.1056/NEJMoa050518.PubMedView ArticleGoogle Scholar
 Spratt JS, Meyer JS, Spratt JA: Rates of growth of human solid neoplasms: part I. J Surg Oncol. 1995, 60: 137146. 10.1002/jso.2930600216.PubMedView ArticleGoogle Scholar
 Spratt JS, Meyer JS, Spratt JA: Rates of growth of human neoplasms: part II. J Surg Oncol. 1996, 61: 6883. 10.1002/10969098(199601)61:1<68::AIDJSO2930610102>3.0.CO;2E.PubMedView ArticleGoogle Scholar
 Peer PG, van Dijck JA, Hendriks JH, Holland R, Verbeek AL: Agedependent growth rate of primary breast cancer. Cancer. 1993, 71: 35473551. 10.1002/10970142(19930601)71:11<3547::AIDCNCR2820711114>3.0.CO;2C.PubMedView ArticleGoogle Scholar
 Prevost TC, Launoy G, Duffy SW, Chen HH: Estimating sensitivity and sojourn time in screening for colorectal cancer: a comparison of statistical approaches. Am J Epidemiol. 1998, 148: 609619.PubMedView ArticleGoogle Scholar
 Karlin S, Taylor H: An Introduction to Stochastic Modeling. 1998, Academic Press, New YorkGoogle Scholar
 Chen HH, Prey MU, Babcock DS, Day NE: Breast carcinoma cell kinetics, morphology, stage, and the preclinical screendetectable phase. Lab Invest. 1997, 2: 923.Google Scholar
 van Oortmarssen GJ, Habbema JD, Maas van der PJ, de Koning HJ, Collette HJ, Verbeek AL, Geerts AT, Lubbe KT: A model for breast cancer screening. Cancer. 1990, 66: 16011612. 10.1002/10970142(19901001)66:7<1601::AIDCNCR2820660727>3.0.CO;2O.PubMedView ArticleGoogle Scholar
 Spratt JA, von Fournier D, Spratt JS, Weber EE: Decelerating growth and human breast cancer. Cancer. 1993, 71: 20132019. 10.1002/10970142(19930315)71:6<2013::AIDCNCR2820710615>3.0.CO;2V.PubMedView ArticleGoogle Scholar
 Wang H, Kåresen R, Hervik A, Thoresen SØ: Mammography screening in Norway; results from the first screening round in four counties and costeffectiveness of a modeled nationwide screening. Cancer Causes Control. 2001, 1: 3945. 10.1023/A:1008999403069.View ArticleGoogle Scholar
 Cancer in Norway 2006. [http://www.kreftregisteret.no]
 Day NE, Walter SD: Simplified models of screening for chronic disease: estimation procedures from mass screening programmes. Biometrics. 1984, 40: 114. 10.2307/2530739.PubMedView ArticleGoogle Scholar
 Duffy SW, Chen HH, Prevost TC, Tabár L: Markov chain models of breast tumour progression and its arrest by screening. Quantitative Methods for the Evaluation of Cancer Screening. 2001, Arnold, London, 4260.Google Scholar
 Duffy SW, Chen HH, Tabar L, Day NE: Estimation of mean sojourn time in breast cancer screening using a Markov chain model of both entry to and exit from the preclinical detectable phase. Stat Med. 1995, 14: 15311543. 10.1002/sim.4780141404.PubMedView ArticleGoogle Scholar
 Habbema JD, van Oortmarssen GJ, Lubbe JT, Maas van der PJ: The MISCAN simulation program for the evaluation of screening for disease. Comput Methods Programs Biomed. 1985, 20: 7993. 10.1016/01692607(85)900483.PubMedView ArticleGoogle Scholar
 Hartveit F, Maehle BO, Pettersen KC: Size of breast carcinomas at operation related to tumour growth rate. Breast Cancer Res Treat. 1987, 10: 4750. 10.1007/BF01806134.PubMedView ArticleGoogle Scholar
 Hart D, Shochat E, Agur Z: The growth law of primary breast cancer as inferred from mammography screening trials data. Br J Cancer. 1998, 78: 382387.PubMedPubMed CentralView ArticleGoogle Scholar
 Norton L: A Gompertzian model of human breast cancer growth. Cancer Res. 1988, 48: 70677071.PubMedGoogle Scholar
 Spratt JA, von Fournier D, Spratt JS, Weber EE: Mammographic assessment of human breast cancer growth and duration. Cancer. 1993, 71: 20202026. 10.1002/10970142(19930315)71:6<2020::AIDCNCR2820710616>3.0.CO;2#.PubMedView ArticleGoogle Scholar
 Meyer JS, Prey MU, Babcock DS, McDivitt RW: Breast carcinoma cell kinetics, morphology, stage, and host characteristics. A thymidine labeling study. Lab Invest. 1986, 54: 4151.PubMedGoogle Scholar
 Moller B, WeedonFekjaer H, Hakulinen T, Tryggvadottir L, Storm HH, Talback M, Haldorsen T: The influence of mammographic screening on national trends in breast cancer incidence. Eur J Cancer Prev. 2005, 14: 117128. 10.1097/0000846920050400000007.PubMedView ArticleGoogle Scholar
 Bakken K, Alsaker E, Eggen AE, Lund E: Hormone replacement therapy and incidence of hormonedependent cancers in the Norwegian Women and Cancer study. Int J Cancer. 2004, 112: 130134. 10.1002/ijc.20389.PubMedView ArticleGoogle Scholar
 Beral V: Breast cancer and hormonereplacement therapy in the Million Women Study. Lancet. 2003, 362: 419427. 10.1016/S01406736(03)145965.PubMedView ArticleGoogle Scholar
 R. [http://www.rproject.org]
 Chernick MR: Bootstrap Methods: A Practitioner's Guide. 1999, WileyInterscience, USAGoogle Scholar
 WeedonFekjaer H, Vatten LJ, Aalen OO, Lindqvist B, Tretli S: Estimating mean sojourn time and screening test sensitivity in breast cancer mammography screening; new results. J Med Screen. 2005, 12: 172178. 10.1258/096914105775220732.PubMedView ArticleGoogle Scholar
 Brown BW, Atkinson EN, Bartoszynski R, Thompson JR, Montague ED: Estimation of human tumor growth rate from distribution of tumor size at detection. J Natl Cancer Inst. 1984, 72: 3138.PubMedGoogle Scholar
 Ma L, Fishell E, Wright B, Hanna W, Allan S, Boyd NF: Casecontrol study of factors associated with failure to detect breast cancer by mammography. J Natl Cancer Inst. 1992, 84: 781785. 10.1093/jnci/84.10.781.PubMedView ArticleGoogle Scholar
 Porter PL, El Bastawissi AY, Mandelson MT, Lin MG, Khalid N, Watney EA, Cousens L, White D, Taplin S, White E: Breast tumor characteristics as predictors of mammographic detection: comparison of interval and screendetected cancers. J Natl Cancer Inst. 1999, 91: 20202028. 10.1093/jnci/91.23.2020.PubMedView ArticleGoogle Scholar
 Kavanagh AM, Mitchell H, Giles GG: Hormone replacement therapy and accuracy of mammographic screening. Lancet. 2000, 355: 270274. 10.1016/S01406736(99)073195.PubMedView ArticleGoogle Scholar
 Litherland JC, Stallard S, Hole D, Cordiner C: The effect of hormone replacement therapy on the sensitivity of screening mammograms. Clin Radiol. 1999, 54: 285288. 10.1016/S00099260(99)90555X.PubMedView ArticleGoogle Scholar
 Tabar L, Fagerberg G, Chen HH, Duffy SW, Smart CR, Gad A, Smith RA: Efficacy of breast cancer screening by age. New results from the Swedish TwoCounty Trial. Cancer. 1995, 75: 25072517. 10.1002/10970142(19950515)75:10<2507::AIDCNCR2820751017>3.0.CO;2H.PubMedView ArticleGoogle Scholar
 Paci E, Duffy SW: Modelling the analysis of breast cancer screening programmes: sensitivity, lead time and predictive value in the Florence District Programme (1975–1986). Int J Epidemiol. 1991, 20: 852858. 10.1093/ije/20.4.852.PubMedView ArticleGoogle Scholar
 Hofvind S, Geller B, Vacek PM, Thoresen S, Skaane P: Using the European guidelines to evaluate the Norwegian Breast Cancer Screening Program. Eur J Epidemiol. 2007, 22: 447455. 10.1007/s106540079137y.PubMedView ArticleGoogle Scholar
 Plevritis SK: A mathematical algorithm that computes breast cancer sizes and doubling times detected by screening. Math Biosci. 2001, 171: 155178. 10.1016/S00255564(01)000542.PubMedView ArticleGoogle Scholar
 Maas van der PJ, de Koning HJ, van Ineveld BM, van Oortmarssen GJ, Habbema JD, Lubbe KT, Geerts AT, Collette HJ, Verbeek AL, Hendriks JH, Rombach JJ: The costeffectiveness of breast cancer screening. Int J Cancer. 1989, 43: 10551060. 10.1002/ijc.2910430617.PubMedView ArticleGoogle Scholar
 Cancer Intervention and Surveillance Modeling Network. [http://cisnet.cancer.gov/]
 World Health Organization: IARC Handbooks of Cancer Prevention: Handbook 7: Breast Cancer Screening. 2001, IARC Press, Lyon, FranceGoogle Scholar
 Ross MB, Buzdar AU, Hortobagyi GN, Lukeman JM: Spontaneous regression of breast carcinoma: followup report and literature review. J Surg Oncol. 1982, 19: 2224. 10.1002/jso.2930190107.PubMedView ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.