# Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England

Matthew James Smith, Aurélien Belot, Matteo Quartagno, Miguel Angel Luque Fernandez, Audrey Bonaventure, Susan Gachau, Sara Benitez Majano,

Bernard Rachet, Edmund Njeru Njagi

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Matthew James Smith, Aurélien Belot, Matteo Quartagno, Miguel Angel Luque Fernandez, Audrey Bonaventure, et al.. Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England. Cancers, 2021, 13 (22), pp.5805. 10.3390/cancers13225805. inserm-04023230

HAL Id: inserm-04023230

https://inserm.hal.science/inserm-04023230

Submitted on 10 Mar 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England

Matthew James Smith <fnref n="1" />,*, Aurélien Belot <fnref n="1" />, Matteo Quartagno <fnref n="2" />, Miguel Angel Luque Fernandez <fnref n="1" />,<fnref n="3" />,<fnref n="4" />, Audrey Bonaventure <fnref n="5" />, Susan Gachau <fnref n="6" />, Sara Benitez Majano <fnref n="1" />, Bernard Rachet <fnref n="1" /> and Edmund Njeru Njagi <fnref n="1" />

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Smith, M.J.; Belot, A.; Quartagno, M.; Luque Fernandez, M.A.; Bonaventure, A.; Gachau, S.; Benitez Majano, S.; Rachet, B.; Njagi, E.N. Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England. *Cancers* **2021**, *13*, 5805. https://doi.org/10.3390/cancers13225805

Academic Editor: Alberto Zamò

Received: 1 October 2021

Accepted: 16 November 2021

Published: 19 November 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Simple Summary:** Diffuse large B-cell (DLBCL) and follicular lymphoma (FL) account for most non-Hodgkin lymphoma diagnoses: around 35% and 20% in England, respectively. Despite the vast contrast in survival between the subtypes, similar socioeconomic inequalities in survival have persisted over the past two decades, possibly due to the presence of comorbidities. The aim of our study was to assess the association between socioeconomic status and survival from DLBCL or FL accounting for patient and health system characteristics. We found that, for both DLBCL and FL, the most deprived patients had a higher excess mortality hazard compared to the least deprived, regardless of the comorbidity status. Our results show the need for the current framework of the National Health Service to improve the survival of DLBCL and FL patients in the most deprived areas of England, and further consideration is needed for patient-tailored management plans amongst patients with comorbidities or multimorbidities.

**Abstract:** (1) Background: Socioeconomic inequalities of survival in patients with lymphoma persist, which may be explained by patients' comorbidities. We aimed to assess the association between comorbidities and the survival of patients diagnosed with diffuse large B-cell (DLBCL) or follicular lymphoma (FL) in England accounting for other socio-demographic characteristics. (2) Methods: Population-based cancer registry data were linked to Hospital Episode Statistics. We used a flexible multilevel excess hazard model to estimate excess mortality and net survival by patient's comorbidity status, adjusted for sociodemographic, economic, and healthcare factors, and accounting for the patient's area of residence. We used the latent normal joint modelling multiple imputation approach for missing data. (3) Results: Overall, 15,516 and 29,898 patients were diagnosed with FL and DLBCL in England between 2005 and 2013, respectively. Amongst DLBCL and FL patients, respectively, those in the most deprived areas showed 1.22 (95% confidence interval (CI): 1.18–1.27) and 1.45 (95% CI: 1.30–1.62) times higher excess mortality hazard compared to those in the least deprived areas, adjusted for comorbidity status, age at diagnosis, sex, ethnicity, and route to diagnosis. (4) Conclusions: Deprivation is consistently associated with poorer survival among patients diagnosed with DLBCL or FL, after adjusting for co/multimorbidities. Comorbidities and multimorbidities need to be considered when planning public health interventions targeting haematological malignancies in England.

<footnote n="1">¹ Inequalities in Cancer Outcomes Network, Department of Non-Communicable Disease Epidemiology, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK; aurelien.belot@lshtm.ac.uk (A.B.); miguel-angel.luque@lshtm.ac.uk (M.A.L.F.); Sara.benitezmajano@lshtm.ac.uk (S.B.M.); bernard.rachet@lshtm.ac.uk (B.R.); edmund.njeru.njagi@lshtm.ac.uk (E.N.N.)</footnote>

<footnote n="2">² MRC Clinical Trials Unit, Institute of Clinical Trials and Methodology, University College London, London WC1V 6LJ, UK; m.quartagno@ucl.ac.uk</footnote>

<footnote n="3">³ Noncommunicable Disease and Cancer Epidemiology Group, Instituto de Investigación Biosanitaria de Granada, Ibs.GRANADA, Andalusian School of Public Health, 18012 Granada, Spain</footnote>

<footnote n="4">⁴ Centro de Investigación Biomédica en Red de Epidemiología y Salud Pública (CIBER of Epidemiology and Public Health, CIBERESP), 28029 Madrid, Spain</footnote>

<footnote n="5">5 Epidemiology of Childhood and Adolescent Cancers Team, Research Centre in Epidemiology and Biostatistics (CRESS), Inserm UMR 1153, Université de Paris, 94801 Villejuif, France; audrey.bonaventure@inserm.fr</footnote>

<footnote n="6">⁶ School of Mathematics, University of Nairobi, Nairobi 30197-00100, Kenya; sgachau06@gmail.com</footnote>

<footnote>* Correspondence: matthew.smith1@lshtm.ac.uk</footnote>
**Keywords:** cancer epidemiology; diffuse large B-cell lymphoma; follicular lymphoma; survival analysis; comorbidity; multimorbidity; socioeconomic status

# 1. Introduction

Non-Hodgkin lymphoma (NHL) is a heterogeneous group of malignancies, and is currently the 6th most commonly diagnosed cancer in England; in 2014, approximately 32 males and 23 females per 100,000 person years were diagnosed [1]. The heterogeneity in morphology leads to variation in survival probability; for instance, 5-year survival of follicular lymphoma (FL) (86.3%) is higher than diffuse large B-cell lymphoma (DLBCL) (54.8%) [2].

The healthcare system in England aims to offer equitable access to care for all patients. However, variability in health outcomes amongst patients with similar cancers and sociodemographic characteristics still occur; [2–4] convincing reasons for this variability remain a topic of interest. In 2001, the National Health Service (NHS) Cancer Plan [5] recognised, and aimed to reduce, the disparities in survival. Since implementation, there is no evidence that the Plan has had an impact on the inequalities [6,7]. The deprivation gap in survival is still apparent, despite the Plan and successive policies [5,8–10], illustrating the incomplete understanding of the mechanisms underlying these inequalities and raising the concern that these policies have missed the relevant targets.

Patients' comorbidity status may impact timely diagnosis, possibly leading to treatment with more adverse effects [11]; comorbidities are, on average, more prevalent and severe amongst more deprived patients [12]. However, recent evidence indicates that comorbidity explains little of the differential cancer survival between socioeconomic groups [13–15]. Variations in healthcare access, such as location of residence, could partly explain the inequalities [16–20].

Since population-based cancer registries rarely hold reliable information on the cause of death, cancer-specific mortality estimates can be estimated with relative survival methods. These methods compare the mortality hazard (i.e., excess mortality hazard) observed in a population of cancer patients to the mortality hazard observed in the general population with identical demographic characteristics. In this context, the survival estimate derived from the excess mortality hazard is termed net survival (or cancer survival), which is interpreted as the survival where death is due directly, or indirectly, to the cancer studied, and death from other causes has been removed [21].

Overall, the association between comorbidity and cancer survival in patients with DLBCL and FL, accounting for other socio-demographic characteristics and the area of residence, remains unclear. We aim to describe the association between comorbidities and cancer survival amongst DLBCL or FL patients, while accounting for sociodemographic and economic factors, hypothesizing that the presence of comorbidities is associated with poorer survival.

# 2. Methods

## 2.1. *Study Design, Participants, and Data Sources*

We developed a population-based multilevel cohort study of adult patients diagnosed with DLBCL or FL between 1 January 2005 and 31 December 2013 in England. Patients were followed up until death or the end of the study on the 31 December 2015, whichever occurred first.

