 Methodology
 Open Access
 Published:
Methodological issues in estimating survival in patients with multiple primary cancers: an application to women with breast cancer as a first tumour
Emerging Themes in Epidemiology volume 6, Article number: 2 (2009)
Abstract
Background
Comparing survival of patients with a single tumour and patients with multiple primaries poses different methodological problems. In population based studies, where we cannot rely on detailed clinical information, the issue is disentangling the share of survival probability from the first and second cancer, and their compounded effect. We examined three hypotheses: A) the survival probability since the first tumour does not change with the occurrence of a second tumour; B) the probability of surviving a tumour does not change with the presence of a previous primary; C) the probabilities of surviving two subsequent primary tumours are independent (additivity hypothesis on mortality rates).
Methods
We studied the survival probabilities modelling mortality rates according to hypotheses A), B) and C). Mortality rates were calculated using AalenJohansen estimators which allowed to discount for the lagtime survival before developing a second tumour. We applied this approach to a cohort of 436 women with breast cancer (BC) and a subsequent tumour in the resident population of Turin, Italy, between 1985 and 2002.
Results
We presented our results in term of a Standardised Mortality Ratio calculated (SMR_{ AJ }) after 10 years of followup. For hypothesis A we observed a significant excess mortality of 2.21 (95% C.I. 1.94 – 2.45). Concerning hypothesis B we found a not significant SMR_{ AJ }of 0.98 (95% C.I. 0.87 – 1.10). The additivity hypothesis (C) was not confirmed as it overestimated the risk of death, in fact SMRs_{ AJ }were all below 1: 0.75 (95% C.I. 0.66 – 0.84) for BC and all subsequent cancers, 0.72 (95% C.I. 0.55 – 0.94) for BC and colonrectum cancer, 0.76 (95% C.I. 0.48 – 1.14) for BC and corpus uteri cancer (not significant).
Conclusion
This method proved to be useful in disentangling the effect of different subsequent cancers on mortality. In our application it shows a worse longterm mortality for women with two cancers than that with BC only. However, the increase in mortality was lower than expected under the additivity assumption.
Introduction
The improvement of patients survival for the vast majority of neoplasms led to a substantial increase in the probability of further developing subsequent primary tumours. However, the study of multiple primary tumours on a population basis posed many additional problems. There is, indeed, a problem of differential diagnosis, when it comes to distinguish between local and distant metastases, recurrences and the onset of a truly new lesion. Classifications may also vary leading to substantial differences in rates. For example Surveillance Epidemiology and End Results (SEER) rules [1] differ substantially from those adopted by International Agency for Research on Cancer (IARC) [2].
Furthermore, survival of patients with multiple tumours has been neglected in populationbased analyses, where they are usually listwise deleted, or analysed for the first tumour occurrence only [3, 4]. Only recently two studies [5, 6] reconsidered this exclusion policy. On the contrary, in clinical series survival of patients with multiple tumours is usually defined clinically and specific cause of death is assessed accordingly. However, in population studies and in series from cancer registries, clinical information on patients followup is often unavailable and assessment of cause of death is based only on death certificates, often liable to gross misclassification. Heinävaara et al. [7] proposed to estimate the differential amount due to first or second tumour with a statistical parametric model. Their application dealt with patients with two primary breast cancers, where the question of disentangling the cancerspecific survival due to the first or the second tumour is more difficult, also from a clinical point of view. In the case of a subsequent primary cancer of a different origin the question is apparently simpler, although not yet investigated on a population basis.
The following questions can be raised: whether the overall survival of patients has decreased because of the interaction between the two cancers, or if it has been left substantially unchanged in comparison to those with one cancer only, or even increased. For example, active surveillance and care due to the first cancer can lead to earlier diagnosis of subsequent cancers and therefore to a longer survival (or a longer lead time). However, before studying the possible effect of surveillance and other prognostic factors (which was not the aim of this study), we should focus on the correct measurement of survival, which is our research objective.
To achieve this, we had to face many complex methodological challenges: first, we had to fix the zero reference time (the time from when we started the followup); second, a person can die only once thus the background death rate is confounded in the follow up information after the diagnosis of the second primary, therefore it is crucial to use models able to suitably describe a situation of competing risks; third, in order to make inferences, for each model we had to define the correct expected survival based on the appropriate comparison group.
We focused our attention on the following questions:

