Next Article in Journal
HPV Type Distribution in Benign, High-Grade Squamous Intraepithelial Lesions and Squamous Cell Cancers of the Anus by HIV Status
Next Article in Special Issue
Pancreatic Cancer Incidence Trends by Race, Ethnicity, Age and Sex in the United States: A Population-Based Study, 2000–2018
Previous Article in Journal
Novel Computed-Tomography-Based Transformer Models for the Noninvasive Prediction of PD-1 in Pre-Operative Settings
Previous Article in Special Issue
Health Burdens and SES in Alabama: Using Geographic Information System to Examine Prostate Cancer Health Disparity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Is the Social Gradient in Net Survival Observed in France the Result of Inequalities in Cancer-Specific Mortality or Inequalities in General Mortality?

1
ANTICIPE U1086 INSERM-UCN, Equipe Labellisée Ligue Contre le Cancer, Centre François Baclesse, Normandie Université UNICAEN, 14000 Caen, France
2
Service de Biostatistique—Bioinformatique, Pôle Santé Publique, Hospices Civils de Lyon, 69000 Lyon, France
3
University of Lyon, 69000 Lyon, France
4
University of Lyon 1, 69100 Villeurbanne, France
5
Équipe Biostatistique-Santé, Laboratoire de Biométrie et Biologie Évolutive, CNRS, UMR 5558, 69100 Villeurbanne, France
6
Inequalities in Cancer Outcomes Network, Department of Non-Communicable Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK
7
French Network of Cancer Registries (FRANCIM), 31000 Toulouse, France
8
Loire-Atlantique-Vendée Cancer Registry, 44000 Nantes, France
9
Centre d’Epidémiologie et de Recherche en santé des POPulations (CERPOP) UMR1295, Université de Toulouse Paul Sabatier, Inserm, 31000 Toulouse, France
10
Research Department, Caen University Hospital Centre, 14000 Caen, France
*
Author to whom correspondence should be addressed.
Collaborators of the Francim Group is provided in the Acknowledgments.
Cancers 2023, 15(3), 659; https://doi.org/10.3390/cancers15030659
Submission received: 24 October 2022 / Revised: 11 January 2023 / Accepted: 17 January 2023 / Published: 20 January 2023
(This article belongs to the Special Issue Advancing Health Equity to Reduce Cancer Health Disparities)

Abstract

:

Simple Summary

That the original French life tables are not stratified in terms of deprivation whilst the background mortality in the general population differs according to socio-economic position, social gradient in the net survival of patients with cancer, as was found in a previous study, could be due, at least partly, to socially-determined co-morbidities. We found that the social gradient in cancer net survival was reduced using simulated deprivation-specific life tables. This study alerts us to the fact of this overestimation in the social gradient in cancer net survival using the original life tables, which, in a few cases, can be so important that conclusions might be wrong (e.g., prostate cancer). As this work relies upon simulated rather than real data, we were not able to precisely quantify the potential bias resulting from the lack of deprivation-specific life tables. This present study points to how important it is to create proper deprivation-specific life tables in order to accurately investigate social inequalities in cancer net survival analyses.

Abstract

Background: In cancer net survival analyses, if life tables (LT) are not stratified based on socio-demographic characteristics, then the social gradient in mortality in the general population is ignored. Consequently, the social gradient estimated on cancer-related excess mortality might be inaccurate. We aimed to evaluate whether the social gradient in cancer net survival observed in France could be attributable to inaccurate LT. Methods: Deprivation-specific LT were simulated, applying the social gradient in the background mortality due to external sources to the original French LT. Cancer registries’ data from a previous French study were re-analyzed using the simulated LT. Deprivation was assessed according to the European Deprivation Index (EDI). Net survival was estimated by the Pohar–Perme method and flexible excess mortality hazard models by using multidimensional penalized splines. Results: A reduction in net survival among patients living in the most-deprived areas was attenuated with simulated LT, but trends in the social gradient remained, except for prostate cancer, for which the social gradient reversed. Flexible modelling additionally showed a loss of effect of EDI upon the excess mortality hazard of esophagus, bladder and kidney cancers in men and bladder cancer in women using simulated LT. Conclusions: For most cancers the results were similar using simulated LT. However, inconsistent results, particularly for prostate cancer, highlight the need for deprivation-specific LT in order to produce accurate results.

1. Introduction

Social inequalities in health are a major concern. Whatever the disease, in order to better understand the underlying mechanisms of the social gradient in survival—and to be able to propose targeted actions to remedy it—it is necessary to know to what extent it can be attributed to general mortality or disease-specific mortality. In the cancer field, it is important to make this distinction in order to know to what extent the observed social gradient in survival is due to differences in management (screening, diagnostic, treatment, and follow-up) associated with the social environment.
Social inequalities in cancer survival have been reported worldwide and represent a major public health burden [1,2,3,4,5,6,7,8,9,10]. As a matter of fact, public health policies, both at the national (French cancer plans and Stratégie Décennale Santé) and the European (Cancer Mission) level are focused on the reduction of social inequalities in cancer screening, incidence, and survival.
In France, a recent study [11] investigated the social gradient in cancer survival for 19 major solid tumors (about 200,000 cancer cases) using data from the French Network of Cancer registries (Francim). The study found lower 5-years survival among patients living in the most-deprived areas, as compared to patients living in the least-deprived ones, for almost all cancer sites. This reduction in survival probability among the most deprived reached 60% for bile ducts in women, and 35% for liver cancer in men. In this study, net survival was estimated using the relative survival setting or excess mortality framework [12]. In this setting, the knowledge of the cause of death of cancer patients is not required, but appropriate life tables (LT) are needed. Indeed, net survival for a given cancer is a conceptual and theoretical measure [13,14] derived from the excess (or cancer-specific) mortality, this latter being the difference between the observed mortality of cancer patients and their expected or background mortality, known from general population LT. However, as in other countries, an important limitation of this previous work was the lack of stratification in terms of deprivation of the French LT produced by the National Institute of Statistics and Economic Studies (INSEE), which prevented accounting for the expected social gradient in general mortality [5,15].
As a matter of fact, mortality in the general population is known to be socially differentiated. Mackenbach and colleagues [16] reported higher rates of mortality among low- versus highly educated people in several European countries. For example, in France, mortality is doubled among the least educated compared with the most educated, and in Lithuania, mortality is up to three times higher among people with low education. The same pattern was recently found in England and Wales [17]. Therefore, in net survival analyses using relative survival setting, if social inequalities in background mortality in the general population are not properly accounted for, they will be included in the social gradient in excess mortality in the cancer registries population: this will mathematically result in an overestimation of the social gradient in excess mortality [18]. In other words, in net survival analyses, without stratification of LT by social situation, if there is a social gradient in the observed mortality, a social gradient in excess mortality (likewise in net survival) will be found, but it might result totally or partly from the unaccounted social gradient in the expected mortality.
As original French LT ignore the social inequalities in mortality in the general population (being stratified only by sex, age, year, and French Département, which is a territorial and administrative division in France) an overestimation of the social gradient in cancer net survival cannot be excluded in this previous French study [11].
Correcting methods for excess mortality models are thus being developed [19,20,21], but to date, they are mainly experimental. One solution is to conduct sensitivity analyses using simulated deprivation-specific LT [5,12,15]. Using this method, Antunes and colleagues [15] have shown that, in Portugal, social inequalities in colorectal cancer net survival were most likely attributable to inequalities in background survival. However, Ito and colleagues [5] still found social inequalities in cancer net survival for several cancers in sensitivity analyses, and showed that the overestimation of the social gradient might be smaller at 1 year than at 5 years of follow-up.
In that context, our aim was to evaluate whether the social deprivation gap in cancer net survival observed in France could be attributable to this lack of appropriate LT. To accomplish this, we simulated deprivation-specific French LT and compared net survival analyses using original LT versus net survival analyses using simulated deprivation-specific LT, thus conducting sensitivity analyses for results from Tron et al. 2019, providing social gradient in survival for almost all cancer sites using French population-based registries data [11].

