- Research article
- Open Access
- Open Peer Review
On the use of flexible excess hazard regression models for describing long-term breast cancer survival: a case-study using population-based cancer registry data
© The Author(s). 2019
- Received: 29 May 2018
- Accepted: 14 January 2019
- Published: 28 January 2019
Breast cancer prognosis has dramatically improved over 40 years. There is, however, no proof of population ‘cure’. This research aimed to examine the pattern of long-term excess mortality due to breast cancer and evaluate its determinants in the context of cancer registry data.
We used data from the Geneva Cancer Registry to identify women younger than 75 years diagnosed with invasive, localised and operated breast cancer between 1995 and 2002. Flexible modelling of excess mortality hazard, including time-dependent (TD) regression parameters, was used to estimate mortality related to breast cancer. We derived a single “final” model using a backward selection procedure and evaluated its stability through sensitivity analyses using a bootstrap technique.
We analysed data from 1574 breast cancer women including 351 deaths (22.3%). The model building strategy retained age at diagnosis (TD), tumour size and grade (TD), chemotherapy and hormonal treatment (TD) as prognostic factors, while the sensitivity analysis on bootstrap samples identified nodes involvement and hormone receptors (TD) as additional long-term prognostic factors but did not identify chemotherapy and hormonal treatment as important prognostic factors.
Two main issues were observed when describing the determinants of long-term survival. First, the modelling strategy presented a lack of robustness, probably due to the limited number of events observed in our study. The second was the misspecification of the model, probably due to confounding by indication. Our results highlight the need for more detailed data and the use of causal inference methods.
- Breast cancer
- Prognostic factors
- Excess hazard model
- Sensitivity analyses
- Cancer registry
Breast cancer is a major disease worldwide. Its prognosis has, however, improved rapidly during the last four decades [1–3]. Accordingly, there are increasing numbers of women who have survived breast cancer. Despite this, there is evidence for a lack of population ‘cure’, that is, the probability of dying as a consequence of the disease persists for many years after diagnosis [4, 5] even for women who were screen-detected .
The estimation of net survival has allowed these trends to be observed [7–9]. Unlike all other metrics, net survival evaluates the mortality arising only from the disease of interest, disregarding the influence of other causes of death . In the context of long-term survival this is fundamental because the likelihood of death from other causes increases with follow-up time (i.e. with ageing of the patients). The use of net survival allows accurate comparisons of patient’s subgroups across space and time, between which mortality from other causes may vary considerably [9, 11].
Although there is a great interest, both clinically and epidemiologically, in the determinants of long-term survival for breast cancer patients, follow-up beyond 5 or 10 years has not been widely considered. The few studies with long-term observations have demonstrated that the associations of some covariables do vary with time since diagnosis in the long term , very few however considered very long-term follow-up and/or multiple covariables [13, 14]. In particular, the influence of treatment represents an interesting line of investigation since it is likely that certain treatments lead to severe long-term side effects [13, 15].
The Geneva Cancer Registry offers an ideal context to study the evaluation of determinants of long-term net survival. The cancer registry, initiated in 1970, allows long follow-up of cancer patients. The availability of detailed information for each woman’s tumour enables multivariable survival analysis.
In this research, we aim to evaluate long-term associations between prognostic factors and the excess mortality hazard for breast cancer patients diagnosed in Geneva, focusing especially on treatment variables. To reach this aim, we focus on early-stage tumours, which were surgically resected. We use flexible excess hazard regression models to account for potential time-varying and non-linear associations. We apply a systematic model selection process to build a “final” regression model, and check the stability of our “final” model by conducting a sensitivity analysis using bootstrap sampling.
The Geneva Cancer Registry collects information on incident cancer cases from various sources, including hospitals, laboratories and private clinics, all of whom report newly diagnosed cancer cases. Trained registrars systematically extract information from the medical records and conduct further investigations in the case of missing data. The registry regularly estimates cancer patient survival, taking as the reference the date the diagnosis was confirmed or, if it preceded the diagnosis and was related to the disease, the date of hospitalisation. In addition to standard examination of death certificates and hospital records, patient’s vital status is assessed annually by matching the Registry’s database with information held by the Cantonal Population Office which maintains a live register of the resident population.
We included all women diagnosed with an invasive primary breast cancer in the Geneva Canton between 1995 and 2002. We restricted the sample to patients diagnosed with pathological TNM stage I and II disease who were treated with surgery (N = 2029). Among those patients, we excluded patients older than 75 years (N = 232). Information on stage was missing among only 60 (2.57%) patients with surgery. All women were followed up until 31st December 2013 (11 years of minimum follow-up).
The Geneva Cancer Registry has a general authorization to collect nominative data, and to analyze the anonymized data. Since the study did not require additional clinical information, approval of the Ethics Committee was not required.
Prognostic factors and treatment
We focused on established prognostic factors and on treatment. Age at diagnosis (years) was included a priori as an irrefutable prognostic factor . We considered tumour size (mm), degree of differentiation (Well vs. Moderately/Not differentiated), nodal involvement (No vs. Yes) and hormone receptor status (Negative vs. Positive) which together reflect the severity of the disease. We included radiotherapy, chemotherapy and hormonal treatment following surgery (each Yes vs. No, within 6 months after diagnosis) in order to examine the long-term associations of these systemic treatments with survival.
Statistical modelling of the excess mortality hazard
We estimated the excess mortality hazard due to cancer for the patient group. The excess mortality hazard corresponds to the mortality hazard related only to the disease of interest (in our case, breast cancer) and is defined as the difference between the mortality observed amongst a cohort of patients and their expected (background) mortality [17, 18]. The association between covariables and excess mortality can vary with time since diagnosis, particularly when considering long-term follow-up. For example, a particular treatment might have a strong influence on excess mortality one year after diagnosis but a weaker influence at ten years (time-dependent, TD, association). Furthermore, continuous variables can display non-linear (NL) associations (for example, excess mortality might increase exponentially with age). We handled such complex associations through the flexible excess hazard model proposed by Charvat et al. , which follows the work of Remontet et al . This excess hazard model is implemented in the “mexhaz” package written for R software [19, 21].
Model building strategy
age at diagnosis (continuous, NL and TD associations included),
tumour size (continuous, log-transformed, NL and TD associations included),
nodal involvement (binary, TD association included, “Yes” as reference category),
grade of the tumour (binary, TD association included, “Moderately/Not differentiated” as reference category),
hormone receptor status (binary, TD association included, “Positive” as reference category),
radiotherapy (binary, TD association included, “No” as reference category),
chemotherapy (binary, TD association included, “No” as reference category) and
hormonal treatment (binary, TD association included, “No” as reference category).
The model building strategy resulted in a single derived model (the “final” model), which included only those variables found to be significant, along with any significant TD and/or NL associations for these variables.
We conducted a sensitivity analysis to examine the stability of the derived “final” model using a bootstrap technique . This involved re-applying the model selection procedure to 300 random samples, drawn, with replacement, from the cohort of cancer patients. This procedure allows the evaluation of the strength of association between a particular covariable and the excess mortality hazard by the calculation of the bootstrap inclusion frequency (BIF) . The BIF is the proportion of times a specific variable was included by the model selection process over the total number of bootstrap samples. We further considered only models where the association was plausible (outliers where estimated values of the Excess Hazard Ratio (EHR) were greater than 100 or less than 0.01 were excluded). We then plotted all the estimated functional forms of each covariable (N ≤ 300 due to convergence issues, see below), along with the averaged functional form calculated on all the retained samples.
Characteristics of the patients diagnosed with breast cancer between 1995 and 2002
Size in mm
Moderately/ poorly differentiated
“Final” derived model
The single derived model resulting from the model building strategy included associations with the excess mortality hazard which were linear and TD for age, linear and time-constant for size of the tumour, TD for grade of the tumour, time-constant for chemotherapy and TD for hormonal treatment. All the other variables (nodal involvement, hormone receptor status and radiotherapy) were considered not associated with the excess mortality hazard and were therefore excluded.
Stability of the “final” derived model
Bootstrap Inclusion Frequency (BIF) for each co-variable and their type of associations following the sensitivity analysis
Size of the tumour
Grade of the tumour
The determinants of long-term survival are currently of particular interest because of the dramatic increase in the number of patients surviving breast cancer matched to the observation that these women are never ‘cured’. Understanding the impact of prognostic factors and of treatment with time since diagnosis is therefore increasingly important. In this context, population-based data are crucial to understand how treatment influences the outcomes for all cancer patients. These aims should be distinguished from those of web-based models, which provide patients with an estimation of his/her survival according to his/her values of prognostic factors.
In order to estimate the long-term associations of prognostic factors and treatment with the risk of dying from breast cancer, we used observational data from the population-based Geneva Cancer Registry. For this purpose, we restricted our cohort to a relatively homogeneous group of younger patients (less than 75) with localised disease (stage I and II) and who had received surgery. The severity of the disease was controlled for through the combination of several covariables, and the analyses accounted for differences in individual characteristics. Furthermore, the estimation of the mortality related to the disease, after controlling for other causes, was based on flexible excess hazard regression models, which enable the assumptions of linear and time constant excess hazard ratios to be relaxed. Both of these assumptions are clinically unlikely in the context of long-term survival. We used a recommended strategy  for selection of covariables and their complex associations, and performed a sensitivity analysis to evaluate the reproducibility of the model .
Despite using a fairly homogeneous group of patients, this optimised and up-to-date modelling strategy, a clear process for variable and complex association selection and a sensitivity analysis, our results demonstrated a lack of stability and model misspecification, associated with unrealistic effects of some treatments (e.g. chemotherapy).
First, our sensitivity analysis demonstrated that the set of covariables included (eventually with NL and/or TD functional forms) in the “final” model for the excess mortality hazard was unstable. Because of this demonstrated instability, results obtained from a single model should be interpreted with caution. This is best illustrated by the fact that 20% of models did not reach convergence during the sensitivity analysis, as well as the fact that several variables selected for the single derived model were rarely retained in the sensitivity analysis (low BIF, e.g. TD for hormonal treatment). Meanwhile others not retained in the derived model were often selected by the sensitivity analysis (high BIF e.g. hormone receptors).
There are a number of possible reasons for this lack of robustness. The first is related to the context in which the study was conducted. Since breast cancer patients present with high survival, the number of events (“excess” death) is relatively low in breast cancer data, even where long-term follow-up is available. This is especially true for the fairly small Geneva population (495,000 inhabitants) and for the study population which was restricted to early-stage cancer patients. It is recommended that at least 10 events per parameter should be included when estimating regression coefficients [25, 26]. Because we considered both time-dependent and non-linear associations for all prognostic variables, the number of parameters included in our model was large relative to the number of deaths. The convergence issues that we encountered are therefore likely to be explained, in part, by a lack of information from the observed data. However, decreasing the number of parameters (either by reducing the number of variables, or excluding some complex associations) would not have been a better strategy, given that our core aim was to try to better understand the long-term associations of prognostic covariables for breast cancer patients. Neither was it practical to increase the number of women in order to increase the number of events since this could only have been done by including women with advanced disease (for which treatment protocols are very different) or by including elderly women (who do not have the opportunity for long-term follow-up, and for whom the excess regression modelling would not make sense on the longer term .
The analysis excluded 12.3% of the cohort because of missing data, thus leading to a loss of information. However this proportion is relatively low for these types of observational data and complete-case analyses have been proved to be sufficiently efficient for such ranges of missing data proportion . Also, our aim was to highlight the difficulties encountered with modelling in the context of observational data. We therefore performed a complete-case analysis in order not to dilute the message with issues related to multiple imputation.
It is possible that the lack of stability may have been a result of the modelling approach. We consider this unlikely, however. The flexible regression model we applied has been purposefully designed to estimate excess mortality hazard  and take into account complex associations. The model selection strategy has previously been shown to be efficient and successful in detecting the correct complex associations as well as eliminating spurious ones .
The second main issue was that our strategy was unable to fully control for confounding by indication leading to model misspecification. This would be an issue even with a perfectly robust model. This confounding is best illustrated by the unexpected results for chemotherapy and hormonal treatment. Women receiving these treatments experienced an increased risk of dying from breast cancer compared to women who did not receive them (Fig. 3). This reflects the fact that the patients in the cohort who received chemotherapy and hormonal treatment were probably those with more advanced disease at diagnosis, among the early-stage cases (Table 1). This represents a limitation of our strategy, which was not able to account for the fact that almost all women who were likely to benefit from these therapies were given them, resulting in a sparse comparison group within the patient cohort (confounding by indication). We performed a stratified analysis to explore this (data not shown). We grouped patients with very similar characteristics together and compared their survival according to receipt of chemotherapy or not. This similarly showed an increased risk in the excess hazard of death associated with chemotherapy. This strongly suggests that additional information about the prognosis of patients not receiving chemotherapy is missing from our dataset, and that this led to misspecification of the model.
In addition, interactions between treatment received and other co-variables might be required. Although we planned to examine the existence of such interactions, they were tricky to implement due to the convergence issues we encountered during the modelling process, and not reasonable to explore in our small sample size dataset.
Other possible strategies
Our results point towards the need for different statistical strategies in addition to our modelling strategy to be better able to examine these effects more than only the associations. Causal inference analyses would be one suitable approach [27–29]. The objective of causal inference is to mimic the randomised trial that would have been set for the research question by using observational data and specific statistical techniques. Propensity score methods could, for example, be implemented within the flexible regression models we have used here [30, 31]. In our work, we assumed that people were treated at the date of diagnosis, which is probably not correct for all patients. Also, some changes in the prognostic factors values for some patients (e.g. growth of the tumor size) may suggest that a treatment needs to be undertaken later on after the diagnosis. In the presence of such time-varying confounding, other approaches such as parametric g-formula , structural nested models or marginal structural models with inverse probability weighting would also be of interest, especially for the long-term treatment effect [33, 34] . All these approaches assume the models to be well specified, which is not so easy to achieve. Various approaches, including using machine learning techniques, have been developed to minimise model misspecification . This would however require much more detailed data, including comorbidities and other factors used to define the treatment choice. Furthermore, software to implement causal inference techniques is not yet available for the excess mortality hazard. Further methodological research is thus required to enable such analyses to be conducted.
Nevertheless, a few cautious clinical interpretations can be drawn from these data. Some co-variables presented high BIFs within the sensitivity analysis and the observed associations appeared stable to the exclusion of outliers suggesting that they are indicative of a robust, underlying associations. Consistent with Jatoi et al.  we found that patients with negative hormone receptors presented a higher excess mortality during the first years after diagnosis compared to those who have positive hormone receptors (BIF 95.4%). Regarding age at diagnosis, our results matched those found by Cluze et al.  which showed the risk of dying from breast cancer was associated with increasing age at 1 and 5 years after diagnosis but that this association reversed at 10 years (BIF 87.1%). In addition to hormone receptor status and age at diagnosis, tumour size, grade and nodal involvement displayed associations which were similar to those described in a previous meta-analysis . Although our results are broadly consistent with previous studies, caution should be exercised in reporting the size of these associations, given that they have been derived from models, which display a lack of robustness. We observed a time-dependent association for radiotherapy: patients treated with radiotherapy exhibited a decreased excess mortality hazard in the first 10 years following their diagnosis but an increased hazard afterwards. This association was, however, sensitive to the inclusion or exclusion of outliers. That said, it could potentially correspond to late side effects of treatment, in particular cardiac complications, which are known as a likely consequence of irradiations given close to the heart [36–38].
Our research aimed to estimate the long-term associations of prognostic factors and treatment for breast cancer using flexible excess hazard-based regression models for patients diagnosed in Geneva between 1995 and 2002. Our study highlights the challenges of interpreting these associations in observational data and as well as the need for high quality and detailed clinical information at a population level so that these associations can be examined in detail. With such data, causal inference methods could be applied to be able to describe an effect rather than an association. However, applying causal inference methods requires further methodological work and the development of specialist software for the use of causal inference in the context of excess hazard modelling.
This collaborative work was conducted with the support of the Geneva Cancer Registry and the Cancer Survival Group at London School of Hygiene and Tropical Medicine. We should like to thank Elisabetta Rapiti, Isabelle Neyroud, Massimo Usel, Christine Bouchardy and all collaborators at the Geneva Cancer Registry as well as Michel Coleman and all collaborators at the London School of Hygiene and Tropical Medicine for their support and encouragement.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. In compliance with data protection regulations, data are stored at the Geneva Cancer Registry, Geneva, Switzerland.
RS planned and carried out the literature review and prepared the draft of the paper. Data analysis was planned and implemented by RS and AB. AB, BR and LW provided input and feedback on the content data analysis and on the paper drafts. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The Geneva Cancer Registry has a general authorization provided by “la commission d’expert du secret professionnel en matière de recherché médicale” (https://www.unige.ch/medecine/rgt/index.php/download_file/view/5/197/) to collect nominative data and analyze the anonymized data. Informed consent to participate from patients is therefore unnecessary. Since the study did not require additional clinical information, approval of the Ethics Committee was not required.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gondos A, Bray F, Hakulinen T, Brenner H, EUNICE survival working group. Trends in cancer survival in 11 European populations from 1990 to 2009: a model-based analysis. Ann Oncol Off J Eur Soc Med Oncol. 2009;20(3):564–73.View ArticleGoogle Scholar
- Gondos A, et al. Recent trends in cancer survival across Europe between 2000 and 2004: a model-based period analysis from 12 cancer registries. Eur J Cancer Oxf Engl 1990. 2008;44(10):1463–75.Google Scholar
- Sant M, Francisci S, Capocaccia R, Verdecchia A, Allemani C, Berrino F. Time trends of breast cancer survival in Europe in relation to incidence and mortality. Int J Cancer. 2006;119(10):2417–22.View ArticleGoogle Scholar
- Woods LM, Rachet B, Lambert PC, Coleman MP. ‘Cure’ from breast cancer among two populations of women followed for 23 years after diagnosis. Ann Oncol. Aug. 2009;20(8):1331–6.View ArticleGoogle Scholar
- Brenner H, Hakulinen T. Are patients diagnosed with breast cancer before age 50 years ever cured? J Clin Oncol Off J Am Soc Clin Oncol. 2004;22(3):432–8.View ArticleGoogle Scholar
- Woods LM, Morris M, Rachet B. No ‘cure’ within 12 years of diagnosis among breast cancer patients who are diagnosed via mammographic screening: women diagnosed in the west midlands region of England 1989–2011. Ann Oncol. 2016;27(11):2025–31.View ArticleGoogle Scholar
- Bossard N, et al. Survival of cancer patients in France: a population-based study from the Association of the French Cancer Registries (FRANCIM). Eur. J. Cancer Oxf. Engl. 1990. 2007;43(1):149–60.Google Scholar
- Crocetti E, De Lisi V, Gafà L, Sechi O, Mangone L. Net survival: comparison between relative and cause-specific survival estimates. Epidemiol Prev. 2001;25(3 Suppl):32–6.PubMedGoogle Scholar
- Uhry Z, Bossard N, Remontet L, Iwaz J, Roche L. New insights into survival trend analyses in cancer population-based studies: the SUDCAN methodology. Eur J Cancer Prev. 2017;26:S9–15.Google Scholar
- Pohar Perme M, Estève J, Rachet B. Analysing population-based cancer survival - settling the controversies. BMC Cancer. 2016;16(1):933.View ArticleGoogle Scholar
- Allemani C, et al. Global surveillance of cancer survival 1995-2009: analysis of individual data for 25,676,887 patients from 279 population-based registries in 67 countries (CONCORD-2). Lancet Lond Engl. 2015;385(9972):977–1010.View ArticleGoogle Scholar
- Soerjomataram I, Louwman MWJ, Ribot JG, Roukema JA, Coebergh JWW. An overview of prognostic factors for long-term survivors of breast cancer. Breast Cancer Res Treat. 2008;107(3):309–30.View ArticleGoogle Scholar
- Jatoi I, et al. Time-varying effects of breast Cancer adjuvant systemic therapy. J Natl Cancer Inst. 2015;108(1):djv304.View ArticleGoogle Scholar
- Jatoi I, Anderson WF, Jeong J-H, Redmond CK. Breast cancer adjuvant therapy: time to consider its time-dependent effects. J Clin Oncol Off J Am Soc Clin Oncol. 2011;29(17):2301–4.View ArticleGoogle Scholar
- Bodai BI, Tuso P. Breast cancer survivorship: a comprehensive review of long-term medical issues and lifestyle recommendations. Perm J. 2015;19(2):48–79.Google Scholar
- Cluze C, et al. Analysis of the effect of age on the prognosis of breast cancer. Breast Cancer Res Treat. 2009;117(1):121–9.View ArticleGoogle Scholar
- Andersen PK, Vaeth M. Simple parametric and nonparametric models for excess and relative mortality. Biometrics. 1989;45(2):523–35.View ArticleGoogle Scholar
- Estève J, Benhamou E, Croasdale M, Raymond L. Relative survival and the estimation of net survival: elements for further discussion. Stat Med. 1990;9(5):529–38.View ArticleGoogle Scholar
- Charvat H, et al. A multilevel excess hazard model to estimate net survival on hierarchical data allowing for non-linear and non-proportional effects of covariates. Stat Med. 2016;35(18):3066–84.View ArticleGoogle Scholar
- Remontet L, Bossard N, Belot A, Estève J, French network of cancer registries FRANCIM. An overall strategy based on regression models to estimate relative survival and model the effects of prognostic factors in cancer survival studies. Stat Med. 2007;26(10):2214–28.View ArticleGoogle Scholar
- R Development Core Team, R. A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2009.Google Scholar
- Wynant W, Abrahamowicz M. Impact of the model-building strategy on inference about nonlinear and time-dependent covariate effects in survival analysis. Stat Med. 2014;33(19):3318–37.View ArticleGoogle Scholar
- Royston P, Sauerbrei W. Multivariable model-building: a pragmatic approach to regression analysis based on fractional polynomials for modelling continuous variables. Chichester, England. Hoboken, NJ: John Wiley; 2008.View ArticleGoogle Scholar
- Marshall A, Altman DG, Royston P, Holder RL. Comparison of techniques for handling missing covariate data within prognostic modelling studies: a simulation study. BMC Med Res Methodol. 2010;10:7.View ArticleGoogle Scholar
- Concato J, Peduzzi P, Holford TR, Feinstein AR. Importance of events per independent variable in proportional hazards analysis. I. Background, goals, and general strategy. J Clin Epidemiol. 1995;48(12):1495–501.View ArticleGoogle Scholar
- Peduzzi P, Concato J, Feinstein AR, Holford TR. Importance of events per independent variable in proportional hazards regression analysis. II. Accuracy and precision of regression estimates. J Clin Epidemiol. 1995;48(12):1503–10.View ArticleGoogle Scholar
- Hernán MA, Robins JM. Using big data to emulate a target trial when a randomized trial is not available: table 1. Am J Epidemiol. 2016;183(8):758–64.View ArticleGoogle Scholar
- Splawa-Neyman J, Dabrowska DM, Speed TP. On the application of probability theory to agricultural experiments. Essay on principles. Section 9. Stat Sci. 1990;5(4):465–72.View ArticleGoogle Scholar
- Wright S. Correlation and causation; 1921.Google Scholar
- Mao H, Li L, Yang W, Shen Y. On the propensity score weighting analysis with survival outcome: Estimands, estimation, and inference. Stat Med. 2018;37(26):3745–63.Google Scholar
- Austin PC. The use of propensity score methods with survival or time-to-event outcomes: reporting measures of effect similar to those used in randomized experiments: propensity scores and survival analysis. Stat Med. 2014;33(7):1242–58.View ArticleGoogle Scholar
- Keil AP, Edwards JK, Richardson DB, Naimi AI, Cole SR. The parametric g-formula for time-to-event data: intuition and a worked example. Epidemiology. 2014;25(6):889–97.View ArticleGoogle Scholar
- Newsome SJ, Keogh RH, Daniel RM. Estimating long-term treatment effects in observational data: a comparison of the performance of different methods under real-world uncertainty. Stat Med. 2018;37(15):2367–90.View ArticleGoogle Scholar
- Vansteelandt S, Joffe M. Structural nested models and G-estimation: the partially realized promise. Stat Sci. 2014;29(4):707–31.View ArticleGoogle Scholar
- van der Laan MJ, Rose S. Targeted learning: causal inference for observational and experimental data. New York: Springer; 2011.View ArticleGoogle Scholar
- Cheng Y-J, Nie X-Y, Ji C-C, Lin X-X, Liu L-J, Chen X-M, et al. Long-Term Cardiovascular Risk After Radiotherapy in Women With Breast Cancer. J Am Heart Assoc. 2017;6(5).Google Scholar
- Boerman LM, et al. Long-term outcome of cardiac function in a population-based cohort of breast cancer survivors: A cross-sectional study. Eur J Cancer Oxf Engl 1990. 2017;81:56–65.Google Scholar
- Rygiel K. Cardiotoxic effects of radiotherapy and strategies to reduce them in patients with breast cancer: an overview. J Cancer Res Ther. Jun. 2017;13(2):186–92.View ArticleGoogle Scholar