A semiparametric bivariate probit model for joint modeling of outcomes in STEMI patients

In this work we analyse the relationship among in-hospital mortality and a treatment effectiveness outcome in patients affected by ST-Elevation myocardial infarction. The main idea is to carry out a joint modeling of the two outcomes applying a Semiparametric Bivariate Probit Model to data arising from a clinical registry called STEMI Archive. A realistic quantification of the relationship between outcomes can be problematic for several reasons. First, latent factors associated with hospitals organization can affect the treatment efficacy and/or interact with patient’s condition at admission time. Moreover, they can also directly influence the mortality outcome. Such factors can be hardly measurable. Thus, the use of classical estimation methods will clearly result in inconsistent or biased parameter estimates. Secondly, covariate-outcomes relationships can exhibit nonlinear patterns. Provided that proper statistical methods for model fitting in such framework are available, it is possible to employ a simultaneous estimation approach to account for unobservable confounders. Such a framework can also provide flexible covariate structures and model the whole conditional distribution of the response

The computational neurology of active vision

In this thesis, we appeal to recent developments in theoretical neurobiology – namely, active inference – to understand the active visual system and its disorders. Chapter 1 reviews the neurobiology of active vision. This introduces some of the key conceptual themes around attention and inference that recur through subsequent chapters. Chapter 2 provides a technical overview of active inference, and its interpretation in terms of message passing between populations of neurons. Chapter 3 applies the material in Chapter 2 to provide a computational characterisation of the oculomotor system. This deals with two key challenges in active vision: deciding where to look, and working out how to look there. The homology between this message passing and the brain networks solving these inference problems provide a basis for in silico lesion experiments, and an account of the aberrant neural computations that give rise to clinical oculomotor signs (including internuclear ophthalmoplegia). Chapter 4 picks up on the role of uncertainty resolution in deciding where to look, and examines the role of beliefs about the quality (or precision) of data in perceptual inference. We illustrate how abnormal prior beliefs influence inferences about uncertainty and give rise to neuromodulatory changes and visual hallucinatory phenomena (of the sort associated with synucleinopathies). We then demonstrate how synthetic pharmacological perturbations that alter these neuromodulatory systems give rise to the oculomotor changes associated with drugs acting upon these systems. Chapter 5 develops a model of visual neglect, using an oculomotor version of a line cancellation task. We then test a prediction of this model using magnetoencephalography and dynamic causal modelling. Chapter 6 concludes by situating the work in this thesis in the context of computational neurology. This illustrates how the variational principles used here to characterise the active visual system may be generalised to other sensorimotor systems and their disorders

Computational Neuropsychology and Bayesian Inference

Computational theories of brain function have become very influential in neuroscience. They have facilitated the growth of formal approaches to disease, particularly in psychiatric research. In this paper, we provide a narrative review of the body of computational research addressing neuropsychological syndromes, and focus on those that employ Bayesian frameworks. Bayesian approaches to understanding brain function formulate perception and action as inferential processes. These inferences combine ‘prior’ beliefs with a generative (predictive) model to explain the causes of sensations. Under this view, neuropsychological deficits can be thought of as false inferences that arise due to aberrant prior beliefs (that are poor fits to the real world). This draws upon the notion of a Bayes optimal pathology – optimal inference with suboptimal priors – and provides a means for computational phenotyping. In principle, any given neuropsychological disorder could be characterized by the set of prior beliefs that would make a patient’s behavior appear Bayes optimal. We start with an overview of some key theoretical constructs and use these to motivate a form of computational neuropsychology that relates anatomical structures in the brain to the computations they perform. Throughout, we draw upon computational accounts of neuropsychological syndromes. These are selected to emphasize the key features of a Bayesian approach, and the possible types of pathological prior that may be present. They range from visual neglect through hallucinations to autism. Through these illustrative examples, we review the use of Bayesian approaches to understand the link between biology and computation that is at the heart of neuropsychology

Causal hierarchy within the thalamo-cortical network in spike and wave discharges

