A Bayesian alternative to mutual information for the hierarchical clustering of dependent random variables

The use of mutual information as a similarity measure in agglomerative hierarchical clustering (AHC) raises an important issue: some correction needs to be applied for the dimensionality of variables. In this work, we formulate the decision of merging dependent multivariate normal variables in an AHC procedure as a Bayesian model comparison. We found that the Bayesian formulation naturally shrinks the empirical covariance matrix towards a matrix set a priori (e.g., the identity), provides an automated stopping rule, and corrects for dimensionality using a term that scales up the measure as a function of the dimensionality of the variables. Also, the resulting log Bayes factor is asymptotically proportional to the plug-in estimate of mutual information, with an additive correction for dimensionality in agreement with the Bayesian information criterion. We investigated the behavior of these Bayesian alternatives (in exact and asymptotic forms) to mutual information on simulated and real data. An encouraging result was first derived on simulations: the hierarchical clustering based on the log Bayes factor outperformed off-the-shelf clustering techniques as well as raw and normalized mutual information in terms of classification accuracy. On a toy example, we found that the Bayesian approaches led to results that were similar to those of mutual information clustering techniques, with the advantage of an automated thresholding. On real functional magnetic resonance imaging (fMRI) datasets measuring brain activity, it identified clusters consistent with the established outcome of standard procedures. On this application, normalized mutual information had a highly atypical behavior, in the sense that it systematically favored very large clusters. These initial experiments suggest that the proposed Bayesian alternatives to mutual information are a useful new tool for hierarchical clustering

fMRI activation detection with EEG priors

The purpose of brain mapping techniques is to advance the understanding of the relationship between structure and function in the human brain in so-called activation studies. In this work, an advanced statistical model for combining functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) recordings is developed to fuse complementary information about the location of neuronal activity. More precisely, a new Bayesian method is proposed for enhancing fMRI activation detection by the use of EEG-based spatial prior information in stimulus based experimental paradigms. I.e., we model and analyse stimulus influence by a spatial Bayesian variable selection scheme, and extend existing high-dimensional regression methods by incorporating prior information on binary selection indicators via a latent probit regression with either a spatially-varying or constant EEG effect. Spatially-varying effects are regularized by intrinsic Markov random field priors. Inference is based on a full Bayesian Markov Chain Monte Carlo (MCMC) approach. Whether the proposed algorithm is able to increase the sensitivity of mere fMRI models is examined in both a real-world application and a simulation study. We observed, that carefully selected EEG--prior information additionally increases sensitivity in activation regions that have been distorted by a low signal-to-noise ratio

Non-standard templates for non-standard populations: optimizing template selection for voxel-based morphometry pre-processing

The human brain is a complex and powerful organ, directing every aspect of life from somatosensory and motor function to visceral responses to higher order cognition. Neurological and psychiatric disorders often disrupt normal functioning. While the clinical symptoms of such disorders are known, their biological underpinnings are not as clearly characterized. Structural neuroimaging is a powerful, non-invasive tool that can play a critical role in finding biomarkers of these illnesses. Currently, variations in pre-processing techniques yield inconsistent and conflicting results. As neuroimaging is a nascent branch of medical research, gold standards in imaging methodologies have not yet been established. Quantitatively validating and optimizing the way these images are preprocessed is the first step towards standardization. Voxel-based morphometry (VBM) is one technique that is commonly used to compare whole-brain structural differences between groups. Statistical tests are used to compare intensities of voxels throughout all brain scans in each group. In order to ensure that comparable voxels are being tested, the images must be fitted into a common space, which is done through image preprocessing. Spatial normalization to templates is an early pre-processing step that is executed unreliably as many options for both templates and normalization algorithms exist. To determine the effect variations in template usage may cause, we utilized a VBM approach to detect simulated lesions. Template performance was analyzed by comparing the accuracy with which the lesion was detected

Graph analysis of functional brain networks: practical issues in translational neuroscience

The brain can be regarded as a network: a connected system where nodes, or units, represent different specialized regions and links, or connections, represent communication pathways. From a functional perspective communication is coded by temporal dependence between the activities of different brain areas. In the last decade, the abstract representation of the brain as a graph has allowed to visualize functional brain networks and describe their non-trivial topological properties in a compact and objective way. Nowadays, the use of graph analysis in translational neuroscience has become essential to quantify brain dysfunctions in terms of aberrant reconfiguration of functional brain networks. Despite its evident impact, graph analysis of functional brain networks is not a simple toolbox that can be blindly applied to brain signals. On the one hand, it requires a know-how of all the methodological steps of the processing pipeline that manipulates the input brain signals and extract the functional network properties. On the other hand, a knowledge of the neural phenomenon under study is required to perform physiological-relevant analysis. The aim of this review is to provide practical indications to make sense of brain network analysis and contrast counterproductive attitudes

Transcranial Focused Ultrasound to the Right Prefrontal Cortex Improves Mood and Alters Functional Connectivity in Humans

Transcranial focused ultrasound (tFUS) is an emerging method for non-invasive neuromodulation akin to transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS). tFUS offers several advantages over electromagnetic methods including high spatial resolution and the ability to reach deep brain targets. Here we describe two experiments assessing whether tFUS could modulate mood in healthy human volunteers by targeting the right inferior frontal gyrus (rIFG), an area implicated in mood and emotional regulation. In a randomized, placebo-controlled, double-blind study, participants received 30 s of 500 kHz tFUS or a placebo control. Visual Analog Mood Scales (VAMS) assessed mood four times within an hour (baseline and three times after tFUS). Participants who received tFUS reported an overall increase in Global Affect (GA), an aggregate score from the VAMS scale, indicating a positive shift in mood. Experiment 2 examined resting-state functional (FC) connectivity using functional magnetic resonance imaging (fMRI) following 2 min of 500 kHz tFUS at the rIFG. As in Experiment 1, tFUS enhanced self-reported mood states and also decreased FC in resting state networks related to emotion and mood regulation. These results suggest that tFUS can be used to modulate mood and emotional regulation networks in the prefrontal cortex