Testing the causal relationships of physical activity and sedentary behaviour with mental health and substance use disorders: a Mendelian randomisation study.

Observational studies suggest that physical activity can reduce the risk of mental health and substance use disorders. However, it is unclear whether this relationship is causal or explained by confounding bias (e.g., common underlying causes or reverse causality). We investigated the bidirectional causal relationship of physical activity (PA) and sedentary behaviour (SB) with ten mental health and substance use disorders, applying two-sample Mendelian Randomisation (MR). Genetic instruments for the exposures and outcomes were derived from the largest available, non-overlapping genome-wide association studies (GWAS). Summary-level data for objectively assessed PA (accelerometer-based average activity, moderate activity, and walking) and SB and self-reported moderate-to-vigorous PA were obtained from the UK Biobank. Data for mental health/substance use disorders were obtained from the Psychiatric Genomics Consortium and the GWAS and Sequencing Consortium of Alcohol and Nicotine Use. MR estimates were combined using inverse variance weighted meta-analysis (IVW). Sensitivity analyses were conducted to assess the robustness of the results. Accelerometer-based average PA was associated with a lower risk of depression (b = -0.043, 95% CI: -0.071 to -0.016, effect size[OR] = 0.957) and cigarette smoking (b = -0.026; 95% CI: -0.035 to -0.017, effect size[ß] = -0.022). Accelerometer-based SB decreased the risk of anorexia (b = -0.341, 95% CI: -0.530 to -0.152, effect size[OR] = 0.711) and schizophrenia (b = -0.230; 95% CI: -0.285 to -0.175, effect size[OR] = 0.795). However, we found evidence of reverse causality in the relationship between SB and schizophrenia. Further, PTSD, bipolar disorder, anorexia, and ADHD were all associated with increased PA. This study provides evidence consistent with a causal protective effect of objectively assessed but not self-reported PA on reduced depression and cigarette smoking. Objectively assessed SB had a protective relationship with anorexia. Enhancing PA may be an effective intervention strategy to reduce depressive symptoms and addictive behaviours, while promoting sedentary or light physical activities may help to reduce the risk of anorexia in at-risk individuals.

Theory and performance of substitution models for estimating relative causal effects in nutritional epidemiology.

BACKGROUND: Estimating relative causal effects (i.e., "substitution effects") is a common aim of nutritional research. In observational data, this is usually attempted using 1 of 2 statistical modeling approaches: the leave-one-out model and the energy partition model. Despite their widespread use, there are concerns that neither approach is well understood in practice. OBJECTIVES: We aimed to explore and illustrate the theory and performance of the leave-one-out and energy partition models for estimating substitution effects in nutritional epidemiology. METHODS: Monte Carlo data simulations were used to illustrate the theory and performance of both the leave-one-out model and energy partition model, by considering 3 broad types of causal effect estimands: 1) direct substitutions of the exposure with a single component, 2) inadvertent substitutions of the exposure with several components, and 3) average relative causal effects of the exposure instead of all other dietary sources. Models containing macronutrients, foods measured in calories, and foods measured in grams were all examined. RESULTS: The leave-one-out and energy partition models both performed equally well when the target estimand involved substituting a single exposure with a single component, provided all variables were measured in the same units. Bias occurred when the substitution involved >1 substituting component. Leave-one-out models that examined foods in mass while adjusting for total energy intake evaluated obscure estimands. CONCLUSIONS: Regardless of the approach, substitution models need to be constructed from clearly defined causal effect estimands. Estimands involving a single exposure and a single substituting component are typically estimated more accurately than estimands involving more complex substitutions. The practice of examining foods measured in grams or portions while adjusting for total energy intake is likely to deliver obscure relative effect estimands with unclear interpretations.

Estimating the effects of lockdown timing on COVID-19 cases and deaths in England: A counterfactual modelling study.