A.
Does the survival probability of a patient with a second primary tumour differ from those with only first type of tumour?

B.
Does the survival probability of a patient with a second primary tumour differ from those with only second type of tumour?
If a difference in survival is found in some of the previous situations, a third more fundamental question arises.

C.
Are the probabilities of surviving two subsequent primary tumours independent?
Studying survival probabilities in terms of the underlying hazard of death, the question can be rephrased as follow:
Is the mortality rate after a second tumour simply the sum of the two intensities (additivity hypothesis), or the way the mortality rates act follows a different functional law?
This paper aims at answering these questions for women with breast cancer and a subsequent primary tumour, paying particular attention to the conditional survival probability due to the time elapsed between the two malignancies.
Methods
Statistical analysis
To correctly defining the probability of surviving conditional to be alive up to the occurrence of a second tumour, we started by writing questions A, B and C as hypothesis in term of mortality hazard. We defined:
λ_{ A }(t): mortality rate for the population with two tumours at time t from the occurrence of the first tumour;
λ_{ B }(t): mortality rate for the population with two tumours at time t from the occurrence of the second tumour;
λ_{C, α}(t): mortality rate at time t from the occurrence of the second tumour for the population with a second tumour given that they already survived a time interval α.
We can break these down as
where, for i = 1, 2, λ_{i0}is the specific mortality rate at time t from the occurrence of tumour i for the population with only that tumour, and λ_{0} is the general mortality.
We assumed that λ_{10}, λ_{20} and λ_{0} were known, by previous studies on mortality and survival in population with the first type of tumour only, with the second type of tumour only, and in the general population, respectively.
We observed that λ_{12}(t) was the possible difference in mortality rate in patients with a tumour of type 1 followed by a tumour of type 2 with respect to that of patients with a tumour of type 1 only, measured from the occurrence of tumour 1; λ_{21}(t) was the possible difference in mortality rate in patients with a tumour of type 1 followed by a tumour of type 2 with respect to that of patients with a tumour of type 2 only, measured from the occurrence of tumour 2.
Questions A, B, and C can be written as follows:
Occurrence probabilities conditioned to different events (occurrence of a second cancer, death) in each time interval can be estimated with the AalenJohansen [8] (AJ) method in the framework of a Markov process, as described later. Once we obtained these conditional probabilities, we calculated the number of expected deaths by sex, age, calendar period and followup time, under the different hypotheses A, B and C. From a practical point of view, we calculated the expected deaths in a similar way to that used to calculate the denominator of relative survival [9]. For example, in the case of a woman diagnosed with breast cancer at 62 who developed a rectal cancer after two years and survived for an additional period of five years, we associated an expected probability of dying with a breast cancer, occurred at the same age, for the two years elapsed with that cancer only. Subsequently, we associated an expected probability of dying with breast and/or with rectal cancer for the following years, taking into consideration the ageing of the patient (i.e. using the annual probability of dying according to the age of the patient, from age 64 to age 69). The way the calculation of the expected number of death (or the expected probability of dying) for the conjoint period when both tumours are present is performed depends on which one of the three hypotheses we are testing. If we consider hypothesis A, we do not add the probability of dying with a colonrectum cancer. If we test hypothesis B, we do not add the probability of dying associated to a breast cancer for the first period. Finally, if we test hypothesis C (additive hypothesis), we sum the two underlying mortality hazards during the second period. Expected probabilities were derived from analyses of the cohort of patients with only one incident cancer included in the cancer registry's data.
For the interested reader, we now explain in details how we calculated expected probabilities. Since different states are concerning, we resorted to the theory of Markov models [8]. In a Markov process individuals can belong to a finite set of states and move to one state to some others with a probability, possibly depending on time. The main hypothesis is that the probability of moving from state i to state j at time t depends on i, j and t only, and not on the previous history of the individual.
We constructed a simple model with three states