DLBCL and FL were defined according to the 10th revision of the International Statistical Classification of Diseases and Related Problems (ICD-10 codes C82.0–C85.9) [22]. Morphology (cell type) and topography (tumour site) were defined using renewed updates of the ICD for Oncology (ICD-O); ICD-O-3 [23] was used for diagnoses up to 2010, and ICD-O-3.1 [24] for diagnoses after 2011. Information on patients with DLBCL or FL was collected from the linkage of English cancer registry data, the Cancer Analysis System [25]
(CAS), and Hospital Episode Statistics [26] (HES) data sets within the National Cancer Registry and Analysis Service (NCRAS). These datasets contained detailed information on patient and tumour characteristics (see details below).

## 2.2. *Outcome, Exposure, and Patients' Sociodemographic Characteristics*

The outcome of the study was the time to death, or censoring, among DLBCL and FL patients 5 years after cancer diagnosis. Net survival was deduced after estimating the excess mortality hazard. Hence, we used England life tables stratified by deprivation, sex, age, and calendar year (2005–2013) to account for the overall mortality rate from the background population [27]. As follow up of patients ended in 2015 and life tables were available until 2013, we assumed that the expected mortality rates plateaued for 2014 and 2015.

Comorbidity status was the main exposure. We defined comorbidity as the existence of other chronic medical disorders, in addition to cancer, the primary disease of interest, which are causally unrelated to the primary disease [28,29]. Records from HES were used to identify patients’ comorbidity status based on a computational algorithm published elsewhere [30]. The algorithm searches for the presence of comorbidities retrospectively and defines a time window of 6 to 24 months prior to cancer diagnosis where comorbidities are recorded to avoid bias due to the presence of comorbidities related to cancer (i.e., cardi- ological comorbidities due to DLBCL or FL cancer treatment). Patient comorbidity status was adapted from the original Charlson comorbidity index [31] (CCI). We used the Royal College of Surgeons (RCS) modified Charlson Score (Table A1) [32]. The score removes patients with a previous malignancy to avoid bias, does not assign different weights to comorbidities, and categorises comorbidities as: no comorbidities, one comorbidity, and two or more comorbidities (multimorbidity).

Socio-demographic and economic characteristics were collected from the HES dataset. Age was specified at time of diagnosis. Sex is recorded as male or female. Ethnicity was recorded as white or other. Area-level deprivation, classified into one of five quintiles, was determined by the Index of Multiple Deprivation [33] (IMD), which was based on the Lower Super Output Area [34] (LSOA) residence of the patient at the time of cancer diagnosis. LSOA is a geographical location with a median of 1500 inhabitants. We also include the information regarding patients’ diagnosis path (route to diagnosis), a UK-specific programme, classified as: accident and emergency room diagnosis, general practitioner referral (routine and urgent referrals where the patient was not referred under two-week-wait), two-week-wait (urgent GP referral with a suspicion of cancer), and secondary care diagnosis (other outpatient and inpatient elective routes) [35].

## 2.3. *Statistical Analysis*

We tabulated the sociodemographic characteristics by DLBCL and FL. To estimate the excess mortality hazard, we used a multilevel excess hazard regression model (EHM) with a cubic B-spline with two knots placed at 1 and 3 years after diagnosis for the baseline hazard $\lambda_0(t)$. We accounted for the hierarchical structure of the data via the inclusion of a random effect [36]. The statistical contribution of the random effect to the overall goodness of fit of the model was tested using a likelihood ratio test statistic with a Chi-square mixture distribution [37]. From the estimated excess hazard, we could deduce the net survival via the classical relationship between hazard and survival [38]. Net survival is the survival associated with the cancer under study, after eliminating the other causes of death.

In the EHM we included the following variables: age, sex, comorbidities (categorical, 3 categories), deprivation (categorical, 5 categories), lymphoma subtype, ethnicity, and route of cancer diagnosis. We included the non-linear effect of age using a regression spline (defined using a truncated power basis) with one knot located at 70 years of age. Furthermore, we assumed a time-dependent effect of age at diagnosis, represented by the interaction between B-spline function of time and age. The parameter estimates for the variables were interpreted conditionally on the random effect, i.e., they have a cluster-
specific interpretation, where a cluster refers to a given LSOA. From the model we derived the excess mortality hazard ratios (EMHR) and their respective 95% confidence intervals (CI) for all the categorical variables, and the variance of the random effect for the LSOA. Empirical Bayes estimates of the random effect were used to explore the between-LSOA variability in the excess mortality hazard from DLBCL or FL. The random effect was tested for using a likelihood ratio test, with the reference distribution being a mixture of chi-squared distributions with 0 and 1 degrees of freedom, to account for the well-known boundary problem for random effect variances [39,40].

## 2.4. *Missing Data Analysis*

We explored the missing data mechanism for the two variables with missing data (ethnicity (FL 24.9%, DLBCL: 22.7%) and route (FL: 7.8%, DLBCL 5.0%)). Due to clustered data and partially observed categorical variables, we used the latent normal joint modelling multiple imputation approach, under a missing at random assumption (MAR) [41]. The imputation model included all fully and partially observed variables, vital status indicator, and the Nelson-Aalen estimate of the cumulative overall hazard, and accounted for clustering of patients within lower super output areas. We generated 10 imputed datasets. The multilevel EHM was fitted to each of these datasets, and results combined using Rubin's rules [42,43]. Overall tests for the effects of age after multiple imputation were carried out using the F-based procedure for the test of multiple parameters after multiple imputation [41].

We used R software (version 4.1.2, R Development Core Team, 2020, R: *A language and environment for statistical computing*, R Foundation for Statistical Computing, Vienna, Austria) for all data analyses; the mexhaz [36] package was used for excess hazard modelling and the jomo [44] package for multiple imputation.

## 3. *Results*

Overall, 15,516 (34.2%) patients were diagnosed with FL and 29,898 (65.8%) diagnosed with DLBCL in England between 2005 and 2013 (Table 1). The prevalence of at least one comorbidity was higher amongst DLBCL (10.7%) compared to FL (7.5%). The average age was lower amongst FL compared to DLBCL, 63.9 compared to 67.4 years, respectively. The prevalence of DLBCL was higher amongst deprived areas (16.0%) than FL (14.4%). 'White' was the most prevalent ethnicity for both FL (94.9%) and DLBCL (94.1%). GP referral was the most common route to diagnosis amongst FL (44.0%), whereas amongst DLBCL, A&E was the most common (33.8%).

<table><caption>**Table 1.** Distribution of cancer subtypes by patient and healthcare system characteristics for patients (n = 45,414) diagnosed with non-Hodgkin lymphoma in England during the period 2005–2013.</caption><thead><tr><th rowspan="2">Patient Characteristics</th><th colspan="4">Subtype of NHL</th></tr><tr><th colspan="2">FL</th><th colspan="2">DLBCL</th></tr></thead><tbody><tr><td>Age (mean, SD)</td><td colspan="2">N = 15,516<br>63.9 (13.6)</td><td colspan="2">N = 29,898<br>67.4 (14.9)</td></tr><tr><td>Sex, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Male</td><td>7318</td><td>(47.2%)</td><td>16,215</td><td>(54.2%)</td></tr><tr><td>Female</td><td>8198</td><td>(52.8%)</td><td>13,683</td><td>(45.8%)</td></tr><tr><td>Deprivation quintiles (Q), n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Least deprived (Q1)</td><td>3547</td><td>(22.9%)</td><td>6340</td><td>(21.2%)</td></tr><tr><td>Q2</td><td>3517</td><td>(22.7%)</td><td>6663</td><td>(22.3%)</td></tr><tr><td>Q3</td><td>3294</td><td>(21.2%)</td><td>6246</td><td>(20.9%)</td></tr><tr><td>Q4</td><td>2925</td><td>(18.9%)</td><td>5863</td><td>(19.6%)</td></tr><tr><td>Most deprived (Q5)</td><td>2233</td><td>(14.4%)</td><td>4786</td><td>(16.0%)</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th rowspan="2">Patient Characteristics</th><th colspan="4">Subtype of NHL</th></tr><tr><th colspan="2">FL</th><th colspan="2">DLBCL</th></tr></thead><tbody><tr><td>Comorbidity status, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>No comorbidity</td><td>14,343</td><td>(92.4%)</td><td>26,718</td><td>(89.4%)</td></tr><tr><td>One comorbidity</td><td>641</td><td>(4.1%)</td><td>1570</td><td>(5.3%)</td></tr><tr><td>Multimorbidity</td><td>532</td><td>(3.4%)</td><td>1610</td><td>(5.4%)</td></tr><tr><td>Route of diagnosis, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>GP referral</td><td>6297</td><td>(44.0%)</td><td>8157</td><td>(28.7%)</td></tr><tr><td>A&E</td><td>1869</td><td>(13.1%)</td><td>9617</td><td>(33.8%)</td></tr><tr><td>Secondary care</td><td>2222</td><td>(15.5%)</td><td>3724</td><td>(13.1%)</td></tr><tr><td>TWW</td><td>3912</td><td>(27.4%)</td><td>6918</td><td>(24.4%)</td></tr><tr><td>Missing *</td><td>1216</td><td>(7.8%)</td><td>1482</td><td>(5.0%)</td></tr><tr><td>Ethnicity, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>White</td><td>11,052</td><td>(94.9%)</td><td>21,739</td><td>(94.1%)</td></tr><tr><td>Others</td><td>600</td><td>(5.2%)</td><td>1369</td><td>(5.9%)</td></tr><tr><td>Missing *</td><td>3864</td><td>(24.9%)</td><td>6790</td><td>(22.7%)</td></tr></tbody></table>