BACKGROUND: During the first wave of the COVID-19 pandemic, the United Kingdom experienced one of the highest per-capita death tolls worldwide. It is debated whether this may partly be explained by the relatively late initiation of voluntary social distancing and mandatory lockdown measures. In this study, we used simulations to estimate the number of cases and deaths that would have occurred in England by 1 June 2020 if these interventions had been implemented one or two weeks earlier, and the impact on the required duration of lockdown. METHODS: Using official reported data on the number of Pillar 1 lab-confirmed cases of COVID-19 and associated deaths occurring in England from 3 March to 1 June, we modelled: the natural (i.e. observed) growth of cases, and the counterfactual (i.e. hypothetical) growth of cases that would have occurred had measures been implemented one or two weeks earlier. Under each counterfactual condition, we estimated the expected number of deaths and the time required to reach the incidence observed under natural growth on 1 June. RESULTS: Introducing measures one week earlier would have reduced by 74% the number of confirmed COVID-19 cases in England by 1 June, resulting in approximately 21,000 fewer hospital deaths and 34,000 fewer total deaths; the required time spent in full lockdown could also have been halved, from 69 to 35 days. Acting two weeks earlier would have reduced cases by 93%, resulting in between 26,000 and 43,000 fewer deaths. CONCLUSIONS: Our modelling supports the claim that the relatively late introduction of social distancing and lockdown measures likely increased the scale, severity, and duration of the first wave of COVID-19 in England. Our results highlight the importance of acting swiftly to minimise the spread of an infectious disease when case numbers are increasing exponentially.

Analyses of 'change scores' do not estimate causal effects in observational data.

BACKGROUND: In longitudinal data, it is common to create 'change scores' by subtracting measurements taken at baseline from those taken at follow-up, and then to analyse the resulting 'change' as the outcome variable. In observational data, this approach can produce misleading causal-effect estimates. The present article uses directed acyclic graphs (DAGs) and simple simulations to provide an accessible explanation for why change scores do not estimate causal effects in observational data. METHODS: Data were simulated to match three general scenarios in which the outcome variable at baseline was a (i) 'competing exposure' (i.e. a cause of the outcome that is neither caused by nor causes the exposure), (ii) confounder or (iii) mediator for the total causal effect of the exposure variable at baseline on the outcome variable at follow-up. Regression coefficients were compared between change-score analyses and the appropriate estimator(s) for the total and/or direct causal effect(s). RESULTS: Change-score analyses do not provide meaningful causal-effect estimates unless the baseline outcome variable is a 'competing exposure' for the effect of the exposure on the outcome at follow-up. Where the baseline outcome is a confounder or mediator, change-score analyses evaluate obscure estimands, which may diverge substantially in magnitude and direction from the total and direct causal effects. CONCLUSION: Future observational studies that seek causal-effect estimates should avoid analysing change scores and adopt alternative analytical strategies.

Demonstration of functional rehabilitation treatment effects in children and young people after severe acquired brain injury.

PURPOSE: To examine relationships between functional outcomes after pediatric acquired brain injury (ABI) and measures of rehabilitation dose. METHODS: An observational study of children receiving residential neurorehabilitation after severe ABI. RESULTS: Basic total rehabilitation dose shows a paradoxical inverse relationship to global outcome. This is due to confounding by both initial injury severity and length of stay, and variation in treatment content for a given total rehabilitation dose. Content-aware rehabilitation dose measures show robust positive correlations between fractions of rehabilitation treatment received and plausibly related aspects of outcome: specifically, between rates of recovery of gross motor function and the fraction of rehabilitation effort directed to active practice and motor learning. This relationship was robust to adjustment for therapists' expectations of recovery. CONCLUSION: Content-aware measures of rehabilitation dose are robustly causally related to pertinent aspects of outcome. These findings are step toward a goal of comparative effectiveness research in pediatric neurorehabilitation.

Adjustment for energy intake in nutritional research: a causal inference perspective.

BACKGROUND: Four models are commonly used to adjust for energy intake when estimating the causal effect of a dietary component on an outcome: 1) the "standard model" adjusts for total energy intake, 2) the "energy partition model" adjusts for remaining energy intake, 3) the "nutrient density model" rescales the exposure as a proportion of total energy, and 4) the "residual model" indirectly adjusts for total energy by using a residual. It remains underappreciated that each approach evaluates a different estimand and only partially accounts for confounding by common dietary causes. OBJECTIVES: We aimed to clarify the implied causal estimand and interpretation of each model and evaluate their performance in reducing dietary confounding. METHODS: Semiparametric directed acyclic graphs and Monte Carlo simulations were used to identify the estimands and interpretations implied by each model and explore their performance in the absence or presence of dietary confounding. RESULTS: The "standard model" and the mathematically identical "residual model" estimate the average relative causal effect (i.e., a "substitution" effect) but provide biased estimates even in the absence of confounding. The "energy partition model" estimates the total causal effect but only provides unbiased estimates in the absence of confounding or when all other nutrients have equal effects on the outcome. The "nutrient density model" has an obscure interpretation but attempts to estimate the average relative causal effect rescaled as a proportion of total energy. Accurate estimates of both the total and average relative causal effects may instead be derived by simultaneously adjusting for all dietary components, an approach we term the "all-components model." CONCLUSIONS: Lack of awareness of the estimand differences and accuracy of the 4 modeling approaches may explain some of the apparent heterogeneity among existing nutritional studies. This raises serious questions regarding the validity of meta-analyses where different estimands have been inappropriately pooled.