1
first tumour

2
second tumour

3
death after a first (but not a second) tumour
where 2, and 3 are absorbing states and the possible moves are: 1 → 2, 1 → 3.
Since our data showed right censoring, transition probabilities P_{ ij }(s, t) from state i to state j, in the time interval (s, t) were calculated using AalenJohansen (AJ) estimators [8].
The procedure we adopted included age standardisation, and precisely:

For each age class k we calculated the AJ estimator P_{ ijk }(s, t). We let N_{ k }be the number of subjects in class k at time 0 and we set a weight ${W}_{k}={\scriptscriptstyle \frac{{N}_{k}}{N}}$, where N equals the sum of the N_{ k }'s.

We defined the standardised estimator as: ${P}_{12}^{stand}(s,t)={\displaystyle {\sum}_{k}{W}_{k}}\cdot {P}_{12k}(s,t)$

It is reasonable to assume that weights are deterministic (fixed) variables; under this assumption we have:
$$var({P}_{12}^{stand}(s,t))={\displaystyle \sum _{k}{W}_{k}^{2}\cdot var({P}_{12k}(s,t))}.$$
Then, from probabilities previously calculated with AJ estimators, it was possible to compare observed mortality with mortality expected in the hypothesis of no interaction between the two tumours; that is the mortality intensities due if the two tumours were independent. As a consequence, the number of expected deaths is the sum of the deaths due to mortality for both tumours acting separately. We calculated the number of expected deaths considering for each patient j the time of occurrence of the first primary malignancy T_{1j}, time of occurrence of the second primary malignancy T_{2j}, and, most important, the time interval between the occurrence of the two tumours α_{ j }= T_{2j} T_{1j}. Each patient, after a time interval t_{2} since the inception of the second tumour, has a probability p_{2j}(t_{2}) of dying for the second tumour or general mortality equal to that of the general population of patients with only that type of tumour, according to her/his age, sex, calendar period of diagnosis and followup time. In addition, that patient has a probability p_{1j}(t_{2} + α_{ j }) of dying at the (t_{2} + α_{ j }) time interval for the first tumour or general mortality again equal to that of the general population of patients with that type of tumour only, according to her/his age, sex, calendar period of diagnosis and followup time.
We set ${\widehat{p}}_{2j}$(t_{2}) = p_{2j}(t_{2})·(1  p_{0j}) where p_{0j}is the general mortality of the subject j according to her/his age, sex, calendar period of diagnosis, taken from the life tables of the general population. Thus, we can say that ${\widehat{p}}_{2j}$(t_{2}) is the specific mortality for the second tumour.
Therefore, within the cohort of patients with two malignancies, at the t_{2} time interval since the second tumour, in the hypothesis of no interaction between the mortality forces of the two tumours, we expect the following number of deaths:
where the probability of dying for the second tumour ${\widehat{p}}_{2j}$(t_{2}) is corrected by the probability of surviving from the first tumour and general mortality 1  p_{1j}(t_{2} + α_{ j }).
Since the output of these calculations was the number of expected deaths, we consequentially compared it with the observed number in a ratio similar to the well known Standardised Mortality Ratio:
We used the term SMR_{ AJ }because it was quite similar to the standard term "SMR"in the sense that it was that ratio of observed to expected deaths; the expected deaths were calculated as a sum over age groups; and finally, it was similar to the indirect method of age standardisation since, as standard, we applied the mortality rates of the cohort of patients with only one tumour.
Patients
We selected all incident breast cancer cases recorded by the Piedmont Cancer Registry in the resident women of Turin from 1985 to 1998. This cohort of patients was followed up for four years until the end of 2002, both for what concerns life status or development of a subsequent tumour (excluding skin carcinoma). Life status of women who emigrated outside the resident population observed by the Piedmont Cancer registry was ascertained with an active followup at the municipality rosters of the new residency. In this analysis we considered women with cancer of corpus uteri or cancer of colonrectum as second primary tumour, since we observed a consistent number of cases (91 for colonrectum and 62 for corpus uteri) for making reasonable stable estimates. In addition, we analysed all types of second primary cancers (escept skin carcinomas) including corpus uteri and colonrectum. For age standardisation, we introduced age at diagnosis in five broad classes: 0 – 44, 45 – 54, 55 – 64, 65 – 74, 75+. Age standardisation for the unconditional survival estimates were calculated using those standards proposed by Corazziari et colleagues [10] for comparisons in international studies.
We calculated λ_{10} from our cohort of 8234 women with breast cancer only. Mortality for the second type of tumour (λ_{20}) was calculated from 1443 women with corpus uteri cancer only, and from 4050 women with colonrectum cancer only. In the case of mortality for all cancers λ_{20} was calculated in two ways: including breast cancers (28737 women), and excluding breast cancers (20082 women). We also used overall mortality including breast cancer, as a reference for comparing available published statistics that usually do not make exclusions for specific type of cancers. Life tables for the general mortality were from the Statistics Office of Turin for the period 1985–2002.
Results
We identified 9233 women with breast cancer in Turin from 1985 to 1998, 563 cases were excluded as they were identified from Death Certificate Only (DCO) or synchronous primary tumours (same day of diagnosis for the two tumours) or they were patients with more than two tumours, leaving 8670 cases for analysis. From this cohort, 436 second (metachronous) primary tumours developed during the prolonged followup period (1985–2002). The most frequent types of subsequent cancers were colonrectum (91 patients) and corpus of uteri (62 patients) (Table 1). The median interval time between the two tumours was five years, but with an highly positive skewed distribution. The completeness of clinical documentation was rather high, considering that registries work on a population basis, with a 94.5% of microscopic confirmation for first tumours that reached 99.3% for second tumours.
At the end of the study period (2002), we observed 285 (65.4%) deaths among women with a second tumour, distributed as in table 1, and 3931 (47.7%) among women with breast cancer only.
Using a traditional approach, starting from the diagnosis of the first tumour (breast cancer), we had, at 1, 5 and 10 years respectively, an age standardised survival of 96.1%, 74.3% and 43.6% for the 436 women with two cancers, while it was 94.0%, 72.9% and 56.9% for women with breast cancer only (Table 2). On the contrary, survival of women with one tumour only (all site including breast cancer) was worse: 69.5%, 47.5%, and 37.1% after 1, 5 and 10 years of followup. Also survival of women with colonrectum or corpus uteri cancer only was worse than that of women with a subsequent cancer, when starting followup since diagnosis of breast cancer. Their higher survival is due to the fact that patients with multiple cancers had survived an extra amount of time (a median of five years) before developing the subsequent cancer in comparison with patients with a single cancer only. A less biased and more appropriate comparison can be obtained moving the zero time point to the date of diagnosis of the second tumour. In this case, the survival of patients with breast cancer and a subsequent tumour was rapidly decreasing since the date of diagnosis of the second tumour, and comparable to that of patients with only one tumour excluding breast cancer. Also survival of patients with breast and colonrectum cancer was rapidly decreasing after the first year, when it was still comparable, and lower (20.8% versus 32.9%), although not statistically significantly, than that of patients with colonrectum cancer only. Survival of patients with breast and corpus cancers was generally quite high, with 70.6% of women still alive after 10 years since the breast cancer and 44.2% alive after 10 years since the corpus uteri cancer diagnosis. However, patients with corpus cancer only showed a higher survival at 10 years (56.6%).
In table 3 we compared the number of observed deaths with the expected deaths under the three hypotheses A, B and C. We calculated the number of expected deaths taking into account the probabilities of surviving the first tumour, developing a second tumour and subsequently dying. We presented the expected deaths broken down according to the three different probabilities. The expected numbers of deaths under the first hypothesis (A) were definitively below the observed numbers since they did not take into account the probability of dying because of the second tumour. On the other hand, the expected number of deaths from the unconditional probabilities of dying for the second tumours (hypothesis B) were overestimated (291.8) in the case of all cancers and slightly underestimated in the case of colonrectum (53.45) and corpus uteri cancers (19.16). The expected numbers of deaths according to the probabilities of dying because of the second tumour, given that the patients survived the first, were 248.30 for all cancers, 44.76 for colonrectum cancer and 15.81 for corpus uteri. Summing up the latter figures with the expected numbers for hypothesis A, we obtained the expected number of deaths according to the additivity hypothesis C. Considering the ratio between observed and expected numbers (SMR_{ AJ }), the additivity hypothesis is not confirmed for all cancers (SMR_{ AJ }0.75 95% CL 0.66–0.84), or for colonrectal cancer (SMR_{ AJ }0.72 95% CL 0.55–0.94) while in the case of corpus uteri cancer we cannot reject the hypothesis that mortality rates are additive. The differences between the three models and observed deaths is visualised in Figure 1 where cumulative numbers of deaths are displayed over years of followup. It can be seen that in all graphs hypothesis C (additivity) tends to overestimate the actual observed trend, hypothesis A strongly underestimates it, while hypothesis B is the closest to observed data, especially during the first time (years) of followup.
Discussion
The dramatic improvement of cancer survival during the last decades in Western countries brought with it a new health threat: the development of second primary cancers in survivors. An editorial in CEBP of David Alberts clearly stated that 'Second cancers are killing us! ' [11]. However, in spite of the fact that several studies on the multiple primary cancer risk were undertaken [12], the rate at which first and second, or higherorder cancers are killing us remains neglected. In clinical studies, when reliable information are available it is often possible to understand if the pathological conditions linked to a specific cancer affected the patient survival and to which extent. However, at a population based level this is often not feasible due to the lack of clinical information or cause of death. Even when cancerspecific causes of death are available, they are subject to various degrees of misclassification, hindering the possibility of a reliable estimate of cancer specific survival. In the main population based statistics on cancer survival worldwide available (Eurocare [3] and SEER [4]) subsequent cancers were excluded: only the first occurring cancer was analysed, or all the subjects with multiple cancers were deleted from analysis. Although, this strategy has recently undergone through a rethinking [5, 6], it was supposed to allow for more comparable results across registries with different back up information, and therefore with a different possibility in identifying those cancers that occurred in prevalent cases. However, we believe that the problem deserves more attention also from its implication in the management and care of such patients. Indeed, a wider availability of effective cancer treatments has prolonged patient survival, so increasing the possibility of developing another cancer. Studying the occurrence of multiple tumours and their association, it emerged as the higher susceptibility to subsequent malignancies can possibly be due to unfavorable genetic pattern or common exogenous risk factors [13, 14]. Multiple cancer survival is also a stimulating topic of study, but received less attention. Recently, an analysis of the SEER data on multiple tumours following breast cancer [15] showed that survival of women 20–29 years old at time of breast cancer diagnosis had a worse 10year survival, compared with women with breast cancer only, while there were no differences in the 5year survival. However, in that analysis the time elapsed until the second cancer occurrence was not taken into account.
Before investigating the reasons influencing survival for patients with multiple tumours, we, indeed, believe that it is essential to have a correct measurement of survival that takes into account the effect of conditional probabilities of surviving given the different timing of primary cancers occurrences. We proposed a method that assigns the correct number of expected events according to the different components of mortality due to each type of cancer. The proposed method is useful only in correctly stating the prediction of mortality probabilities while cannot explain the causes of the different mortality probabilities.
The expected number of deaths was calculated taking into account the exact time spent at risk of dying for one or another cancer by age classes and calendar period, using conditional probabilities estimated by the AJ estimator from a simple oneway Markov process with two absorbing states. Such approach was recommended since it allowed a better control of probabilities of events arising from different states. In the model referring to hypothesis A, we calculated the expected number of deaths due to the first occurring cancer starting since its time of occurrence. This model is similar to model 2 proposed by Heinävaara and colleagues [7] in the absence of cancer specific cause of death. We wrote the model's parameters in terms of risk excess (hazard rate), rather than estimating the specific mortality rates. While survival of patients with a second primary tumour was comparable or higher with that of those patients with breast cancer only during the first years, it was rapidly declining at a higher rate than the reference group after five years of followup. This effect was explained by the fact those patients had survived an extra amount of time (a median of five years) before developing the subsequent cancer. Indeed, results from hypothesis A showed an increased cumulative mortality only at ten years for women with two cancers when compared to those with breast cancer only, as found in the study of Raymond and Hogue [15].
The second model (hypothesis B) was built with the same structure as model A, calculating the expected number of deaths due to the second occurring cancer starting since its time of occurrence. However, the change in the baseline population and the shift in the time zero reference made the hazard rates not comparable. Indeed, for a proper comparison with those patients with the second type of cancer only, we set the starting time at the diagnosis of the second cancer. In this case, the survival was comparable at 1 and 5 years of followup, than that of patients with one type of cancer only, while it was slightly shorter at 10 years. In summary, results from hypothesis B showed no extra mortality compared to patients with only one cancer of the same type, and observed and expected number of deaths closely get on during the years of observation.
We then addressed the question of evaluating the eventual extra mortality due to the combination of effects of the two primary neoplasms, checking the hypothesis if the mortality of women with two cancers was due to the sum of the baseline mortality rates of breast and other cancers (additivity hypothesis C). It clearly emerged how observed cumulative mortality was lower than expected under the additivity assumption, with a statistically significant difference in the case of all cancers and colonrectum after 10 years of followup. The agreement of a specific model to observed data was therefore useful for having further hints of the underlying mechanisms. In our study, the less than expected results can be explained by the fact that the second cancer can have a less advanced stage and therefore a better prognosis, since a subsequent cancer is usually diagnosed because of a deeper clinical surveillance due to the first cancer. It is clear that women with breast cancer and a subsequent cancer survive less than women with breast cancer only, but their survival is not always decreased simply as it would be if the forces of mortality work together in an additive way.
The study has some possible limitations. First of all, the method of correction is based on observed rates (mortality rates measured in the cohort of patients with only one tumour) that, when based on small numbers, can be unstable. Then, this method, being inherently nonparametric, does not give information on the underlying incidence/mortality competing laws. In calculating expected number of deaths a possible bias could have been introduced, depending on the numbers of patients who emigrated outside the Cancer Registry's area. In this case, information on life status were still available and collected, but we did not know if the patient had developed a subsequent cancer when resident in another area. During the study period, we observed about 8% of women who emigrated among those classified with breast cancer only. Their median time of emigration was 6.5 years since the breast cancer diagnosis. As a consequence, considering that the median time for developing a second primary cancer was about 5 years, the detection bias should be very limited. Finally, the method was presently tested only on a limited set of data: patients with breast cancer as a first primary tumour. As few studies are still available on this topic, more research is needed, with larger samples and including clinical data (e.g. stage at presentation, hormone receptor status), therapies (e.g. tamoxifen), information on followup circumstances, and modality of diagnosis. In conclusion, we showed that the presented approach for calculating conditional probabilities was correct when dealing with situations, as with multiple tumours, where competing causes of death can bias the results of survival probabilities. We also pointed out how shifted reference times can be considered in correctly comparing survival. In addition, departure from the expected additive model can give hints towards which direction to further investigate.
References
 1.
