Modeling the multi-state natural history of rare diseases with heterogeneous individual patient data: A simulation study.

Multi-state survival models are used to represent the natural history of a disease, forming the basis of a health technology assessment comparing a novel treatment to current practice. Constructing such models for rare diseases is problematic, since evidence sources are typically much sparser and more heterogeneous. This simulation study investigated different one-stage and two-stage approaches to meta-analyzing individual patient data (IPD) in a multi-state survival setting when the number and size of studies being meta-analyzed are small. The objective was to assess methods of different complexity to see when they are accurate, when they are inaccurate and when they struggle to converge due to the sparsity of data. Biologically plausible multi-state IPD were simulated from study- and transition-specific hazard functions. One-stage frailty and two-stage stratified models were estimated, and compared to a base case model that did not account for study heterogeneity. Convergence and the bias/coverage of population-level transition probabilities to, and lengths of stay in, each state were used to assess model performance. A real-world application to Duchenne Muscular Dystrophy, a neuromuscular rare disease, was conducted, and a software demonstration is provided. Models not accounting for study heterogeneity were consistently out-performed by two-stage models. Frailty models struggled to converge, particularly in scenarios of low heterogeneity, and predictions from models that did converge were also subject to bias. Stratified models may be better suited to meta-analyzing disparate sources of IPD in rare disease natural history/economic modeling, as they converge more consistently and produce less biased predictions of lengths of stay.

Evaluation of Flexible Parametric Relative Survival Approaches for Enforcing Long-Term Constraints When Extrapolating All-Cause Survival.

OBJECTIVES: Parametric models are used to estimate the lifetime benefit of an intervention beyond the range of trial follow-up. Recent recommendations have suggested more flexible survival approaches and the use of external data when extrapolating. Both of these can be realized by using flexible parametric relative survival modeling. The overall aim of this article is to introduce and contrast various approaches for applying constraints on the long-term disease-related (excess) mortality including cure models and evaluate the consequent implications for extrapolation. METHODS: We describe flexible parametric relative survival modeling approaches. We then introduce various options for constraining the long-term excess mortality and compare the performance of each method in simulated data. These methods include fitting a standard flexible parametric relative survival model, enforcing statistical cure, and forcing the long-term excess mortality to converge to a constant. We simulate various scenarios, including where statistical cure is reasonable and where the long-term excess mortality persists. RESULTS: The compared approaches showed similar survival fits within the follow-up period. However, when extrapolating the all-cause survival beyond trial follow-up, there is variation depending on the assumption made about the long-term excess mortality. Altering the time point from which the excess mortality is constrained enables further flexibility. CONCLUSIONS: The various constraints can lead to applying explicit assumptions when extrapolating, which could lead to more plausible survival extrapolations. The inclusion of general population mortality directly into the model-building process, which is possible for all considered approaches, should be adopted more widely in survival extrapolation in health technology assessment.

Perils of Randomized Controlled Trial Survival Extrapolation Assuming Treatment Effect Waning: Why the Distinction Between Marginal and Conditional Estimates Matters.

OBJECTIVES: A long-term, constant, protective treatment effect is a strong assumption when extrapolating survival beyond clinical trial follow-up; hence, sensitivity to treatment effect waning is commonly assessed for economic evaluations. Forcing a hazard ratio (HR) to 1 does not necessarily estimate loss of individual-level treatment effect accurately because of HR selection bias. A simulation study was designed to explore the behavior of marginal HRs under a waning conditional (individual-level) treatment effect and demonstrate bias in forcing a marginal HR to 1 when the estimand is "survival difference with individual-level waning". METHODS: Data were simulated under 4 parameter combinations (varying prognostic strength of heterogeneity and treatment effect). Time-varying marginal HRs were estimated in scenarios where the true conditional HR attenuated to 1. Restricted mean survival time differences, estimated having constrained the marginal HR to 1, were compared with true values to assess bias induced by marginal constraints. RESULTS: Under loss of conditional treatment effect, the marginal HR took a value >1 because of covariate imbalances. Constraining this value to 1 lead to restricted mean survival time difference bias of up to 0.8 years (57% increase). Inflation of effect size estimates also increased with the magnitude of initial protective treatment effect. CONCLUSIONS: Important differences exist between survival extrapolations assuming marginal versus conditional treatment effect waning. When a marginal HR is constrained to 1 to assess efficacy under individual-level treatment effect waning, the survival benefits associated with the new treatment will be overestimated, and incremental cost-effectiveness ratios will be underestimated.