GP: general practitioner referral, A&E: accident and emergency room, TWW: two-week-wait. Complete case analysis: missing ethnicity 23.5%; missing route to diagnosis 5.9%. * Proportions are of the total number of patients.

In the multivariable analysis (Table 2), amongst DLBCL, and after multiple imputation, patients with comorbidity and multimorbidity showed 23% and 40% increased excess mortality compared to patients without comorbidity (i.e., EMHR: 1.23; 95% CI: 1.14–1.32, and EMHR: 1.40; CI: 1.01–1.94, respectively). Patients living in the most deprived areas had 1.22 (95% CI: 1.18–1.27) times higher excess mortality than those living in the least deprived areas. Patients diagnosed through A&E had nearly three times a higher excess mortality compared to GP referral (i.e., EMHR: 2.75; 95% CI: 2.54–2.98). Females had a significantly lower excess mortality compared to males (i.e., EMHR 0.93; 95% CI: 0.90–0.96). There was, however, no evidence of a difference in excess mortality by ethnicity (Table 2). Using a likelihood ratio test (a mixture of chi-square distributions) there was strong evidence (*p* < 0.001) that including the random effect improved the fit of the model.

<table><caption>Table 2. Adjusted excess mortality hazard ratios for age, sex, deprivation, comorbidity, cancer subtype, route of diagnosis, ethnicity, and LSOA as random intercept for (i) complete case analysis, and (ii) after multiple imputation for patients (n = 29,898) diagnosed with diffuse large B-cell lymphoma in England during the period 2005–2013.</caption><thead><tr><th rowspan="2">Patient Characteristics</th><th colspan="3">Model (i): Complete Case</th><th colspan="3">Model (ii): After Imputation</th></tr><tr><th>HR</th><th>CI</th><th>p-Value</th><th>HR</th><th>CI</th><th>p-Value</th></tr></thead><tbody><tr><td>Sex</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Male</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>Female</td><td>0.93</td><td>0.89–0.98</td><td>0.003</td><td>0.93</td><td>0.90–0.96</td><td>&lt;0.001</td></tr><tr><td>Ethnicity</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>White</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>Other</td><td>0.97</td><td>0.87–1.08</td><td>0.556</td><td>0.99</td><td>0.91–1.08</td><td>0.809</td></tr><tr><td>Deprivation quintiles (Q)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Least deprived Q1</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>Q2</td><td>1.03</td><td>0.96–1.11</td><td>0.372</td><td>1.00</td><td>0.93–1.08</td><td>0.922</td></tr><tr><td>Q3</td><td>1.08</td><td>1.00–1.16</td><td>0.045</td><td>1.07</td><td>1.00–1.14</td><td>0.045</td></tr><tr><td>Q4</td><td>1.17</td><td>1.08–1.26</td><td>&lt;0.001</td><td>1.13</td><td>1.04–1.23</td><td>0.003</td></tr><tr><td>Most deprived Q5</td><td>1.26</td><td>1.16–1.37</td><td>&lt;0.001</td><td>1.22</td><td>1.18–1.27</td><td>&lt;0.001</td></tr><tr><td>Comorbidity status</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No comorbidity</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>One comorbidity</td><td>1.26</td><td>1.15–1.38</td><td>&lt;0.001</td><td>1.23</td><td>1.14–1.32</td><td>&lt;0.001</td></tr><tr><td>Multimorbidity</td><td>1.50</td><td>1.38–1.64</td><td>&lt;0.001</td><td>1.40</td><td>1.01–1.94</td><td>0.043</td></tr></tbody></table>
<table><caption>**Table 2.** *Cont.*</caption><thead><tr><th rowspan="2">Patient Characteristics</th><th colspan="3">Model (i): Complete Case</th><th colspan="3">Model (ii): After Imputation</th></tr><tr><th>HR</th><th>CI</th><th>p-Value</th><th>HR</th><th>CI</th><th>p-Value</th></tr></thead><tbody><tr><td>Route of diagnosis</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>GP referral</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>A&E</td><td>2.75</td><td>2.60-2.91</td><td>&lt;0.001</td><td>2.75</td><td>2.54-2.98</td><td>&lt;0.001</td></tr><tr><td>Secondary Care</td><td>1.43</td><td>1.22-1.67</td><td>&lt;0.001</td><td>1.23</td><td>1.11-1.36</td><td>&lt;0.001</td></tr><tr><td>TWW</td><td>1.33</td><td>1.23-1.45</td><td>&lt;0.001</td><td>0.83</td><td>0.56-1.24</td><td>0.362</td></tr><tr><td>Random Effect</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SD (SE)</td><td>0.48 (0.08)</td><td>-</td><td>-</td><td>0.39 (0.04)</td><td>-</td><td>-</td></tr></tbody></table>

GP: general practitioner referral. A&E: accident and emergency room. TWW: two-week-wait.

In the multivariable analysis (Table 3), amongst FL, patients with comorbidity and multimorbidity showed 1.52 and 2.19 times the excess mortality compared to patients without comorbidity (i.e., EMHR: 1.52; 95% CI: 1.25–1.84, and EMHR: 2.19; CI: 1.45–3.31, respectively). Patients living in the most deprived areas had 1.45 (95% CI: 1.30–1.62) times higher excess mortality than those living in the least deprived areas. Patients diagnosed through A&E had nearly three times a higher excess mortality compared to GP referral (i.e., EMHR: 3.32; 95% CI: 2.49–4.43). Females had a significantly lower excess mortality compared to males (i.e., EMHR 0.89; 95% CI: 0.81–0.97). There was, however, no evidence of a difference in excess mortality by ethnicity (Table 3). Using a likelihood ratio test (a mixture of chi-square distributions), there was strong evidence (*p* < 0.001) that including the random effect improved the fit of the model.