2. Materials and Methods

2.1. Cancer Registries Data

Data were collected from French cancer registries (Francim) as described previously [11]. Cancer cases diagnosed between 1 January 2006 and 31 December 2009 in patients over 15 years old were included and followed-up until the date of death or 30 June 2013. The quality of the data derived from the registries was regularly assessed by the National Committee of Registries (CER) at the national level and by the International Agency for Research on Cancer (IARC) at the international level. The study was approved by the Consultative Committee for the Processing of Health Research Data (CCTIRS) and the French Data Protection Authority (CNIL, authorization n°913013). Available variables were survival time (i.e., difference between the date of death or last information on vital status and the date of cancer diagnosis), sex, age at diagnosis, and social environment, assessed by the 2007 version of the European Deprivation Index (EDI) [22,23]. The EDI is an aggregate ecological measure of relative poverty (deprivation), available in France at the “IRIS” (Ilots Regroupés pour l’Information Statistique) level and based on census data and information from the European Union Statistics on Income and Living Conditions survey (EU-SILC). IRIS are small (around 2000 inhabitants) geographic areas created by INSEE, for which census data are available. For each cancer case recorded in the registries, the patient’s address was geolocalized using Geographic Information Systems (ArcGIS 10.5, ESRI, Redlands, CA, USA) in order to be allocated to an IRIS and thus the corresponding EDI score. The EDI score is a continuous variable ranging from −17 to 51 (median value: −0.75) at the national level—the higher the score, the greater the deprivation in the IRIS. It has also been categorized according to quintiles of distribution at the national level (Q1 referring to the least- and Q5 to the most-deprived environment). Missing data for EDI accounted for less than 1%; therefore, we performed complete cases analyses.

2.2. Simulation of Deprivation-Specific French LT

Following the same methodology used in previous studies [5,12,15,24,25], two sets of simulated deprivation-specific French LT were created by applying, to the original French LT, the social gradient in mortality observed in two external sources. The first source was the deprivation-specific LT available for the general population of England (Eng LT) [26], which provide the mortality rate ratios from 2006 to 2011, by sex, age, and national quintiles of the income domain score of the Indice of Multiple Deprivation (IMD)—an aggregated ecological index at the contextual level [27]. England deprivation-specific LT were chosen, because they are fully available, and because England is known to have large mortality inequalities, as is France [16]. The second source was a large-scale socio-demographic panel established in France in the general population, the Echantillon Démographique Permanent (EDP, Permanent Demographic Sample) [28] which provides the mortality rate ratios (mortality quotient per 100,000 people) for the period 2012–2016, by sex, age and by educational level, socio-professionnal category or net income per consumption unit of household. We chose to use the net income per consumption unit of household —a single indicator at the individual level for better comparability. In the EDP panel, 20-quantiles of net income per consumption unit of household were available, and we calculated the means of 20-quantile 1 to 4, 5 to 8, 9 to 12, 13 to 16, and 17 to 20 in order to obtain a variable with 5 categories for the analyses.
In both cases (Eng LT and EDP), we extracted the smoothed expected mortality rates from the external source according to five levels of deprivation. Next, we estimated the mortality rate ratios between each of these five categories and the overall general mortality provided by the external source (see Supplementary Material, Figure S1). Finally, for each stratum of sex, age, year, and French Département in the French LT, we applied these rate ratios to the original mortality rate in order to obtain new ones for 5 levels of deprivation. Thus, we applied the social gradient in mortality observed in the general population from external sources (Eng LT and EDP) to the original French LT. This resulted in simulated French LT stratified by sex, age, year, French Département, and deprivation. The detailed calculation process is presented in Supplementary Information.

2.3. Statistical Analyses

In the relative survival setting, net survival can be studied using non-parametric Pohar–Perme method [14] or by the flexible modelling of excess mortality hazard using multidimensional penalized splines [29,30]. These two methods have been used in the previous studies using the original French LT; results from the Pohar–Perme method have been published in Tron et al., 2019 [11], whereas results from flexible modelling have been published in Poiseuil et al., 2022 [31] and Tron et al., 2021 [32].
In the present study (which corresponds to a sensitivity analysis of the previous ones), we then compared these initial results with those obtained using the same exact methodology (Pohar–Perme and modelling) but using the two sets of simulated LT. We recall below some characteristics concerning the methods.
(i)
The 5-year age-standardized net survival probabilities were estimated using the Pohar–Perme method [14]. For each cancer site and sex, we calculated the 5-year deprivation gaps (and their 95% confidence intervals [CI]), which are the difference in 5-year age-standardized net survival probabilities between patients from the least- and most-deprived environments defined by the 1st and 5th national quintile of EDI, respectively (see Tron et al., 2019 [11] method section for further calculation details).
(ii)
In flexible modelling, at given values of time (t), age at diagnosis, and EDI, the observed mortality hazard λ of a given patient is decomposed as follows: λ(t,age,EDI,z) = λE(t,age,EDI) + λP(age + t,year + t,z), where λE is its excess mortality hazard (EMH), which is the mortality directly or indirectly due to cancer, and λP is its expected mortality, i.e., the all-cause mortality hazard of the general French population at age at diagnosis + t, and year of diagnosis + t, given demographic characteristics z of that individual. Here z is composed of the variables sex, year of death, residence French Département, and deprivation in the simulated LT. The EMH was modeled using (multidimensional) penalized splines, which allows for the modelling of flexible baseline hazard, the non-linear and non-proportional (i.e., time-dependent) effects of covariates, as well as interactions [29,30]. More precisely, four models based on penalized splines were adjusted, and the best one was selected according to the corrected Akaike Information Criterion (AIC) [33] indicating the overall effect of EDI on cancer net survival and its form, either 1—no effect; 2—proportional (i.e., not time-dependent) effect; 3—time-dependent effect; or 4—time- and/or age-dependent effect (i.e., interaction EDI*t and/or EDI*age). Then, excess mortality hazard ratios (EHR) by EDI were calculated based on the selected model. See the methods section of Poiseuil et al., 2022 [31] and Tron et al., 2021 [32] for further details about the modelling strategy.
For each cancer and sex, three analyses were compared: main analyses based on the original French LT, sensitivity analyses using English LT, and external source and sensitivity analyses using EDP external source. The following elements were compared between the three analyses:
-
value of the 5-year deprivation gaps (and 95% CI)
-
the selected flexible model (indicating the effect of EDI on survival and its form)
-
the curves of EHR as a function of EDI (using the 10th percentile of EDI as a reference, i.e., EDI score = −3.9)
-
and the curves of the EHR of the 90th percentile (p90 = 4.4) as compared to the 10th percentile (p10 = −3.9) of EDI as a function of time since cancer diagnosis.
All of the described models were fitted using R software (R Core Team, Vienna, Austria, version 3.5.1) with package ‘relsurv’ (2.2.3) and ‘survPen’ (1.0.1) [34].

2.4. Data Availability Statement

The data generated in this study are available upon request from the corresponding author.

3. Results