Shifting incidence and survival of epithelial ovarian cancer (1995-2014): A SurvMark-2 study.

The aim of the study is to provide a comprehensive assessment of incidence and survival trends of epithelial ovarian cancer (EOC) by histological subtype across seven high income countries (Australia, Canada, Denmark, Ireland, New Zealand, Norway and the United Kingdom). Data on invasive EOC diagnosed in women aged 15 to 99 years during 1995 to 2014 were obtained from 20 cancer registries. Age standardized incidence rates and average annual percentage change were calculated by subtype for all ages and age groups (15-64 and 65-99 years). Net survival (NS) was estimated by subtype, age group and 5-year period using Pohar-Perme estimator. Our findings showed marked increase in serous carcinoma incidence was observed between 1995 and 2014 among women aged 65 to 99 years with average annual increase ranging between 2.2% and 5.8%. We documented a marked decrease in the incidence of adenocarcinoma "not otherwise specified" with estimates ranging between 4.4% and 7.4% in women aged 15 to 64 years and between 2.0% and 3.7% among the older age group. Improved survival, combining all EOC subtypes, was observed for all ages combined over the 20-year study period in all countries with 5-year NS absolute percent change ranging between 5.0 in Canada and 12.6 in Denmark. Several factors such as changes in guidelines and advancement in diagnostic tools may potentially influence the observed shift in histological subtypes and temporal trends. Progress in clinical management and treatment over the past decades potentially plays a role in the observed improvements in EOC survival.

Improving communication of cancer survival statistics-feasibility of implementing model-based algorithms in routine publications.

BACKGROUND: Routine reporting of cancer patient survival is important, both to monitor the effectiveness of health care and to inform about prognosis following a cancer diagnosis. A range of different survival measures exist, each serving different purposes and targeting different audiences. It is important that routine publications expand on current practice and provide estimates on a wider range of survival measures. We examine the feasibility of automated production of such statistics. METHODS: We used data on 23 cancer sites obtained from the Cancer Registry of Norway (CRN). We propose an automated way of estimating flexible parametric relative survival models and calculating estimates of net survival, crude probabilities, and loss in life expectancy across many cancer sites and subgroups of patients. RESULTS: For 21 of 23 cancer sites, we were able to estimate survival models without assuming proportional hazards. Reliable estimates of all desired measures were obtained for all cancer sites. DISCUSSION: It may be challenging to implement new survival measures in routine publications as it can require the application of modeling techniques. We propose a way of automating the production of such statistics and show that we can obtain reliable estimates across a range of measures and subgroups of patients.

Using temporal recalibration to improve the calibration of risk prediction models in competing risk settings when there are trends in survival over time.

We have previously proposed temporal recalibration to account for trends in survival over time to improve the calibration of predictions from prognostic models for new patients. This involves first estimating the predictor effects using data from all individuals (full dataset) and then re-estimating the baseline using a subset of the most recent data whilst constraining the predictor effects to remain the same. In this article, we demonstrate how temporal recalibration can be applied in competing risk settings by recalibrating each cause-specific (or subdistribution) hazard model separately. We illustrate this using an example of colon cancer survival with data from the Surveillance Epidemiology and End Results (SEER) program. Data from patients diagnosed in 1995-2004 were used to fit two models for deaths due to colon cancer and other causes respectively. We discuss considerations that need to be made in order to apply temporal recalibration such as the choice of data used in the recalibration step. We also demonstrate how to assess the calibration of these models in new data for patients diagnosed subsequently in 2005. Comparison was made to a standard analysis (when improvements over time are not taken into account) and a period analysis which is similar to temporal recalibration but differs in the data used to estimate the predictor effects. The 10-year calibration plots demonstrated that using the standard approach over-estimated the risk of death due to colon cancer and the total risk of death and that calibration was improved using temporal recalibration or period analysis.

Flexible parametric methods for calculating life expectancy in small populations.