<table><caption>**Table 3.** Adjusted excess mortality hazard ratios for age, sex, deprivation, comorbidity, cancer subtype, route of diagnosis, ethnicity, and LSOA as random intercept for (i) complete case analysis, and (ii) after multiple imputation for patients (*n* = 15,516) diagnosed with follicular lymphoma in England during the period 2005–2013.</caption><thead><tr><th rowspan="2">Characteristics</th><th colspan="3">Model (i): Complete Case</th><th colspan="3">Model (ii): After Imputation</th></tr><tr><th>HR</th><th>CI</th><th>p-Value</th><th>HR</th><th>CI</th><th>p-Value</th></tr></thead><tbody><tr><td>Sex</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Male</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>Female</td><td>0.86</td><td>0.76-0.96</td><td>0.010</td><td>0.89</td><td>0.81-0.97</td><td>0.009</td></tr><tr><td>Ethnicity</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>White</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>Other</td><td>0.59</td><td>0.41-0.83</td><td>0.003</td><td>0.76</td><td>0.60-0.96</td><td>0.019</td></tr><tr><td>Deprivation quintiles (Q)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Least deprived Q1</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>Q2</td><td>1.09</td><td>0.91-1.31</td><td>0.364</td><td>1.10</td><td>0.92-1.32</td><td>0.309</td></tr><tr><td>Q3</td><td>1.23</td><td>1.02-1.48</td><td>0.030</td><td>1.11</td><td>0.96-1.29</td><td>0.166</td></tr><tr><td>Q4</td><td>1.37</td><td>1.13-1.65</td><td>0.001</td><td>1.34</td><td>1.06-1.69</td><td>0.015</td></tr><tr><td>Most deprived Q5</td><td>1.69</td><td>1.38-2.06</td><td>&lt;0.001</td><td>1.45</td><td>1.30-1.62</td><td>&lt;0.001</td></tr><tr><td>Comorbidity status</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No comorbidity</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>One comorbidity</td><td>1.51</td><td>1.19-1.91</td><td>&lt;0.001</td><td>1.52</td><td>1.25-1.84</td><td>&lt;0.001</td></tr><tr><td>Multimorbidity</td><td>2.38</td><td>1.90-3.00</td><td>&lt;0.001</td><td>2.19</td><td>1.45-3.31</td><td>&lt;0.001</td></tr><tr><td>Route of diagnosis</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>GP referral</td><td>Ref</td><td>Ref</td><td></td><td>Ref</td><td>Ref</td><td></td></tr><tr><td>A&E</td><td>3.18</td><td>2.69-3.76</td><td>&lt;0.001</td><td>3.32</td><td>2.49-4.43</td><td>&lt;0.001</td></tr><tr><td>Secondary Care</td><td>1.27</td><td>0.86-1.90</td><td>0.233</td><td>1.22</td><td>0.96-1.55</td><td>0.107</td></tr><tr><td>TWW</td><td>1.17</td><td>0.98-1.40</td><td>0.084</td><td>1.06</td><td>0.63-1.78</td><td>0.830</td></tr><tr><td>Random Effect</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SD (SE)</td><td>0.87 (0.14)</td><td>-</td><td>-</td><td>0.69 (0.16)</td><td>-</td><td>-</td></tr></tbody></table>

GP: general practitioner referral. A&E: accident and emergency room. TWW: two-week-wait.
Figures 1 and 2 show the EMHR for patients with DLBCL and FL, respectively, according to age at diagnosis at different time since diagnosis (Figures 1A and 2A), and according to time since diagnosis for different age at diagnosis (Figures 1B and 2B). The excess mortality hazard for DLBCL and FL patients for different values of age at diagnosis is shown in the Appendix A (Figures A1 and A2, respectively). These plots were obtained from the three-dimensional plots of EMHR, as shown in the Appendix A (Figures A3 and A4, respectively). For DLBCL (Figure 1), the EMHR was higher for older patients whatever the follow-up time (Figure 1A). For those of older or younger ages, in comparison to 70-year-olds, the EMHR was markedly different immediately after, or at 5 years since, diagnosis, but was most similar around 18 months after diagnosis (Figure 1B).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Excess mortality hazard ratios according to (A) age at diagnosis at different time since diagnosis (3 months, 1 and 5 years), and (B) time since diagnosis for different age groups, amongst patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.</figcaption></figure>

For FL (Figure 2), the non-linear effect of age was almost similar whatever the time since diagnosis; being older was associated with a higher excess mortality hazard (Figure 2A). For those of older or younger ages, in comparison to 70-year-olds, the EMHR was markedly different immediately after or at 5 years since diagnosis but was most similar around 18 months after diagnosis (Figure 2B).
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Excess mortality hazard ratios according to (A) age at diagnosis at different time since diagnosis (3 months, 1 and 5 years), and (B) time since diagnosis for different age groups, amongst patients diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.</figcaption></figure>

Figures 3 and 4 show the net survival probability as predicted from the regression model amongst patients with DLBCL and FL, respectively. Amongst DLBCL patients (Figure 3), those living in more deprived areas experienced approximately 7% lower 5-year survival compared to patients in the least deprived areas (e.g., 5-year net survival, amongst those without comorbidities, was 56% for the least deprived compared to 49% for the most deprived). Amongst FL patients (Figure 4), those living in more deprived areas experienced approximately 4% lower 5-year survival compared to those living in the least deprived areas (e.g., 5-year net survival, amongst those without comorbidities, was 86% for the least deprived compared to 82% for the most deprived). For DLBCL only (Figure 3), the deprivation gap in survival was apparent from approximately 6 months after diagnosis, regardless of the comorbidity status.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Net survival model-based prediction for diffuse large B-cell lymphoma for each comorbidity status by deprivation level (n = 29,898) in England between 2005 and 2013. The values here are presented for 70-year-old white males diagnosed via a general practitioner referral. Values will change for other covariates, but the pattern observed here will remain.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Net survival model-based prediction for follicular lymphoma for each comorbidity status by deprivation level (n = 15,516) in England between 2005 and 2013. The values here are presented for 70-year-old white males diagnosed via a general practitioner referral. Values will change for other covariates, but the pattern observed here will remain.</figcaption></figure>
In the Appendix A, we graphically illustrate the empirical Bayes (EB) estimates of the LSOA random effect for the excess mortality hazard from DLBCL and FL (Figures A5 and A6, respectively). A positive EB estimate indicated a higher excess mortality hazard for a patient from that LSOA in comparison with a patient who has similar observed characteristics but from an LSOA with either a less positive, or negative EB estimate. The EB estimates were grouped by deprivation level, to which the LSOA contributed. For both DLBCL and FL (Figures A5 and A6, respectively), the results show there were no outliers and approximately equal distribution of the EB estimates for each deprivation level.

**4. Discussion**

We found strong evidence of a higher excess mortality amongst DLBCL and FL patients diagnosed with comorbidities compared to patients without comorbidities after adjusting for age, deprivation level, ethnicity, and route to diagnosis and accounting for the patient's area of residence; we also found a noticeable deprivation gap in cancer survival.

Differences in access to treatments, or risk of adverse effects, may explain some of the disparities in survival among DLBCL patients. Immunotherapy (rituximab) for the treatment of aggressive lymphomas (e.g., DLBCL) is known to be effective for those of an advanced age [45–47]. Rituximab is often used in combination with doxorubicin, an increase in dosage of which is associated with an increased incidence of adverse effects (cardiotoxicity), such as congestive heart failure [48]. Guidelines based on National Institute for Health and Care Excellence (NICE) recommend that patients at risk of cardiotoxicity, or low tolerance of intensive therapy, consider a less-intensive treatment regimen [49–51]. This less-intensive treatment allocation may partly explain the comorbidity inequalities in survival from DLBCL. For patients with FL, the standard management is ‘watch-and-wait’; thus, in the absence of a treatment, the comorbidity inequalities in survival may be largely explained by the presence of a comorbidity itself rather than being explained by the effect of comorbidity on treatment. Novel treatment strategies are being explored to ascertain the survival and quality of life benefits in comparison to current standards of care [52]; however, access to these treatments will depend on (i) the specialist centre available to the patient, and (ii) the patients’ influence on the decision of a treatment allocation. While the results of clinical trials provide insights into the efficacy of a treatment, they may lack external validity when the treatment is administered within the healthcare system. Pragmatic trials could be developed to understand the real-world benefit of NHL treatments, and how they are delivered, to patients of older ages or with underlying health conditions who would often not be eligible for clinical trials [53]. In addition, real-world data are becoming increasingly more available and would be of great utility to evaluate the effectiveness of treatment on the whole population [54,55].

For FL patients, we showed that the excess mortality hazard among older patients compared to younger patients is highest after just 4 years since diagnosis (Figure 2B). Since we accounted for background population mortality, and adjusted for comorbidity, the higher excess hazard could be because of histological transformation from lower to higher grades of FL. Studies suggest the risk of histological transformation increases by 3% per year since cancer onset [56]. Thus, the increased excess hazard amongst older patients may be because histological transformation, which complicates the treatment and management of FL. The effects of time-varying variables (e.g., progression/relapse, treatment response, transplant, etc.) on the excess mortality hazard would be strong predictive factors, but this information was not available from the data. This represents an interesting topic for further investigation, where we could either use an extension of the flexible parametric models for time-dependent variables, or a landmark approach [57].

