This article has Open Peer Review reports available.
Different survival analysis methods for measuring long-term outcomes of Indigenous and non-Indigenous Australian cancer patients in the presence and absence of competing risks
© The Author(s). 2017
Received: 12 February 2016
Accepted: 9 December 2016
Published: 17 January 2017
Net survival is the most common measure of cancer prognosis and has been used to study differentials in cancer survival between ethnic or racial population subgroups. However, net survival ignores competing risks of deaths and so provides incomplete prognostic information for cancer patients, and when comparing survival between populations with different all-cause mortality. Another prognosis measure, “crude probability of death”, which takes competing risk of death into account, overcomes this limitation. Similar to net survival, it can be calculated using either life tables (using Cronin-Feuer method) or cause of death data (using Fine-Gray method). The aim of this study is two-fold: (1) to compare the multivariable results produced by different survival analysis methods; and (2) to compare the Cronin-Feuer with the Fine-Gray methods, in estimating the cancer and non-cancer death probability of both Indigenous and non-Indigenous cancer patients and the Indigenous cancer disparities.
Cancer survival was investigated for 9,595 people (18.5% Indigenous) diagnosed with cancer in the Northern Territory of Australia between 1991 and 2009. The Cox proportional hazard model along with Poisson and Fine-Gray regression were used in the multivariable analysis. The crude probabilities of cancer and non-cancer methods were estimated in two ways: first, using cause of death data with the Fine-Gray method, and second, using life tables with the Cronin-Feuer method.
Multivariable regression using the relative survival, cause-specific survival, and competing risk analysis produced similar results. In the presence of competing risks, the Cronin-Feuer method produced similar results to Fine-Gray in the estimation of cancer death probability (higher Indigenous cancer death probabilities for all cancers) and non-cancer death probabilities (higher Indigenous non-cancer death probabilities for all cancers except lung cancer and head and neck cancers). Cronin-Feuer estimated much lower non-cancer death probabilities than Fine-Gray for non-Indigenous patients with head and neck cancers and lung cancers (both smoking-related cancers).
Despite the limitations of the Cronin-Feuer method, it is a reasonable alternative to the Fine-Gray method for assessing the Indigenous survival differential in the presence of competing risks when valid and reliable subgroup-specific life tables are available and cause of death data are unavailable or unreliable.
In Australia, Indigenous cancer patients have lower net survival than non-Indigenous patients , but there is limited information about whether Indigenous cancer patients also have higher probability of non-cancer death, given the Indigenous population’s higher prevalence of chronic diseases and much higher all-cause mortality rates .
Net survival is the measure most commonly used to compare cancer prognosis in different populations (e.g., between nations, regions, or racial/ethnic groups) [1, 3–9] because net survival removes the effect of non-cancer deaths (which may vary between populations because of different all-cause mortality rates ) and measures the hypothetical scenario in which patients are only able to die of their cancer. However, this is a disadvantage of net survival for patients and clinicians who, when considering the treatment options and weighing up the benefits, drawbacks, and toxicities of cancer therapy, want to know about actual prognosis: what is the chance of dying (from cancer or any cause) compared with the chance of surviving [5, 9, 10]? This is particularly important for Indigenous Australians, and other populations who also have higher probability of death from “competing” (i.e., non-cancer) causes [3, 5, 10].
For this purpose, “crude probability of death” is a more appropriate measure of prognosis than net survival because it takes competing risk into account and enables us to identify whether survival disparities are caused by cancer deaths, non-cancer deaths, or both. Crude probability of cancer death (hereafter referred to as “cancer death probability”) is the probability of death from cancer in the presence of other causes; while crude probability of non-cancer death (hereafter referred to as “non-cancer death probability”) is the probability of death from other causes in the presence of cancer.
A summary of different cancer prognosis statistics and survival analysis methods used in this study
Using cause of death data
Using life tables
Net survival (excludes competing risks)
Divide observed survival by expected survival
Multivariable analysis method
Cox proportional hazard model
Crude probability (includes competing risks)
Crude probability of death
Crude probability of death using life table
Multivariable analysis method
Life table-based methods are commonly used for comparing cancer survival between population groups because of deficiencies in cause of death data; such deficiencies may be greater for disadvantaged subgroups. However, life tables are often not available for such subgroups because of uncertainty about population estimates and absent or incomplete identification of subgroup members in death registrations. Australia has high-quality cause of death data but incomplete life tables for the Indigenous population, except for the Northern Territory (NT) where Indigenous Australians are very high prevalence (30% of the NT population compared to 2–4% in other states) and consequently are accurately identified in deaths data and have reliable population estimates . Indigenous status is also recorded with a high degree of accuracy in the NT Cancer Register (NTCR) . The availability of high-quality data for the Indigenous population and Indigenous cancer patients puts the NT in the best position to compare the implications of using the Cronin-Feuer and Fine-Gray methods to assess the Indigenous disparities in cancer and non-cancer deaths.
The aim of this study is two-fold: (1) to compare the multivariable results (i.e., Indigenous differentials) produced by relative survival, cause-specific survival, and competing risk (Fine-Gray) survival analysis methods; and (2) to compare the Cronin-Feuer with the Fine-Gray survival analysis methods, in estimating the cancer and non-cancer death probability of both Indigenous and non-Indigenous cancer patients and the Indigenous cancer disparities in the presence of competing causes of death.
Cancer registrations data for all NT residents diagnosed with cancer between 1 January 1991 and 31 December 2009 was obtained from the NTCR for the following data items: sex; date of birth; Indigenous status; remoteness of residence (urban or remote); date of diagnosis; cancer site, coded according to the International Classification of Diseases 10 (ICD-10); date of death; and underlying cause of death. Vital status was verified by matching the NTCR dataset to the National Death Index for deaths occurring up to 31 December 2011 (at the time of analysis, cause of death data were only available up until 31 December 2011). For the survival analysis of site-specific cancers, female breast (C50), colorectal (C18-20), lung (C33-34), and head and neck (C1-C14 & C30-32) cancers were chosen because there were sufficient numbers of both Indigenous and non-Indigenous cases for the analysis, and they are Australia’s designated national “priority” cancers  (except head and neck cancers). Head and neck cancer was chosen as it is a major cause of morbidity and mortality associated with smoking.
The censoring date for the survival analysis was 31 December 2011. The survival time for cancer patients who died on the day of diagnosis was counted as half a day.
Relative survival was calculated using the “strs” procedures  of Stata, in which the expected survival was estimated using the Ederer II method . The NT non-Indigenous population has very high migration to and from other parts of Australia and similar health status and mortality to Australians generally, so life tables for the relevant years (1991–2009) for the total Australian population were used for non-Indigenous cancer patients. Life tables for the NT Indigenous population (for which deaths and population data are more reliable than for Indigenous people in other states and territories ) were used for the Indigenous cancer patients. Cancer and non-cancer death probabilities were calculated in two ways: the Fine-Gray model used cause of death information from the cancer cohort, and was derived using the “stcrreg” command in Stata, while the Cronin-Feuer method used life table data  derived using the “strs” command in Stata .
For multivariable analysis, the Cox proportional hazard regression, Poisson regression (generalized linear model with Poisson error structure), and Fine-Gray regression were used to investigate the Indigenous survival disparity; multivariable analysis was not available for the Cronin-Feuer method. The same independent variables were included in all final models: Indigenous status (Indigenous compared to non-Indigenous); age at diagnosis (per year); gender (female compared with male); and cancer site. Since the effect of age at diagnosis was found to be different for Indigenous compared with non-Indigenous patients, an interaction term for Indigenous status by age at diagnosis (base age 55 years) was added to the models. For the Cox proportional hazard regression, scaled Schoenfeld residuals were used to check if the proportional hazards assumptions of each variable were satisfied. As the proportional hazards assumption was not met for Indigenous status, a step function of Indigenous status with follow-up time (as annual intervals) was also included in the model.
In the multivariable analysis for time trends, as all subjects had at least 2 years of potential follow-up, the follow-up was limited to the first 2 years after diagnosis so that the shorter follow-up time for subjects diagnosed late in the study period did not bias time trends. Regression models included the same terms as above, plus year of diagnosis.
All analyses were conducted using Stata/SE 13 (StataCorp, College Station, TX). Ethics approval was obtained from the Human Research Ethics Committee of the Menzies School of Health Research and NT Department of Health. Approval to use the cancer registrations data was obtained from the NT Cancer Registry.
Demographic characteristics of people diagnosed with cancer, Australia NT, 1991–2009
Indigenous (n = 1789)
Non-Indigenous (n = 7806)
Age at diagnosis (%)
0 to 49 years
50 to 59 years
60 to 69 years
70 years and over
Median age (years)
Vital Status at 31 Dec 2011 (%)
One-year and 5-year cumulative relative survival (%) and 95% confidence interval by Indigenous status and cancer site (age and sex adjusted), Australia NT, 2001–2009
Head and necka
Regression analysis of cause-specific mortality (Cox proportional hazard regression), relative survival (Poisson regression), and competing risk (Fine-Gray regression), all cancers combineda, Australia NT, 1991–2009
Competing (due to cancers)
Competing (other death)
HR (95% CI)b
HR (95% CI)
SHR (95% CI)
SHR (95% CI)
1st year after diagnosis
2nd year after diagnosis
3rd year after diagnosis
4th year after diagnosis
5th year after diagnosis
Female vs male
Age at diagnosisd
Regression analysis of time trend after diagnosis using cause-specific mortality (Cox proportional hazard regression), relative survival (Poisson regression), and competing risk analysis (Fine-Gray regression), all cancers combineda, Australia NT, 1991–2009
Competing (due to cancers)
Competing (other death)
HR (95% CI)b
HR (95% CI)
SHR (95% CI)
SHR (95% CI)
Female vs male
Age at diagnosisd
Year of diagnosis
Five-year cumulative probabilities of cancer and non-cancer death (%), age and sex adjusted
Cancer death probabilities
Head & neckc
Head & neckc
This study is the first to examine the implications of various survival methods to assess the survival inequalities faced by Indigenous cancer patients in the NT, in particular the impact of competing causes of death. Multivariable regression using relative survival (Poisson model), cause-specific survival (Cox regression) and competing risk analysis (Fine-Gray model) produced similar results. For all cancers combined, death rates were higher for Indigenous, male and older cancer patients for both cancer and non-cancer death. In the absence of competing risk, Indigenous people have lower 5-year net survival than non-Indigenous people for all three cancers (breast, colorectal, and head and neck cancers).
Crude probability of death potentially offers a more understandable way than net survival to communicate cancer prognosis to patients and clinicians, showing the different mortality risk profiles (cancer and non-cancer death probabilities) among the Indigenous and non-Indigenous populations, while acknowledging that cancer patients can also die of non-cancer causes . It enables us to identify whether survival disparities are caused by cancer, non-cancer, or both cancer and non-cancer deaths, as well as identify the disparities in both cancer and non-cancer deaths across time. The multivariable competing risk analysis, Fine-Gray model, found that Indigenous patients were more likely than non-Indigenous to die of non-cancer causes and that, in contrast to cancer deaths, this differential increased rather than decreased with time after diagnosis. Five years after diagnosis, Indigenous cancer patients had higher rates of non-cancer death than non-Indigenous patients and this differential was increasing with time since diagnosis, while there was no longer a differential in cancer deaths. This finding suggests the need for integrated chronic disease management strategies for cancer patients, particularly for the Indigenous population with higher rate of chronic diseases comorbidities. The finding also demonstrates an important advantage of the competing risk approach over conventional relative survival in analyzing long-term outcomes for Indigenous cancer patients.
In the presence of competing risks, the Cronin-Feuer method produced similar results to the Fine-Gray method in the estimation of cancer death probability (higher Indigenous cancer death probabilities for all cancers) and non-cancer death probabilities (higher Indigenous non-cancer death probabilities for all cancers except lung cancer and head and neck cancers) 5 years after cancer diagnosis. The Cronin-Feuer method estimated much lower non-cancer death probabilities than the Fine-Gray method for non-Indigenous patients with head and neck cancers and lung cancer, which are both smoking-related cancers.
Comparison of death probabilities estimated using the Cronin-Feuer method and Fine-Gray models provides insights into “external factors” affecting measures of survival for various cancer types such as smoking-related cancers. The death probabilities estimated using the Cronin-Feuer method can be seen as the death probabilities that the cancer patients were “expected to have”, in the hypothetical scenario where cancer patients have the same death risks as the general population, unaffected by “external factors” that cancer patients experienced such as unhealthy behaviors, access to health care, or screening effects. The death probabilities estimated using cause of death data and the Fine-Gray method can be seen as the “observed death probabilities” in which cancer patients might have different risks of death to the general population due to various “external factors”. Therefore, higher (or lower) non-cancer death probabilities estimated by the Fine-Gray model imply higher (or lower) non-cancer death risks in cancer patients than the general population.
When analyzing smoking-related cancers, caution is required for prognosis measures calculated using general population life tables such as relative survival and the Cronin-Feuer method [18–20], which assume that cancer patients have the same death risks as the general population. This assumption is probably violated for smoking-related cancers, because smoking increases mortality from other causes of death such as respiratory and cardiovascular diseases, resulting in higher non-cancer death risks in cancer patients than the general population. Therefore, the Cronin-Feuer method will underestimate the non-cancer death probabilities of the patients with smoking-related cancers if it uses the general population life tables for the cancer patients with higher non-cancer death risks. The results from our study support the hypothesis that methods using life tables (relative survival and Cronin-Feuer) underestimate non-cancer deaths and overestimate cancer deaths for smoking-related cancers for non-Indigenous patients, suggesting that Cronin-Feuer method should be used with caution in estimating cancer and non-cancer death probability for non-Indigenous patients with smoking related cancers and quantifying the Indigenous disparity for smoking-related cancers. Our study has shown that the Indigenous disparity in smoking-related cancer death probabilities estimated by Cronin-Feuer method is lower than that estimated by the Fine-Gray method, which suggests that survival analysis methods that use the life table such as relative survival might underestimate the Indigenous disparity in cancer survival for smoking-related cancers.
Previous studies suggested that relative survival tends to underestimate social inequalities in cancer survival [21–23]. Our study suggests that the Cronin-Feuer method does underestimate Indigenous to non-Indigenous disparities in cancer death probabilities (except for breast cancer) while overestimating disparities in non-cancer death probabilities (except for lung, head and neck cancers, which are both smoking-related cancers) after considering competing risks. Underestimation of cancer death disparity by the Cronin-Feuer method was largest for colorectal cancer and head and neck cancers. For colorectal cancer, this underestimation was mainly due to underestimation of the probability of cancer death for Indigenous patients. For head and neck cancers, underestimation of disparity is due to a combination of underestimation of probability of cancer death for Indigenous patients and overestimation of cancer death probabilities for non-Indigenous patients. However, for breast cancer patients the disparity in cancer deaths was similar for the two methods; the Fine-Gray method produced lower estimates of both cancer and non-cancer death probability than the Cronin-Feuer method for both Indigenous and non-Indigenous patients.
The lower estimates of cancer death probabilities of the Fine-Gray method might indicate that cancer survival for breast cancer patients is better than expected (calculated using life table). The lower estimates of non-cancer death probabilities of the Fine-Gray method might indicate lower non-cancer 5-year death risk in breast cancer patients than the background population, supported by the findings of another study  (not specifically for Indigenous cancer patients) that breast cancer patients have reduced risks for death (compared to the general population) from cardiovascular disease (SMR = 83.9), which accounts for the majority (55%) of non-cancer deaths among cancer patients generally in Australia; higher risks for non-cancer deaths were observed for all cancer types except breast cancer and melanoma. The lower than expected cancer and non-cancer death probability for breast cancer patients may be due to these patients having higher socioeconomic status, which might be a protective factor of cancer and non-cancer deaths, as indicated by the finding that breast cancer is more common in areas of higher socioeconomic status in Australia  (breast cancer incidence rate was 122 per 100,000 in the highest socioeconomic status group, compared to 103 per 100,000 in the lowest socioeconomic status group).
While net survival is useful for reporting trends in cancer survival, comparing different groups of cancer patients and investigating the impact of various factors on cancer treatment, crude probability of death is important when communicating risks to patients  during clinical decision making [9, 26]. Previous studies have suggested that relative survival tends to underestimate social inequalities in cancer survival; our study suggests that the Cronin-Feuer method does underestimate Indigenous to non-Indigenous disparities in cancer death probabilities (except for breast cancer) in the presence of competing risks. Our study also suggests that when analyzing smoking-related cancers, caution is required when measuring cancer survival disparity using population life tables. Despite the limitations of the Cronin-Feuer method , it is a reasonable alternative to the Fine-Gray method for assessing the Indigenous survival differential in the presence of competing risks when valid and reliable subgroup-specific life tables are available and cause of death data are unavailable or unreliable.
We thank Karen Dempsey (the former NT Cancer Register data-custodian at NT DoH Health Gains Planning Branch) for extracting the cancer register data. Part of this manuscript was written during a writing retreat funded by the Centre of Research Excellence in Discovering Indigenous Strategies to improve Cancer Outcomes via Engagement, Research Translation and Training (funded by the National Health and Medical Research Council #1041111). VH was supported by a University Postgraduate Research Scholarship from Charles Darwin University.
Availability of data and materials
Data are not available to the public. As the NT Cancer register data are government health administrative data, the ethics protocol requires (1) the data to be stored on the secure information system at the NT Department of Health (2) access to the data to be restricted (by password protection) to the study investigators and members of the research team and (3) no data on individual patients will be published and no further copies will be made or distributed to any person or organisation.
VH conceived and designed the experiment, applied for and obtained the ethics and data-custodian approval (cancer register data), performed the data management and statistical analysis, conducted the literature review, and wrote the manuscript draft. JC, YZ, PB, and XZ contributed to the review and revision of the related journal paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Ethics approval was obtained from the Human Research Ethics Committee of the Menzies School of Health Research and NT Department of Health (HREC 2013–2090). Approval to use the cancer registrations data was obtained from the NT Cancer Registry.
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.
- Condon J, Zhang X, Baade P, Griffiths K, Cunningham J, Roder D, Coory M, Jelfs P, Threlfall T. Cancer survival for Aboriginal and Torres Strait Islander Australians: a national study of survival rates and excess mortality. Popul Health Metrics. 2014;12:1.View ArticleGoogle Scholar
- Australian Institute of Health and Welfare. The health and welfare of Australia’s Aboriginal and Torres Strait Islander peoples 2015. Cat. no. IHW 147. Canberra: AIHW; 2015.Google Scholar
- Mariotto AB, Noone A-M, Howlader N, Cho H, Keel GE, Garshell J, Woloshin S, Schwartz LM. Cancer survival: an overview of measures, uses, and interpretation. JNCI Monogr. 2014;2014:145–86.View ArticleGoogle Scholar
- Perme MP, Stare J, Estève J. On estimation in relative survival. Biometrics. 2012;68:113–20.View ArticlePubMedGoogle Scholar
- Howlader N, Mariotto AB, Woloshin S, Schwartz LM. Providing clinicians and patients with actual prognosis: cancer in the context of competing causes of death. JNCI Monogr. 2014;2014:255–64.View ArticleGoogle Scholar
- Withrow DR, Racey CS, Jamal S. A critical review of methods for assessing cancer survival disparities in indigenous population. Ann Epidemiol. 2016;26:579–91.View ArticlePubMedGoogle Scholar
- Howlader N, Ries LA, Mariotto AB, Reichman ME, Ruhl J, Cronin KA. Improved estimates of cancer-specific survival rates from population-based data. J Natl Cancer Inst. 2010;102:1584–98.View ArticlePubMedPubMed CentralGoogle Scholar
- Baade PD, Dasgupta P, Dickman PW, Cramb S, Williamson JD, Condon JR, Garvey G. Quantifying the changes in survival inequality for Indigenous people diagnosed with cancer in Queensland, Australia. Cancer Epidemiol. 2016;43:1–8.View ArticlePubMedGoogle Scholar
- Zakeri K, MacEwan I, Vazirnia A, Cohen EEW, Spiotto MT, Haraf DJ, Vokes EE, Weichselbaum RR, Mell LK. Race and competing mortality in advanced head and neck cancer. Oral Oncol. 2014;50:40–4.View ArticlePubMedGoogle Scholar
- Eloranta S, Adolfsson J, Lambert PC, Stattin P, Akre O, Andersson TM, Dickman PW. How can we make cancer survival statistics more useful for patients and clinicians: An illustration using localized prostate cancer in Sweden. Cancer Causes Control. 2013;24:505–15.View ArticlePubMedGoogle Scholar
- Condon JR, Barnes T, Cunningham J, Smith L. Demographic characteristics and trends of the Northern Territory Indigenous population, 1966 to 2001. Darwin: Cooperative Research Centre for Aboriginal Health; 2004.Google Scholar
- Condon J, Zhao Y, Armstrong B, Barnes T. Northern Territory Cancer Register data quality, 1981–2001. Darwin: NT Cancer Registry, the Cooperative Research Centre for Aboriginal Health, Charles Darwin University and the Menzies School of Health Research; 2004.Google Scholar
- AIHW, Australasian Association of Cancer Registries. Cancer survival in Australia 1992–1997: geographic categories and socioeconomic status. Cancer Series no. 22. Cat. no. CAN 17. Canberra: AIHW; 2003. Viewed 25 November 2015. http://www.aihw.gov.au/publication-detail/?id=6442467451.
- Dickman PW, Coviello E, Hills M. Estimating and modelling relative survival. Stata J. 2009;15:186–215.Google Scholar
- Dickman PW, Sloggett A, Hills M, Hakulinen T. Regression models for relative survival. Stat Med. 2004;23:51–64.View ArticlePubMedGoogle Scholar
- Cronin KA, Feuer EJ. Cumulative cause-specific mortality for cancer patients in the presence of other causes: a crude analogue of relative survival. Stat Med. 2000;19:1729–40.View ArticlePubMedGoogle Scholar
- Eloranta S, Lambert PC, Andersson TM, Czene K, Hall P, Björkholm M, Dickman PW. Partitioning of excess mortality in population-based cancer patient survival studies using flexible parametric survival models. BMC Med Res Methodol. 2012;12:1.View ArticleGoogle Scholar
- Hinchliffe SR, Rutherford MJ, Crowther MJ, Nelson CP, Lambert PC. Should relative survival be used with lung cancer data? Br J Cancer. 2012;106:1854–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Blakely T, Soeberg M, Carter K, Costilla R, Atkinson J, Sarfati D. Bias in relative survival methods when using incorrect life‐tables: Lung and bladder cancer by smoking status and ethnicity in New Zealand. Int J Cancer. 2012;131:E974–82.View ArticlePubMedGoogle Scholar
- Ellis L, Coleman M, Rachet B. The impact of life tables adjusted for smoking on the socio-economic difference in net survival for laryngeal and lung cancer. Br J Cancer. 2014;111:195–202.View ArticlePubMedPubMed CentralGoogle Scholar
- Baade PD, Turrell G, Aitken JF. A multilevel study of the determinants of area-level inequalities in colorectal cancer survival. BMC Cancer. 2010;10:24.View ArticlePubMedPubMed CentralGoogle Scholar
- Dejardin O, Remontet L, Bouvier A, Danzon A, Tretarre B, Delafosse P, Molinié F, Maarouf N, Velten M, Sauleau E. Socioeconomic and geographic determinants of survival of patients with digestive cancer in France. Br J Cancer. 2006;95:944–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Auvinen A, Karjalainen S. Possible explanations for social class differences in cancer patient survival. In: Kogenivas MPN, Susser M, Boffetta P, editors. Social inequalities and cancer. Lyon: IARC Scientific Publications; 1997. p. 377–97.Google Scholar
- Baade PD, Fritschi L, Eakin EG. Non-cancer mortality among people diagnosed with cancer (Australia). Cancer Causes Control. 2006;17:287–97.View ArticlePubMedGoogle Scholar
- Australian Institute of Health and Welfare & Cancer Australia. Breast cancer in Australia: an overview. Cancer series no. 71. Cat. no. CAN 67. Canberra: AIHW; 2012.Google Scholar
- Charvat H, Bossard N, Daubisse L, Binder F, Belot A, Remontet L. Probabilities of dying from cancer and other causes in French cancer patients based on an unbiased estimator of net survival: a study of five common cancers. Cancer Epidemiol. 2013;37:857–63.View ArticlePubMedGoogle Scholar