BACKGROUND: Life expectancy is a simple measure of assessing health differences between two or more populations but current life expectancy calculations are not reliable for small populations. A potential solution to this is to borrow strength from larger populations from the same source, but this has not formally been investigated. METHODS: Using data on 451,222 individuals from the Clinical Practice Research Datalink on the presence/absence of intellectual disability and type 2 diabetes mellitus, we compared stratified and combined flexible parametric models, and Chiang's methods, for calculating life expectancy. Confidence intervals were calculated using the Delta method, Chiang's adjusted life table approach and bootstrapping. RESULTS: The flexible parametric models allowed calculation of life expectancy by exact age and beyond traditional life expectancy age thresholds. The combined model that fit age interaction effects as a spline term provided less bias and greater statistical precision for small covariate subgroups by borrowing strength from the larger subgroups. However, careful consideration of the distribution of events in the smallest group was needed. CONCLUSIONS: Life expectancy is a simple measure to compare health differences between populations. The use of combined flexible parametric methods to calculate life expectancy in small samples has shown promising results by allowing life expectancy to be modelled by exact age, greater statistical precision, less bias and prediction of different covariate patterns without stratification. We recommend further investigation of their application for both policymakers and researchers.

Protection of cochlear synapses from noise-induced excitotoxic trauma by blockade of Ca2+-permeable AMPA receptors.

Exposure to loud sound damages the postsynaptic terminals of spiral ganglion neurons (SGNs) on cochlear inner hair cells (IHCs), resulting in loss of synapses, a process termed synaptopathy. Glutamatergic neurotransmission via α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA)-type receptors is required for synaptopathy, and here we identify a possible involvement of GluA2-lacking Ca2+-permeable AMPA receptors (CP-AMPARs) using IEM-1460, which has been shown to block GluA2-lacking AMPARs. In CBA/CaJ mice, a 2-h exposure to 100-dB sound pressure level octave band (8 to 16 kHz) noise results in no permanent threshold shift but does cause significant synaptopathy and a reduction in auditory brainstem response (ABR) wave-I amplitude. Chronic intracochlear perfusion of IEM-1460 in artificial perilymph (AP) into adult CBA/CaJ mice prevented the decrease in ABR wave-I amplitude and the synaptopathy relative to intracochlear perfusion of AP alone. Interestingly, IEM-1460 itself did not affect the ABR threshold, presumably because GluA2-containing AMPARs can sustain sufficient synaptic transmission to evoke low-threshold responses during blockade of GluA2-lacking AMPARs. On individual postsynaptic densities, we observed GluA2-lacking nanodomains alongside regions with robust GluA2 expression, consistent with the idea that individual synapses have both CP-AMPARs and Ca2+-impermeable AMPARs. SGNs innervating the same IHC differ in their relative vulnerability to noise. We found local heterogeneity among synapses in the relative abundance of GluA2 subunits that may underlie such differences in vulnerability. We propose a role for GluA2-lacking CP-AMPARs in noise-induced cochlear synaptopathy whereby differences among synapses account for differences in excitotoxic susceptibility. These data suggest a means of maintaining normal hearing thresholds while protecting against noise-induced synaptopathy, via selective blockade of CP-AMPARs.

International variation in oesophageal and gastric cancer survival 2012-2014: differences by histological subtype and stage at diagnosis (an ICBP SURVMARK-2 population-based study).

OBJECTIVE: To provide the first international comparison of oesophageal and gastric cancer survival by stage at diagnosis and histological subtype across high-income countries with similar access to healthcare. METHODS: As part of the ICBP SURVMARK-2 project, data from 28 923 patients with oesophageal cancer and 25 946 patients with gastric cancer diagnosed during 2012-2014 from 14 cancer registries in seven countries (Australia, Canada, Denmark, Ireland, New Zealand, Norway and the UK) were included. 1-year and 3-year age-standardised net survival were estimated by stage at diagnosis, histological subtype (oesophageal adenocarcinoma (OAC) and oesophageal squamous cell carcinoma (OSCC)) and country. RESULTS: Oesophageal cancer survival was highest in Ireland and lowest in Canada at 1 (50.3% vs 41.3%, respectively) and 3 years (27.0% vs 19.2%) postdiagnosis. Survival from gastric cancer was highest in Australia and lowest in the UK, for both 1-year (55.2% vs 44.8%, respectively) and 3-year survival (33.7% vs 22.3%). Most patients with oesophageal and gastric cancer had regional or distant disease, with proportions ranging between 56% and 90% across countries. Stage-specific analyses showed that variation between countries was greatest for localised disease, where survival ranged between 66.6% in Australia and 83.2% in the UK for oesophageal cancer and between 75.5% in Australia and 94.3% in New Zealand for gastric cancer at 1-year postdiagnosis. While survival for OAC was generally higher than that for OSCC, disparities across countries were similar for both histological subtypes. CONCLUSION: Survival from oesophageal and gastric cancer varies across high-income countries including within stage groups, particularly for localised disease. Disparities can partly be explained by earlier diagnosis resulting in more favourable stage distributions, and distributions of histological subtypes of oesophageal cancer across countries. Yet, differences in treatment, and also in cancer registration practice and the use of different staging methods and systems, across countries may have impacted the comparisons. While primary prevention remains key, advancements in early detection research are promising and will likely allow for additional risk stratification and survival improvements in the future.