Jonhson C: The SEER coding and staging manual 2004. NIH Pub No 045581. Fourth Edition Bethesda, MD: National Cancer Institute; 2004.
 2.
IARC/IACR: International rules for multiple primary cancer (ICDO Third edition), Volume 2004/02 of Internal Reports. 2004, Lyon: International Agency for Research on Cancer; 2004.
 3.
Berrino F, Sant M, Verdecchia A, Capocaccia R, Hakulinen T, Estéve J: Survival of cancer patients in Europe. The EUROCARE Study Volume 132. IARC Scientific Publications. Lyon: International Agency for Research on Cancer; 1995.
 4.
Ries L, Harkins D, Krapcho M, Mariotto A, Miller B, Feuer E, Clegg L, Eisner M, Horner M, Howlander N, Hayat M, Hankey B, Edwards B: SEER Cancer Statistics Review, 1975–2003. November 2005 SEER data submission, posted to the SEER web site 2006. Bethesda, MD: National Cancer Institute; 2006.
 5.
Brenner H, Hakulinen T: Patients with previous cancer should not be excluded in international comparative cancer survival studies. Int J Cancer. 2007, 121: 22742278. 10.1002/ijc.22932
 6.
Rosso S, De Angelis R, Ciccolallo L, Carrani E, Soerjomataran I, Grande E, Zigon G, Brenner H, : Multiple tumours in survival estimates. Eur J Cancer. 2009.
 7.