Hepatectomy risk assessment with functional magnetic resonance imaging (HEPARIM).

BACKGROUND: Post hepatectomy liver failure (PHLF) remains a significant risk in patients undergoing curative liver resection for cancer, however currently available PHLF risk prediction investigations are not sufficiently accurate. The Hepatectomy risk assessment with functional magnetic resonance imaging trial (HEPARIM) aims to establish if quantitative MRI biomarkers of liver function & perfusion can be used to more accurately predict PHLF risk and FLR function, measured against indocyanine green (ICG) liver function test. METHODS: HEPARIM is an observational cohort study recruiting patients undergoing liver resection of 2 segments or more, prior to surgery patients will have both Dynamic Gadoxetate-enhanced (DGE) liver MRI and ICG testing. Day one post op ICG testing is repeated and R15 compared to the Gadoxetate Clearance (GC) of the future liver remnant (FLR-GC) as measure by preoperative DGE- MRI which is the primary outcome, and preoperative ICG R15 compared to GC of whole liver (WL-GC) as a secondary outcome. Data will be collected from medical records, biochemistry, pathology and radiology reports and used in a multi-variate analysis to the value of functional MRI and derive multivariant prediction models for future validation. DISCUSSION: If successful, this test will potentially provide an efficient means to quantitatively assess FLR function and PHLF risk enabling surgeons to push boundaries of liver surgery further while maintaining safe practice and thereby offering chance of cure to patients who would previously been deemed inoperable. MRI has the added benefit of already being part of the routine diagnostic pathway and as such would have limited additional burden on patients time or cost to health care systems. (Hepatectomy Risk Assessment With Functional Magnetic Resonance Imaging - Full Text View - ClinicalTrials.gov , n.d.) TRIAL REGISTRATION: ClinicalTrials.gov, ClinicalTrials.gov NCT04705194 - Registered 12th January 2021 - Retrospectively registered.

Early childhood weight gain: Latent patterns and body composition outcomes.

BACKGROUND: Despite early childhood weight gain being a key indicator of obesity risk, we do not have a good understanding of the different patterns that exist. OBJECTIVES: To identify and characterise distinct groups of children displaying similar early-life weight trajectories. METHODS: A growth mixture model captured heterogeneity in weight trajectories between 0 and 60 months in 1390 children in the Avon Longitudinal Study of Parents and Children. Differences between the classes in characteristics and body size/composition at 9 years were investigated. RESULTS: The best model had five classes. The "Normal" (45%) and "Normal after initial catch-down" (24%) classes were close to the 50th centile of a growth standard between 24 and 60 months. The "High-decreasing" (21%) and "Stable-high" (7%) classes peaked at the ~91st centile at 12-18 months, but while the former declined to the ~75th centile and comprised constitutionally big children, the latter did not. The "Rapidly increasing" (3%) class gained weight from below the 50th centile at 4 months to above the 91st centile at 60 months. By 9 years, their mean body mass index (BMI) placed them at the 98th centile. This class was characterised by the highest maternal BMI; highest parity; highest levels of gestational hypertension and diabetes; and the lowest socio-economic position. At 9 years, the "Rapidly increasing" class was estimated to have 68.2% (95% confidence interval [CI] 48.3, 88.1) more fat mass than the "Normal" class, but only 14.0% (95% CI 9.1, 18.9) more lean mass. CONCLUSIONS: Criteria used in growth monitoring practice are unlikely to consistently distinguish between the different patterns of weight gain reported here.

Latent class regression improves the predictive acuity and clinical utility of survival prognostication amongst chronic heart failure patients.