A total of 185,169 solid tumors (18 different cancer sites, 14 in men and 17 in women) were analyzed as in Tron et al., 2019 [11] (see Table 2 of this previous article for the description of the studied population and cancer sites).

3.1. Comparison of Main and Sensitivity Analyses Based on Non-Parametric Method

Overall, 5-year deprivation gaps were reduced in both sensitivity analyses as compared to main analyses.
In men (Figure 1a), clear 5-year deprivation gaps were observed in all three analyses for ear, neck, throat (ENT) cancers and liver cancer. The results regarding kidney, stomach, pancreas, central nervous system (CNS), and bile ducts cancers as well as melanoma and sarcomas were identical according to the three analyses, with negligible 5-year deprivation gaps (DG) in any of them. For some cancer sites such as colon–rectum, bladder, and lung, the 5-year deprivation gap became negligible in both sensitivity analyses, but a similar trend remained as compared to main analyses. Inconsistent results were observed for prostate cancer, for which the simulated LT not only reduced the social gradient in net survival but also inversed it (DGmain analyses: 3.0, 95% CI: 1.1;4.9, DGSA_Eng LT: −3.5, 95% CI: −5.4;−1.6, DGSA_EDP: −4.9, 95% CI: −6.8;−3.0) and for esophagus for which a social gradient was found only in sensitivity analyses using England LT (DGmain analyses: 0.4, 95% CI: −4.2;5.0, DGSA_Eng LT: 7.0, 95% CI: 1.4;12.6, DGSA_EDP: −0.7, 95% CI: −5.3;3.8).
In women (Figure 1b), clear 5-year deprivation gaps were observed in all three analyses for ENT cancers, bile ducts, esophagus, and cervix uteri. The results showed a trend towards a social gradient in net survival, which was moderate in the three analyses for bladder, corpus uteri, lung, kidney, liver, pancreas, and melanoma, with close values in the 5-year deprivation gaps. Results regarding ovarian cancer, stomach cancer, sarcoma, and CNS cancer were identical according to the three analyses, with negligible 5-year deprivation gaps. The social gradient found in the main analyses was attenuated, but it remained clear with English LT, but not in sensitivity analyses using EDP for breast cancer (DGmain analyses: 5.1, 95% CI: 3.0;7.3, DGSA_Eng LT: 2.5, 95% CI: 0.3;4.6, DGSA_EDP: 1.4, 95% CI: −0.7;3.6) and colorectal cancer (DGmain analyses: 5.4, 95% CI: 2.2;8.6, DGSA_Eng LT: 3.5, 95% CI: 0.4;6.7, DGSA_EDP: 2.8, 95% CI: −0.4;6.0).
For cancers with a worse prognosis (i.e., lung, liver, pancreas, esophagus, bile ducts, and stomach), the results were stable in all three analyses and for both men and women when considering 1-year survival.

3.2. Comparison of Main and Sensitivity Analyses Based on Flexible Modeling

(a)
Same model selected in the three analyses.
As presented in Table 1, the same flexible model was selected by the three analyses for 8/14 cancers in men and 11/17 cancers in women, meaning that the effect of EDI on EMH was identical according to the three analyses. Among these consistent cases, results showed that for four cancers in men (bile ducts, CNS, sarcoma, and stomach) and three cancers in women (kidney, lung, and sarcoma), there was no effect of EDI on EMH according to the three analyses. For four cancers in men (ENT, liver, lung, and melanoma) and five cancers in women (breast, CNS, ENT, esophagus, and liver), EDI had an overall proportional effect on EMH in the three analyses according to model selection (Table 1). For these nine cases, although EHR were reduced in both sensitivity analyses (Figure 2), the curves of EHR as a function of EDI according to the three analyses overlapped for liver, esophagus, and CNS cancers in women, or followed the same pattern with close values between main and sensitivity analyses (ENT in men and women, breast in women, lung in men), though sometimes with loss of EHR significance (according to 95% CI) in sensitivity analyses (liver and melanoma in men). The effect of EDI on EMH was found to be time-dependent for cervical and stomach cancer in woman and time- and/or age-dependent for corpus uteri according to the three analyses (Table 1). Curves of EHR of p90 versus p10 of EDI (EHRp90/p10) as a function of time (Figure 3) according to the three analyses overlapped for cervical and stomach cancer and followed the same pattern with close values for corpus uteri (with lower EHRp90/p10 in sensitivity analyses as compared to main analyses).
(b)
Different models selected in the three analyses.
On the other hand, the selected model was different for 6/14 cancers in men and 6/17 cancers in women (Table 1). Among these cases we can distinguish consistent (i), mitigated (ii), inconsistent (iii), and inconsistent and inversed (iv) results.
(i)
For cancer of the bile ducts in women (Figure 4), an effect of EDI on EMH was found in all three analyses, which was time-dependent in main analyses (with EHRp90/p10 reaching a maximum of 1.96, 95% CI: 1.12;3.43 at 3.8 years of follow-up) and proportional in sensitivity analyses (EHRp90/p10 (SA_Eng LT): 1.38, 95% CI: 1.13;1.69; EHRp90/p10 (SA_EDP): 1.37, 95% CI:1.12;1.68).
(ii)
Mitigated results were observed for three cancers in men (colon–rectum, kidney, and pancreas) and four in women (colon–rectum, pancreas, melanoma, and ovarian).
For colon–rectum cancers in men and melanoma and ovarian cancers in women, sensitivity analysis using EDP showed no effect of EDI based on model selection as opposed to the two other analyses. Regarding colon–rectum in men (Figure 4), main analyses showed a time-dependent effect of EDI on EMH (EHRp90/p10 reaching a maximum of 1.37, 95% CI: 1.12;1.66 at 5 years of follow-up), while sensitivity analysis using English LT presented an overall proportional (according to model selection) but moderate effect (EHRp90/p10: 1.05, 95% CI: 0.99;1.12). Regarding melanoma in women (Figure 4), in the main analysis, the effect of EDI on EMH was proportional (according to model selection) and moderate (herp90/p10: 1.21, 95% CI: 0.86;1.70), while in the sensitivity analysis using English LT, it was time- and/or age-dependent (according to model selection) but moderate according to the 95% CI. Regarding ovarian cancer (see Supplementary Materials, Figure S2), the effect of EDI was found proportional (according to model selection) according to main and sensitivity analysis using English LT, with similar values of EHR as a function of EDI (curves almost overlapped) but a moderate effect of EDI on EMH according to 95% CI in both analyses.
For colon–rectum cancers in women, kidney cancer in men, and pancreas cancer both in men and women, the three analyses found an overall effect of EDI according to the model selection but with partial loss of significance of this effect (according to 95% CI) during the follow-up in at least one of the three analyses. Regarding colon–rectum cancer in women (Figure 4), very similar results were found according to main and sensitivity analysis using English LT, with an overall proportional effect of EDI on EMH (EHRp90/p10 (main analyses): 1.15, 95% CI: 1.06;1.24; EHRp90/p10 (SA_Eng LT): 1.09, 95% CI: 1.02;1.17), whereas a more mitigated situation occurred in sensitivity analyses using EDP, where the effect of EDI on EMH was observed only at the very beginning of follow-up. A similar situation was observed in women regarding pancreas cancers (Figure 4), with a more moderate effect of EDI on the EMH (EHRp90/p10 (main analyses): 1.08, 95% CI: 0.99;1.17; EHRp90/p10 (SA_Eng LT): 1.07, 95% CI: 0.98;1.16). Regarding pancreas cancer in men, the effect of EDI on EMH was observed around the date of diagnosis but not thereafter (time-dependent effect) in main analysis, and an overall proportional (according to model selection) but moderate effect was found in both sensitivity analyses (EHRp90/p10 (SA_Eng LT): 1.06, 95% CI: 0.99;1.15; EHRp90/p10 (SA_EDP): 1.06, 95% CI: 0.98;1.15). Regarding kidney cancer in men, results were contrasted between main and sensitivity analyses (Figure 4), since main analyses reported a proportional effect of EDI on EMH (EHRp90/p10: 1.24, 95% CI: 1.07;1.43), while in both sensitivity analyses a more complex, time-dependent effect was found, with a social gradient at the time of diagnosis, an inverse social gradient at the end of follow-up, and no effect of EDI on EMH in between.
(iii)
Regarding bladder (both in men and women) and esophagus cancers in men, an overall proportional effect of EDI on EMH was found in the main analyses according to model selection (with EHRp90/p10 (bladder, men): 1.24, 95% CI: 1.09;1.4; EHRp90/p10 (bladder, women): 1.14, 95% CI: 0.97;1.35; EHRp90/p10 (esophagus, men): 1.13, 95% CI: 1.01;1.26), but no effect was found in any sensitivity analysis for these cancer sites.
(iv)
Finally, regarding prostate cancer (Figure 5), a time- and/or age-dependent effect of EDI on EMH was found in the main analysis, with an EHR of p90 versus p10 of EDI around 2 for 60- and 70-year-olds, and around 1.5 for 90-year-olds, and no effect at the end of follow-up. In both sensitivity analyses, a time-dependent effect of EDI on EMH was found, with an inverse social gradient reaching a maximum at 5 years of follow-up (EHRp90/p10 (SA_Eng LT): 0.37, 95% CI: 0.19;0.72; EHRp90/p10 (SA_EDP): 0.27, 95% CI: 0.13;0.55).

