What can hospital emergency admissions prior to cancer diagnosis tell us about socio-economic inequalities in cancer diagnosis? Evidence from population-based data in England.

BACKGROUND: More deprived cancer patients are at higher risk of Emergency Presentation (EP) with most studies pointing to lower symptom awareness and increased comorbidities to explain those patterns. With the example of colon cancer, we examine patterns of hospital emergency admissions (HEAs) history in the most and least deprived patients as a potential precursor of EP. METHODS: We analysed the rates of hospital admissions and their admission codes (retrieved from Hospital Episode Statistics) in the two years preceding cancer diagnosis by sex, deprivation and route to diagnosis (EP, non-EP). To select the conditions (grouped admission codes) that best predict emergency admission, we adapted the purposeful variable selection to mixed-effects logistic regression. RESULTS: Colon cancer patients diagnosed through EP had the highest number of HEAs than all the other routes to diagnosis, especially in the last 7 months before diagnosis. Most deprived patients had an overall higher rate and higher probability of HEA but fewer conditions associated with it. CONCLUSIONS: Our findings point to higher use of emergency services for non-specific symptoms and conditions in the most deprived patients, preceding colon cancer diagnosis. Health system barriers may be a shared factor of socio-economic inequalities in EP and HEAs.

Immortal-time bias in older vs younger age groups: a simulation study with application to a population-based cohort of patients with colon cancer.

BACKGROUND: In observational studies, the risk of immortal-time bias (ITB) increases with the likelihood of early death, itself increasing with age. We investigated how age impacts the magnitude of ITB when estimating the effect of surgery on 1-year overall survival (OS) in patients with Stage IV colon cancer aged 50-74 and 75-84 in England. METHODS: Using simulations, we compared estimates from a time-fixed exposure model to three statistical methods addressing ITB: time-varying exposure, delayed entry and landmark methods. We then estimated the effect of surgery on OS using a population-based cohort of patients from the CORECT-R resource and conducted the analysis using the emulated target trial framework. RESULTS: In simulations, the magnitude of ITB was larger among older patients when their probability of early death increased or treatment was delayed. The bias was corrected using the methods addressing ITB. When applied to CORECT-R data, these methods yielded a smaller effect of surgery than the time-fixed exposure approach but effects were similar in both age groups. CONCLUSION: ITB must be addressed in all longitudinal studies, particularly, when investigating the effect of exposure on an outcome in different groups of people (e.g., age groups) with different distributions of exposure and outcomes.

Target trial emulation with multi-state model analysis to assess treatment effectiveness using clinical COVID-19 data.

BACKGROUND: Real-world observational data are an important source of evidence on the treatment effectiveness for patients hospitalized with coronavirus disease 2019 (COVID-19). However, observational studies evaluating treatment effectiveness based on longitudinal data are often prone to methodological biases such as immortal time bias, confounding bias, and competing risks. METHODS: For exemplary target trial emulation, we used a cohort of patients hospitalized with COVID-19 (n = 501) in a single centre. We described the methodology for evaluating the effectiveness of a single-dose treatment, emulated a trial using real-world data, and drafted a hypothetical study protocol describing the main components. To avoid immortal time and time-fixed confounding biases, we applied the clone-censor-weight technique. We set a 5-day grace period as a period of time when treatment could be initiated. We used the inverse probability of censoring weights to account for the selection bias introduced by artificial censoring. To estimate the treatment effects, we took the multi-state model approach. We considered a multi-state model with five states. The primary endpoint was defined as clinical severity status, assessed by a 5-point ordinal scale on day 30. Differences between the treatment group and standard of care treatment group were calculated using a proportional odds model and shown as odds ratios. Additionally, the weighted cause-specific hazards and transition probabilities for each treatment arm were presented. RESULTS: Our study demonstrates that trial emulation with a multi-state model analysis is a suitable approach to address observational data limitations, evaluate treatment effects on clinically heterogeneous in-hospital death and discharge alive endpoints, and consider the intermediate state of admission to ICU. The multi-state model analysis allows us to summarize results using stacked probability plots that make it easier to interpret results. CONCLUSIONS: Extending the emulated target trial approach to multi-state model analysis complements treatment effectiveness analysis by gaining information on competing events. Combining two methodologies offers an option to address immortal time bias, confounding bias, and competing risk events. This methodological approach can provide additional insight for decision-making, particularly when data from randomized controlled trials (RCTs) are unavailable.