Obtaining long-term stage-specific relative survival estimates in the presence of incomplete historical stage information.

BACKGROUND: Completeness of recording for cancer stage at diagnosis is often historically poor in cancer registries, making it challenging to provide long-term stage-specific survival estimates. Stage-specific survival differences are driven by differences in short-term prognosis, meaning estimated survival metrics using period analysis are unlikely to be sensitive to imputed historical stage data. METHODS: We used data from the Surveillance, Epidemiology, and End Results (SEER) Program for lung, colon and breast cancer. To represent missing data patterns in less complete registry data, we artificially inflated the proportion of missing stage information conditional on stage at diagnosis and calendar year of diagnosis. Period analysis was applied and missing stage at diagnosis information was imputed under four different conditions to emulate extreme imputed stage distributions. RESULTS: We fit a flexible parametric model for each cancer stage on the excess hazard scale and the differences in stage-specific marginal relative survival were assessed. Estimates were also obtained from non-parametric approaches for validation. There was little difference between the 10-year stage-specific marginal relative survival estimates, regardless of the assumed historical stage distribution. CONCLUSIONS: When conducting a period analysis, multiple imputation can be used to obtain stage-specific long-term estimates of relative survival, even when the historical stage information is largely incomplete.

Five ways to improve international comparisons of cancer survival: lessons learned from ICBP SURVMARK-2.

BACKGROUND: Comparisons of population-based cancer survival between countries are important to benchmark the overall effectiveness of cancer management. The International Cancer Benchmarking Partnership (ICBP) Survmark-2 study aims to compare survival in seven high-income countries across eight cancer sites and explore reasons for the observed differences. A critical aspect in ensuring comparability in the reported survival estimates are similarities in practice across cancer registries. While ICBP Survmark-2 has shown these differences are unlikely to explain the observed differences in cancer-specific survival between countries, it is important to keep in mind potential biases linked to registry practice and understand their likely impact. METHODS: Based on experiences gained within ICBP Survmark-2, we have developed a set of recommendations that seek to optimally harmonise cancer registry datasets to improve future benchmarking exercises. RESULTS: Our recommendations stem from considering the impact on cancer survival estimates in five key areas: (1) the completeness of the registry and the availability of registration sources; (2) the inclusion of death certification as a source of identifying cases; (3) the specification of the date of incidence; (4) the approach to handling multiple primary tumours and (5) the quality of linkage of cases to the deaths register. CONCLUSION: These recommendations seek to improve comparability whilst maintaining the opportunity to understand and act upon international variations in outcomes among cancer patients.

Pancreatic cancer survival by stage and age in seven high-income countries (ICBP SURVMARK-2): a population-based study.

BACKGROUND: The global burden of pancreatic cancer has steadily increased, while the prognosis after pancreatic cancer diagnosis remains poor. This study aims to compare the stage- and age-specific pancreatic cancer net survival (NS) for seven high-income countries: Australia, Canada, Denmark, Ireland, New Zealand, Norway, and United Kingdom. METHODS: The study included over 35,000 pancreatic cancer cases diagnosed during 2012-2014, followed through 31 December 2015. The stage- and age-specific NS were calculated using the Pohar-Perme estimator. RESULTS: Pancreatic cancer survival estimates were low across all 7 countries, with 1-year NS ranging from 21.1% in New Zealand to 30.9% in Australia, and 3-year NS from 6.6% in the UK to 10.9% in Australia. Most pancreatic cancers were diagnosed with distant stage, ranging from 53.9% in Ireland to 83.3% in New Zealand. While survival differences were evident between countries across all stage categories at one year after diagnosis, this survival advantage diminished, particularly in cases with distant stage. CONCLUSION: This study demonstrated the importance of stage and age at diagnosis in pancreatic cancer survival. Although progress has been made in improving pancreatic cancer prognosis, the disease is highly fatal and will remain so without major breakthroughs in the early diagnosis and management.