4. Discussion

In summary, in non-parametric analyses, the results were consistent across the three analyses for 9/14 cancers in men and 15/17 cancers in women. They were inconsistent but with similar trends remaining for 3/14 cancers (bladder, colon–rectum, and lung) in men and for 2/17 cancers in women (breast and colon–rectum), and they were fully inconsistent only for 2 cancers in men (prostate and esophagus). According to flexible modelling, the results were consistent for 8/14 cancers in men and 12/17 cancers in women. They were mitigated for 2/14 cancers in men (colon–rectum and pancreas) and 4/17 cancers in women (colon–rectum, pancreas, ovarian, and melanoma), for which the social gradient was attenuated for all or part of the follow-up, or only in one of the sensitivity analyses, and final conclusions were not contradictory. However, the results were inconsistent for 4/14 cancers in men (esophagus, bladder, kidney, and especially prostate) and 1/17 cancer in women (bladder).
As expected, the social gradient in cancer net survival was reduced in sensitivity analyses using simulated deprivation-specific LT, suggesting an overestimation of the social gradient in cancer net survival in main analyses (using the original French LT not stratified in deprivation). This might result from the absence of consideration of the social gradient in general population mortality due to the social determination of other causes of mortality. Overall, our main results were not contradicted, since for most cancers, a significant social gradient in survival remained in sensitivity analyses or at least a tendency was still observed. Nevertheless, this study alerts us to the fact that the overestimation, in few cases, can be so important that conclusions might be questioned, such as those for prostate cancer. The study should therefore warn researchers about this important problem of accounting for the social gradient of mortality in the general population as well as the importance of developing national deprivation-specific LT to produce accurate results in cancer net survival analyses.
While the overall results were consistent or attenuated in sensitivity analyses, for some cancers the results were surprising and not fully understood. For prostate cancer, which is a cancer with very good survival, an explanation could be that the least-deprived patients do not die from other causes as often as those which are the most deprived. Consequently, they are more likely to die from prostate cancer. Conversely, least-deprived patients are more likely to die from other causes and less likely from prostate cancer, resulting in a reversed social gradient. A study in the United States [35] comparing net survival analyses using national vs. state-specific life tables pointed out a possible bias resulting from men with prostate cancer diagnosed at early stages having a better health status than the US general population. We can also make that hypothesis in our own study, and, also assuming that men with prostate cancer diagnosed at early stages are more likely to be less deprived, it could contribute to explaining the inconsistencies found for prostate cancer. The ecologic bias of the EDI could also contribute to this surprising result. Also, one can suppose that the correction brought by these simulated LT is in fact too strong compared to reality (strong overestimation of the expected social gradient of general mortality). Consequently, when the excess mortality is very small (as in the case of prostate cancer) this could lead to false corrected results. Regarding kidney, bladder, or esophageal cancers, the hypotheses to explain such inconsistent results are less obvious but our findings at least suggest that the social gradient in cancer net survival observed in the main analyses might be in a large part explained by the social gradient in background survival, which is also important to consider in the cancer care management and follow-up of patients as part of their general healthcare. We could also assume that socially-stratified risk factors and comorbidities play a part. Surprisingly, previous research [6,36,37,38] relying on appropriate deprivation-specific LT still found a social gradient in cancer net survival for bladder, kidney, esophagus and prostate (with no inverse gradient). One study [39] (using appropriate deprivation-specific LT) did not find either an association between deprivation or cancer net survival for bladder (in women), kidney and esophagus (in men). However, the issues of social inequalities in health and cancer care management may substantially differ across countries, and the statistical methods used in these studies were not always comparable, which could explain such conflicting results. The overestimation of the correction could also explain the discrepancies with some results of previous studies.
Results from the main and sensitivity analyses were even more consistent among women, possibly resulting from lower social inequalities in mortality in the general population than for men [40,41,42].
Although studies of this kind are scarce, similar sensitivity analyses have been conducted [5,15]. In particular, Antunes et al. [15] also found a disappearance of the social gradient of survival for colorectal cancer after the stratification of LT on the socio-economic level in Portugal, but no other cancer site was studied. Ito et al. [5] also found a systematic reduction in deprivation gaps when using simulated deprivation-specific LT for Japan in sensitivity analyses, but they did not found any inversed social gradient or inconsistent results as compared to us. However, this study was based only on a non-parametric method. Our research brings additional information about the impact of using inappropriate LT in cancer net survival analyses, and it relies on a wide variety of cancers, and large-sample, high-quality population-based data as well as very robust statistical analyses.
One strength of this study was the conduction of comparisons between the three analyses using both a non-parametric method and flexible modelling. This allowed us to obtain a confirmation of the results through two different and validated methods, including a very sophisticated one (flexible modelling of excess mortality hazard using multidimensional penalized splines) which permitted adjustments and the consideration of possible non-linear and non-proportional effects as well as interactions [29,30]. The results from the two methods were consistent, both showing inconsistent results for esophageal and prostate cancer in men. Flexible modelling additionally highlighted the overestimation of the social gradient in the main analyses for kidney cancer in men and bladder cancer in both men and women.
Even though this work allowed us to point out possible cancer sites for which the social gradient in cancer survival was overestimated in our previous study, it is important to note that, the present analyses being based on simulated and not real data, we were not able to precisely quantify the potential bias resulting from the lack of stratification on deprivation of the original French LT nor to shed light on the true extent of the overestimation. Moreover, social inequalities in the background mortality from external sources and thus in the simulated LT were not based on the same variable as in the net survival analyses since the index used in England LT was not the EDI but the IMD (i.e., aggregated ecological index at the contextual level) and the indicator used in EDP was the income (i.e., single indicator at the individual level). In the EDP, the net income per consumption unit of a given household as well as the education level and socio-professional categories were available. All three variables are known to be good indicators of socio-economic position. We considered net income per consumption unit of the household to be the most relevant for use in our analyses, because it represents a continuum and can be classified into five categories according to its natural distribution in the population without the need to create arbitrary categories. Although EDI and IMD are different in terms of construction (e.g., geographic scale) and information provided, they are probably more comparable to each other than income. Indeed, using an individual measure of deprivation such as income to stratify the LT was not the most appropriate way in a study based on an ecological aggregated index such as EDI. The correction of the LT by the income indicator was probably too extreme to study deprivation with EDI. This would explain why the results from the sensitivity analyses using English LT were generally intermediate between results from main analyses (strongest social gradient) and sensitivity analyses using EDP (smallest social gradient). On the other hand, while English LT provided a social gradient of the background mortality based on a quite close index as compared to EDI, this cannot be considered as a perfect reflection of the social gradient in background mortality in France, since the social inequalities in mortality are not the same between these two countries [16]. EDP data were the only source providing information on the social gradient in background mortality for the French population. It therefore seemed interesting to also consider French data, so we chose to use both external sources and produced two sets of simulated LT. Additionally, using deprivation-specific LT from a European country with substantial social inequalities in background mortality (such as Lithuania, according to Mackenbach et al. [16]) would have allowed us to consider a “worst case scenario” for these sensitivity analyses. We also acknowledge some limitations regarding the implementation of the simulated LT. In the set simulated using English LT, we applied the social gradient of mortality of England, but this also changes over time (between 2006 and 2013, duplicating the data from 2011 in the years 2012 and 2013, as described in the Methods section and in the Supplementary Materials) without knowing whether this trend was comparable to that of France, whereas in the LT simulated using EDP data, we assumed that the social gradient of mortality from the period 2012–2016 was stable over time in the French general population.
Regarding other possible weaknesses of the study, we cannot be certain of the absence of noise for some results due to the lack of power after stratifying the analyses by cancer site and by sex. Moreover, data on the stage at cancer diagnosis were not available, and as discussed in Tron et al. 2019 [11], this lack of information limits the interpretation of the results regarding the social gradient in cancer net survival. However, this is unlikely to have a substantial impact on the results of the sensitivity analyses per se.
Given all of these elements, results from neither the main nor the sensitivity analyses can be considered as right or “true” results. These can only be determined when validated deprivation-specific French LT are available, which should be a future research and development prospect, as it has been done recently in Portugal [43]. Also, the problem investigated in this paper only occurs when the estimation of net survival is based on excess mortality rates, and it does not concern studies relying on the causes of death (which have their own limitations).
In conclusion, except for prostate cancer and to a lesser extent bladder (among both men and women), esophageal (among men), and kidney (among men) cancers, this study confirms previously published findings about the disparities in cancer net survival according to the socioeconomic environment for several cancer sites in France, despite the absence of consideration given to the social gradient in background mortality in the LT, which was a main limitation to this previous work. Nevertheless, the present study indicates how important it is to create proper deprivation-specific LT in France in order to accurately investigate the social gradient in cancer net survival, which should be done in the next coming years.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15030659/s1, Supplementary information: Detailed calculation process for the construction of the simulated deprivation-specific French life tables used in the analyses; Figure S1: Mortality rate ratios (RR) according to (a) quintiles of IMD in the England LT (examples for 1990 and 2009 since all the years were available) and (b) the 5 categories of income in the EDP (for the period 2012–2016); Figure S2: Excess mortality hazard ratio (EHR) as a function of deprivation (EDI) according to main analyses and sensitivity analyses using England life tables, for ovary cancer (discordant results between main and sensitivity analyses).

