Estimating the prevalence of breast cancer using a disease model: data problems and trends
© Kruijshaar et al; licensee BioMed Central Ltd. 2003
Received: 28 March 2003
Accepted: 14 April 2003
Published: 14 April 2003
Health policy and planning depend on quantitative data of disease epidemiology. However, empirical data are often incomplete or are of questionable validity. Disease models describing the relationship between incidence, prevalence and mortality are used to detect data problems or supplement missing data. Because time trends in the data affect their outcome, we compared the extent to which trends and known data problems affected model outcome for breast cancer.
We calculated breast cancer prevalence from Dutch incidence and mortality data (the Netherlands Cancer Registry and Statistics Netherlands) and compared this to regionally available prevalence data (Eindhoven Cancer Registry, IKZ). Subsequently, we recalculated the model adjusting for 1) limitations of the prevalence data, 2) a trend in incidence, 3) secondary primaries, and 4) excess mortality due to non-breast cancer deaths.
There was a large discrepancy between calculated and IKZ prevalence, which could be explained for 60% by the limitations of the prevalence data plus the trend in incidence. Secondary primaries and excess mortality had relatively small effects only (explaining 17% and 6%, respectively), leaving a smaller part of the difference unexplained.
IPM models can be useful both for checking data inconsistencies and for supplementing incomplete data, but their results should be interpreted with caution. Unknown data problems and trends may affect the outcome and in the absence of additional data, expert opinion is the only available judge.
Estimates of disease-specific incidence, prevalence and mortality, specified by age and sex, are important information to health care policy and planning. They are essential inputs to cost-effectiveness analyses and burden of disease calculations. Empirical data, however, are often difficult to obtain or are of questionable validity. To remedy some of these data problems, disease models have been developed that describe the relationship between the epidemiological parameters, by exploiting the causal structure of a disease. Incidence, prevalence, mortality models (IPM models) formalise the relationship between the three parameters, using the fact that incidence has to precede prevalence, and that cause-specific mortality can only follow disease. IPM models have been used frequently both to supplement missing data and to study the agreement between different epidemiological data [1–4]. Our previous study supported the formal validity of IPM models, but when the modelling was applied to empirical data on four types of cancer, the model calculations differed to a large extent from the empirical data . For breast cancer the difference was particularly large. It was argued that these discrepancies may indicate inconsistencies in the data, but that they may also be caused by time trends.
When the data for one disease are not in accordance with each other, they are internally inconsistent. Inconsistencies may be caused by differences in the completeness of the data. For example, when more incident cases are missed than deaths, incidence and mortality are inconsistent. Also, inconsistencies may arise when the data were derived from different contexts (e.g. a different region) or measured differently (e.g. varying case-definitions). Applying inconsistent incidence and mortality to an IPM model will result in under- or overestimating prevalence, and thus in discrepancies between model estimations and empirical prevalence data. Time trends, on the other hand, may cause the data to appear inconsistent in a steady state model, while in fact they are not. Because prevalence is the resultant of incident cases from the past, it cannot react instantaneously to changes in incidence and case-fatality, but only with a certain delay. It is possible to account for the effects of time trends in a dynamic model, but this requires additional input data on the nature and size of the trends, which are not available for most diseases. Often, we do not even know whether a trend is present or not, and the researcher faces a dilemma what to do with the discrepancies. Adjusting observed data for apparent inconsistencies that are in fact the consequence of past trends would rather defeat the purpose of IPM models.
For breast cancer in the Netherlands the discrepancy between observed prevalence and prevalence calculated from incidence and mortality was particularly large . Fortunately, for breast cancer several data problems are known, and, in addition, there is a tentative estimate of the trend in incidence. This allowed us to quantify the relative contribution of the trend and several known data problems on the discrepancy and to throw some light upon the researcher's dilemma. Even though prevalence of breast cancer is not a very useful epidemiological measure, it does allow us to illustrate the difficulties in the use of IPM modelling because of the relative abundance of data.
We calculated the point-prevalence of breast cancer in the Netherlands and its 95 % confidence interval from national incidence rates and cause-specific mortality rates using the IPM model described by Barendregt et al. . We compared it to regional prevalence data. To determine the separate effects of trends and known data problems we next recalculated the model:
a) incorporating the incompleteness in the prevalence data (see below),
b) adjusting for a trend in incidence,
c) adjusting for double counting of incident women with secondary primaries, and
d) taking excess mortality from non-breast cancer deaths into account.
For point c a single cohort model was used, while for a, b, and d it was necessary to use a multi-cohort model that takes into account both age and calendar time. We estimated the proportion of the difference that was explained by each of the recalculations from the overall differences in the number of prevalent cases using 1993 population figures and summing over all age groups.
First, we calculated prevalence from national incidence and mortality rates of female breast cancer (ICD-9 code 174) for 1991–1995, averaged by five-year age groups up to 85+. Incidence data obtained from the Netherlands Cancer Registry (NKR) are based on pathology and hospital admission data, and the mortality data from Statistics Netherlands (CBS) on death certificates. To enable comparison with the prevalence data we used incidence rates excluding in-situ tumours. Incidence and mortality rates were first interpolated to one-year age groups using the cubic-spline method and then converted into probabilities (see the Appendix Additional file: 1for formulas). A 95% confidence interval of the calculated prevalence was obtained by parametric bootstrapping assuming a Poisson distribution for numbers of incident cases and deaths in each age group. The @risk software programme  was used to simulate 10000 iterations by Monte Carlo sampling.
The Regional Cancer centre South (IKZ) was the only cancer registry in the Netherlands that has estimated prevalence. The IKZ determined prevalence for a specific part (the core region) of the region South at 1-1-1993 by checking the vital status of all incident cases registered in the region since 1970 against the municipal population administration and the National Death Index. If a person had moved to another municipality in the Netherlands, data from that municipal population administration was used (under Dutch law, registration with the municipal population administration is obligatory within five days of changing address). In case of migration to another country, a case was lost to follow up. In-situ tumours are not included. For more information on these data we refer to Coebergh et al. . A 95% confidence interval around these data (by five-year age group up to 85+) was calculated assuming a binomial distribution (see formula 2 in the Appendix Additional file: 1). Point estimates and confidence intervals were interpolated to one-year age groups as described in the Appendix Additional file: 1.
Limitations of the prevalence data
Prevalence registered by the IKZ does not include patients diagnosed before 1970 and is consequently underestimated. To quantify the effect of this, we created a dynamic model, incorporating a parameter y that represents the number of years prior to incidence. Prevalence in the reference year y0, was calculated from prevalence in the 23 years prior to y0 (see Appendix Additional file: 1).
Also, prevalence data refer to a specific region only, whereas the model calculations are based on national input. We showed in our previous study , that differences between regional and national incidence and mortality rates hardly affected the calculated prevalence. Therefore, we did not further examine this here.
Trend in incidence
Next, we estimated the effect of a trend of increasing incidence on the calculated prevalence. While incidence increased, the population mortality rate for breast cancer has remained approximately the same in the Netherlands [8–12], thus case-fatality must have declined. Coebergh et al. have estimated the yearly rise in incidence in the region South between 1975 and 1986, to be approximately one percent . The effects of screening are not included in this estimate, as it was not introduced until later, but effects of other types of increased case finding were included – if present.
We estimated the effect using the dynamic model described above, decreasing incidence by one percent for each additional year y prior to the year of reference (see Appendix Additional file: 1). We assumed the trend was present up to 95 years before the year of reference. Sensitivity to the trend was inspected, by applying a 50% higher and lower trend (1.5% and 0.5%) as input values.
Double counting incident women with a secondary primary
The NKR registers the number of incident tumours, whereas prevalence refers to women. Women are thus counted twice, if they have a second primary tumour in their breast (SP). The percentage of SPs by age group were provided for 1991–1995 by the IKZ. We estimated the incidence of women with breast cancer, by subtracting the age-specific proportion SP from the reported incidence rate and recalculated our model.
Excess mortality from non-breast cancer deaths
Breast cancer patients have been shown to be at increased risk of dying from other causes of death , although with longer survival (>20 years) they may experience a decreased risk . Cause-specific mortality data thus underestimate the total excess mortality from breast cancer. We adjusted for this estimating excess mortality from duration-specific relative survival data assuming a lognormal survival distribution from breast cancer and a proportion not dying from breast cancer (proportion cured) , as described in the Appendix Additional file: 1(formula 5). Relative survival data, i.e. excess mortality over and above the background mortality, were reported by the IKZ registry for three cohorts (1970–1979, 1980–1986 and 1987–1992) for ages < 70, and > 70 . Because survival improves over time, we used survival probabilities at one, three and five years after incidence of the youngest 1987–1992 cohort. As longer follow-up was not available for this cohort, survival at 10 and 20 years was estimated using the conditional 10 and 20-year survival in the older cohorts. Fitting the lognormal model allowed us to estimate the cumulative probability of excess mortality with time after incidence, from which we calculated the yearly mortality probability (see Appendix Additional file: 1).
We extended our baseline model to include duration (see Appendix Additional file: 1 for the mathematical description): prevalence at age n, d years after incidence was calculated from the prevalence at the previous age and year. Summation across all years d provided age-specific prevalence.
Comparison of the total number of prevalent cases of breast cancer, estimated in different ways, and 1993 prevalence data.
Total number of cases estimated
% difference from 1993 data
% of the gap explained
Figure 2 and Table 1 also show the effects of the known data problems and trend in incidence (coloured lines). First, excluding cases that became incident more than 23 years before the year of reference had the largest effect and altered the age pattern. Prevalence increased with age more slowly, levelling off to a maximum of 4.4 % at age 78, and declining thereafter. The total number of prevalent cases was 53 % higher than the IKZ prevalence, explaining 39% of the difference between calculated and IKZ prevalence. Second, incorporating a secular trend in incidence in the calculation of prevalence resulted in a lesser increase with age, levelling off after age 70, without reaching a plateau. The secular trend of one percent decreased the discrepancy with 44% (ranging between 49% and 18% for 50% higher and lower estimates of the trend). Third, adjusting for double counting of incident women with a second primary tumour of the breast (SP) had a smaller effect. The percentage of SP did not exceed nine percent for any age group. Subtracting this percentage from the incidence rate explained 17% of the difference. Fourth, taking excess mortality from non-breast cancer deaths into account made a difference of six percent.
We quantified the separate effects of known data problems and a trend in incidence on IPM model calculations of breast cancer prevalence, and inspected their influence on the discrepancy between model estimations and prevalence data. Two factors had a major effect on the estimated prevalence: the limitation of the IKZ prevalence data including only incident cases since 1970, and the trend in incidence. Together, they accounted for a major part (60%) of the discrepancy. Still, their combined effect leaves a part of the discrepancy unexplained. The effects of adjusting for double counting of incident women with secondary primaries of the breast and of taking excess mortality from non-breast cancer deaths into account were small. Consequently, the combined effect of all four factors does not explain the entire difference. Either unknown data problems are present contributing to the remaining difference, or we underestimated some of the effects.
The two factors with the largest effects could influence the calculations because prevalence is a stock variable: it contains cases that became incident in the past. The long survival time of breast cancer explains why these two factors can exert such a strong influence. The time lag after which prevalence fully reflects a change in incidence is determined by the rate at which the pool of prevalent cases is replaced by new cases, which, in turn, depends on the survival time. The underestimation of the prevalence data due to its limitation to incident cases after 1970 may be as large as 44 % at age 85, but is lower on average, as the effect increases with age. When the effect is calculated additional to adjusting for the trend in incidence, it is also smaller: 24 % at age 85. In a small part of the IKZ region prevalence can be based upon incidence registration since 1958 . These data are indeed higher than the IKZ data for the larger region indicating the data will improve with longer follow-up in the future.
The estimated trend in incidence and the effect of adjusting for it are subject to some uncertainty for several reasons. First, the one percent trend is an average figure; its magnitude differs with age and calendar year. Second, the regional estimate may deviate somewhat from the national trend. Furthermore, it is a rounded figure. Nevertheless, as the results of the sensitivity analysis showed that increasing or decreasing the trend by 50% altered the percentage explained by a combined effect (of limiting prevalence and the secular trend) by only 10%-points, we expect these three effects to be small. The estimated one percent may also be too low because it does not include the effect of the introduction of screening around 1990 in the Netherlands. We estimated this effect in an additional analysis, increasing the estimated one percent to five percent for the last four years. The additional effect of screening was very small. Finally, we may have overestimated the effect assuming that a one percent trend was present many years before the period for which it was estimated (1975 to 1986). Assuming no trend before 1975 in an extra analysis increased the calculated prevalence by almost 10 % at age 85, but inspection of the incidence rates since 1958 for the smaller part of the IKZ region showed that incidence has increased since 1958 . The effect of assuming no trend before 1960 was negligible.
The effects of adjusting for secondary primaries and excess mortality were only small, but this is not surprising. The proportion of women with a secondary primary in the breast did not exceed nine percent, preventing a much larger change in the calculated prevalence. Also, the relative risk of breast cancer patients to die from other causes of death is not that high, and may even inverse with longer follow-up [14, 15].
Thus, although the effect of the trend may be somewhat uncertain, other unknown data problems are likely to cause the remaining difference between calculated and IKZ prevalence. One possible explanation for the remaining deviation could be that when the cancer registries started completeness of incidence was not as high as it is now. As a result the prevalence is underestimated by the IKZ. On the other hand, IKZ prevalence was determined by checking population administrations, which involves matching of the registrations, which is never 100% accurate. Therefore, some deaths may have been missed and individuals inaccurately assumed to be alive, resulting in an overestimation of the IKZ prevalence. How large both counteracting effects are, is difficult to determine. An additional explanation may be found in variation between national and regional all-cause mortality, as the equations in our model are based on the assumption of similar background mortality. A higher background mortality in the IKZ region would decrease the discrepancy, but we expect only a minimal effect as the differences in all-cause mortality will be minor. Additional reasons for overestimating incidence or underestimating mortality (both resulting in overestimating prevalence) are difficult to think of. Cancer incidence data in the Netherlands are reliable with a completeness of 96.2 % around 1990 , and probably even higher in the years thereafter. Furthermore, mortality by cause-of-death statistics of the Netherlands are assumed to be reliable and, compared to other European countries, the detection fraction for deaths from cancer is high . The causes of the remaining discrepancy thus remain uncertain.
Even when data are regarded as relatively reliable, as was the case for breast cancer in the Netherlands, data problems may be present. 1993 prevalence data for breast cancer in the IKZ region are underestimated, as they do not include incident cases before 1970. Our analyses show the importance of using IPM modelling to detect data problems. However, we also showed the trend in incidence to have a large effect on the model estimations for breast cancer, complicating the use of IPM models.
IPM models can be useful both for checking for data inconsistencies and for supplementing incomplete data, but in both cases there remains the need for careful interpretation of the results. In the all too common situation where, unlike for breast cancer, no data on the size and nature of trends are available, the effects of trends cannot be estimated. Furthermore, unknown data problems may affect the model estimations in unknown directions. In the absence of additional data the researcher is faced with the dilemma of how to interpret model discrepancies, and expert opinion is the only available judge.
MK carried out the analyses and drafted the manuscript. JB participated in the design of the analyses. LvdP-F provided additional points for analyses. All authors read and approved the final manuscript.
We would like to thank Dr. Jan-Willem Coebergh for his constructive suggestions and comments on this manuscript. This study was sponsored by the Netherlands Institute of Health Sciences.
- Murray CJ, Lopez AD: Quantifying disability: data, methods and results. Bull World Health Organ 1994, 72: 481-494.PubMedPubMed CentralGoogle Scholar
- Murray CJL, Lopez AD: The global burden of disease: a comprehensive assessment of mortality and disability from diseases, injuries and risk factors in 1990 and projected to 2020. Cambridge: Harvard University Press 1996.Google Scholar
- Barendregt JJ, Baan CA, Bonneux L: An indirect estimate of the incidence of non-insulin-dependent diabetes mellitus. Epidemiology 2000, 11: 274-279. 10.1097/00001648-200005000-00008View ArticlePubMedGoogle Scholar
- Mathers CD, Vos ET, Stevenson CE, Begg SJ: The Australian Burden of Disease Study: measuring the loss of health from diseases, injuries and risk factors. Med J Aust 2000, 172: 592-596.PubMedGoogle Scholar
- Kruijshaar ME, Barendregt JJ, Hoeymans N: The use of models in the estimation of disease epidemiology. Bull World Health Organ 2002, 80: 622-628.PubMedPubMed CentralGoogle Scholar
- Barendregt JJ, Oortmarssen GJ, van BA Hout, van JM Bosch, van den, Bonneux L: Coping with multiple morbidity in a life table. Mathematical Population Studies 1998, 7: 29-49.View ArticlePubMedGoogle Scholar
- Palisade : @RISK, advanced risk analysis for spreadsheets. Newfield, NY, USA: Corporation Palisade 4 Edition 2000.Google Scholar
- Coebergh JWW, van der Heijden LH, Janssen-Heijnen MLG: Cancer Incidence and Survival in the Southeast of the Netherlands 1955 – 1994: a report from the Eindhoven Cancer Registry. Eindhoven, the Netherlands: Integraal Kankercentrum Zuid 1995.Google Scholar
- Coebergh JWW, Janssen-Heijnen MLG, Louwman WJ, AC V: Cancer. Incidence, care and survival in the Southeast of the Netherlands 1955 – 1999: a report from the Eindhoven Cancer Registry (IKZ) with cross-border implications. Eindhoven, the Netherlands: Integraal Kankercentrum Zuid 2001.Google Scholar
- Maas IAM, Gijsen R, Lobbezo IE, Poos MJJC: Volksgezondheid Toekomst Verkenning 1997. I De gezondheidstoestand: een actualisering. Bilthoven: Rijksinstituut voor Volksgezondheid en Milieu 1997.Google Scholar
- Visser O, Coebergh JWH, Schouten LJ, van Dijck JAAM: Incidence of Cancer in the Netherlands 1995. Utrecht: Vereniging van Integrale Kankercentra 1998.Google Scholar
- Visser O, Coebergh JWW, Schouten LJ, van Dijck JAAM: Incidence of cancer in the Netherlands 1997. Utrecht: Vereniging van Integrale Kankercentra 2001.Google Scholar
- Coebergh JW, Crommelin MA, Kluck HM, van Beek M, van der Horst F, Verhagen-Teulings FT: Breast cancer in southeast North Brabant and in North Limburg; trends in incidence and earlier diagnosis in an unscreened female population, 1975–1986 [in Dutch]. Ned Tijdschr Geneeskd 1990, 134: 760-765.PubMedGoogle Scholar
- Brown BW, Brauner C, Minnotte MC: Noncancer deaths in white adult cancer patients. J Natl Cancer Inst 1993, 85: 979-987.View ArticlePubMedGoogle Scholar
- Louwman WJ, Klokman WJ, Coebergh JW: Excess mortality from breast cancer 20 years after diagnosis when life expectancy is normal. Br J Cancer 2001, 84: 700-703. 10.1054/bjoc.2000.1616View ArticlePubMedPubMed CentralGoogle Scholar
- Rutqvist LE: On the utility of the lognormal model for analysis of breast cancer survival in Sweden 1961–1973. Br J Cancer 1985, 52: 875-883.View ArticlePubMedPubMed CentralGoogle Scholar
- Schouten LJ, Hoppener P, van den Brandt PA, Knottnerus JA, Jager JJ: Completeness of cancer registration in Limburg, The Netherlands. Int J Epidemiol 1993, 22: 369-376.View ArticlePubMedGoogle Scholar
- Mackenbach JP, Van Duyne WM, Kelson MC: Certification and coding of two underlying causes of death in The Netherlands and other countries of the European Community. J Epidemiol Community Health 1987, 41: 156-160.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.