The importance of understanding the association of comorbid conditions with cancer patients' outcomes has been well documented [58]. To our knowledge, this is the first study of England cancer registry data investigating survival by comorbidity status among DLBCL and FL patients. Our results are consistent with previous findings from a Danish study, which showed that the hazard of death increased with severity of comorbidity status [59];
however, the study did not account for missing data and the association with comorbidities was potentially overestimated. The EHR associated with comorbidity decreased after accounting for missing data. The deprivation gap in survival persists even after accounting for prognostic factors such as comorbidity [59–61]. Smith et al. [3] reported no deprivation gap in survival; however, their study may have lacked power, and their study used the relative survival ratio, which can be biased over longer-term follow up [62].

Consistent with previous studies [4], survival after GP referral (non-emergency) diagnosis is significantly better compared with A&E. However, our study also finds that patients diagnosed through TWW, who would be expected to have worse symptoms and survival, showed no evidence of a difference in survival compared with GP referral. There are two possible reasons for the absence of a difference in the associations. Firstly, GPs could advocate for a prompt referral even though the patient is not on the TWW pathway, resulting in patients with similar access to healthcare facilities. Secondly, on the other hand, patients referred through the TWW pathway have more severe symptoms and are expected to have a higher excess hazard. Our results show no difference in the excess mortality, indicating that the TWW pathway prevents patients with more severe symptoms from having a higher excess hazard. This suggests that the performance of TWW pathway is at least as beneficial to a patient's survival as GP referral. Other studies have suggested ways to improve outcomes for patients diagnosed with comorbidities, which include: novel treatment strategies [63], the inclusion of elderly patients in clinical trials [64,65], and the investigation of dose allocation amongst those with higher comorbidity scores [66]. However, further factors associated with the interactions between comorbidities and health care systems leading to poorer survival among DLBCL and FL cancer patients need to be studied.

The strengths of this study are that, firstly, we used a large population-based sample size obtained from cancer registry databases linked to HES, which encompasses all patients in England with a diagnosis of DLBCL and FL between 2005 and 2013. HES data encapsulate a national coverage of comorbidities diagnosed during hospital admission and may have missed comorbidities diagnosed during primary care (e.g., diabetes diagnosed during a GP consultation). However, the addition of information provided from comorbidity records captured during primary care may not improve the prediction of cancer patient survival beyond what is captured in HES data [67]. For example, information on comorbidities, such as diabetes, diagnosed outside of hospital admission are likely to have a minimal impact on the prediction of survival beyond information captured in HES. Secondly, we used the Royal College of Surgeons' adaptation [32] of the Charlson comorbidity score, which provides a more valid measure of the patient's comorbidity status, because it was developed within the England population healthcare data setting. Thirdly, we used a latent normal joint modelling multiple imputation to treat missing data in ethnicity and diagnostic route. This approach allows the imputation of a mix of variable types, while accounting for multilevel structures arising from the clustering of patients within LSOAs [41,68,69]. We assumed that missing data on partially observed variables were missing at random, given the observed variables; further analysis could explore the violation of this assumption and impute under a not-missing-at-random assumption.

This study has its limitations. Firstly, individual-level socioeconomic measures are recommended in addition to area-level measures [70]. Information on individual-level socioeconomic status was unavailable, but using area-level measures captures the multidimensional composition of a patient's deprivation level, in addition to the contextual level [33,71]. Furthermore, using area-level measures, there is greater consistency in the measurement of deprivation between time periods, because deprivation scores have a high concordance amongst updates [33]. Secondly, due to data availability, we did not include tumour stage, which may have partly explained the socioeconomic inequalities in survival. However, even though reliable estimates can be obtained after the multiple imputation of partially observed variables with high proportions of missing data [72], the inclusion of tumour stage may not have provided further information for the prediction of
survival beyond that of diagnostic route, because late cancer stage is strongly associated
with delayed diagnostic route [73].

Survival at 1 and 5 years since diagnosis of DLBCL and FL in England trails that of other European countries [74]; however, restricting estimates to those surviving at least 1 year after diagnosis (conditional survival) shows a comparable 5-year survival [75]. This indicates that long-term survival differences are largely explained by the increased short-term mortality. Understanding long-term survival from FL is more complex due to the histological transformation of indolent lymphomas, which would require an adaptation of the treatment, support, and management from healthcare facilities. This adaptation could be compounded by the patient's susceptibility to cardiotoxic treatments. Further studies could focus on the mechanisms and inequalities of short-term mortality, the long-term survival of patients with transformed lymphomas, and the survival of patients at risk of cardiotoxicity.

**5. Conclusions**

After accounting for sociodemographic factors, healthcare factors, socioeconomic
deprivation, and the patient's area of residence, comorbidities were consistently associated
with poorer survival and an increased excess mortality amongst patients with DLBCL
or FL in England. Furthermore, survival inequalities between socioeconomic levels in
patients with DLBCL or FL persisted after accounting for the presence of comorbidities and
multimorbidities. These results show the need for the current framework of the National
Health Service to improve the survival of DLBCL and FL patients in the most deprived
areas of England, and further consideration is needed for patient-tailored management
plans amongst patients with comorbidities or multimorbidities.

**Author Contributions:** M.J.S., E.N.N., A.B. (Aurélien Belot) and B.R. contributed to the conception of the study and designed the study. E.N.N., A.B. (Aurélien Belot), B.R., S.G. and M.Q. provided advice on statistical methods. M.J.S. conducted the analyses of the data and prepared the draft of the manuscript, tables, and figures. E.N.N., A.B. (Aurélien Belot), B.R. and M.A.L.F. supervised the study and provided comments on the manuscript draft. E.N.N., B.R., M.Q., M.A.L.F., S.G., S.B.M., A.B. (Aurélien Belot) and A.B. (Audrey Bonaventure) provided comments on the final draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Cancer Research UK (grant number C7923/A18525) and the Cancer Research UK Studentship in Cancer Survival (grant number C7923/A20987). The analyses, and the writing of the manuscript were solely the responsibility of the authors. The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the views of Cancer Research UK.

**Institutional Review Board Statement:** We obtained the statutory approvals required for this re-
search from the Confidentiality Advisory Group (CAG) of the Health Research Authority (HRA):
PIAG 1–05(c) 2007. Ethical approval was obtained from the Research Ethics Committee (REC) of the
Health Research Authority (HRA): 07/MRE01/52.

**Informed Consent Statement:** Informed consent from participants was waived by the ethics committee. We used anonymised National Cancer Registry and Hospital Episode Statistics data. All methods were carried out in accordance with relevant guidelines and regulations.

**Data Availability Statement:** The data that support the findings of this study are available via application to the Public Health England Office for Data Release, but restrictions apply to the availability of these data.

**Acknowledgments:** We would like to thank Adrian Turculet, Data Manager of the LSHTM Inequalities in Cancer Outcomes Network, for his support and assistance with the data linkage.

**Conflicts of Interest:** The authors declare no potential conflicts of interest. The authors declare no support from any organisations for the submitted work. The design of the study.
Appendix A

<figure><img src="image_9.png" /><figcaption>**Figure A1.** Excess mortality hazard (i.e., white males, least deprived, no comorbidities, diagnosed through general practitioner referral within an average LSOA (random effect of zero)) over time since diagnosis, for different ages, amongst those diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure A2.** Excess mortality hazard (i.e., white males, least deprived, no comorbidities, diagnosed through general practitioner referral within an average LSOA (random effect of zero)) over time since diagnosis, for different ages, amongst those diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure A3.** Excess mortality hazard ratio according to age at diagnosis and time since diagnosis for patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England between 2005 and 2013.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure A4.** Excess mortality hazard ratio according to age at diagnosis and time since diagnosis for patients diagnosed with follicular lymphoma (n = 15,516) in England between 2005 and 2013.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure A5.** Empirical Bayes estimates of the random effect of LSOA from the excess mortality hazard model for patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure A6.** Empirical Bayes estimates of the random effect of LSOA from the excess mortality hazard model for patients diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.</figcaption></figure>
<table><caption>**Table A1.** Comorbidities and their diagnostic ICD-10 codes.</caption><thead><tr><th>Comorbidity</th><th>ICD-10</th></tr></thead><tbody><tr><td>Myocardial infarction</td><td>I21.x, I22.x, I25.2</td></tr><tr><td>Congestive heart failure</td><td>I11.0, I13.0, I13.2, I25.5, I42.0, I42.5-I42.9, I43.x, I50.x, P29.0</td></tr><tr><td>Peripheral vascular disease</td><td>I70.x, I71.x, I73.1, I73.8, I73.9, I77.1, I79.0, I79.2, K55.1, K55.8, K55.9, Z95.8, Z95.9</td></tr><tr><td>Cerebrovascular disease</td><td>G45.x, G46.x, H34.0, I60.x-I69.x</td></tr><tr><td>Dementia</td><td>F00.x-F03.x, F05.1, G30.x, G31.1</td></tr><tr><td>Chronic obstructive pulmonary disease</td><td>I27.9, J40.x-J47.x, J60.x-J67.x, J68.4, J70.1, J70.3</td></tr><tr><td>Rheumatic disease</td><td>M05.x, M06.x, M31.5, M32.x-M34.x, M35.1, M35.3, M36.0</td></tr><tr><td>Liver disease</td><td>B18.x, K70.0-K70.3, K70.9, K71.3-K71.5, K71.7, K73.x, K74.x, K76.0, K76.2-K76.4, K76.8, K76.9, Z94.4, K71.1, K72.1, K72.9, K76.5, K76.6, K76.7, I85.0, I85.9, I86.4, I98.2, K70.4</td></tr><tr><td>Diabetes without chronic complication</td><td>E10.0, E10.1, E10.6, E10.8, E10.9, E11.0, E11.1, E11.6, E11.8, E11.9, E12.0, E12.1, E12.6, E12.8, E12.9, E13.0, E13.1, E13.6, E13.8, E13.9, E14.0, E14.1, E14.6, E14.8, E14.9</td></tr><tr><td>Diabetes with chronic complication</td><td>E10.7, E11.2-E11.5, E11.7, E12.2-E12.5, E12.7, E13.2-E13.5, E13.7, E14.2-E14.5, E14.7</td></tr><tr><td>Hemiplegia or paraplegia</td><td>G04.1, G11.4, G80.1, G80.2, G81.x, G82.x, G83.0-G83.4, G83.9</td></tr><tr><td>Renal disease</td><td>I12.0, I13.1, N03.2-N03.7, N05.2-N05.7, N18.x, N19.x, N25.0, Z49.0-Z49.2, Z94.0, Z99.2</td></tr><tr><td>AIDS/HIV</td><td>B20.x-B22.x, B24.x</td></tr></tbody></table>

ICD-10: International Classification of Diseases, 10th Revision. Diabetes with/without chronic complication is combined in the RCS Charlson Comorbidity Score.

## References

1. Smittenaar, C.R.; Petersen, K.A.; Stewart, K.; Moitt, N. Cancer Incidence and Mortality Projections in the UK until 2035. *Br. J. Cancer* **2016**, *115*, 1147–1155. [CrossRef] [PubMed]

2. Smith, A.; Crouch, S.; Lax, S.; Li, J.; Painter, D.; Howell, D.; Patmore, R.; Jack, A.; Roman, E. Lymphoma Incidence, Survival and Prevalence 2004–2014: Sub-Type Analyses from the UK’s Haematological Malignancy Research Network. *Br. J. Cancer* **2015**, *112*, 1575–1584. [CrossRef] [PubMed]

3. Smith, A.; Crouch, S.; Howell, D.; Burton, C.; Patmore, R.; Roman, E. Impact of Age and Socioeconomic Status on Treatment and Survival from Aggressive Lymphoma: A UK Population-Based Study of Diffuse Large B-Cell Lymphoma. *Cancer Epidemiol. 2015*, *39*, 1103–1112. [CrossRef]

4. Kane, E.; Howell, D.; Smith, A.; Crouch, S.; Burton, C.; Roman, E.; Patmore, R. Emergency Admission and Survival from Aggressive Non-Hodgkin Lymphoma: A Report from the UK’s Population-Based Haematological Malignancy Research Network. *Eur. J. Cancer* **2017**, *78*, 53–60. [CrossRef] [PubMed]

5. The NHS Cancer Plan: A Plan for Investment, A Plan for Reform. Available online: https://www.thh.nhs.uk/documents/_Departments/Cancer/NHSCancerPlan.pdf (accessed on 24 March 2021).

6. 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. [CrossRef] [PubMed]

7. Maringe, C.; Li, R.; Mangtani, P.; Coleman, M.P.; Rachet, B. Cancer Survival Differences between South Asians and Non-South Asians of England in 1986–2004, Accounting for Age at Diagnosis and Deprivation. *Br. J. Cancer* **2015**, *113*, 173. [CrossRef] [PubMed]

8. Improving Outcomes: A Strategy for Cancer. Available online: https://www.gov.uk/government/publications/the-national-cancer-strategy (accessed on 24 March 2021).

9. Improving Outcomes in Haematological Cancers: The Manual. Available online: https://www.nice.org.uk/guidance/NG47/documents/supporting-evidence (accessed on 29 March 2020).

10. Haematological Cancers: Improving Outcomes. Available online: https://www.nice.org.uk/guidance/ng47 (accessed on 29 March 2020).

11. Renzi, C.; Lyratzopoulos, G.; Hamilton, W.; Maringe, C.; Rachet, B. Contrasting Effects of Comorbidities on Emergency Colon Cancer Diagnosis: A Longitudinal Data-Linkage Study in England. *BMC Health Serv. Res.* **2019**, *19*, 311. [CrossRef] [PubMed]

12. Fowler, H.; Belot, A.; Ellis, L.; Maringe, C.; Luque-Fernandez, M.A.; Njagi, E.N.; Navani, N.; Sarfati, D.; Rachet, B. Comorbidity Prevalence among Cancer Patients: A Population-Based Cohort Study of Four Cancers. *BMC Cancer* **2020**, *20*, 2. [CrossRef] [PubMed]

13. Li, R.; Daniel, R.; Rachet, B. How Much Do Tumor Stage and Treatment Explain Socioeconomic Inequalities in Breast Cancer Survival? Applying Causal Mediation Analysis to Population-Based Data. *Eur. J. Epidemiol.* **2016**, *31*, 603–611. [CrossRef]
14. Fowler, H.; Belot, A.; Njagi, E.N.; Luque-Fernandez, M.A.; Maringe, C.; Quaresma, M.; Kajiwara, M.; Rachet, B. Persistent Inequalities in 90-Day Colon Cancer Mortality: An English Cohort Study. *Br. J. Cancer* **2017**, *117*, 1396. [CrossRef] [PubMed]

15. Belot, A.; Fowler, H.; Njagi, E.N.; Luque-Fernandez, M.-A.; Maringe, C.; Magadi, W.; Exarchakou, A.; Quaresma, M.; Turculet, A.; Peake, M.D.; et al. Association between Age, Deprivation and Specific Comorbid Conditions and the Receipt of Major Surgery in Patients with Non-Small Cell Lung Cancer in England: A Population-Based Study. *Thorax* **2019**, *74*, 51–59. [CrossRef] [PubMed]

16. 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. [CrossRef] [PubMed]

17. Comber, H.; De Camargo Cancela, M.; Haase, T.; Johnson, H.; Sharp, L.; Pratschke, J. Affluence and Private Health Insurance Influence Treatment and Survival in Non-Hodgkin's Lymphoma. *PLoS ONE* **2016**, *11*, e0168684. [CrossRef] [PubMed]

18. Woods, L.M.; Rachet, B.; Coleman, M.P. Origins of Socio-Economic Inequalities in Cancer Survival: A Review. *Ann. Oncol.* **2005**, *17*, 5–19. [CrossRef] [PubMed]

19. Quaglia, A.; Vercelli, M.; Lillini, R.; Mugno, E.; Coebergh, J.W.; Quinn, M.; Martinez-Garcia, C.; Capocaccia, R.; Micheli, A. Socio-Economic Factors and Health Care System Characteristics Related to Cancer Survival in the Elderly: A Population-Based Analysis in 16 European Countries (ELDCARE Project). *Crit. Rev. Oncol. Hematol.* **2005**, *54*, 117–128. [CrossRef] [PubMed]