Background: Generalised spike wave (GSW) discharges are the electroencephalographic (EEG) hallmark of absence seizures, clinically characterised by a transitory interruption of ongoing activities and impaired consciousness, occurring during states of reduced awareness. Several theories have been proposed to explain the pathophysiology of GSW discharges and the role of thalamus and cortex as generators. In this work we extend the existing theories by hypothesizing a role for the precuneus, a brain region neglected in previous works on GSW generation but already known to be linked to consciousness and awareness. We analysed fMRI data using dynamic causal modelling (DCM) to investigate the effective connectivity between precuneus, thalamus and prefrontal cortex in patients with GSW discharges. Methodology and Principal Findings: We analysed fMRI data from seven patients affected by Idiopathic Generalized Epilepsy (IGE) with frequent GSW discharges and significant GSW-correlated haemodynamic signal changes in the thalamus, the prefrontal cortex and the precuneus. Using DCM we assessed their effective connectivity, i.e. which region drives another region. Three dynamic causal models were constructed: GSW was modelled as autonomous input to the thalamus (model A), ventromedial prefrontal cortex (model B), and precuneus (model C). Bayesian model comparison revealed Model C (GSW as autonomous input to precuneus), to be the best in 5 patients while model A prevailed in two cases. At the group level model C dominated and at the population-level the p value of model C was ∼1. Conclusion: Our results provide strong evidence that activity in the precuneus gates GSW discharges in the thalamo-(fronto) cortical network. This study is the first demonstration of a causal link between haemodynamic changes in the precuneus - an index of awareness - and the occurrence of pathological discharges in epilepsy. © 2009 Vaudano et al

Development, Calibration and Performance of an HIV Transmission Model Incorporating Natural History and Behavioral Patterns: Application in South Africa

Understanding HIV transmission dynamics is critical to estimating the potential population-wide impact of HIV prevention and treatment interventions. We developed an individual-based simulation model of the heterosexual HIV epidemic in South Africa and linked it to the previously published Cost-Effectiveness of Preventing AIDS Complications (CEPAC) International Model, which simulates the natural history and treatment of HIV. In this new model, the CEPAC Dynamic Model (CDM), the probability of HIV transmission per sexual encounter between short-term, long-term and commercial sex worker partners depends upon the HIV RNA and disease stage of the infected partner, condom use, and the circumcision status of the uninfected male partner. We included behavioral, demographic and biological values in the CDM and calibrated to HIV prevalence in South Africa pre-antiretroviral therapy. Using a multi-step fitting procedure based on Bayesian melding methodology, we performed 264,225 simulations of the HIV epidemic in South Africa and identified 3,750 parameter sets that created an epidemic and had behavioral characteristics representative of a South African population pre-ART. Of these parameter sets, 564 contributed 90% of the likelihood weight to the fit, and closely reproduced the UNAIDS HIV prevalence curve in South Africa from 1990–2002. The calibration was sensitive to changes in the rate of formation of short-duration partnerships and to the partnership acquisition rate among high-risk individuals, both of which impacted concurrency. Runs that closely fit to historical HIV prevalence reflect diverse ranges for individual parameter values and predict a wide range of possible steady-state prevalence in the absence of interventions, illustrating the value of the calibration procedure and utility of the model for evaluating interventions. This model, which includes detailed behavioral patterns and HIV natural history, closely fits HIV prevalence estimates

Avoiding the misuse of BLUP in behavioral ecology

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record.Having recognized that variation around the population-level “Golden Mean” of labile traits contains biologically meaningful information, behavioural ecologists have focused increasingly on exploring the causes and consequences of individual variation in behaviour. These are exciting new directions for the field, assisted in no small part by the adoption of mixed-effects modelling techniques that enable the partitioning of among- and within-individual behavioural variation. It has become commonplace to extract predictions of individual random effects from such models for use in subsequent analyses (for example, between a personality trait and other individual traits such as cognition, physiology, or fitness-related measures). However, these predictions are made with large amounts of error that is not carried forward, rendering further tests susceptible to spurious P values from these individual-level point estimates. We briefly summarize the problems with such statistical methods that are used regularly by behavioural ecologists, and highlight the robust solutions that exist within the mixed model framework, providing tutorials to aid in their implementation.This work was supported by a Biotechnology and Biological Sciences Research Council grant (BB/L022656/1

GABAergic cortical network physiology in frontotemporal lobar degeneration.

The clinical syndromes caused by frontotemporal lobar degeneration are heterogeneous, including the behavioural variant frontotemporal dementia (bvFTD) and progressive supranuclear palsy. Although pathologically distinct, they share many behavioural, cognitive and physiological features, which may in part arise from common deficits of major neurotransmitters such as γ-aminobutyric acid (GABA). Here, we quantify the GABAergic impairment and its restoration with dynamic causal modelling of a double-blind placebo-controlled crossover pharmaco-magnetoencephalography study. We analysed 17 patients with bvFTD, 15 patients with progressive supranuclear palsy, and 20 healthy age- and gender-matched controls. In addition to neuropsychological assessment and structural MRI, participants undertook two magnetoencephalography sessions using a roving auditory oddball paradigm: once on placebo and once on 10 mg of the oral GABA reuptake inhibitor tiagabine. A subgroup underwent ultrahigh-field magnetic resonance spectroscopy measurement of GABA concentration, which was reduced among patients. We identified deficits in frontotemporal processing using conductance-based biophysical models of local and global neuronal networks. The clinical relevance of this physiological deficit is indicated by the correlation between top-down connectivity from frontal to temporal cortex and clinical measures of cognitive and behavioural change. A critical validation of the biophysical modelling approach was evidence from parametric empirical Bayes analysis that GABA levels in patients, measured by spectroscopy, were related to posterior estimates of patients' GABAergic synaptic connectivity. Further evidence for the role of GABA in frontotemporal lobar degeneration came from confirmation that the effects of tiagabine on local circuits depended not only on participant group, but also on individual baseline GABA levels. Specifically, the phasic inhibition of deep cortico-cortical pyramidal neurons following tiagabine, but not placebo, was a function of GABA concentration. The study provides proof-of-concept for the potential of dynamic causal modelling to elucidate mechanisms of human neurodegenerative disease, and explains the variation in response to candidate therapies among patients. The laminar- and neurotransmitter-specific features of the modelling framework, can be used to study other treatment approaches and disorders. In the context of frontotemporal lobar degeneration, we suggest that neurophysiological restoration in selected patients, by targeting neurotransmitter deficits, could be used to bridge between clinical and preclinical models of disease, and inform the personalized selection of drugs and stratification of patients for future clinical trials

GABAergic cortical network physiology in frontotemporal lobar degeneration.

The clinical syndromes caused by frontotemporal lobar degeneration are heterogeneous, including the behavioural variant frontotemporal dementia (bvFTD) and progressive supranuclear palsy. Although pathologically distinct, they share many behavioural, cognitive and physiological features, which may in part arise from common deficits of major neurotransmitters such as γ-aminobutyric acid (GABA). Here, we quantify the GABAergic impairment and its restoration with dynamic causal modelling of a double-blind placebo-controlled crossover pharmaco-magnetoencephalography study. We analysed 17 patients with bvFTD, 15 patients with progressive supranuclear palsy, and 20 healthy age- and gender-matched controls. In addition to neuropsychological assessment and structural MRI, participants undertook two magnetoencephalography sessions using a roving auditory oddball paradigm: once on placebo and once on 10 mg of the oral GABA reuptake inhibitor tiagabine. A subgroup underwent ultrahigh-field magnetic resonance spectroscopy measurement of GABA concentration, which was reduced among patients. We identified deficits in frontotemporal processing using conductance-based biophysical models of local and global neuronal networks. The clinical relevance of this physiological deficit is indicated by the correlation between top-down connectivity from frontal to temporal cortex and clinical measures of cognitive and behavioural change. A critical validation of the biophysical modelling approach was evidence from parametric empirical Bayes analysis that GABA levels in patients, measured by spectroscopy, were related to posterior estimates of patients' GABAergic synaptic connectivity. Further evidence for the role of GABA in frontotemporal lobar degeneration came from confirmation that the effects of tiagabine on local circuits depended not only on participant group, but also on individual baseline GABA levels. Specifically, the phasic inhibition of deep cortico-cortical pyramidal neurons following tiagabine, but not placebo, was a function of GABA concentration. The study provides proof-of-concept for the potential of dynamic causal modelling to elucidate mechanisms of human neurodegenerative disease, and explains the variation in response to candidate therapies among patients. The laminar- and neurotransmitter-specific features of the modelling framework, can be used to study other treatment approaches and disorders. In the context of frontotemporal lobar degeneration, we suggest that neurophysiological restoration in selected patients, by targeting neurotransmitter deficits, could be used to bridge between clinical and preclinical models of disease, and inform the personalized selection of drugs and stratification of patients for future clinical trials