Heinävaara S, Teppo L, Hakulinen T: Cancerspecific survival of patients with multiple cancers: an application to patients with multiple breast cancers. Statistics in Medicine. 2002, 21: 31833195. 10.1002/sim.1247
 8.
Andersen P, Borgan O, Gill R, Keiding N: Statistical models based on counting processes. In Springer Series in Statistics SpringerVerlag; 1993
 9.
Hakulinen T: Cancer survival corrected for heterogeneity in patients withdrawal. Biometrics. 1982, 38: 993942. 10.2307/2529873. 10.2307/2529873
 10.
Corazziari I, Quinn M, Capocaccia R: Standard cancer patient population for age standardising survival rates. Eur J Cancer. 2004, 40: 23072316. 10.1016/j.ejca.2004.07.002
 11.
Alberts D: Second cancers are killing us! Cancer Epidemiol. Biomarkers Prev. 2006, 15: 201910.1158/10559965.EPI060417. 10.1158/10559965.EPI060417
 12.
Travis L: The epidemiology of second primary cancers. Cancer Epidemiol Biomarkers Prev. 2006, 15: 20202026. 10.1158/10559965.EPI060414
 13.
Crocetti E, Buiatti E, Falini P, the Italian Multiple Primary Cancer Working Group: Multiple primary cancer incidence in Italy. Eur J Cancer. 2001, 37: 24492456. 10.1016/S09598049(01)003148
 14.