Extracorporeal Membrane Oxygenation for Severe Acute Respiratory Distress Syndrome Associated with COVID-19: An Emulated Target Trial Analysis.

Rationale: Whether patients with coronavirus disease (COVID-19) may benefit from extracorporeal membrane oxygenation (ECMO) compared with conventional invasive mechanical ventilation (IMV) remains unknown. Objectives: To estimate the effect of ECMO on 90-day mortality versus IMV only. Methods: Among 4,244 critically ill adult patients with COVID-19 included in a multicenter cohort study, we emulated a target trial comparing the treatment strategies of initiating ECMO versus no ECMO within 7 days of IMV in patients with severe acute respiratory distress syndrome (PaO2/FiO2 < 80 or PaCO2 ⩾ 60 mm Hg). We controlled for confounding using a multivariable Cox model on the basis of predefined variables. Measurements and Main Results: A total of 1,235 patients met the full eligibility criteria for the emulated trial, among whom 164 patients initiated ECMO. The ECMO strategy had a higher survival probability on Day 7 from the onset of eligibility criteria (87% vs. 83%; risk difference, 4%; 95% confidence interval, 0-9%), which decreased during follow-up (survival on Day 90: 63% vs. 65%; risk difference, -2%; 95% confidence interval, -10 to 5%). However, ECMO was associated with higher survival when performed in high-volume ECMO centers or in regions where a specific ECMO network organization was set up to handle high demand and when initiated within the first 4 days of IMV and in patients who are profoundly hypoxemic. Conclusions: In an emulated trial on the basis of a nationwide COVID-19 cohort, we found differential survival over time of an ECMO compared with a no-ECMO strategy. However, ECMO was consistently associated with better outcomes when performed in high-volume centers and regions with ECMO capacities specifically organized to handle high demand.

On models for the estimation of the excess mortality hazard in case of insufficiently stratified life tables.

In cancer epidemiology using population-based data, regression models for the excess mortality hazard is a useful method to estimate cancer survival and to describe the association between prognosis factors and excess mortality. This method requires expected mortality rates from general population life tables: each cancer patient is assigned an expected (background) mortality rate obtained from the life tables, typically at least according to their age and sex, from the population they belong to. However, those life tables may be insufficiently stratified, as some characteristics such as deprivation, ethnicity, and comorbidities, are not available in the life tables for a number of countries. This may affect the background mortality rate allocated to each patient, and it has been shown that not including relevant information for assigning an expected mortality rate to each patient induces a bias in the estimation of the regression parameters of the excess hazard model. We propose two parametric corrections in excess hazard regression models, including a single-parameter or a random effect (frailty), to account for possible mismatches in the life table and thus misspecification of the background mortality rate. In an extensive simulation study, the good statistical performance of the proposed approach is demonstrated, and we illustrate their use on real population-based data of lung cancer patients. We present conditions and limitations of these methods and provide some recommendations for their use in practice.

Introduction to computational causal inference using reproducible Stata, R, and Python code: A tutorial.

The main purpose of many medical studies is to estimate the effects of a treatment or exposure on an outcome. However, it is not always possible to randomize the study participants to a particular treatment, therefore observational study designs may be used. There are major challenges with observational studies; one of which is confounding. Controlling for confounding is commonly performed by direct adjustment of measured confounders; although, sometimes this approach is suboptimal due to modeling assumptions and misspecification. Recent advances in the field of causal inference have dealt with confounding by building on classical standardization methods. However, these recent advances have progressed quickly with a relative paucity of computational-oriented applied tutorials contributing to some confusion in the use of these methods among applied researchers. In this tutorial, we show the computational implementation of different causal inference estimators from a historical perspective where new estimators were developed to overcome the limitations of the previous estimators (ie, nonparametric and parametric g-formula, inverse probability weighting, double-robust, and data-adaptive estimators). We illustrate the implementation of different methods using an empirical example from the Connors study based on intensive care medicine, and most importantly, we provide reproducible and commented code in Stata, R, and Python for researchers to adapt in their own observational study. The code can be accessed at https://github.com/migariane/Tutorial_Computational_Causal_Inference_Estimators.