Author Contributions

L.T., L.R., M.F., B.R., A.B., L.L., O.M., F.M., O.D., F.G., and G.L. contributed to the study conception and design. Material preparation and data collection were provided by F.G. Analyses were performed by L.T. The first draft of the manuscript was written by L.T., L.R., M.F., B.R., A.B., L.L., O.M., F.M., O.D., F.G., and G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the ‘Institut de Recherche en Santé Publique’ (IReSP) (no grant number), the Institut National du Cancer (INCa) (no grant number), the ‘Ligue Nationale Contre le Cancer’ (no grant number), and Santé Publique France (SpF) (no grant number).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Consultative Committee for the Processing of Health Research Data (CCTIRS) and the French Data Protection Authority (CNIL, authorization n°913013).

Informed Consent Statement

Consent to participate is not applicable, since patients are not identifiable (anonymous data). Consent for publication is not applicable, since patients are not identifiable (anonymous data).

Data Availability Statement

Data and code are available upon reasonable request.

Acknowledgments

The authors thank Camille Maringe and Libby Ellis for their methodological support in implementing the sensitivity analyses and the Francim working group on “deprivation-specific French life tables simulation”, who contributed to the reflection on this project (authors of the paper and Laurent Roche, Nadine Bossard, Zoé Uhry, Valérie Jooste, Pascale Grosclaude, Véronique Bouvier, Sébastien Lamy, Michael Mounie, Anne-Sophie Woronoff). Members of the French Network of Cancer Registries (FRANCIM): Françoise GALATEAU-SALLE (Registre multicentrique du mésothéliome à vocation nationale (Mesonat)), Anne-Marie BOUVIER (Registre Bourguignon des cancers digestifs), Simona BARA (Registre des cancers de la Manche), Clarisse JOACHIM-CONTARET (Registre des cancers de la Martinique), Bénédicte LAPOTRE-LEDOUX (Registre général des cancers de la Somme), Claire SCHVARTZ (Registre des cancers de la thyroïde Marne Ardennes), Sandrine PLOUVIER (Registre des cancers de Lille et de sa région), Guy LAUNOY (Registre des cancers digestifs du Calvados), Karima HAMMAS (Registre des cancers du Haut-Rhin), Sandrine DABAKUYO (Registre des cancers du sein et des cancers gynécologiques de Côte d’Or), Pascale GROSCLAUDE (Registre des cancers généraux du Tarn), Xavier TROUSSARD (Registre des hémopathies malignes de Basse-Normandie), Marc MAYNA-DIE (Registre des hémopathies malignes de Côte d’Or), Alain MONNEREAU (Registre des hémopathies malignes de la Gironde), Brigitte TRETARRE (Registre général des tumeurs de l’Hérault), Florence MOLINIE and Anne COWPPLI-BONY (Registre des tumeurs de Loire-Atlantique/Vendée), Anne-Sophie WORONOFF (Registre des tumeurs du Doubs et du Territoire de Belfort), Isabelle BALDI (Registre des tumeurs primitives du système nerveux en Gironde), Jean-Baptiste NOUSBBAUM (Registre Finistérien des tumeurs digestives), Gaëlle COUREAU (Registre général des cancers de la Gironde), Jacqueline DELOUMEAUX (Registre général des cancers de la Guadeloupe), Marc COLONNA (Registre général des cancers de l’Isère), Michel VELTEN (Registre général des cancers du Bas-Rhin), Tania D’ALMEIDA (Registre général des cancers en région Limousin), Anne-Valérie GUIZARD (Registre général des tumeurs du Calvados), Jacqueline CLAVEL (Registre national des hémopathies malignes de l’enfant (RNHME)), Brigitte LACOUR (Registre national des tumeurs solides de l’enfant (RNTSE)), Françoise BORSON-CHAZOT (Registre Rhône Alpin des cancers thyroïdiens), Gautier DEFOSSEZ (Registre des cancers de Poitou-Charentes), Sylvie LAUMOND (Cancers généraux-Nouvelle Calédonie), Emmanuel CHIRPAZ (Registre des Cancers de la Réunion), Laure-Manuella DESROZIERS-IMOUNGA (Registre des cancers généraux de Guyane).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dalton, S.O.; Olsen, M.H.; Johansen, C.; Olsen, J.H.; Andersen, K.K. Socioeconomic inequality in cancer survival—Changes over time. A population-based study, Denmark, 1987–2013. Acta Oncol. 2019, 58, 737–744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Exarchakou, A.; Rachet, B.; Belot, A.; Maringe, C.; Coleman, M.P. Impact of national cancer policies on cancer survival trends and socioeconomic inequalities in England, 1996–2013: Population based study. BMJ 2018, 360, k764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Arias-Ortiz, N.E.; de Vries, E. Health inequities and cancer survival in Manizales, Colombia: A population-based study. Colomb. Med. 2018, 49, 63–72. [Google Scholar] [CrossRef] [PubMed]
  4. Hagedoorn, P.; Vandenheede, H.; Vanthomme, K.; Gadeyne, S. Socioeconomic position, population density and site-specific cancer mortality: A multilevel analysis of Belgian adults, 2001–2011. Int. J. Cancer 2018, 142, 23–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ito, Y.; Nakaya, T.; Nakayama, T.; Miyashiro, I.; Ioka, A.; Tsukuma, H.; Rachet, B. Socioeconomic inequalities in cancer survival: A population-based study of adult patients diagnosed in Osaka, Japan, during the period 1993–2004. Acta Oncol. 2014, 53, 1423–1433. [Google Scholar] [CrossRef]
  6. Jansen, L.; Eberle, A.; Emrich, K.; Gondos, A.; Holleczek, B.; Kajüter, H.; Maier, W.; Nennecke, A.; Pritzkuleit, R.; Brenner, H. Socioeconomic deprivation and cancer survival in Germany: An ecological analysis in 200 districts in Germany. Int. J. Cancer 2014, 134, 2951–2960. [Google Scholar] [CrossRef]
  7. Kwak, M.; Kim, C. Disparities by Age, Sex, Tumor Stage, Diagnosis Path, and Area-level Socioeconomic Status in Survival Time for Major Cancers: Results from the Busan Cancer Registry. J. Korean Med. Sci. 2017, 32, 1974–1983. [Google Scholar] [CrossRef]
  8. Fleisch Marcus, A.; Illescas, A.H.; Hohl, B.C.; Llanos, A.A. Relationships between social isolation, neighborhood poverty, and cancer mortality in a population-based study of US adults. PLoS ONE 2017, 12, e0173370. [Google Scholar] [CrossRef] [Green Version]
  9. Skyrud, K.D.; Bray, F.; Eriksen, M.T.; Nilssen, Y.; Møller, B. Regional variations in cancer survival: Impact of tumour stage, socioeconomic status, comorbidity and type of treatment in Norway. Int. J. Cancer 2016, 138, 2190–2200. [Google Scholar] [CrossRef] [Green Version]
  10. Tervonen, H.E.; Aranda, S.; Roder, D.; You, H.; Walton, R.; Morrell, S.; Baker, D.; Currow, D.C. Cancer survival disparities worsening by socio-economic disadvantage over the last 3 decades in new South Wales, Australia. BMC Public Health 2017, 17, 691. [Google Scholar] [CrossRef]
  11. Tron, L.; Belot, A.; Fauvernier, M.; Remontet, L.; Bossard, N.; Launay, L.; Bryere, J.; Monnereau, A.; Dejardin, O.; Launoy, G.; et al. Socioeconomic environment and disparities in cancer survival for 19 solid tumor sites: An analysis of the French Network of Cancer Registries (FRANCIM) data. Int. J. Cancer 2019, 144, 1262–1274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Schaffar, R.; Rachet, B.; Belot, A.; Woods, L.M. Estimation of net survival for cancer patients: Relative survival setting more robust to some assumption violations than cause-specific setting, a sensitivity analysis on empirical data. Eur. J. Cancer 2017, 72, 78–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. 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, 529–538. [Google Scholar] [CrossRef]
  14. Perme, M.P.; Stare, J.; Esteve, J. On estimation in relative survival. Biometrics 2012, 68, 113–120. [Google Scholar] [CrossRef]
  15. Antunes, L.; Mendonça, D.; Bento, M.J.; Rachet, B. No inequalities in survival from colorectal cancer by education and socioeconomic deprivation—A population-based study in the North Region of Portugal, 2000–2002. BMC Cancer 2016, 16, 608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Mackenbach, J.P.; Bopp, M.; Deboosere, P.; Kovacs, K.; Leinsalu, M.; Martikainen, P.; Menvielle, G.; Regidor, E.; Geldera, R. Determinants of the magnitude of socioeconomic inequalities in mortality: A study of 17 European countries. Health Place 2017, 47, 44–53. [Google Scholar] [CrossRef]
  17. Ingleby, F.C.; Woods, L.M.; Atherton, I.M.; Baker, M.; Elliss-Brookes, L.; Belot, A. Describing socio-economic variation in life expectancy according to an individual’s education, occupation and wage in England and Wales: An analysis of the ONS Longitudinal Study. SSM Popul. Health 2021, 14, 100815. [Google Scholar] [CrossRef]
  18. Grafféo, N.; Jooste, V.; Giorgi, R. The impact of additional life-table variables on excess mortality estimates. Stat. Med. 2012, 31, 4219–4230. [Google Scholar] [CrossRef]
  19. Rubio, F.J.; Rachet, B.; Giorgi, R.; Maringe, C.; Belot, A. On models for the estimation of the excess mortality hazard in case of insufficiently stratified life tables. Biostatistics 2019, 22, 51–67. [Google Scholar] [CrossRef] [Green Version]
  20. Goungounga, J.A.; Touraine, C.; Grafféo, N.; Giorgi, R. Correcting for misclassification and selection effects in estimating net survival in clinical trials. BMC Med. Res. Methodol. 2019, 19, 104. [Google Scholar] [CrossRef]
  21. Touraine, C.; Grafféo, N.; Giorgi, R.; Group CWS. More accurate cancer-related excess mortality through correcting background mortality for extra variables. Stat. Methods Med. Res. 2020, 29, 122–136. [Google Scholar] [CrossRef]
  22. Guillaume, E.; Pornet, C.; Dejardin, O.; Launay, L.; Lillini, R.; Vercelli, M.; Marí-Dell, M.; Fontelo, A.F.; Borrell, C.; Ribeiro, A.I.; et al. Development of a cross-cultural deprivation index in five European countries. J. Epidemiol. Community Health 2016, 70, 493–499. [Google Scholar] [CrossRef] [Green Version]
  23. Pornet, C.; Delpierre, C.; Dejardin, O.; Grosclaude, P.; Launay, L.; Guittet, L.; Lang, T.; Launoy, G. Construction of an adaptable European transnational ecological deprivation index: The French version. J. Epidemiol. Community Health 2012, 66, 982–989. [Google Scholar] [CrossRef]
  24. 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–E982. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, Y.C.; Graubard, B.I.; Rosenberg, M.A.; Kuntz, K.M.; Zauber, A.G.; Kahle, L.; Schechter, C.B.; Feuer, E.J. Derivation of background mortality by smoking and obesity in cancer simulation models. Med. Decis. Making 2013, 33, 176–197. [Google Scholar] [CrossRef] [Green Version]
  26. Di Carlo, V.R.B.; Bannon, F.; Woods, L.M.; Maringe, C.; Bonaventure, A.; Coleman, M.P.; Allemani, C. Life Tables for the CONCORD Programme. Available online: http://csg.lshtm.ac.uk/life-tables (accessed on 19 January 2023).
  27. Jordan, H.; Roderick, P.; Martin, D. The Index of Multiple Deprivation 2000 and accessibility effects on health. J. Epidemiol. Community Health 2004, 58, 250–257. [Google Scholar] [CrossRef]
  28. Blanpain, N. “L’espérance de Vie par Niveau de Vie-Méthode et Principaux Résultats” Documents de Travail, n° F1801. INSEE, Paris. 2018. Available online: https://www.insee.fr/fr/statistiques/3322051 (accessed on 19 January 2023).
  29. Fauvernier, M.; Roche, L.; Uhry, Z.; Tron, L.; Bossard, N.; Remontet, L. Multi-dimensional penalized hazard model with continuous covariates: Applications for studying trends and social inequalities in cancer survival. J. R. Stat. Soc. Ser. C Appl. Stat. 2019, 68, 1233–1257. [Google Scholar] [CrossRef]
  30. Remontet, L.; Uhry, Z.; Bossard, N.; Iwaz, J.; Belot, A.; Danieli, C.; Charvat, H.; Roche, L.; The CENSUR Working Survival Group. Flexible and structured survival model for a simultaneous estimation of non-linear and non-proportional effects and complex interactions between continuous variables: Performance of this multidimensional penalized spline approach in net survival trend analysis. Stat. Methods Med. Res. 2019, 28, 2368–2384. [Google Scholar] [CrossRef]
  31. Poiseuil, M.; Tron, L.; Woronoff, A.S.; Trétarre, B.; Dabakuyo-Yonli, T.S.; Fauvernier, M.; Roche, L.; Dejardin, O.; Molinié, F.; Launoy, G.; et al. How do age and social environment affect the dynamics of death hazard and survival in patients with breast or gynecological cancer in France? Int. J. Cancer 2022, 150, 253–262. [Google Scholar] [CrossRef]
  32. Tron, L.; Fauvernier, M.; Bouvier, A.M.; Robaszkiewicz, M.; Bouvier, V.; Cariou, M.; Jooste, V.; Dejardin, O.; Remontet, L.; Alves, A.; et al. Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study. Cancers 2021, 13, 5156. [Google Scholar] [CrossRef]
  33. Wood, S. Generalized Additive Models: An Introduction with R, 2nd ed.; Chapman & Hall/CRC: London, UK, 2017. [Google Scholar]
  34. Fauvernier, M.; Remontet, L.; Uhry, Z.; Bossard, N.; Roche, L. SurvPen: An R package for hazard and excess hazard modelling with multidimensional penalized splines. J. Open Source Softw. 2019, 4, 1434. [Google Scholar] [CrossRef]
  35. Stroup, A.M.; Cho, H.; Scoppa, S.M.; Weir, H.K.; Mariotto, A.B. The impact of state-specific life tables on relative survival. J Natl Cancer Inst. Monogr. 2014, 49, 218–227. [Google Scholar] [CrossRef] [PubMed]
  36. Rachet, B.; Ellis, L.; Maringe, C.; Chu, T.; Nur, U.; Quaresma, M.; Shah, A.; Walters, S.; Woods, L.; Forman, D.; et al. Socioeconomic inequalities in cancer survival in England after the NHS cancer plan. Br. J. Cancer 2010, 103, 446–453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Shafique, K.; Morrison, D.S. Socio-economic inequalities in survival of patients with prostate cancer: Role of age and Gleason grade at diagnosis. PLoS ONE 2013, 8, e56184. [Google Scholar] [CrossRef]
  38. Stanbury, J.F.; Baade, P.D.; Yu, Y.; Yu, X.Q. Cancer survival in New South Wales, Australia: Socioeconomic disparities remain despite overall improvements. BMC Cancer 2016, 16, 48. [Google Scholar] [CrossRef] [Green Version]
  39. Shack, L.G.; Rachet, B.; Brewster, D.H.; Coleman, M.P. Socioeconomic inequalities in cancer survival in Scotland 1986–2000. Br. J. Cancer 2007, 97, 999–1004. [Google Scholar] [CrossRef] [Green Version]
  40. Saurel-Cubizolles, M.J.; Chastang, J.F.; Menvielle, G.; Leclerc, A.; Luce, D. Social inequalities in mortality by cause among men and women in France. J. Epidemiol. Community Health 2009, 63, 197–202. [Google Scholar] [CrossRef] [Green Version]
  41. Layte, R.; Banks, J.; Walsh, C.; McKnight, G. Trends in socio-economic inequalities in mortality by sex in Ireland from the 1980s to the 2000s. Ir. J. Med. Sci. 2015, 184, 613–621. [Google Scholar] [CrossRef]
  42. Mustard, C.A.; Etches, J. Gender differences in socioeconomic inequality in mortality. J. Epidemiol. Community Health 2003, 57, 974–980. [Google Scholar] [CrossRef] [Green Version]
  43. Antunes, L.; Mendonça, D.; Ribeiro, A.I.; Maringe, C.; Rachet, B. Deprivation-specific life tables using multivariable flexible modelling—Trends from 2000–2002 to 2010–2012, Portugal. BMC Public Health 2019, 19, 276. [Google Scholar] [CrossRef]