Evans H, Lewis C, Robinson D, Bell C, Møller H, Hodgson SV: Incidence of multiple primary cancers in a cohort of women diagnosed with breast cancer in southeast England. Br J Cancer. 2001, 84: 435440. 10.1054/bjoc.2000.1603
 15.
Raymond J, Hogue C: Multiple primary tumours in women following breast cancer, 1973–2000. Br J Cancer. 2006, 94: 17451750.
Acknowledgements
This research was partly founded by Regione Piemonte – Ricerca Sanitaria Finalizzata (years 2004 and 2006).
Preliminary results were presented by Fulvio Ricceri at the GRELL meeting 2006 in Palma de Mallorca awarding the "Enrico Anglesio Prize "offered by the "Anglesio/Moroni Foundation "of Turin, Italy. We thank the researchers and professors of the Me.Ri.Ma. group of the University of Turin (Department of Mathemathics) who shared their ideas with us and gave us their time and comments.
We also thank Federica and Simona Gallo for their editorial assistance.
The authors declare that they have no competing economic or financial interests.
Author information
Affiliations
Corresponding author
Additional information
Authors' contributions
SR conceived the idea for the study. SR, FR and LT planned and designed the research. FR developed the statistical models. LT revised the mathematical foundations of the proposed model. FR and LT analysed the data. SR and RZ wrote the first draft of the manuscript. RZ coordinated this project. All authors edited and approved the final version of the manuscript. The corresponding author has final responsibility to submit for publication.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
Rosso, S., Ricceri, F., Terracini, L. et al. Methodological issues in estimating survival in patients with multiple primary cancers: an application to women with breast cancer as a first tumour. Emerg Themes Epidemiol 6, 2 (2009). https://doi.org/10.1186/1742762262
Received:
Accepted:
Published:
Keywords
 Breast Cancer
 Breast Cancer
 Standardise Mortality Ratio
 Calendar Period
 Multiple Primary Cancer