The impact of the COVID-19 pandemic on cancer deaths due to delays in diagnosis in England, UK: a national, population-based, modelling study.

BACKGROUND: Since a national lockdown was introduced across the UK in March, 2020, in response to the COVID-19 pandemic, cancer screening has been suspended, routine diagnostic work deferred, and only urgent symptomatic cases prioritised for diagnostic intervention. In this study, we estimated the impact of delays in diagnosis on cancer survival outcomes in four major tumour types. METHODS: In this national population-based modelling study, we used linked English National Health Service (NHS) cancer registration and hospital administrative datasets for patients aged 15-84 years, diagnosed with breast, colorectal, and oesophageal cancer between Jan 1, 2010, and Dec 31, 2010, with follow-up data until Dec 31, 2014, and diagnosed with lung cancer between Jan 1, 2012, and Dec 31, 2012, with follow-up data until Dec 31, 2015. We use a routes-to-diagnosis framework to estimate the impact of diagnostic delays over a 12-month period from the commencement of physical distancing measures, on March 16, 2020, up to 1, 3, and 5 years after diagnosis. To model the subsequent impact of diagnostic delays on survival, we reallocated patients who were on screening and routine referral pathways to urgent and emergency pathways that are associated with more advanced stage of disease at diagnosis. We considered three reallocation scenarios representing the best to worst case scenarios and reflect actual changes in the diagnostic pathway being seen in the NHS, as of March 16, 2020, and estimated the impact on net survival at 1, 3, and 5 years after diagnosis to calculate the additional deaths that can be attributed to cancer, and the total years of life lost (YLLs) compared with pre-pandemic data. FINDINGS: We collected data for 32 583 patients with breast cancer, 24 975 with colorectal cancer, 6744 with oesophageal cancer, and 29 305 with lung cancer. Across the three different scenarios, compared with pre-pandemic figures, we estimate a 7·9-9·6% increase in the number of deaths due to breast cancer up to year 5 after diagnosis, corresponding to between 281 (95% CI 266-295) and 344 (329-358) additional deaths. For colorectal cancer, we estimate 1445 (1392-1591) to 1563 (1534-1592) additional deaths, a 15·3-16·6% increase; for lung cancer, 1235 (1220-1254) to 1372 (1343-1401) additional deaths, a 4·8-5·3% increase; and for oesophageal cancer, 330 (324-335) to 342 (336-348) additional deaths, 5·8-6·0% increase up to 5 years after diagnosis. For these four tumour types, these data correspond with 3291-3621 additional deaths across the scenarios within 5 years. The total additional YLLs across these cancers is estimated to be 59 204-63 229 years. INTERPRETATION: Substantial increases in the number of avoidable cancer deaths in England are to be expected as a result of diagnostic delays due to the COVID-19 pandemic in the UK. Urgent policy interventions are necessary, particularly the need to manage the backlog within routine diagnostic services to mitigate the expected impact of the COVID-19 pandemic on patients with cancer. FUNDING: UK Research and Innovation Economic and Social Research Council.

Comorbidity prevalence among cancer patients: a population-based cohort study of four cancers.

BACKGROUND: The presence of comorbidity affects the care of cancer patients, many of whom are living with multiple comorbidities. The prevalence of cancer comorbidity, beyond summary metrics, is not well known. This study aims to estimate the prevalence of comorbid conditions among cancer patients in England, and describe the association between cancer comorbidity and socio-economic position, using population-based electronic health records. METHODS: We linked England cancer registry records of patients diagnosed with cancer of the colon, rectum, lung or Hodgkin lymphoma between 2009 and 2013, with hospital admissions records. A comorbidity was any one of fourteen specific conditions, diagnosed during hospital admission up to 6 years prior to cancer diagnosis. We calculated the crude and age-sex adjusted prevalence of each condition, the frequency of multiple comorbidity combinations, and used logistic regression and multinomial logistic regression to estimate the adjusted odds of having each condition and the probability of having each condition as a single or one of multiple comorbidities, respectively, by cancer type. RESULTS: Comorbidity was most prevalent in patients with lung cancer and least prevalent in Hodgkin lymphoma patients. Up to two-thirds of patients within each of the four cancer patient cohorts we studied had at least one comorbidity, and around half of the comorbid patients had multiple comorbidities. Our study highlighted common comorbid conditions among the cancer patient cohorts. In all four cohorts, the odds of having a comorbidity and the probability of multiple comorbidity were consistently highest in the most deprived cancer patients. CONCLUSIONS: Cancer healthcare guidelines may need to consider prominent comorbid conditions, particularly to benefit the prognosis of the most deprived patients who carry the greater burden of comorbidity. Insight into patterns of cancer comorbidity may inform further research into the influence of specific comorbidities on socio-economic inequalities in receipt of cancer treatment and in short-term mortality.