International differences in lung cancer survival by sex, histological type and stage at diagnosis: an ICBP SURVMARK-2 Study.

INTRODUCTION: Lung cancer has a poor prognosis that varies internationally when assessed by the two major histological subgroups (non-small cell (NSCLC) and small cell (SCLC)). METHOD: 236 114 NSCLC and 43 167 SCLC cases diagnosed during 2010-2014 in Australia, Canada, Denmark, Ireland, New Zealand, Norway and the UK were included in the analyses. One-year and 3-year age-standardised net survival (NS) was estimated by sex, histological type, stage and country. RESULTS: One-year and 3-year NS was consistently higher for Canada and Norway, and lower for the UK, New Zealand and Ireland, irrespective of stage at diagnosis. Three-year NS for NSCLC ranged from 19.7% for the UK to 27.1% for Canada for men and was consistently higher for women (25.3% in the UK; 35.0% in Canada) partly because men were diagnosed at more advanced stages. International differences in survival for NSCLC were largest for regional stage and smallest at the advanced stage. For SCLC, 3-year NS also showed a clear female advantage with the highest being for Canada (13.8% for women; 9.1% for men) and Norway (12.8% for women; 9.7% for men). CONCLUSION: Distribution of stage at diagnosis among lung cancer cases differed by sex, histological subtype and country, which may partly explain observed survival differences. Yet, survival differences were also observed within stages, suggesting that quality of treatment, healthcare system factors and prevalence of comorbid conditions may also influence survival. Other possible explanations include differences in data collection practice, as well as differences in histological verification, staging and coding across jurisdictions.

Minimum sample size calculations for external validation of a clinical prediction model with a time-to-event outcome.

Previous articles in Statistics in Medicine describe how to calculate the sample size required for external validation of prediction models with continuous and binary outcomes. The minimum sample size criteria aim to ensure precise estimation of key measures of a model's predictive performance, including measures of calibration, discrimination, and net benefit. Here, we extend the sample size guidance to prediction models with a time-to-event (survival) outcome, to cover external validation in datasets containing censoring. A simulation-based framework is proposed, which calculates the sample size required to target a particular confidence interval width for the calibration slope measuring the agreement between predicted risks (from the model) and observed risks (derived using pseudo-observations to account for censoring) on the log cumulative hazard scale. Precise estimation of calibration curves, discrimination, and net-benefit can also be checked in this framework. The process requires assumptions about the validation population in terms of the (i) distribution of the model's linear predictor and (ii) event and censoring distributions. Existing information can inform this; in particular, the linear predictor distribution can be approximated using the C-index or Royston's D statistic from the model development article, together with the overall event risk. We demonstrate how the approach can be used to calculate the sample size required to validate a prediction model for recurrent venous thromboembolism. Ideally the sample size should ensure precise calibration across the entire range of predicted risks, but must at least ensure adequate precision in regions important for clinical decision-making. Stata and R code are provided.

Generating high-fidelity synthetic time-to-event datasets to improve data transparency and accessibility.

BACKGROUND: A lack of available data and statistical code being published alongside journal articles provides a significant barrier to open scientific discourse, and reproducibility of research. Information governance restrictions inhibit the active dissemination of individual level data to accompany published manuscripts. Realistic, high-fidelity time-to-event synthetic data can aid in the acceleration of methodological developments in survival analysis and beyond by enabling researchers to access and test published methods using data similar to that which they were developed on. METHODS: We present methods to accurately emulate the covariate patterns and survival times found in real-world datasets using synthetic data techniques, without compromising patient privacy. We model the joint covariate distribution of the original data using covariate specific sequential conditional regression models, then fit a complex flexible parametric survival model from which to generate survival times conditional on individual covariate patterns. We recreate the administrative censoring mechanism using the last observed follow-up date information from the initial dataset. Metrics for evaluating the accuracy of the synthetic data, and the non-identifiability of individuals from the original dataset, are presented. RESULTS: We successfully create a synthetic version of an example colon cancer dataset consisting of 9064 patients which aims to show good similarity to both covariate distributions and survival times from the original data, without containing any exact information from the original data, therefore allowing them to be published openly alongside research. CONCLUSIONS: We evaluate the effectiveness of the methods for constructing synthetic data, as well as providing evidence that there is minimal risk that a given patient from the original data could be identified from their individual unique patient information. Synthetic datasets using this methodology could be made available alongside published research without breaching data privacy protocols, and allow for data and code to be made available alongside methodological or applied manuscripts to greatly improve the transparency and accessibility of medical research.