The present study aimed to compare the predictive acuity of latent class regression (LCR) modelling with: standard generalised linear modelling (GLM); and GLMs that include the membership of subgroups/classes (identified through prior latent class analysis; LCA) as alternative or additional candidate predictors. Using real world demographic and clinical data from 1,802 heart failure patients enrolled in the UK-HEART2 cohort, the study found that univariable GLMs using LCA-generated subgroup/class membership as the sole candidate predictor of survival were inferior to standard multivariable GLMs using the same four covariates as those used in the LCA. The inclusion of the LCA subgroup/class membership together with these four covariates as candidate predictors in a multivariable GLM showed no improvement in predictive acuity. In contrast, LCR modelling resulted in a 18-22% improvement in predictive acuity and provided a range of alternative models from which it would be possible to balance predictive acuity against entropy to select models that were optimally suited to improve the efficient allocation of clinical resources to address the differential risk of the outcome (in this instance, survival). These findings provide proof-of-principle that LCR modelling can improve the predictive acuity of GLMs and enhance the clinical utility of their predictions. These improvements warrant further attention and exploration, including the use of alternative techniques (including machine learning algorithms) that are also capable of generating latent class structure while determining outcome predictions, particularly for use with large and routinely collected clinical datasets, and with binary, count and continuous variables.

Distinct Body Mass Index Trajectories to Young-Adulthood Obesity and Their Different Cardiometabolic Consequences.

[Figure: see text].

Reflection on modern methods: generalized linear models for prognosis and intervention-theory, practice and implications for machine learning.

Prediction and causal explanation are fundamentally distinct tasks of data analysis. In health applications, this difference can be understood in terms of the difference between prognosis (prediction) and prevention/treatment (causal explanation). Nevertheless, these two concepts are often conflated in practice. We use the framework of generalized linear models (GLMs) to illustrate that predictive and causal queries require distinct processes for their application and subsequent interpretation of results. In particular, we identify five primary ways in which GLMs for prediction differ from GLMs for causal inference: (i) the covariates that should be considered for inclusion in (and possibly exclusion from) the model; (ii) how a suitable set of covariates to include in the model is determined; (iii) which covariates are ultimately selected and what functional form (i.e. parameterization) they take; (iv) how the model is evaluated; and (v) how the model is interpreted. We outline some of the potential consequences of failing to acknowledge and respect these differences, and additionally consider the implications for machine learning (ML) methods. We then conclude with three recommendations that we hope will help ensure that both prediction and causal modelling are used appropriately and to greatest effect in health research.

Use of directed acyclic graphs (DAGs) to identify confounders in applied health research: review and recommendations.

BACKGROUND: Directed acyclic graphs (DAGs) are an increasingly popular approach for identifying confounding variables that require conditioning when estimating causal effects. This review examined the use of DAGs in applied health research to inform recommendations for improving their transparency and utility in future research. METHODS: Original health research articles published during 1999-2017 mentioning 'directed acyclic graphs' (or similar) or citing DAGitty were identified from Scopus, Web of Science, Medline and Embase. Data were extracted on the reporting of: estimands, DAGs and adjustment sets, alongside the characteristics of each article's largest DAG. RESULTS: A total of 234 articles were identified that reported using DAGs. A fifth (n = 48, 21%) reported their target estimand(s) and half (n = 115, 48%) reported the adjustment set(s) implied by their DAG(s). Two-thirds of the articles (n = 144, 62%) made at least one DAG available. DAGs varied in size but averaged 12 nodes [interquartile range (IQR): 9-16, range: 3-28] and 29 arcs (IQR: 19-42, range: 3-99). The median saturation (i.e. percentage of total possible arcs) was 46% (IQR: 31-67, range: 12-100). 37% (n = 53) of the DAGs included unobserved variables, 17% (n = 25) included 'super-nodes' (i.e. nodes containing more than one variable) and 34% (n = 49) were visually arranged so that the constituent arcs flowed in the same direction (e.g. top-to-bottom). CONCLUSION: There is substantial variation in the use and reporting of DAGs in applied health research. Although this partly reflects their flexibility, it also highlights some potential areas for improvement. This review hence offers several recommendations to improve the reporting and use of DAGs in future research.

Time to reality check the promises of machine learning-powered precision medicine.

Machine learning methods, combined with large electronic health databases, could enable a personalised approach to medicine through improved diagnosis and prediction of individual responses to therapies. If successful, this strategy would represent a revolution in clinical research and practice. However, although the vision of individually tailored medicine is alluring, there is a need to distinguish genuine potential from hype. We argue that the goal of personalised medical care faces serious challenges, many of which cannot be addressed through algorithmic complexity, and call for collaboration between traditional methodologists and experts in medical machine learning to avoid extensive research waste.

A causal inference perspective on the analysis of compositional data.