Surgical treatment and survival from colorectal cancer in Denmark, England, Norway, and Sweden: a population-based study.

BACKGROUND: Survival from colorectal cancer has been shown to be lower in Denmark and England than in comparable high-income countries. We used data from national colorectal cancer registries to assess whether differences in the proportion of patients receiving resectional surgery could contribute to international differences in colorectal cancer survival. METHODS: In this population-based study, we collected data from all patients aged 18-99 years diagnosed with primary, invasive, colorectal adenocarcinoma from Jan 1, 2010, to Dec 31, 2012, in Denmark, England, Norway, and Sweden, from national colorectal cancer registries. We estimated age-standardised net survival using multivariable modelling, and we compared the proportion of patients receiving resectional surgery by stage and age. We used logistic regression to predict the resectional surgery status patients would have had if they had been treated as in the best performing country, given their individual characteristics. FINDINGS: We extracted registry data for 139 457 adult patients with invasive colorectal adenocarcinoma: 12 958 patients in Denmark, 97 466 in England, 11 450 in Norway, and 17 583 in Sweden. 3-year colon cancer survival was lower in England (63·9%, 95% CI 63·5-64·3) and Denmark (65·7%, 64·7-66·8) than in Norway (69·5%, 68·4-70·5) and Sweden (72·1%, 71·2-73·0). Rectal cancer survival was lower in England (69·7%, 69·1-70·3) than in the other three countries (Denmark 72·5%, 71·1-74·0; Sweden 74·1%, 72·7-75·4; and Norway 75·0%, 73·1-76·8). We found no significant differences in survival for patients with stage I disease in any of the four countries. 3-year survival after stage II or III rectal cancer and stage IV colon cancer was consistently lower in England (stage II rectal cancer 86·4%, 95% CI 85·0-87·6; stage III rectal cancer 75·5%, 74·2-76·7; and stage IV colon cancer 20·5%, 19·9-21·1) than in Norway (94·1%, 91·5-96·0; 83·4%, 80·1-86·1; and 33·0%, 31·0-35·1) and Sweden (92·9%, 90·8-94·6; 80·6%, 78·2-82·7; and 23·7%, 22·0-25·3). 3-year survival after stage II rectal cancer and stage IV colon cancer was also lower in England than in Denmark (stage II rectal cancer 91·2%, 88·8-93·1; and stage IV colon cancer 23·5%, 21·9-25·1). The total proportion of patients treated with resectional surgery ranged from 47 803 (68·4%) of 69 867 patients in England to 9582 (81·3%) of 11 786 in Sweden for colon cancer, and from 16 544 (59·9%) of 27 599 in England to 4106 (70·8%) of 5797 in Sweden for rectal cancer. This range was widest for patients older than 75 years (colon cancer 19 078 [59·7%] of 31 946 patients in England to 4429 [80·9%] of 5474 in Sweden; rectal cancer 4663 [45·7%] of 10 195 in England to 1342 [61·9%] of 2169 in Sweden), and the proportion of patients treated with resectional surgery was consistently lowest in England. The age gradient of the decline in the proportion of patients treated with resectional surgery was steeper in England than in the other three countries in all stage categories. In the hypothetical scenario where all patients were treated as in Sweden, given their age, sex, and disease stage, the largest increase in resectional surgery would be for patients with stage III rectal cancer in England (increasing from 70·3% to 88·2%). INTERPRETATION: Survival from colon cancer and rectal cancer in England and colon cancer in Denmark was lower than in Norway and Sweden. Survival paralleled the relative provision of resectional surgery in these countries. Differences in patient selection for surgery, especially in patients older than 75 years or individuals with advanced disease, might partly explain these differences in international colorectal cancer survival. FUNDING: Early Diagnosis Policy Research Grant from Cancer Research UK (C7923/A18348).

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.

INTRODUCTION: We investigated socioeconomic disparities and the role of the main prognostic factors in receiving major surgical treatment in patients with lung cancer in England. METHODS: Our study comprised 31 351 patients diagnosed with non-small cell lung cancer in England in 2012. Data from the national population-based cancer registry were linked to Hospital Episode Statistics and National Lung Cancer Audit data to obtain information on stage, performance status and comorbidities, and to identify patients receiving major surgical treatment. To describe the association between prognostic factors and surgery, we performed two different analyses: one using multivariable logistic regression and one estimating cause-specific hazards for death and surgery. In both analyses, we used multiple imputation to deal with missing data. RESULTS: We showed strong evidence that the comorbidities 'congestive heart failure', 'cerebrovascular disease' and 'chronic obstructive pulmonary disease' reduced the receipt of surgery in early stage patients. We also observed gender differences and substantial age differences in the receipt of surgery. Despite accounting for sex, age at diagnosis, comorbidities, stage at diagnosis, performance status and indication of having had a PET-CT scan, the socioeconomic differences persisted in both analyses: more deprived people had lower odds and lower rates of receiving surgery in early stage lung cancer. DISCUSSION: Comorbidities play an important role in whether patients undergo surgery, but do not completely explain the socioeconomic difference observed in early stage patients. Future work investigating access to and distance from specialist hospitals, as well as patient perceptions and patient choice in receiving surgery, could help disentangle these persistent socioeconomic inequalities.

Comparison of model-building strategies for excess hazard regression models in the context of cancer epidemiology.

BACKGROUND: Large and complex population-based cancer data are becoming broadly available, thanks to purposeful linkage between cancer registry data and health electronic records. Aiming at understanding the explanatory power of factors on cancer survival, the modelling and selection of variables need to be understood and exploited properly for improving model-based estimates of cancer survival. METHOD: We assess the performances of well-known model selection strategies developed by Royston and Sauerbrei and Wynant and Abrahamowicz that we adapt to the relative survival data setting and to test for interaction terms. RESULTS: We apply these to all male patients diagnosed with lung cancer in England in 2012 (N = 15,688), and followed-up until 31/12/2015. We model the effects of age at diagnosis, tumour stage, deprivation, comorbidity and emergency presentation, as well as interactions between age and all of the above. Given the size of the dataset, all model selection strategies favoured virtually the same model, except for a non-linear effect of age at diagnosis selected by the backward-based selection strategies (versus a linear effect selected otherwise). CONCLUSION: The results from extensive simulations evaluating varying model complexity and sample sizes provide guidelines on a model selection strategy in the context of excess hazard modelling.

Deprivation-specific life tables using multivariable flexible modelling - trends from 2000-2002 to 2010-2012, Portugal.

BACKGROUND: Completing mortality data by information on possible socioeconomic inequalities in mortality is crucial for policy planning. The aim of this study was to build deprivation-specific life tables using the Portuguese version of the European Deprivation Index (EDI) as a measure of area-level socioeconomic deprivation, and to evaluate mortality trends between the periods 2000-2002 and 2010-2012. METHODS: Statistics Portugal provided the counts of deaths and population by sex, age group, calendar year and area of residence (parish). A socioeconomic deprivation level was assigned to each parish according to the quintile of their national EDI distribution. Death counts were modelled within the generalised linear model framework as a function of age, deprivation level and calendar period. Mortality Rate Ratios (MRR) were estimated to evaluate variations in mortality between deprivation groups and periods. RESULTS: Life expectancy at birth increased from 74.0 and 80.9 years in 2000-2002, for men and women, respectively, and to 77.6 and 83.8 years in 2010-2012. Yet, life expectancy at birth differed by deprivation, with, compared to least deprived population, a deficit of about 2 (men) and 1 (women) years among most deprived in the whole study period. The higher mortality experienced by most deprived groups at birth (in 2010-2012, mortality rate ratios of 1.74 and 1.29 in men and women, respectively) progressively disappeared with increasing age. CONCLUSIONS: Persistent differences in mortality and life expectancy were observed according to ecological socioeconomic deprivation. These differences were larger among men and mostly marked at birth for both sexes.

Contrasting effects of comorbidities on emergency colon cancer diagnosis: a longitudinal data-linkage study in England.

BACKGROUND: One in three colon cancers are diagnosed as an emergency, which is associated with worse cancer outcomes. Chronic conditions (comorbidities) affect large proportions of adults and they might influence the risk of emergency presentations (EP). METHODS: We aimed to evaluate the effect of specific pre-existing comorbidities on the risk of colon cancer being diagnosed following an EP rather than through non-emergency routes. The cohort study included 5745 colon cancer patients diagnosed in England 2005-2010, with individually-linked cancer registry, primary and secondary care data. In addition to multivariable analyses we also used potential-outcomes methods. RESULTS: Colon cancer patients with comorbidities consulted their GP more frequently with cancer symptoms during the pre-diagnostic year, compared with non-comorbid cancer patients. EP occurred more frequently in patients with 'serious' or complex comorbidities (diabetes, cardiac and respiratory diseases) diagnosed/treated in hospital during the years pre-cancer diagnosis (43% EP in comorbid versus 27% in non-comorbid individuals; multivariable analysis Odds Ratio (OR), controlling for socio-demographic factors and symptoms: men OR = 2.40; 95% CI 2.0-2.9 and women OR = 1.98; 95% CI 1.6-2.4. Among women younger than 60, gynaecological (OR = 3.41; 95% CI 1.2-9.9) or recent onset gastro-intestinal conditions (OR = 2.84; 95% CI 1.1-7.7) increased the risk of EP. In contrast, primary care visits for hypertension monitoring decreased EPs for both genders. CONCLUSIONS: Patients with comorbidities have a greater risk of being diagnosed with cancer as an emergency, although they consult more frequently with cancer symptoms during the year pre-cancer diagnosis. This suggests that comorbidities may interfere with diagnostic reasoning or investigations due to 'competing demands' or because they provide 'alternative explanations'. In contrast, the management of chronic risk factors such as hypertension may offer opportunities for earlier diagnosis. Interventions are needed to support the diagnostic process in comorbid patients. Appropriate guidelines and diagnostic services to support the evaluation of new or changing symptoms in comorbid patients may be useful.

Data-Adaptive Estimation for Double-Robust Methods in Population-Based Cancer Epidemiology: Risk Differences for Lung Cancer Mortality by Emergency Presentation.

In this paper, we propose a structural framework for population-based cancer epidemiology and evaluate the performance of double-robust estimators for a binary exposure in cancer mortality. We conduct numerical analyses to study the bias and efficiency of these estimators. Furthermore, we compare 2 different model selection strategies based on 1) Akaike's Information Criterion and the Bayesian Information Criterion and 2) machine learning algorithms, and we illustrate double-robust estimators' performance in a real-world setting. In simulations with correctly specified models and near-positivity violations, all but the naive estimators had relatively good performance. However, the augmented inverse-probability-of-treatment weighting estimator showed the largest relative bias. Under dual model misspecification and near-positivity violations, all double-robust estimators were biased. Nevertheless, the targeted maximum likelihood estimator showed the best bias-variance trade-off, more precise estimates, and appropriate 95% confidence interval coverage, supporting the use of the data-adaptive model selection strategies based on machine learning algorithms. We applied these methods to estimate adjusted 1-year mortality risk differences in 183,426 lung cancer patients diagnosed after admittance to an emergency department versus persons with a nonemergency cancer diagnosis in England (2006-2013). The adjusted mortality risk (for patients diagnosed with lung cancer after admittance to an emergency department) was 16% higher in men and 18% higher in women, suggesting the importance of interventions targeting early detection of lung cancer signs and symptoms.

Persistent inequalities in unplanned hospitalisation among colon cancer patients across critical phases of their care pathway, England, 2011-13.

BACKGROUND: Reducing hospital emergency admissions is a key target for all modern health systems. METHODS: We analysed colon cancer patients diagnosed in 2011-13 in England. We screened their individual Hospital Episode Statistics records in the 90 days pre-diagnosis, the 90 days post-diagnosis, and the 90 days pre-death (in the year following diagnosis), for the occurrence of hospital emergency admissions (HEAs). RESULTS: Between a quarter and two thirds of patients experience HEA in the three 90-day periods examined: pre-diagnosis, post-diagnosis and before death. Patients with tumour stage I-III from more deprived backgrounds had higher proportions of HEAs than less deprived patients during all studied periods. This remains even after adjusting for differing distributions of risk factors such as age, sex, comorbidity and stage at diagnosis. CONCLUSIONS: Although in some cases HEAs might be unavoidable or even appropriate, the proportion of HEAs varies by socioeconomic status, even after controlling for the usual patient factors, suggestive of remediable causes of excess emergency healthcare utilisation in patients belonging to higher deprivation groups. Future inquiries should address the potential role of clinical complications, sub-optimal healthcare administration, premature discharge or a lack of social support as potential explanations for these patterns of inequality.

Trends in lung cancer emergency presentation in England, 2006-2013: is there a pattern by general practice?

BACKGROUND: Emergency presentations (EP) represent over a third of all lung cancer admissions in England. Such presentations usually reflect late stage disease and are associated with poor survival. General practitioners (GPs) act as gate-keepers to secondary care and so we sought to understand the association between GP practice characteristics and lung cancer EP. METHODS: Data on general practice characteristics were extracted for all practices in England from the Quality Outcomes Framework, the Health and Social Care Information Centre, the GP Patient Survey, the Cancer Commissioning Toolkit and the area deprivation score for each practice. After linking these data to lung cancer patient registrations in 2006-2013, we explored trends in three types of EP, patient-led, GP-led and 'other', by general practice characteristics and by socio-demographic characteristics of patients. RESULTS: Overall proportions of lung cancer EP decreased from 37.9% in 2006 to 34.3% in 2013. Proportions of GP-led EP nearly halved during this period, from 28.3 to 16.3%, whilst patient-led emergency presentations rose from 62.1 to 66.7%. When focusing on practice-specific levels of EP, 14% of general practices had higher than expected proportions of EP at least once in 2006-13, but there was no evidence of clustering of patients within practice, meaning that none of the practice characteristics examined explained differing proportions of EP by practice. CONCLUSION: We found that the high proportion of lung cancer EP is not the result of a few practices with very abnormal patterns of EP, but of a large number of practices susceptible to reaching high proportions of EP. This suggests a system-wide issue, rather than problems with specific practices. High proportions of lung cancer EP are mainly the result of patient-initiated attendances in A&E. Our results demonstrate that interventions to encourage patients not to bypass primary care must be system wide rather than targeted at specific practices.

Explained variation of excess hazard models.

The availability of longstanding collection of detailed cancer patient information makes multivariable modelling of cancer-specific hazard of death appealing. We propose to report variation in survival explained by each variable that constitutes these models. We adapted the ranks explained (RE) measure to the relative survival data setting, ie, when competing risks of death are accounted for through life tables from the general population. RE is calculated at each event time. We introduce weights for each death reflecting its probability to be a cancer death. RE varies between -1 and +1 and can be reported at given times in the follow-up and as a time-varying measure from diagnosis onward. We present an application for patients diagnosed with colon or lung cancer in England. The RE measure shows reasonable properties and is comparable in both relative and cause-specific settings. One year after diagnosis, RE for the most complex excess hazard models reaches 0.56, 95% CI: 0.54 to 0.58 (0.58 95% CI: 0.56-0.60) and 0.69, 95% CI: 0.68 to 0.70 (0.67, 95% CI: 0.66-0.69) for lung and colon cancer men (women), respectively. Stage at diagnosis accounts for 12.4% (10.8%) of the overall variation in survival among lung cancer patients whereas it carries 61.8% (53.5%) of the survival variation in colon cancer patients. Variables other than performance status for lung cancer (10%) contribute very little to the overall explained variation. The proportion of the variation in survival explained by key prognostic factors is a crucial information toward understanding the mechanisms underpinning cancer survival. The time-varying RE provides insights into patterns of influence for strong predictors.

Adjusting for overdispersion in piecewise exponential regression models to estimate excess mortality rate in population-based research.