Figure 1. Difference in 5-year age-standardized net survival (ASNS) between patients living in least-deprived environment (Q1) and patients living in most-deprived environment (Q5) (i.e., 5-Year Deprivation Gap), estimated by non-parametric Pohar–Perme method, for men (a) and women (b), for each cancer site, according to main analyses, sensitivity analyses using English life tables (SA_Eng LT), and sensitivity analyses using EDP data (SA_EDP). ASNS: age-standardized net survival; CNS: central nervous system; EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; ENT: ear, neck, throat (head and neck cancers); Q1: 1st quintile of the national distribution of EDI; Q5: 5th quintile of the national distribution of EDI; SA_EDP: sensitivity analyses based on EDP data; SA_Eng LT: sensitivity analyses based on English life tables.
Figure 1. Difference in 5-year age-standardized net survival (ASNS) between patients living in least-deprived environment (Q1) and patients living in most-deprived environment (Q5) (i.e., 5-Year Deprivation Gap), estimated by non-parametric Pohar–Perme method, for men (a) and women (b), for each cancer site, according to main analyses, sensitivity analyses using English life tables (SA_Eng LT), and sensitivity analyses using EDP data (SA_EDP). ASNS: age-standardized net survival; CNS: central nervous system; EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; ENT: ear, neck, throat (head and neck cancers); Q1: 1st quintile of the national distribution of EDI; Q5: 5th quintile of the national distribution of EDI; SA_EDP: sensitivity analyses based on EDP data; SA_Eng LT: sensitivity analyses based on English life tables.
Cancers 15 00659 g001
Figure 2. Excess mortality hazard ratio (EHR) as a function of deprivation (EDI) according to main analyses, sensitivity analyses using English life tables, and sensitivity analyses using EDP data, for cancer sites for which all three analyses were concordant (same flexible model selected), and model selection showed an overall proportional effect of EDI on excess mortality hazard (independent of age and time of follow-up). EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; ENT: ear, neck, throat (head and neck cancers); LT: life tables.
Figure 2. Excess mortality hazard ratio (EHR) as a function of deprivation (EDI) according to main analyses, sensitivity analyses using English life tables, and sensitivity analyses using EDP data, for cancer sites for which all three analyses were concordant (same flexible model selected), and model selection showed an overall proportional effect of EDI on excess mortality hazard (independent of age and time of follow-up). EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; ENT: ear, neck, throat (head and neck cancers); LT: life tables.
Cancers 15 00659 g002
Figure 3. Excess mortality hazard ratio (EHR) between the 90th and the 10th percentile of the national distribution of EDI as a function of the time since cancer diagnosis, according to main analyses, sensitivity analyses using English life tables and sensitivity analyses using EDP data, for cancer sites for which all three analyses were concordant (same flexible model selected) and model showed an overall time-dependent effect of EDI on excess mortality hazard. For corpus uteri cancer, the effect of EDI on excess mortality was also age-dependent, therefore the results are presented for three values of age corresponding to the 25th, 50th, and 75th percentiles of the age distribution in the studied population with corpus uteri cancer. EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; LT: life tables; p90, p10: 90th and 10th percentile of the national distribution of EDI.
Figure 3. Excess mortality hazard ratio (EHR) between the 90th and the 10th percentile of the national distribution of EDI as a function of the time since cancer diagnosis, according to main analyses, sensitivity analyses using English life tables and sensitivity analyses using EDP data, for cancer sites for which all three analyses were concordant (same flexible model selected) and model showed an overall time-dependent effect of EDI on excess mortality hazard. For corpus uteri cancer, the effect of EDI on excess mortality was also age-dependent, therefore the results are presented for three values of age corresponding to the 25th, 50th, and 75th percentiles of the age distribution in the studied population with corpus uteri cancer. EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; LT: life tables; p90, p10: 90th and 10th percentile of the national distribution of EDI.
Cancers 15 00659 g003
Figure 4. Excess mortality hazard ratio (EHR) between the 90th and the 10th percentile of the national distribution of EDI, as a function of the time since cancer diagnosis according to main and sensitivity analyses, for cancer sites for which the effect of EDI was different according to the analysis (i.e., different models selected in the three analyses) and time-dependent. For melanoma in women, the effect of EDI on excess mortality was also age-dependent in the sensitivity analysis using English life tables, therefore results are presented for 3 values of age corresponding to the 25th, 50th and 75th percentiles of the age distribution in the studied female population with melanoma. EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; LT: life tables; p90, p10: 90th and 10th percentile of the national distribution of EDI.
Figure 4. Excess mortality hazard ratio (EHR) between the 90th and the 10th percentile of the national distribution of EDI, as a function of the time since cancer diagnosis according to main and sensitivity analyses, for cancer sites for which the effect of EDI was different according to the analysis (i.e., different models selected in the three analyses) and time-dependent. For melanoma in women, the effect of EDI on excess mortality was also age-dependent in the sensitivity analysis using English life tables, therefore results are presented for 3 values of age corresponding to the 25th, 50th and 75th percentiles of the age distribution in the studied female population with melanoma. EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; LT: life tables; p90, p10: 90th and 10th percentile of the national distribution of EDI.
Cancers 15 00659 g004aCancers 15 00659 g004b
Figure 5. Excess mortality hazard ratio (EHR) between the 90th and the 10th percentile of the national distribution of EDI, as a function of time since cancer diagnosis and according to the main and sensitivity analyses, for prostate cancer (discordant results between main and sensitivity analyses). EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; LT: life tables; p90, p10: 90th and 10th percentile of the national distribution of EDI.
Figure 5. Excess mortality hazard ratio (EHR) between the 90th and the 10th percentile of the national distribution of EDI, as a function of time since cancer diagnosis and according to the main and sensitivity analyses, for prostate cancer (discordant results between main and sensitivity analyses). EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; EHR: excess mortality hazard ratio; LT: life tables; p90, p10: 90th and 10th percentile of the national distribution of EDI.
Cancers 15 00659 g005
Table 1. The effect of deprivation (EDI) on cancer net survival for each cancer site, from the selected flexible model, according to main analyses, sensitivity analyses using English life tables (SA_Eng LT), and sensitivity analyses using EDP data (SA_EDP).
Table 1. The effect of deprivation (EDI) on cancer net survival for each cancer site, from the selected flexible model, according to main analyses, sensitivity analyses using English life tables (SA_Eng LT), and sensitivity analyses using EDP data (SA_EDP).
MenWomen
Main AnalysesSA_Eng LTSA_EDPMain AnalysesSA_Eng LTSA_EDP
Bile ductNo effectNo effectNo effectTime-dependentProportionalProportional
BladderProportionalNo effectNo effectProportionalNo effectNo effect
Breast---ProportionalProportionalProportional
Cervix uteri---Time-dependentTime-dependentTime-dependent
CNSNo effectNo effectNo effectProportionalProportionalProportional
Colon-rectumS and time-dependentProportionalNo effectProportionalProportionalTime-dependent
Corpus uteri---Time- and/or age-dependentTime- and/or age-dependentTime- and/or age-dependent
ENTProportionalProportionalProportionalProportionalProportionalProportional
EsophagusProportionalNo effectNo effectProportionalProportionalProportional
KidneyProportionalTime-dependentTime-dependentNo effectNo effectNo effect
LiverProportionalProportionalProportionalProportionalProportionalProportional
LungProportionalProportionalProportionalNo effectNo effectNo effect
MelanomaProportionalProportionalProportionalProportionalTime- and/or age-dependentNo effect
Ovary---ProportionalProportionalNo effect
PancreasTime-dependentProportionalProportionalProportionalProportionalTime-dependent
ProstateTime- and/or age-dependentTime-dependentTime-dependent---
SarcomaNo effectNo effectNo effectNo effectNo effectNo effect
StomachNo effectNo effectNo effectTime-dependentTime-dependentTime-dependent
CNS: central nervous system; EDI: European Deprivation Index; EDP: French permanent demographic sample [Echantillon Démographique Permanent]; ENT: ear, neck, throat (head and neck cancers); SA_EDP: sensitivity analyses based on EDP data; SA_Eng LT: sensitivity analyses based on English life tables.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tron, L.; Remontet, L.; Fauvernier, M.; Rachet, B.; Belot, A.; Launay, L.; Merville, O.; Molinié, F.; Dejardin, O.; Francim Group; et al. Is the Social Gradient in Net Survival Observed in France the Result of Inequalities in Cancer-Specific Mortality or Inequalities in General Mortality? Cancers 2023, 15, 659. https://doi.org/10.3390/cancers15030659

AMA Style

Tron L, Remontet L, Fauvernier M, Rachet B, Belot A, Launay L, Merville O, Molinié F, Dejardin O, Francim Group, et al. Is the Social Gradient in Net Survival Observed in France the Result of Inequalities in Cancer-Specific Mortality or Inequalities in General Mortality? Cancers. 2023; 15(3):659. https://doi.org/10.3390/cancers15030659

Chicago/Turabian Style

Tron, Laure, Laurent Remontet, Mathieu Fauvernier, Bernard Rachet, Aurélien Belot, Ludivine Launay, Ophélie Merville, Florence Molinié, Olivier Dejardin, Francim Group, and et al. 2023. "Is the Social Gradient in Net Survival Observed in France the Result of Inequalities in Cancer-Specific Mortality or Inequalities in General Mortality?" Cancers 15, no. 3: 659. https://doi.org/10.3390/cancers15030659

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop