A data-adaptive method for investigating effect heterogeneity with high-dimensional covariates in Mendelian randomization.

BACKGROUND: Mendelian randomization is a popular method for causal inference with observational data that uses genetic variants as instrumental variables. Similarly to a randomized trial, a standard Mendelian randomization analysis estimates the population-averaged effect of an exposure on an outcome. Dividing the population into subgroups can reveal effect heterogeneity to inform who would most benefit from intervention on the exposure. However, as covariates are measured post-"randomization", naive stratification typically induces collider bias in stratum-specific estimates. METHOD: We extend a previously proposed stratification method (the "doubly-ranked method") to form strata based on a single covariate, and introduce a data-adaptive random forest method to calculate stratum-specific estimates that are robust to collider bias based on a high-dimensional covariate set. We also propose measures based on the Q statistic to assess heterogeneity between stratum-specific estimates (to understand whether estimates are more variable than expected due to chance alone) and variable importance (to identify the key drivers of effect heterogeneity). RESULT: We show that the effect of body mass index (BMI) on lung function is heterogeneous, depending most strongly on hip circumference and weight. While for most individuals, the predicted effect of increasing BMI on lung function is negative, it is positive for some individuals and strongly negative for others. CONCLUSION: Our data-adaptive approach allows for the exploration of effect heterogeneity in the relationship between an exposure and an outcome within a Mendelian randomization framework. This can yield valuable insights into disease aetiology and help identify specific groups of individuals who would derive the greatest benefit from targeted interventions on the exposure.

Trivariate mover-stayer counting process models for investigating joint damage in psoriatic arthritis.

In psoriatic arthritis, many patients do not develop permanent joint damage even after a prolonged follow-up. This has led several authors to consider the possibility of a subpopulation of stayers (those who do not have the propensity to experience the event of interest), as opposed to assuming the entire population consist of movers (those who have the propensity to experience the event of interest). In addition, it is recognised that the damaged joints process may act very differently across different joint areas, particularly the hands, feet and large joints. From a clinical perspective, interest lies in identifying possible relationships between the damaged joints processes in these joint areas for the movers and estimating the proportion of stayers in these joint areas, if they exist. For this purpose, this paper proposes a novel trivariate mover-stayer model consisting of mover-stayer truncated negative binomial margins, and patient-level dynamic covariates and random effects in the models for the movers and stayers, respectively. The model is then extended to have a two-level mover-stayer structure for its margins so that the nature of the stayer property can be investigated. A particularly attractive feature of the proposed models is that only an optimisation routine is required in their model fitting procedures. © 2016 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.

The versatility of multi-state models for the analysis of longitudinal data with unobservable features.

Multi-state models provide a convenient statistical framework for a wide variety of medical applications characterized by multiple events and longitudinal data. We illustrate this through four examples. The potential value of the incorporation of unobserved or partially observed states is highlighted. In addition, joint modelling of multiple processes is illustrated with application to potentially informative loss to follow-up, mis-measured or missclassified data and causal inference.

Bayesian profile regression for clustering analysis involving a longitudinal response and explanatory variables.

The identification of sets of co-regulated genes that share a common function is a key question of modern genomics. Bayesian profile regression is a semi-supervised mixture modelling approach that makes use of a response to guide inference toward relevant clusterings. Previous applications of profile regression have considered univariate continuous, categorical, and count outcomes. In this work, we extend Bayesian profile regression to cases where the outcome is longitudinal (or multivariate continuous) and provide PReMiuMlongi, an updated version of PReMiuM, the R package for profile regression. We consider multivariate normal and Gaussian process regression response models and provide proof of principle applications to four simulation studies. The model is applied on budding yeast data to identify groups of genes co-regulated during the Saccharomyces cerevisiae cell cycle. We identify 4 distinct groups of genes associated with specific patterns of gene expression trajectories, along with the bound transcriptional factors, likely involved in their co-regulation process.

Dynamics of cognitive variability with age and its genetic underpinning in NIHR BioResource Genes and Cognition cohort participants.

A leading explanation for translational failure in neurodegenerative disease is that new drugs are evaluated late in the disease course when clinical features have become irreversible. Here, to address this gap, we cognitively profiled 21,051 people aged 17-85 years as part of the Genes and Cognition cohort within the National Institute for Health and Care Research BioResource across England. We describe the cohort, present cognitive trajectories and show the potential utility. Surprisingly, when studied at scale, the APOE genotype had negligible impact on cognitive performance. Different cognitive domains had distinct genetic architectures, with one indicating brain region-specific activation of microglia and another with glycogen metabolism. Thus, the molecular and cellular mechanisms underpinning cognition are distinct from dementia risk loci, presenting different targets to slow down age-related cognitive decline. Participants can now be recalled stratified by genotype and cognitive phenotype for natural history and interventional studies of neurodegenerative and other disorders.

Mixture distributions in multi-state modelling: some considerations in a study of psoriatic arthritis.

In many studies, interest lies in determining whether members of the study population will undergo a particular event of interest. Such scenarios are often termed 'mover-stayer' scenarios, and interest lies in modelling two sub-populations of 'movers' (those who have a propensity to undergo the event of interest) and 'stayers' (those who do not). In general, mover-stayer scenarios within data sets are accounted for through the use of mixture distributions, and in this paper, we investigate the use of various random effects distributions for this purpose. Using data from the University of Toronto psoriatic arthritis clinic, we present a multi-state model to describe the progression of clinical damage in hand joints of patients with psoriatic arthritis. We consider the use of mover-stayer gamma, inverse Gaussian and compound Poisson distributions to account for both the correlation amongst joint locations and the possible mover-stayer situation with regard to clinical hand joint damage. We compare the fits obtained from these models and discuss the extent to which a mover-stayer scenario exists in these data. Furthermore, we fit a mover-stayer model that allows a dependence of the probability of a patient being a stayer on a patient-level explanatory variable.

A randomized placebo-controlled trial of methotrexate in psoriatic arthritis.

OBJECTIVE: MTX is widely used to treat synovitis in PsA without supporting trial evidence. The aim of our study was to test the value of MTX in the first large randomized placebo-controlled trial (RCT) in PsA. METHODS: A 6-month double-blind RCT compared MTX (15 mg/week) with placebo in active PsA. The primary outcome was PsA response criteria (PsARC). Other outcomes included ACR20, DAS-28 and their individual components. Missing data were imputed using multiple imputation methods. Treatments were compared using logistic regression analysis (adjusted for age, sex, disease duration and, where appropriate, individual baseline scores). RESULTS: Four hundred and sixty-two patients were screened and 221 recruited. One hundred and nine patients received MTX and 112 received placebo. Forty-four patients were lost to follow-up (21 MTX, 23 placebo). Twenty-six patients discontinued treatment (14 MTX, 12 placebo). Comparing MTX with placebo in all randomized patients at 6 months showed no significant effect on PsARC [odds ratio (OR) 1.77, 95% CI 0.97, 3.23], ACR20 (OR 2.00, 95% CI 0.65, 6.22) or DAS-28 (OR 1.70, 95% CI 0.90, 3.17). There were also no significant treatment effects on tender and swollen joint counts, ESR, CRP, HAQ and pain. The only benefits of MTX were reductions in patient and assessor global scores and skin scores at 6 months (P = 0.03, P < 0.001 and P = 0.02, respectively). There were no unexpected adverse events. CONCLUSIONS: This trial of active PsA found no evidence for MTX improving synovitis and consequently raises questions about its classification as a disease-modifying drug in PsA. Trial registration. Current Controlled Trials, www.controlled-trials.com, ISRCTN:54376151.

Characterization of missing data patterns and mechanisms in longitudinal composite outcome trial in rheumatoid arthritis.

Background: Composite measures, like the Disease Activity Score for 28 joints (DAS28), are key primary outcomes in rheumatoid arthritis (RA) trials. DAS28 combines four different components in a continuous measure. When one or more of these components are missing the overall composite score is also missing at intermediate or trial endpoint assessments. Objectives: This study examined missing data patterns and mechanisms in a longitudinal RA trial to evaluate how best to handle missingness when analysing composite outcomes. Design: The Tumour-Necrosis-Factor Inhibitors against Combination Intensive Therapy (TACIT) trial was an open label, pragmatic randomized multicentre two arm non-inferiority study. Patients were followed up for 12 months, with monthly measurement of the composite outcome and its components. Active RA patients were randomized to conventional disease modifying drugs (cDMARDs) or Tumour Necrosis Factor-α inhibitors (TNFis). Methods: The TACIT trial was used to explore the extent of missing data in the composite outcome, DAS28. Patterns of missing data in components and the composite outcome were examined graphically. Longitudinal multivariable logistic regression analysis assessed missing data mechanisms during follow-up. Results: Two hundred and five patients were randomized: at 12 months 59/205 (29%) had unobserved composite outcome and 146/205 (71%) had an observed DAS28 outcome; however, 34/146 had one or more intermediate assessments missing. We observed mixed missing data patterns, especially for the missing composite outcome due to one component missing rather than patient not attending thier visit. Age and gender predicted missingness components, providing strong evidence the missing observations were unlikely to be Missing Completely at Random (MCAR). Conclusion: Researchers should undertake detailed evaluations of missing data patterns and mechanisms at the final and intermediate time points, whether or not the outcome variable is a composite outcome. In addition, the impact on treatment estimates in patients who only provide data at milestone assessments need to be assessed. Trial Registration ISRCTN Number: 37438295.

Development and validation of a dynamic 48-hour in-hospital mortality risk stratification for COVID-19 in a UK teaching hospital: a retrospective cohort study.

OBJECTIVES: To develop a disease stratification model for COVID-19 that updates according to changes in a patient's condition while in hospital to facilitate patient management and resource allocation. DESIGN: In this retrospective cohort study, we adopted a landmarking approach to dynamic prediction of all-cause in-hospital mortality over the next 48 hours. We accounted for informative predictor missingness and selected predictors using penalised regression. SETTING: All data used in this study were obtained from a single UK teaching hospital. PARTICIPANTS: We developed the model using 473 consecutive patients with COVID-19 presenting to a UK hospital between 1 March 2020 and 12 September 2020; and temporally validated using data on 1119 patients presenting between 13 September 2020 and 17 March 2021. PRIMARY AND SECONDARY OUTCOME MEASURES: The primary outcome is all-cause in-hospital mortality within 48 hours of the prediction time. We accounted for the competing risks of discharge from hospital alive and transfer to a tertiary intensive care unit for extracorporeal membrane oxygenation. RESULTS: Our final model includes age, Clinical Frailty Scale score, heart rate, respiratory rate, oxygen saturation/fractional inspired oxygen ratio, white cell count, presence of acidosis (pH <7.35) and interleukin-6. Internal validation achieved an area under the receiver operating characteristic (AUROC) of 0.90 (95% CI 0.87 to 0.93) and temporal validation gave an AUROC of 0.86 (95% CI 0.83 to 0.88). CONCLUSIONS: Our model incorporates both static risk factors (eg, age) and evolving clinical and laboratory data, to provide a dynamic risk prediction model that adapts to both sudden and gradual changes in an individual patient's clinical condition. On successful external validation, the model has the potential to be a powerful clinical risk assessment tool. TRIAL REGISTRATION: The study is registered as 'researchregistry5464' on the Research Registry (www.researchregistry.com).

Intermittent observation of time-dependent explanatory variables: a multistate modelling approach.

Motivated by investigations of factors related to various patient-reported outcome measures in psoriatic arthritis patients, after controlling for the effect of disease activity on these outcomes, we outline an approach for dealing with a rapidly fluctuating explanatory variable in a multistate model. On the basis of a representation of this variable as an ordinal classification, we suggest the use of an expanded multistate model. We examine the bias in estimating effects associated with other variables via simulation for different modelling choices. We present an analysis of a motivating data set on physical functional disability in psoriatic arthritis patients.

Disease Modelling of Cognitive Outcomes and Biomarkers in the European Prevention of Alzheimer's Dementia Longitudinal Cohort.

A key challenge for the secondary prevention of Alzheimer's dementia is the need to identify individuals early on in the disease process through sensitive cognitive tests and biomarkers. The European Prevention of Alzheimer's Dementia (EPAD) consortium recruited participants into a longitudinal cohort study with the aim of building a readiness cohort for a proof-of-concept clinical trial and also to generate a rich longitudinal data-set for disease modelling. Data have been collected on a wide range of measurements including cognitive outcomes, neuroimaging, cerebrospinal fluid biomarkers, genetics and other clinical and environmental risk factors, and are available for 1,828 eligible participants at baseline, 1,567 at 6 months, 1,188 at one-year follow-up, 383 at 2 years, and 89 participants at three-year follow-up visit. We novelly apply state-of-the-art longitudinal modelling and risk stratification approaches to these data in order to characterise disease progression and biological heterogeneity within the cohort. Specifically, we use longitudinal class-specific mixed effects models to characterise the different clinical disease trajectories and a semi-supervised Bayesian clustering approach to explore whether participants can be stratified into homogeneous subgroups that have different patterns of cognitive functioning evolution, while also having subgroup-specific profiles in terms of baseline biomarkers and longitudinal rate of change in biomarkers.

Baseline predictors of remission, pain and fatigue in rheumatoid arthritis: the TITRATE trial.

BACKGROUND: Clinical trials show intensive treatment to induce remission is effective in patients with highly active rheumatoid arthritis (RA). The TITRATE trial showed that the benefits of intensive treatment also extend to moderately active RA. However, many patients failed to achieve remission or show improvements in pain and fatigue. We investigated whether baseline predictors could identify treatment non-responders. METHODS: The impact of obesity, depression, anxiety and illness perception on RA outcomes, including disease activity, remission, pain and fatigue were determined using a pre-planned secondary analysis of the TITRATE trial data. RESULTS: Body mass index was associated with disease activity levels and remission: obese patients had a higher overall disease activity and fewer obese patients achieved remission. Intensive management was not associated with increased remission in these patients. Obesity was also associated with increased overall pain and fatigue. Anxiety, depression and health perceptions had no discernible impact on disease activity but were associated with high levels of pain and fatigue. There was a strong association between anxiety and high pain scores; and between depression and high fatigue scores; and health perception was strongly related to both. None of the predictors had an important impact on pain and fatigue reduction in cross-sectional analysis. CONCLUSIONS: Disease activity is higher in obese patients and they have fewer remissions over 12 months. Anxiety, depression and health perceptions were associated with higher pain and fatigue scores. Intensive management strategies need to account for these baseline features as they impact significantly on clinical and psychological outcomes. TRIAL REGISTRATION: ISRCTN 70160382 ; date registered 16 January 2014.

Bias in 2-part mixed models for longitudinal semicontinuous data.

Semicontinuous data in the form of a mixture of zeros and continuously distributed positive values frequently arise in biomedical research. Two-part mixed models with correlated random effects are an attractive approach to characterize the complex structure of longitudinal semicontinuous data. In practice, however, an independence assumption about random effects in these models may often be made for convenience and computational feasibility. In this article, we show that bias can be induced for regression coefficients when random effects are truly correlated but misspecified as independent in a 2-part mixed model. Paralleling work on bias under nonignorable missingness within a shared parameter model, we derive and investigate the asymptotic bias in selected settings for misspecified 2-part mixed models. The performance of these models in practice is further evaluated using Monte Carlo simulations. Additionally, the potential bias is investigated when artificial zeros, due to left censoring from some detection or measuring limit, are incorporated. To illustrate, we fit different 2-part mixed models to the data from the University of Toronto Psoriatic Arthritis Clinic, the aim being to examine whether there are differential effects of disease activity and damage on physical functioning as measured by the health assessment questionnaire scores over the course of psoriatic arthritis. Some practical issues on variance component estimation revealed through this data analysis are considered.

Efficient Real-Time Monitoring of an Emerging Influenza Pandemic: How Feasible?

A prompt public health response to a new epidemic relies on the ability to monitor and predict its evolution in real time as data accumulate. The 2009 A/H1N1 outbreak in the UK revealed pandemic data as noisy, contaminated, potentially biased and originating from multiple sources. This seriously challenges the capacity for real-time monitoring. Here, we assess the feasibility of real-time inference based on such data by constructing an analytic tool combining an age-stratified SEIR transmission model with various observation models describing the data generation mechanisms. As batches of data become available, a sequential Monte Carlo (SMC) algorithm is developed to synthesise multiple imperfect data streams, iterate epidemic inferences and assess model adequacy amidst a rapidly evolving epidemic environment, substantially reducing computation time in comparison to standard MCMC, to ensure timely delivery of real-time epidemic assessments. In application to simulated data designed to mimic the 2009 A/H1N1 epidemic, SMC is shown to have additional benefits in terms of assessing predictive performance and coping with parameter nonidentifiability.

Prescreening for European Prevention of Alzheimer Dementia (EPAD) trial-ready cohort: impact of AD risk factors and recruitment settings.

BACKGROUND: Recruitment is often a bottleneck in secondary prevention trials in Alzheimer disease (AD). Furthermore, screen-failure rates in these trials are typically high due to relatively low prevalence of AD pathology in individuals without dementia, especially among cognitively unimpaired. Prescreening on AD risk factors may facilitate recruitment, but the efficiency will depend on how these factors link to participation rates and AD pathology. We investigated whether common AD-related factors predict trial-ready cohort participation and amyloid status across different prescreen settings. METHODS: We monitored the prescreening in four cohorts linked to the European Prevention of Alzheimer Dementia (EPAD) Registry (n = 16,877; mean ± SD age = 64 ± 8 years). These included a clinical cohort, a research in-person cohort, a research online cohort, and a population-based cohort. Individuals were asked to participate in the EPAD longitudinal cohort study (EPAD-LCS), which serves as a trial-ready cohort for secondary prevention trials. Amyloid positivity was measured in cerebrospinal fluid as part of the EPAD-LCS assessment. We calculated participation rates and numbers needed to prescreen (NNPS) per participant that was amyloid-positive. We tested if age, sex, education level, APOE status, family history for dementia, memory complaints or memory scores, previously collected in these cohorts, could predict participation and amyloid status. RESULTS: A total of 2595 participants were contacted for participation in the EPAD-LCS. Participation rates varied by setting between 3 and 59%. The NNPS were 6.9 (clinical cohort), 7.5 (research in-person cohort), 8.4 (research online cohort), and 88.5 (population-based cohort). Participation in the EPAD-LCS (n = 413 (16%)) was associated with lower age (odds ratio (OR) age = 0.97 [0.95-0.99]), high education (OR = 1.64 [1.23-2.17]), male sex (OR = 1.56 [1.19-2.04]), and positive family history of dementia (OR = 1.66 [1.19-2.31]). Among participants in the EPAD-LCS, amyloid positivity (33%) was associated with higher age (OR = 1.06 [1.02-1.10]) and APOE ɛ4 allele carriership (OR = 2.99 [1.81-4.94]). These results were similar across prescreen settings. CONCLUSIONS: Numbers needed to prescreen varied greatly between settings. Understanding how common AD risk factors link to study participation and amyloid positivity is informative for recruitment strategy of studies on secondary prevention of AD.

Clustered multistate models with observation level random effects, mover-stayer effects and dynamic covariates: modelling transition intensities and sojourn times in a study of psoriatic arthritis.

In psoriatic arthritis, it is important to understand the joint activity (represented by swelling and pain) and damage processes because both are related to severe physical disability. The paper aims to provide a comprehensive investigation into both processes occurring over time, in particular their relationship, by specifying a joint multistate model at the individual hand joint level, which also accounts for many of their important features. As there are multiple hand joints, such an analysis will be based on the use of clustered multistate models. Here we consider an observation level random-effects structure with dynamic covariates and allow for the possibility that a subpopulation of patients is at minimal risk of damage. Such an analysis is found to provide further understanding of the activity-damage relationship beyond that provided by previous analyses. Consideration is also given to the modelling of mean sojourn times and jump probabilities. In particular, a novel model parameterization which allows easily interpretable covariate effects to act on these quantities is proposed.

Event-biased referral can distort estimation of hepatitis C virus progression rate to cirrhosis, and of prognostic influences.

OBJECTIVE: To investigate how cirrhosis-biased referral to liver clinics can explain the wide variation in progression rates for differently recruited cohorts and, in particular, for liver clinic cohorts compared to community-based studies of the natural history of hepatitis C virus (HCV). STUDY DESIGN AND SETTING: A simulation was designed to illustrate the sort of referral bias pattern that is capable of converting a 20-year progression rate to cirrhosis of around 5% in the community of HCV-infected individuals into a 20% progression rate for patients who have been selectively referred to a liver clinic. RESULTS: We show that event-biased recruitment, such as occurs if referral to liver clinics is increasingly likely the closer a patient is to cirrhosis, can produce severely upwardly biased estimates of progression rates, can dampen the influence of "poor prognostic" factors (such as history of excessive alcohol consumption), but overrepresents the proportion of patients in the community of HCV-infected individuals who have poor prognosis. CONCLUSION: When attempting to establish the natural history of new diseases with long incubation periods, researchers should be on the look out for potential biases that result from the way patients are referred into clinical cohorts.

Power analysis to detect treatment effects in longitudinal clinical trials for Alzheimer's disease.

INTRODUCTION: Assessing cognitive and functional changes at the early stage of Alzheimer's disease (AD) and detecting treatment effects in clinical trials for early AD are challenging. METHODS: Under the assumption that transformed versions of the Mini-Mental State Examination, the Clinical Dementia Rating Scale-Sum of Boxes, and the Alzheimer's Disease Assessment Scale-Cognitive Subscale tests'/components' scores are from a multivariate linear mixed-effects model, we calculated the sample sizes required to detect treatment effects on the annual rates of change in these three components in clinical trials for participants with mild cognitive impairment. RESULTS: Our results suggest that a large number of participants would be required to detect a clinically meaningful treatment effect in a population with preclinical or prodromal Alzheimer's disease. We found that the transformed Mini-Mental State Examination is more sensitive for detecting treatment effects in early AD than the transformed Clinical Dementia Rating Scale-Sum of Boxes and Alzheimer's Disease Assessment Scale-Cognitive Subscale. The use of optimal weights to construct powerful test statistics or sensitive composite scores/endpoints can reduce the required sample sizes needed for clinical trials. CONCLUSION: Consideration of the multivariate/joint distribution of components' scores rather than the distribution of a single composite score when designing clinical trials can lead to an increase in power and reduced sample sizes for detecting treatment effects in clinical trials for early AD.

Exploring the existence of a stayer population with mover-stayer counting process models: application to joint damage in psoriatic arthritis.

Many psoriatic arthritis patients do not progress to permanent joint damage in any of the 28 hand joints, even under prolonged follow-up. This has led several researchers to fit models that estimate the proportion of stayers (those who do not have the propensity to experience the event of interest) and to characterize the rate of developing damaged joints in the movers (those who have the propensity to experience the event of interest). However, when fitted to the same data, the paper demonstrates that the choice of model for the movers can lead to widely varying conclusions on a stayer population, thus implying that, if interest lies in a stayer population, a single analysis should not generally be adopted. The aim of the paper is to provide greater understanding regarding estimation of a stayer population by comparing the inferences, performance and features of multiple fitted models to real and simulated data sets. The models for the movers are based on Poisson processes with patient level random effects and/or dynamic covariates, which are used to induce within-patient correlation, and observation level random effects are used to account for time varying unobserved heterogeneity. The gamma, inverse Gaussian and compound Poisson distributions are considered for the random effects.