BACKGROUND: In population-based cancer research, piecewise exponential regression models are used to derive adjusted estimates of excess mortality due to cancer using the Poisson generalized linear modelling framework. However, the assumption that the conditional mean and variance of the rate parameter given the set of covariates x i are equal is strong and may fail to account for overdispersion given the variability of the rate parameter (the variance exceeds the mean). Using an empirical example, we aimed to describe simple methods to test and correct for overdispersion. METHODS: We used a regression-based score test for overdispersion under the relative survival framework and proposed different approaches to correct for overdispersion including a quasi-likelihood, robust standard errors estimation, negative binomial regression and flexible piecewise modelling. RESULTS: All piecewise exponential regression models showed the presence of significant inherent overdispersion (p-value <0.001). However, the flexible piecewise exponential model showed the smallest overdispersion parameter (3.2 versus 21.3) for non-flexible piecewise exponential models. CONCLUSION: We showed that there were no major differences between methods. However, using a flexible piecewise regression modelling, with either a quasi-likelihood or robust standard errors, was the best approach as it deals with both, overdispersion due to model misspecification and true or inherent overdispersion.

Multivariable flexible modelling for estimating complete, smoothed life tables for sub-national populations.

BACKGROUND: The methods currently available to estimate age- and sex-specific mortality rates for sub-populations are subject to a number of important limitations. We propose two alternative multivariable approaches: a relational model and a Poisson model both using restricted cubic splines. METHODS: We evaluated a flexible Poisson and flexible relational model against the Elandt-Johnson approach in a simulation study using 100 random samples of population and death counts, with different sampling proportions and data arrangements. Estimated rates were compared to the original mortality rates using goodness-of-fit measures and life expectancy. We further investigated an approach for determining optimal knot locations in the Poisson model. RESULTS: The flexible Poisson model outperformed the flexible relational and Elandt-Johnson methods with the smallest sample of data (1%). With the largest sample of data (20%), the flexible Poisson and flexible relational models performed comparably, though the flexible Poisson model displayed a slight advantage. Both approaches tended to underestimate infant mortality and thereby overestimate life expectancy at birth. The flexible Poisson model performed much better at young ages when knots were fixed a priori. For ages 30 and above, results were similar to the model with no fixed knots. CONCLUSIONS: The flexible Poisson model is recommended because it derives robust and unbiased estimates for sub-populations without making strong assumptions about age-specific mortality profiles. Fixing knots a priori in the final model greatly improves fit at the young ages.

Effectiveness of post-COVID-19 primary care attendance in improving survival in very old patients with multimorbidity: a territory-wide target trial emulation.

BACKGROUND: Older adults with multimorbidity are at high risk of mortality following COVID-19 hospitalisation. However, the potential benefit of timely primary care follow-up on severe outcomes post-COVID-19 has not been well established. AIM: To examine the effectiveness of attending general outpatient within 30 days after discharge from COVID-19 on 1-year survival among older adults aged ≥85 years, with multimorbidity. METHOD: We emulated a target trial using a comprehensive public healthcare database in Hong Kong. The cloning-censoring-weighting technique was used to minimise immortal time bias and confounding bias by adjusting for demographics, hospitalisation duration and ICU admission, baseline chronic conditions, and medication history. The outcome included all-cause and cause-specific mortality. RESULTS: Of 6183 eligible COVID-19 survivors, the all-cause mortality rate following COVID-19 hospitalisation was lower in general out-patient clinics (GOPC) group compared to non-GOPC group (17.1 versus 42.8 deaths per 100 person-year). After adjustment, primary care consultations within 30 days after discharge were associated with a significantly greater 1-year survival (difference in 1-year survival: 11.2%, 95% CI = 8.1% to 14.4%). We also observed better survival from respiratory diseases in the GOPC group. In a sensitivity analysis for different grace period lengths, we found that the earlier participants had a GOPC visit after COVID-19 discharge, the better the survival. CONCLUSION: Timely primary care consultations after discharge may improve survival following COVID-19 hospitalisation among older adults aged ≥85 years, with multimorbidity. Expanding primary care services and implementing follow-up mechanisms are crucial to support this vulnerable population's recovery and well-being.