20. Afshar, N.; English, D.R.; Milne, R.L. Rural–Urban Residence and Cancer Survival in High-Income Countries: A Systematic Review. *Cancer* **2019**, *125*, 2172–2184. [CrossRef] [PubMed]

21. Belot, A.; Pohar-Perme, M. Social Disparities in Cancer Survival: Methodological Considerations. In *Social Environment and Cancer in Europe: Towards an Evidence-Based Public Health Policy*; Launoy, G., Zadnik, V., Coleman, M.P., Eds.; Springer International Publishing: Cham, Switzerland, 2021; pp. 39–54, ISBN 978-3-030-69329-9.

22. International Agency for Research on Cancer International Classification of Diseases for Oncology. Available online: http://codes.iarc.fr/ (accessed on 4 October 2019).

23. Fritz, A.; Percy, C.; Jack, A.; Shanmugaratnam, K.; Sobin, L.H.; Parkin, D.M.; Whelan, S.L. *International Classification of Diseases for Oncology*, 3rd ed.; World Health Organisation: Geneva, Switzerland, 2000.

24. Campo, E.; Swerdlow, S.H.; Harris, N.L.; Pileri, S.; Stein, H.; Jaffe, E.S. The 2008 WHO Classification of Lymphoid Neoplasms and beyond: Evolving Concepts and Practical Applications. *Blood* **2011**, *117*, 5019–5032. [CrossRef] [PubMed]

25. gov.uk. National Cancer Registry and Analysis Service. Available online: https://www.gov.uk/guidance/national-cancer-registration-and-analysis-service-ncras (accessed on 4 October 2019).

26. NHS Digital Hospital Episode Statistics. Available online: https://digital.nhs.uk/data-and-information/data-tools-and-services/data-services/hospital-episode-statistics (accessed on 4 October 2019).

27. Rachet, B.; Maringe, C.; Woods, L.M.; Ellis, L.; Spika, D.; Allemani, C. Multivariable Flexible Modelling for Estimating Complete Smoothed Life Tables for Sub-National Populations. *BMC Public Health* **2015**, *15*, 1240. [CrossRef] [PubMed]

28. Porta, M. *A Dictionary of Epidemiology*; Oxford University Press: Oxford, UK, 2014; ISBN 9780195314496.

29. Feinstein, A.R. The Pre-Therapeutic Classification of Co-Morbidity in Chronic Disease. *J. Chronic. Dis.* **1970**, *23*, 455–468. [CrossRef]

30. Maringe, C.; Fowler, H.; Rachet, B.; Luque-Fernandez, M.A. Reproducibility, Reliability and Validity of Population-Based Administrative Health Data for the Assessment of Cancer Non-Related Comorbidities. *PLoS ONE* **2017**, *12*, e0172814. [CrossRef] [PubMed]

31. Charlson, M.E.; Pompei, P.; Ales, K.L.; MacKenzie, C.R. A New Method of Classifying Prognostic Comorbidity in Longitudinal Studies: Development and Validation. *J. Chronic. Dis.* **1987**, *40*, 373–383. [CrossRef]

32. Armitage, J.N.; van der Meulen, J.H. Identifying Co-Morbidity in Surgical Patients Using Administrative Data with the Royal College of Surgeons Charlson Score. *Br. J. Surg.* **2010**, *97*, 772–781. [CrossRef]

33. gov.uk. Indices of Multiple Deprivation. Available online: https://www.gov.uk/government/statistics/english-indices-of-deprivation-2015 (accessed on 4 October 2019).

34. National Health Service: Data Dictionary Lower Super Output Area. Available online: https://www.datadictionary.nhs.uk/data_dictionary (accessed on 4 October 2019).

35. Elliss-Brookes, L.; McPhail, S.; Ives, A.; Greenslade, M.; Shelton, J.; Hiom, S.; Richards, M. Routes to Diagnosis for Cancer—Determining the Patient Journey Using Multiple Routine Data Sets. *Br. J. Cancer* **2012**, *107*, 1220–1226. [CrossRef] [PubMed]

36. Charvat, H.; Remontet, L.; Bossard, N.; Roche, L.; Dejardin, O.; Rachet, B.; Launoy, G.; Belot, A. A Multilevel Excess Hazard Model to Estimate Net Survival on Hierarchical Data Allowing for Non-Linear and Non-Proportional Effects of Covariates. *Stat. Med.* **2016**, *35*, 3066–3084. [CrossRef]

37. Verbeke, G.; Molenberghs, G. *Linear Mixed Models for Longitudinal Data*, 1st ed.; Springer: New York, NY, USA, 2000.

38. Belot, A.; Ndiaye, A.; Luque-Fernandez, M.-A.; Kipourou, D.-K.; Maringe, C.; Rubio, F.J.; Rachet, B. Summarizing and Communicating on Survival Data According to the Audience: A Tutorial on Different Measures Illustrated with Population-Based Cancer Registry Data. *Clin. Epidemiol.* **2019**, *11*, 53–65. [CrossRef] [PubMed]

39. Molenberghs, G.; Verbeke, G. *Models for Discrete Longitudinal Data*, 1st ed.; Springer: New York, NY, USA, 2005.

40. Agresti, A. *Categorical Data Analysis*, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2002.

41. Carpenter, J.R.; Kenward, M.G. *Multiple Imputation and Its Application*, 1st ed.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2013.

42. Little, R.J.; Rubin, D.B. *Statistical Analysis with Missing Data*; John Wiley & Sons, Inc.: New York, NY, USA, 1987.

43. Rubin, D.B. *Multiple Imputation for Nonresponse in Surveys*; Wiley: New York, NY, USA, 1987; ISBN 978-0-471-65574-9.
44. Quartagno, M.; Carpenter, J.R. Jomo: A Package for Multilevel Joint Modeling Multiple Imputation. *Stat. Med.* **2016**. [CrossRef]

45. Coiffier, B.; Lepage, E.; Brière, J.; Herbrecht, R.; Tilly, H.; Bouabdallah, R.; Morel, P.; Van Den Neste, E.; Salles, G.; Gaulard, P.; et al. CHOP Chemotherapy plus Rituximab Compared with CHOP Alone in Elderly Patients with Diffuse Large-B-Cell Lymphoma. *N. Engl. J. Med.* **2002**, *346*, 235-242. [CrossRef] [PubMed]

46. Coiffier, B. Rituximab in Combination with CHOP Improves Survival in Elderly Patients with Aggressive Non-Hodgkin's Lymphoma. *Semin. Oncol.* **2002**, *29*, 18-22. [CrossRef] [PubMed]

47. Delarue, R.; Tilly, H.; Mounier, N.; Petrella, T.; Salles, G.; Thieblemont, C.; Bologna, S.; Ghesquieres, H.; Hacini, M.; Fruchart, C.; et al. Dose-Dense Rituximab-CHOP Compared with Standard Rituximab-CHOP in Elderly Patients with Diffuse Large B-Cell Lymphoma (the LNH03-6B Study): A Randomised Phase 3 Trial. *Lancet Oncol.* **2013**, *14*, 525-533. [CrossRef]

48. McGowan, J.V.; Chung, R.; Maulik, A.; Piotrowska, I.; Walker, J.M.; Yellon, D.M. Anthracycline Chemotherapy and Cardiotoxicity. *Cardiovasc. Drugs Ther.* **2017**, *31*, 63-75. [CrossRef] [PubMed]

49. National Institute for Health and Care Excellence. *Non-Hodgkin's Lymphoma: Diagnosis and Management*; National Institute for Health and Care Excellence: London, UK, 2016.

50. Tilly, H.; Gomes da Silva, M.; Vitolo, U.; Jack, A.; Meignan, M.; Lopez-Guillermo, A.; Walewski, J.; André, M.; Johnson, P.W.; Pfreundschuh, M.; et al. Diffuse Large B-Cell Lymphoma (DLBCL): ESMO Clinical Practice Guidelines for Diagnosis, Treatment and Follow-Up. *Ann. Oncol.* **2015**, *26*, v116-v125. [CrossRef] [PubMed]