Immortal time bias for life-long conditions in retrospective observational studies using electronic health records.

BACKGROUND: Immortal time bias is common in observational studies but is typically described for pharmacoepidemiology studies where there is a delay between cohort entry and treatment initiation. METHODS: This study used the Clinical Practice Research Datalink (CPRD) and linked national mortality data in England from 2000 to 2019 to investigate immortal time bias for a specific life-long condition, intellectual disability. Life expectancy (Chiang's abridged life table approach) was compared for 33,867 exposed and 980,586 unexposed individuals aged 10+ years using five methods: (1) treating immortal time as observation time; (2) excluding time before date of first exposure diagnosis; (3) matching cohort entry to first exposure diagnosis; (4) excluding time before proxy date of inputting first exposure diagnosis (by the physician); and (5) treating exposure as a time-dependent measure. RESULTS: When not considered in the design or analysis (Method 1), immortal time bias led to disproportionately high life expectancy for the exposed population during the first calendar period (additional years expected to live: 2000-2004: 65.6 [95% CI: 63.6,67.6]) compared to the later calendar periods (2005-2009: 59.9 [58.8,60.9]; 2010-2014: 58.0 [57.1,58.9]; 2015-2019: 58.2 [56.8,59.7]). Date of entry of diagnosis (Method 4) was unreliable in this CPRD cohort. The final methods (Method 2, 3 and 5) appeared to solve the main theoretical problem but residual bias may have remained. CONCLUSIONS: We conclude that immortal time bias is a significant issue for studies of life-long conditions that use electronic health record data and requires careful consideration of how clinical diagnoses are entered onto electronic health record systems.

Estimating causal effects in the presence of competing events using regression standardisation with the Stata command standsurv.

BACKGROUND: When interested in a time-to-event outcome, competing events that prevent the occurrence of the event of interest may be present. In the presence of competing events, various estimands have been suggested for defining the causal effect of treatment on the event of interest. Depending on the estimand, the competing events are either accommodated or eliminated, resulting in causal effects with different interpretations. The former approach captures the total effect of treatment on the event of interest while the latter approach captures the direct effect of treatment on the event of interest that is not mediated by the competing event. Separable effects have also been defined for settings where the treatment can be partitioned into two components that affect the event of interest and the competing event through different causal pathways. METHODS: We outline various causal effects that may be of interest in the presence of competing events, including total, direct and separable effects, and describe how to obtain estimates using regression standardisation with the Stata command standsurv. Regression standardisation is applied by obtaining the average of individual estimates across all individuals in a study population after fitting a survival model. RESULTS: With standsurv several contrasts of interest can be calculated including differences, ratios and other user-defined functions. Confidence intervals can also be obtained using the delta method. Throughout we use an example analysing a publicly available dataset on prostate cancer to allow the reader to replicate the analysis and further explore the different effects of interest. CONCLUSIONS: Several causal effects can be defined in the presence of competing events and, under assumptions, estimates of those can be obtained using regression standardisation with the Stata command standsurv. The choice of which causal effect to define should be given careful consideration based on the research question and the audience to which the findings will be communicated.

Non-parametric estimation of reference adjusted, standardised probabilities of all-cause death and death due to cancer for population group comparisons.

BACKGROUND: Ensuring fair comparisons of cancer survival statistics across population groups requires careful consideration of differential competing mortality due to other causes, and adjusting for imbalances over groups in other prognostic covariates (e.g. age). This has typically been achieved using comparisons of age-standardised net survival, with age standardisation addressing covariate imbalance, and the net estimates removing differences in competing mortality from other causes. However, these estimates lack ease of interpretability. In this paper, we motivate an alternative non-parametric approach that uses a common rate of other cause mortality across groups to give reference-adjusted estimates of the all-cause and cause-specific crude probability of death in contrast to solely reporting net survival estimates. METHODS: We develop the methodology for a non-parametric equivalent of standardised and reference adjusted crude probabilities of death, building on the estimation of non-parametric crude probabilities of death. We illustrate the approach using regional comparisons of survival following a diagnosis of rectal cancer for men in England. We standardise to the covariate distribution and other cause mortality of England as a whole to offer comparability, but with close approximation to the observed all-cause region-specific mortality. RESULTS: The approach gives comparable estimates to observed crude probabilities of death, but allows direct comparison across population groups with different covariate profiles and competing mortality patterns. In our illustrative example, we show that regional variations in survival following a diagnosis of rectal cancer persist even after accounting for the variation in deprivation, age at diagnosis and other cause mortality. CONCLUSIONS: The methodological approach of using standardised and reference adjusted metrics offers an appealing approach for future cancer survival comparison studies and routinely published cancer statistics. Our non-parametric estimation approach through the use of weighting offers the ability to estimate comparable survival estimates without the need for statistical modelling.

Comparing the number and length of primary care consultations in people with and without intellectual disabilities and health needs: observational cohort study using electronic health records.

BACKGROUND: In the United Kingdom, 15-min appointments with the general practitioner (GP) are recommended for people with complex health conditions, including intellectual disabilities and health needs, but we do not know whether this happens. AIMS: We compared number and length of primary care consultations (GP, nurse, other allied health, other) for people with and without intellectual disabilities and health needs. METHODS: Linked primary care data from the Clinical Practice Research Datalink (CPRD) in England were used to investigate face-to-face and telephone primary care consultations in 2017-2019. Health needs investigated were: epilepsy; incontinence; severe visual/hearing impairments; severe mobility difficulties; cerebral palsy; and percutaneous endoscopic gastrostomy feeding. Age and gender-standardized consultation rates per year (Poisson), duration of consultations, and the proportion of "long consultations" (≥15 min) were reported. RESULTS: People with intellectual disabilities (n = 7,794) had 1.9 times as many GP consultations per year as those without (n = 176,807; consultation rate ratio = 1.87 [95% confidence interval 1.86-1.89]). Consultation rates with nurses and allied healthcare professionals were also twice as high. Mean GP consultation time was 9-10 min regardless of intellectual disability/health need status. Long GP consultations were less common in people with intellectual disabilities (18.2% [17.8-18.7] vs. 20.9% [20.8-21.0]). Long consultations with practice nurses were more common in people with health needs, particularly severe visual loss. CONCLUSIONS: People with intellectual disabilities and/or health needs tend to have more, rather than longer, GP consultations compared with the rest of the population. We recommend further investigation into the role of practice nurses to support people with intellectual disabilities and health needs.

LRRC52 regulates BK channel function and localization in mouse cochlear inner hair cells.

The perception of sound relies on sensory hair cells in the cochlea that convert the mechanical energy of sound into release of glutamate onto postsynaptic auditory nerve fibers. The hair cell receptor potential regulates the strength of synaptic transmission and is shaped by a variety of voltage-dependent conductances. Among these conductances, the Ca2+- and voltage-activated large conductance Ca2+-activated K+ channel (BK) current is prominent, and in mammalian inner hair cells (IHCs) displays unusual properties. First, BK currents activate at unprecedentedly negative membrane potentials (-60 mV) even in the absence of intracellular Ca2+ elevations. Second, BK channels are positioned in clusters away from the voltage-dependent Ca2+ channels that mediate glutamate release from IHCs. Here, we test the contributions of two recently identified leucine-rich-repeat-containing (LRRC) regulatory γ subunits, LRRC26 and LRRC52, to BK channel function and localization in mouse IHCs. Whereas BK currents and channel localization were unaltered in IHCs from Lrrc26 knockout (KO) mice, BK current activation was shifted more than +200 mV in IHCs from Lrrc52 KO mice. Furthermore, the absence of LRRC52 disrupted BK channel localization in the IHCs. Given that heterologous coexpression of LRRC52 with BK α subunits shifts BK current gating about -90 mV, to account for the profound change in BK activation range caused by removal of LRRC52, we suggest that additional factors may help define the IHC BK gating range. LRRC52, through stabilization of a macromolecular complex, may help retain some other components essential both for activation of BK currents at negative membrane potentials and for appropriate BK channel positioning.