BACKGROUND: Compositional data comprise the parts of some whole, for which all parts sum to that whole. They are prevalent in many epidemiological contexts. Although many of the challenges associated with analysing compositional data have been discussed previously, we do so within a formal causal framework by utilizing directed acyclic graphs (DAGs). METHODS: We depict compositional data using DAGs and identify two distinct effect estimands in the generic case: (i) the total effect, and (ii) the relative effect. We consider each in the context of three specific example scenarios involving compositional data: (1) the relationship between the economically active population and area-level gross domestic product; (2) the relationship between fat consumption and body weight; and (3) the relationship between time spent sedentary and body weight. For each, we consider the distinct interpretation of each effect, and the resulting implications for related analyses. RESULTS: For scenarios (1) and (2), both the total and relative effects may be identifiable and causally meaningful, depending upon the specific question of interest. For scenario (3), only the relative effect is identifiable. In all scenarios, the relative effect represents a joint effect, and thus requires careful interpretation. CONCLUSIONS: DAGs are useful for considering causal effects for compositional data. In all analyses involving compositional data, researchers should explicitly consider and declare which causal effect is sought and how it should be interpreted.

Response to: Simpson's Paradox is suppression, but Lord's Paradox is neither: clarification of and correction to Tu, Gunnell, and Gilthorpe (2008) by Nickerson CA & Brown NJL (https://doi.org/10.1186/1742-7622-5-2).

We commend Nickerson and Brown on their insightful exposition of the mathematical algebra behind Simpson's paradox, suppression and Lord's paradox; we also acknowledge there can be differences in how Lord's paradox is approached analytically, compared to Simpson's paradox and suppression, though not in every example of Lord's paradox. Furthermore, Simpson's paradox, suppression and Lord's paradox ask the same contextual questions, seeking to understand if statistical adjustment is valid and meaningful, identifying which analytical option is correct. In our exposition of this, we focus on the perspective of context, which must invoke causal thinking. From a causal thinking perspective, Simpson's paradox, suppression and Lord's paradox present very similar analytical challenges.

Analysing trajectories of a longitudinal exposure: A causal perspective on common methods in lifecourse research.

Longitudinal data is commonly analysed to inform prevention policies for diseases that may develop throughout life. Commonly methods interpret the longitudinal data as a series of discrete measurements or as continuous patterns. Some of the latter methods condition on the outcome, aiming to capture 'average' patterns within outcome groups, while others capture individual-level pattern features before relating these to the outcome. Conditioning on the outcome may prevent meaningful interpretation. Repeated measurements of a longitudinal exposure (weight) and later outcome (glycated haemoglobin levels) were simulated to match three scenarios: one with no causal relationship between growth rate and glycated haemoglobin; two with a positive causal effect of growth rate on glycated haemoglobin. Two methods that condition on the outcome and one that did not were applied to the data in 1000 simulations. The interpretation of the two-step method matched the simulation in all causal scenarios, but that of the methods conditioning on the outcome did not. Methods that condition on the outcome do not accurately represent a causal relationship between a longitudinal pattern and outcome. Researchers considering longitudinal data should carefully determine if they wish to analyse longitudinal data as a series of discrete time points or by extracting pattern features.

PLATOON: Premature Loss of bAby Teeth and its impact On Orthodontic Need - protocol.

OBJECTIVE: To investigate the impact of premature extraction of primary teeth (PEPT) on orthodontic treatment need in a cohort of children participating in the Born in Bradford (BiB) longitudinal birth cohort. DESIGN: Observational, cross-sectional cohort. PARTICIPANTS: We aim to recruit 1000 children aged 7-11 years: 500 with a history of PEPT and 500 matched non-PEPT controls. METHODS: After informed consent/assent, orthodontic records will be collected, including extra and intra-oral photographs and alginate impressions for study models. Participants will also complete a measure of oral health-related quality of life (COHIP-SF 19). The records will be used to quantify space loss, identify other occlusal anomalies and assess orthodontic treatment need using the Index of Orthodontic Treatment Need. For each outcome, summary statistics will be calculated and the data for children with and without PEPT compared. The records of the children identified to be in need of orthodontic treatment will be examined by an expert orthodontic panel to judge if this treatment should be undertaken at the time of the records or delayed until the early permanent dentition. Collecting robust records in the mixed dentition provides the clinical basis to link each stage of the causal chain and enable the impact of PEPT on orthodontic need to be characterised. This study is the first to provide the foundations for future longitudinal data collection allowing the long-term impact of PEPT to be studied.