51. Bröckelmann, P.J.; McMullen, S.; Wilson, J.B.; Mueller, K.; Goring, S.; Stamatoullas, A.; Zagadailov, E.; Gautam, A.; Huebner, D.; Dalal, M.; et al. Patient and Physician Preferences for First-Line Treatment of Classical Hodgkin Lymphoma in Germany, France and the United Kingdom. *Br. J. Haematol.* **2019**, *184*, 202-214. [CrossRef]

52. Solimando, A.G.; Ribatti, D.; Vacca, A.; Einsele, H. Targeting B-Cell Non Hodgkin Lymphoma: New and Old Tricks. *Leuk. Res.* **2016**, *42*, 93-104. [CrossRef] [PubMed]

53. Merali, Z.; Wilson, J.R. Explanatory Versus Pragmatic Trials: An Essential Concept in Study Design and Interpretation. *Clin. Spine Surg.* **2017**, *30*, 404-406. [CrossRef]

54. Ghesquières, H.; Rossi, C.; Cherblanc, F.; Le Guyader-Peyrou, S.; Bijou, F.; Sujobert, P.; Fabbro-Peray, P.; Bernier, A.; Belot, A.; Chartier, L.; et al. A French Multicentric Prospective Prognostic Cohort with Epidemiological, Clinical, Biological and Treatment Information to Improve Knowledge on Lymphoma Patients: Study Protocol of the "REal World DAta in LYmphoma and Survival in Adults" (REALYSA) Cohort. *BMC Public Health* **2021**, *21*, 432. [CrossRef]

55. El-Galaly, T.C.; Cheah, C.Y.; Villa, D. Real World Data as a Key Element in Precision Medicine for Lymphoid Malignancies: Potentials and Pitfalls. *Br. J. Haematol.* **2019**, *186*, 409-419. [CrossRef] [PubMed]

56. Lossos, I.S.; Gascoyne, R.D. Transformation of Follicular Lymphoma. *Best Pract. Res. Clin. Haematol.* **2011**, *24*, 147-163. [CrossRef] [PubMed]

57. Putter, H.; van Houwelingen, H.C. Understanding Landmarking and Its Relation with Time-Dependent Cox Regression. *Stat. Biosci.* **2017**, *9*, 489-503. [CrossRef]

58. Sogaard, M.; Thomsen, R.W.; Bossen, K.S.; Sorensen, H.T.; Norgaard, M. The Impact of Comorbidity on Cancer Survival: A Review. *Clin. Epidemiol.* **2013**, *5*, 3-29. [CrossRef]

59. Frederiksen, B.L.; Dalton, S.O.; Osler, M.; Steding-Jessen, M.; de Nully Brown, P. Socioeconomic Position, Treatment, and Survival of Non-Hodgkin Lymphoma in Denmark-a Nationwide Study. *Br. J. Cancer* **2012**, *106*, 988-995. [CrossRef]

60. Rachet, B.; Mitry, E.; Shah, A.; Cooper, N.; Coleman, M.P. Survival from Non-Hodgkin Lymphoma in England and Wales up to 2001. *Br. J. Cancer* **2008**, *99*, S104-S106. [CrossRef] [PubMed]

61. Bray, C.; Morrison, D.S.; McKay, P. Socio-Economic Deprivation and Survival of Non-Hodgkin Lymphoma in Scotland. *Leuk. Lymphoma* **2008**, *49*, 917-923. [CrossRef] [PubMed]

62. Pohar Perme, M.; Estève, J.; Rachet, B. Analysing Population-Based Cancer Survival—Settling the Controversies. *BMC Cancer* **2016**, *16*, 933. [CrossRef]

63. Kobayashi, Y.; Miura, K.; Hojo, A.; Hatta, Y.; Tanaka, T.; Kurita, D.; Iriyama, N.; Kobayashi, S.; Takeuchi, J. Charlson Comorbidity Index Is an Independent Prognostic Factor among Elderly Patients with Diffuse Large B-Cell Lymphoma. *J. Cancer Res. Clin. Oncol.* **2011**, *137*, 1079-1084. [CrossRef]

64. Saygin, C.; Jia, X.; Hill, B.; Dean, R.; Pohlman, B.; Smith, M.R.; Jagadeesh, D. Impact of Comorbidities on Outcomes of Elderly Patients with Diffuse Large B-cell Lymphoma. *Am. J. Hematol.* **2017**, *92*, 989-996. [CrossRef] [PubMed]

65. Chihara, D.; Westin, J.R.; Oki, Y.; Ahmed, M.A.; Do, B.; Fayad, L.E.; Hagemeister, F.B.; Romaguera, J.E.; Fanale, M.A.; Lee, H.J.; et al. Management Strategies and Outcomes for Very Elderly Patients with Diffuse Large B-Cell Lymphoma. *Cancer* **2016**, *122*, 3145-3151. [CrossRef]

66. Janssen-Heijnen, M.L.; van Spronsen, D.J.; Lemmens, V.E.; Houterman, S.; Verheij, K.D.; Coebergh, J.W. A Population-Based Study of Severity of Comorbidity among Patients with Non-Hodgkin's Lymphoma: Prognostic Impact Independent of International Prognostic Index. *Br. J. Haematol.* **2005**, *129*, 597-606. [CrossRef] [PubMed]

67. Crooks, C.J.; West, J.; Card, T.R. A Comparison of the Recording of Comorbidity in Primary and Secondary Care by Using the Charlson Index to Predict Short-Term and Long-Term Survival in a Routine Linked Data Cohort. *BMJ Open* **2015**, *5*, e007974. [CrossRef] [PubMed]

68. Carpenter, J.; Goldstein, H.; Kenward, M. REALCOM-IMPUTE Software for Multilevel Multiple Imputation with Mixed Response Types. *J. Stat. Softw.* **2011**, *45*, 1-20. [CrossRef]
69. Quartagno, M.; Carpenter, J.R. Multiple Imputation for Discrete Data: Evaluation of the Joint Latent Normal Model. *Biom. J.* **Biom. Z. 2019**, *61*, 1003–1019. [CrossRef] [PubMed]

70. Ingleby, F.C.; Belot, A.; Atherton, I.; Baker, M.; Elliss-Brookes, L.; Woods, L.M. Assessment of the Concordance between Individual-Level and Area-Level Measures of Socio-Economic Deprivation in a Cancer Patient Cohort in England and Wales. *BMJ Open* **2020**, *10*, e041714. [CrossRef]

71. Belot, A.; Remontet, L.; Rachet, B.; Dejardin, O.; Charvat, H.; Bara, S.; Guizard, A.-V.; Roche, L.; Launoy, G.; Bossard, N. Describing the Association between Socioeconomic Inequalities and Cancer Survival: Methodological Guidelines and Illustration with Population-Based Data. *Clin. Epidemiol.* **2018**, *10*, 561–573. [CrossRef] [PubMed]

72. Madley-Dowd, P.; Hughes, R.; Tilling, K.; Heron, J. The Proportion of Missing Data Should Not Be Used to Guide Decisions on Multiple Imputation. *J. Clin. Epidemiol.* **2019**, *110*, 63–73. [CrossRef] [PubMed]

73. Carnerio, I.; Fry, A.; Ironmonger, L.; Connor, K.; Ormiston-Smith, N.; Henson, K.; Johnson, S.; McPhail, S.; Elliss-Brookes, L. Variation in the Routes to Cancer Diagnosis and Stage for Ten Cancer Sites. Available online: https://www.researchgate.net/publication/311602679_Variation_in_the_Routes_to_Cancer_Diagnosis_and_Stage_for_Ten_Cancer_Sites (accessed on 15 August 2021).

74. Allemani, C.; Matsuda, T.; Di Carlo, V.; Harewood, R.; Matz, M.; Nikšić, M.; Bonaventure, A.; Valkov, M.; Johnson, C.J.; Estève, J.; et al. Global Surveillance of Trends in Cancer Survival 2000–14 (CONCORD-3): Analysis of Individual Records for 37,513,025 Patients Diagnosed with One of 18 Cancers from 322 Population-Based Registries in 71 Countries. *Lancet* **2018**, *391*, 1023–1075. [CrossRef]

75. Thomson, C.S.; Forman, D. Cancer Survival in England and the Influence of Early Diagnosis: What Can We Learn from Recent EUROCARE Results? *Br. J. Cancer* **2009**, *101*, S102–S109. [CrossRef] [PubMed]