Bayesian Methods in Functional Magnetic Resonance Imaging
|
|
- Jade Melton
- 5 years ago
- Views:
Transcription
1 Bayesian Methods in Functional Magnetic Resonance Imaging Galin L. Jones Kuo-Jung Lee Brian S. Caffo Susan Spear Bassett Abstract: One of the major objectives of functional magnetic resonance imaging studies is to infer regional neuronal activity in response to a stimulus or task. We consider a Bayesian approach which incorporates spatial and temporal dependence and allows for the task-related response to change dynamically over the scanning session. AMS 2000 subject classifications: Primary 62-02; secondary 62P10. Keywords and phrases: BOLD response, functional MRI, neuroimaging, Variable selection. 1. Introduction Most functional neuroimaging experiments aim to either uncover localized regions where the brain activates during a task or describe the networks required for a particular brain function. For example, is there a common area of the brain which is active when presented with auditory stimuli? Or are there remote regions of the brain which interact to accomplish this task? Our focus is on functional magnetic resonance imaging (fmri) techniques to study brain activation. Neuronal activation occurs in milliseconds and is not observed directly in fmri experiments. However, the well-established principle of neurovascular coupling (that is, local neuronal activation is related to changes in cerebral blood flow) allows us to observe this phenomenon via the blood oxygenation level dependent (BOLD) signal contrast. The structure of a typical fmri experiment is straightforward. A subject in an MRI scanner performs a task in response to a stimulus. Three-dimensional images of the subject s brain are captured every 2-3 seconds. The experiment may be repeated on the same subject and is often performed on several subjects. The subjects may be divided into groups. Even if we limit attention to a single subject in a single session, the data collected from an fmri experiment has a complicated structure. Imagine that the subject s brain is divided into a regular grid of 100, ,000 volume elements, or voxels. The BOLD signal is observed at each voxel at each of say time points. Thus there is an enormous amount of data which exhibits both spatial and temporal dependence. Moreover, the data tend to be noisy and have a weak signal. There are few, if any, off-the-shelf methods for constructing 1
2 Jones et al./bayes for fmri 2 sensible, computationally feasible statistical models in this situation. Indeed, Lindquist (2008) considered only non-bayesian approaches and concluded The size and complexity of the data make it difficult to create a full statistical model for describing its behavior, and a number of shortcuts are required to balance computational feasibility with model efficiency. In reality, the situation is much more delicate than we have indicated. Raw fmri data require preprocessing to remove artifacts accrued during acquisition. Throughout we will use standard fmri preprocessing techniques so that we can assume the data have the basic characteristics described in the previous paragraph. The interested reader should consult Friston et al. (2007) and Lindquist (2008) for accessible introductions to this important topic. Now we consider a specific fmri experiment which motivated the development of the general Bayesian approach described in Section An example fmri experiment The following experiment is fully analyzed by Lee et al. (2011) and is part of a longitudinal study of Alzheimer s disease. (We will consider a partial analysis of the data arising from this experiment in Section 2.1.) However, none of the subjects we considered have clinically diagnosed neurologic disorders, including Alzheimer s disease. The goal is to develop useful models of single-subject taskrelated activation in order to develop effective imaging biomarkers. Each subject completed a Stroop exam in an MRI scanner. That is, the subject is shown a word and subsequently asked to identify the color of the ink. There are three possible characteristics of each word: Ink only; for example, when presented with, XXXX, the subject should answer red. Congruence; for example, when presented with BLUE, the subject should answer blue. Interference; for example, when presented with BLUE, the subject should answer red. The Stroop exam was administered in a block design (see e.g. Friston et al., 2007) with a two second scanning time repetition. The battery of Ink, Congruence and Interference tasks were repeated in sequence 3 times with observations taken at a total of 465 time points. The Stroop task is an attempt to assess directed attention and hence it may be unsurprising that there is a well-documented age effect. That is, younger adults will typically experience a smaller interference effect than older adults. Also, interference effects may decrease with practice (Davidson, Zacks and Williams, 2003). In this study all subjects are older, generally well-educated, and healthy so we will not model an age effect, although it could be easily incorporated into our model below, but we should try to account for the possibility of a temporal drift in the response over the scanning session.
3 Jones et al./bayes for fmri 3 2. A Bayesian approach to detecting activation Our goal is detecting activation of a voxel when a stimulus is presented. For voxel v = 1,..., N let y v,i be the BOLD image intensity at time i = 1,..., T v and set y v = (y v,1,..., y v,tv ) T. Let X v be a known T v p matrix of full column rank and β v = (β v,1,..., β v,p ) T be a p 1 vector. (Typically, X v is a stimulus function convolved with the assumed haemodynamic response function (Friston et al., 2007). ) We assume a linear model for the conditional mean so that E[y v X v ] = X v β v. Detecting activation is equivalent to identifying the nonzero β v, a variable selection problem. We build on the work of Smith and Fahrmeir (2007) who were the first to use a variable selection approach in fmri settings, but our model differs substantively from theirs in that we model both spatial and temporal correlation, while allowing for the possibility of temporal drift in the BOLD signal over the scanning session. Let γ v = (γ v,1,..., γ v,p ) T be binary random variables used to indicate whether the voxel is activated by a sequence of input stimuli. That is, the coefficient β v,j is equal to zero if γ v,j = 0 and β v,j is nonzero if γ v,j = 1. The zero of γ v,j implies no effect on voxel v is caused by the corresponding experimental task j. Therefore, the model in can be expressed as y v = X v (γ v )β v (γ v ) + ε v, ε v N Tv (0, σ 2 vλ v ), where β v (γ v ) is the vector of nonzero regression coefficients and X v (γ v ) is the corresponding matrix for a given indicator variable γ v. An appropriate choice for the structure of the matrix Λ v will allow us to account for the temporal dependence between observations on a given voxel. We now turn our attention to specification of prior distributions. First consider the prior on β v (γ v ) given γ v for a particular voxel, v. We will use Zellner s g-prior (Zellner, 1996) given by where β v (γ v ) y v, σ 2 v, Λ v, γ v N( ˆβ v (γ v ), T v σ 2 v[x T v (γ v )Λ 1 v X v (γ v )] 1 ), (1) ˆβ v (γ v ) = [Xv T (γ v )Λ 1 v X v (γ v )] 1 Xv T (γ v )Λ 1 v y v. (2) We assume the σ 2 v are independent with prior π(σ 2 v) 1 σ 2 v. (3) Lee et al. (2011) demonstrate that the choice of prior on Λ v is critically important to the resulting posterior inference. It is natural to think of using a Wishart distribution or to assume the elements of Λ v (i, j) = ρ i j j, that is an AR(1) structure, and further that ρ v Uniform( 1, 1). Other autoregressive structures (e.g. AR(2)) are possible and popular. However, all of these possibilities require substantial computational effort. Thus some have suggested that Λ v
4 Jones et al./bayes for fmri 4 be taken to be the identity matrix. Note that this assumes independence in the time series on each voxel. Lee et al. (2011) show that little inferential efficacy is lost, compared to using autoregressive structures, by taking the prior on ρ v to be a point mass at the maximum likelihood estimator while the computational effort is comparable to assuming Λ v to be the identity matrix. We will use this point mass prior in the rest of this paper. We assume a binary spatial Ising prior for the indicator variables. Let γ (j) = {γ 1,j,..., γ N,j } be the vector of indicator variables for regression j over locations {1, 2,..., N}. We specify a three-dimensional neighborhood structure where each interior voxel has 6 neighbors. Let v k denote that voxels v and k are neighbors. Next, let the ω v,k be prespecified constants that allow us to weigh the interaction between neighboring locations on lattices v and k and let θ j > 0 represent the strength of the interaction between any two voxels. The prior on γ is then π(γ θ) = p j=1 π(γ (j) θ j ) with { N } π(γ (j) θ j ) exp α v,j γ v,j + θ j ω v,k I(γ v,j = γ k,j ). (4) The term v=1 N α v,j γ v,j. v=1 is known as the external field and is where we incorporate anatomical prior information. Finally, we assume a uniform prior on θ = (θ 1,..., θ p ), i.e., π(θ) p j=1 I(0 < θ j < θmax), where θmax is a user-specified hyperparameter. Notice that we assume θ, ρ, and σ 2 are a priori independent, γ conditionally independent, and independence across voxels. Thus the posterior density is characterized by q(β(γ), γ, ρ, θ, σ 2 y) v k N [p(y v β v (γ v ), γ v, σv, 2 Λ v )π(β v (γ v ) y v, γ v, σv, 2 Λ v ) v=1 π(σ 2 v)]π(γ θ)π(ρ)π(θ). (5) The posterior quantities of interest in detecting activation are the activation probability and its magnitude for each voxel, specifically, q(γ v,j = 1 y) and E(β v y). These quantities are analytically intractable and hence sophisticated Markov chain Monte Carlo methods are required to perform posterior inference; see Lee et al. (2011) for the details Application to example fmri data In the fmri experiment described in Section 1.1 the Stroop tasks are repeated three times. Here we consider how to study possible changes to the activation pattern over time. Let x i,j be the transformed input function of task j in ith trial, z i is a vector used to remove low-frequency, stimulus independent effects,
5 Jones et al./bayes for fmri 5 β i,j is the parameter of interest corresponding to the jth task in ith trial, the α i are nonzero nuisance parameters to model the baseline of a brain signal. Set Y v =α 0 z 0 + α 1 x 1 + β 1,1 x 1,1 + β 1,2 x 1,2 + β 1,3 x 1,3 + β 2,1 x 2,1 + β 2,2 x 2,2 + β 2,3 x 2,3 + β 3,1 x 3,1 + β 3,2 x 3,2 + β 3,3 x 3,3 + ε v, (6) where ε v N(0, σ 2 vλ v ). A sequence of binary variables is used to indicate if the corresponding parameter is zero or not so that for voxel v we have γ v = {1, 1, γ 3,v, γ 4,v, γ 5,v, γ 6,v, γ 7,v, γ 8,v, γ 9,v, γ 10,v, γ 11,v }. The rest of the model specification is as described above. The results for the interference task for each of the 3 trials are presented in Figure 1. Notice that there does appear to be a difference in the activation patterns over the trials. 3. Bayesian modeling versus voxel-level analyses A standard non-bayesian approach to analyzing fmri data is statistical parametric mapping (SPM). A normal-theory general linear model is fit for the time series of observations for each voxel, often assuming an AR correlation structure for the time series and a working assumption that all voxels are independent. Detecting activation of a voxel is done with a voxel-specific hypothesis test, producing a p-value (equivalently a test statistic which is usually a t- or F -test). The collection of p-values forms the SPM. In order to obtain an activation map of the brain, the next step is to threshold the spatially dependent p-values at a given overall error rate, leading to a problem of multiplicity. The most common way of addressing this problem is to use random field theory in which one must assume that the SPM is a good approximation to an underlying continuous random field. However, even with this assumption, the thresholding is done indirectly as it is based on the so-called expected Euler characteristic. For a full description of the regularity conditions and the underlying theory the reader can consult Friston et al. (1994), Worsley et al. (1992), Worsley (1994), Worsley et al. (1996) and Worsley (2003). Attempts to avoid random field assumptions employ multivariate permutation tests (Nichols and Holmes, 2002). Even if we accept the SPM approach it does not allow us to employ a principled approach to including prior information. In some fmri experiments we have a lot of scientific prior information about which regions of the brain should activate. In the SPM approach we would just restrict attention to that region plus a little bit more to be safe. Of course, incorporating prior information are conceptually straightforward in the Bayesian framework. Permutation tests, error rate thresholding (Chen et al., 2009; Genovese, Lazar and Nichols, 2002) and SPM random field corrections all have a similar theme of not producing a full model for the neuroimaging data. As such, activities such as a jointly modeling activation and connectivity are not possible. In contrast,
6 Jones et al./bayes for fmri 6 Fig 1. Activation maps for the time-varying model for the Interference task. The top panel contains activation maps for the first trial, the middle panel for the second trial and the bottom panel for the third trial.
7 Jones et al./bayes for fmri 7 Bayesian efforts are capable of producing a far more complete, thorough and robust approach towards investigating fmri data. However, more research is required to ensure that Bayesian modeling approaches are as computationally feasible as voxel-level approaches. References Chen, S., Wang, C., Eberly, L. E., Caffo, B. S. and Schwartz, B. S. (2009). Adaptive control of the false discovery rate in voxel-based morphometry. Human Brain Mapping Davidson, D. J., Zacks, R. T. and Williams, C. C. (2003). Stroop interference, practice, and aging. Aging, Neuropsychology, and Cognition Friston, K. J., Worsley, K. J., Frackowiak, R. S. J., Mazziotta, J. C. and Evans, A. C. (1994). Assessing the significance of focal activations using their spatial extent. Human Brain Mapping Friston, K. J., Ashburner, J. T., Kiebel, S. J., Nichols, T. E. and Penny, W. D. (2007). Statistical Parametric Mapping: The analysis of functional brain images. Academic Press, Amsterdam. Genovese, C. R., Lazar, N. A. and Nichols, T. (2002). Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage Lee, K.-J., Jones, G. L., Caffo, B. S. and Bassett, S. S. (2011). Spatial Bayesian variable selection models for functional magnetic resonance imaging time-series data. Preprint. Lindquist, M. A. (2008). The statistical analysis of fmri data. Statistical Science Nichols, T. E. and Holmes, A. P. (2002). Nonparametric permutation tests for functional neuroimaging: a primer with examples. Human Brain Mapping Smith, M. and Fahrmeir, L. (2007). Spatial Bayesian variable selection with application to functional magnetic resonance imaging. Journal of the American Statistical Association Worsley, K. J. (1994). Local maxima and the expected Euler characteristic of excursion sets of χ 2, F, and t fields. Advances in Applied Probability Worsley, K. J. (2003). Detecting activation in fmri data. Statistical Methods in Medical Research Worsley, K. J., Marrett, S., Neelin, P. and Evans, A. C. (1992). A three-dimensional statistical analysis for CBF activation studies in human brain. Journal of Cerebral Blood Flow and Metabolism Worsley, K. J., Marrett, S., Neelin, P., Vandal, A., Friston, K. J. and Evans, A. C. (1996). A unified statistical approach for determining significant signals in images of cerebral activation. Human Brain Mapping Zellner, A. (1996). On assessing prior distributions and Bayesian regression analysis with g-prior distributions. In Bayesian Inference and Decision Tech-
8 Jones et al./bayes for fmri 8 niques: Essays in Honor of Bruno de Finetti North-Holland/Elsevier
Bayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri
Bayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri Galin L. Jones 1 School of Statistics University of Minnesota March 2015 1 Joint with Martin Bezener and John Hughes Experiment
More informationBias in Resampling-Based Thresholding of Statistical Maps in fmri
Bias in Resampling-Based Thresholding of Statistical Maps in fmri Ola Friman and Carl-Fredrik Westin Laboratory of Mathematics in Imaging, Department of Radiology Brigham and Women s Hospital, Harvard
More informationIntroduction to Neuroimaging Janaina Mourao-Miranda
Introduction to Neuroimaging Janaina Mourao-Miranda Neuroimaging techniques have changed the way neuroscientists address questions about functional anatomy, especially in relation to behavior and clinical
More informationSpatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields
Spatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields Wei Liu 1, Peihong Zhu 1, Jeffrey S. Anderson 2, Deborah Yurgelun-Todd 3, and P. Thomas Fletcher 1 1 Scientific
More informationThe Effect of Correlation and Error Rate Specification on Thresholding Methods in fmri Analysis
The Effect of Correlation and Error Rate Specification on Thresholding Methods in fmri Analysis Brent R. Logan and Daniel B. Rowe, Division of Biostatistics and Department of Biophysics Division of Biostatistics
More informationEMPIRICALLY INVESTIGATING THE STATISTICAL VALIDITY OF SPM, FSL AND AFNI FOR SINGLE SUBJECT FMRI ANALYSIS
EMPIRICALLY INVESTIGATING THE STATISTICAL VALIDITY OF SPM, FSL AND AFNI FOR SINGLE SUBJECT FMRI ANALYSIS Anders Eklund a,b,c, Thomas Nichols d, Mats Andersson a,c, Hans Knutsson a,c a Department of Biomedical
More informationParametric Response Surface Models for Analysis of Multi-Site fmri Data
Parametric Response Surface Models for Analysis of Multi-Site fmri Data Seyoung Kim 1, Padhraic Smyth 1, Hal Stern 1, Jessica Turner 2, and FIRST BIRN 1 Bren School of Information and Computer Sciences,
More informationGraphical Models, Bayesian Method, Sampling, and Variational Inference
Graphical Models, Bayesian Method, Sampling, and Variational Inference With Application in Function MRI Analysis and Other Imaging Problems Wei Liu Scientific Computing and Imaging Institute University
More informationBayesian Inference in fmri Will Penny
Bayesian Inference in fmri Will Penny Bayesian Approaches in Neuroscience Karolinska Institutet, Stockholm February 2016 Overview Posterior Probability Maps Hemodynamic Response Functions Population
More informationBayesian Spatiotemporal Modeling for Detecting Neuronal Activation via Functional Magnetic Resonance Imaging
Bayesian Spatiotemporal Modeling for Detecting Neuronal Activation via Functional Magnetic Resonance Imaging Martin Bezener Stat-Ease, Inc. martin.bezener@gmail.com Lynn E. Eberly Division of Biostatistics
More informationA Nonparametric Bayesian Approach to Detecting Spatial Activation Patterns in fmri Data
A Nonparametric Bayesian Approach to Detecting Spatial Activation Patterns in fmri Data Seyoung Kim, Padhraic Smyth, and Hal Stern Bren School of Information and Computer Sciences University of California,
More informationDetecting Changes In Non-Isotropic Images
Detecting Changes In Non-Isotropic Images K.J. Worsley 1, M. Andermann 1, T. Koulis 1, D. MacDonald, 2 and A.C. Evans 2 August 4, 1999 1 Department of Mathematics and Statistics, 2 Montreal Neurological
More informationLinear Models in Medical Imaging. John Kornak MI square February 19, 2013
Linear Models in Medical Imaging John Kornak MI square February 19, 2013 Acknowledgement / Disclaimer Many of the slides in this lecture have been adapted from slides available in talks available on the
More informationPreprocessing II: Between Subjects John Ashburner
Preprocessing II: Between Subjects John Ashburner Pre-processing Overview Statistics or whatever fmri time-series Anatomical MRI Template Smoothed Estimate Spatial Norm Motion Correct Smooth Coregister
More informationNeuroimaging and mathematical modelling Lesson 2: Voxel Based Morphometry
Neuroimaging and mathematical modelling Lesson 2: Voxel Based Morphometry Nivedita Agarwal, MD Nivedita.agarwal@apss.tn.it Nivedita.agarwal@unitn.it Volume and surface morphometry Brain volume White matter
More informationAdvances in FDR for fmri -p.1
Advances in FDR for fmri Ruth Heller Department of Statistics, University of Pennsylvania Joint work with: Yoav Benjamini, Nava Rubin, Damian Stanley, Daniel Yekutieli, Yulia Golland, Rafael Malach Advances
More informationImproving CCA based fmri Analysis by Covariance Pooling - Using the GPU for Statistical Inference
Improving CCA based fmri Analysis by Covariance Pooling - Using the GPU for Statistical Inference Anders Eklund, Mats Andersson and Hans Knutsson Linköping University Post Print N.B.: When citing this
More informationLinear Models in Medical Imaging. John Kornak MI square February 22, 2011
Linear Models in Medical Imaging John Kornak MI square February 22, 2011 Acknowledgement / Disclaimer Many of the slides in this lecture have been adapted from slides available in talks available on the
More informationEffect of age and dementia on topology of brain functional networks. Paul McCarthy, Luba Benuskova, Liz Franz University of Otago, New Zealand
Effect of age and dementia on topology of brain functional networks Paul McCarthy, Luba Benuskova, Liz Franz University of Otago, New Zealand 1 Structural changes in aging brain Age-related changes in
More informationNew and best-practice approaches to thresholding
New and best-practice approaches to thresholding Thomas Nichols, Ph.D. Department of Statistics & Warwick Manufacturing Group University of Warwick FIL SPM Course 17 May, 2012 1 Overview Why threshold?
More informationStatistical Modeling of Neuroimaging Data: Targeting Activation, Task-Related Connectivity, and Prediction
Statistical Modeling of Neuroimaging Data: Targeting Activation, Task-Related Connectivity, and Prediction F. DuBois Bowman Department of Biostatistics and Bioinformatics Emory University, Atlanta, GA,
More informationStatistical Analysis of Neuroimaging Data. Phebe Kemmer BIOS 516 Sept 24, 2015
Statistical Analysis of Neuroimaging Data Phebe Kemmer BIOS 516 Sept 24, 2015 Review from last time Structural Imaging modalities MRI, CAT, DTI (diffusion tensor imaging) Functional Imaging modalities
More informationBHMSMAfMRI User Manual
BHMSMAfMRI User Manual Nilotpal Sanyal Bayesian and Interdisciplinary Research Unit Indian Statistical Institute, Calcutta, India nsanyal v@isical.ac.in Welcome to BHMSMAfMRI, an R package to analyze functional
More informationAnalysis of Functional MRI Timeseries Data Using Signal Processing Techniques
Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October
More informationCS/NEUR125 Brains, Minds, and Machines. Due: Wednesday, April 5
CS/NEUR125 Brains, Minds, and Machines Lab 8: Using fmri to Discover Language Areas in the Brain Due: Wednesday, April 5 In this lab, you will analyze fmri data from an experiment that was designed to
More informationAn improved theoretical P-value for SPMs based on discrete local maxima
An improved theoretical P-value for SPMs based on discrete local maxima K.J. Worsley August 9, 25 Department of Mathematics and Statistics, McGill University, 85 Sherbrooke St. West, Montreal, Québec,
More informationFunctional and dynamic Magnetic Resonance Imaging using vector adaptive weights smoothing
Functional and dynamic Magnetic Resonance Imaging using vector adaptive weights smoothing Jörg Polzehl and Vladimir G. Spokoiny Weierstrass Institute for Applied Analysis and Stochastics, Germany Summary.
More information7/15/2016 ARE YOUR ANALYSES TOO WHY IS YOUR ANALYSIS PARAMETRIC? PARAMETRIC? That s not Normal!
ARE YOUR ANALYSES TOO PARAMETRIC? That s not Normal! Martin M Monti http://montilab.psych.ucla.edu WHY IS YOUR ANALYSIS PARAMETRIC? i. Optimal power (defined as the probability to detect a real difference)
More informationLinear Models in Medical Imaging. John Kornak MI square February 21, 2012
Linear Models in Medical Imaging John Kornak MI square February 21, 2012 Acknowledgement / Disclaimer Many of the slides in this lecture have been adapted from slides available in talks available on the
More informationFUNCTIONAL magnetic resonance imaging (fmri) is
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. X, NO. X, MARCH 2 A Bayesian Mixture Approach to Modeling Spatial Activation Patterns in Multi-site fmri Data Seyoung Kim, Padhraic Smyth, Member, IEEE, and Hal
More informationMultiVariate Bayesian (MVB) decoding of brain images
MultiVariate Bayesian (MVB) decoding of brain images Alexa Morcom Edinburgh SPM course 2015 With thanks to J. Daunizeau, K. Brodersen for slides stimulus behaviour encoding of sensorial or cognitive state?
More informationTime-Frequency Method Based Activation Detection in Functional MRI Time-Series
Time-Frequency Method Based Activation Detection in Functional MRI Time-Series Arun Kumar 1,2 and Jagath C. Rajapakse 1,3 1 School of Computing, NTU, Singapore 2 School of EEE,Singapore Polytechnic, Singapore
More informationMethods for data preprocessing
Methods for data preprocessing John Ashburner Wellcome Trust Centre for Neuroimaging, 12 Queen Square, London, UK. Overview Voxel-Based Morphometry Morphometry in general Volumetrics VBM preprocessing
More informationSupplementary Figure 1. Decoding results broken down for different ROIs
Supplementary Figure 1 Decoding results broken down for different ROIs Decoding results for areas V1, V2, V3, and V1 V3 combined. (a) Decoded and presented orientations are strongly correlated in areas
More informationStatistical Analysis of MRI Data
Statistical Analysis of MRI Data Shelby Cummings August 1, 2012 Abstract Every day, numerous people around the country go under medical testing with the use of MRI technology. Developed in the late twentieth
More informationA Functional Connectivity Inspired Approach to Non-Local fmri Analysis
A Functional Connectivity Inspired Approach to Non-Local fmri Analysis Anders Eklund, Mats Andersson and Hans Knutsson Linköping University Post Print N.B.: When citing this work, cite the original article.
More informationMultivariate analyses & decoding
Multivariate analyses & decoding Kay Henning Brodersen Computational Neuroeconomics Group Department of Economics, University of Zurich Machine Learning and Pattern Recognition Group Department of Computer
More informationIntroductory Concepts for Voxel-Based Statistical Analysis
Introductory Concepts for Voxel-Based Statistical Analysis John Kornak University of California, San Francisco Department of Radiology and Biomedical Imaging Department of Epidemiology and Biostatistics
More informationAn Empirical Comparison of SPM Preprocessing Parameters to the Analysis of fmri Data
NeuroImage 17, 19 28 (2002) doi:10.1006/nimg.2002.1113 An Empirical Comparison of SPM Preprocessing Parameters to the Analysis of fmri Data Valeria Della-Maggiore, Wilkin Chau, Pedro R. Peres-Neto,* and
More informationSession scaling; global mean scaling; block effect; mean intensity scaling
Types of Scaling Session scaling; global mean scaling; block effect; mean intensity scaling Purpose remove intensity differences between runs (i.e., the mean of the whole time series). whole time series
More informationCorrection of Partial Volume Effects in Arterial Spin Labeling MRI
Correction of Partial Volume Effects in Arterial Spin Labeling MRI By: Tracy Ssali Supervisors: Dr. Keith St. Lawrence and Udunna Anazodo Medical Biophysics 3970Z Six Week Project April 13 th 2012 Introduction
More informationA User s Guide to Graphical-Model-based Multivariate Analysis
A User s Guide to Graphical-Model-based Multivariate Analysis 1. Introduction Rong Chen December 2006 Graphical-Model-based Multivariate Analysis (GAMMA) is a Bayesian data mining software for structural
More informationAn Introduction To Automatic Tissue Classification Of Brain MRI. Colm Elliott Mar 2014
An Introduction To Automatic Tissue Classification Of Brain MRI Colm Elliott Mar 2014 Tissue Classification Tissue classification is part of many processing pipelines. We often want to classify each voxel
More informationLinear Models in Medical Imaging. John Kornak MI square February 23, 2010
Linear Models in Medical Imaging John Kornak MI square February 23, 2010 Acknowledgement / Disclaimer Many of the slides in this lecture have been adapted from slides available in talks available on the
More informationVoxel selection algorithms for fmri
Voxel selection algorithms for fmri Henryk Blasinski December 14, 2012 1 Introduction Functional Magnetic Resonance Imaging (fmri) is a technique to measure and image the Blood- Oxygen Level Dependent
More informationCHAPTER 2. Morphometry on rodent brains. A.E.H. Scheenstra J. Dijkstra L. van der Weerd
CHAPTER 2 Morphometry on rodent brains A.E.H. Scheenstra J. Dijkstra L. van der Weerd This chapter was adapted from: Volumetry and other quantitative measurements to assess the rodent brain, In vivo NMR
More informationExtending the GLM. Outline. Mixed effects motivation Evaluating mixed effects methods Three methods. Conclusions. Overview
Extending the GLM So far, we have considered the GLM for one run in one subject The same logic can be applied to multiple runs and multiple subjects GLM Stats For any given region, we can evaluate the
More informationCS 229 Final Project Report Learning to Decode Cognitive States of Rat using Functional Magnetic Resonance Imaging Time Series
CS 229 Final Project Report Learning to Decode Cognitive States of Rat using Functional Magnetic Resonance Imaging Time Series Jingyuan Chen //Department of Electrical Engineering, cjy2010@stanford.edu//
More informationSupervised Learning for Image Segmentation
Supervised Learning for Image Segmentation Raphael Meier 06.10.2016 Raphael Meier MIA 2016 06.10.2016 1 / 52 References A. Ng, Machine Learning lecture, Stanford University. A. Criminisi, J. Shotton, E.
More informationFactor Topographic Latent Source Analysis: Factor Analysis for Brain Images
Factor Topographic Latent Source Analysis: Factor Analysis for Brain Images Jeremy R. Manning 1,2, Samuel J. Gershman 3, Kenneth A. Norman 1,3, and David M. Blei 2 1 Princeton Neuroscience Institute, Princeton
More informationStatistical Matching using Fractional Imputation
Statistical Matching using Fractional Imputation Jae-Kwang Kim 1 Iowa State University 1 Joint work with Emily Berg and Taesung Park 1 Introduction 2 Classical Approaches 3 Proposed method 4 Application:
More informationfor Images A Bayesian Deformation Model
Statistics in Imaging Workshop July 8, 2004 A Bayesian Deformation Model for Images Sining Chen Postdoctoral Fellow Biostatistics Division, Dept. of Oncology, School of Medicine Outline 1. Introducing
More informationCognitive States Detection in fmri Data Analysis using incremental PCA
Department of Computer Engineering Cognitive States Detection in fmri Data Analysis using incremental PCA Hoang Trong Minh Tuan, Yonggwan Won*, Hyung-Jeong Yang International Conference on Computational
More informationNIH Public Access Author Manuscript Med Image Comput Comput Assist Interv. Author manuscript; available in PMC 2014 August 15.
NIH Public Access Author Manuscript Med Image Comput Comput Assist Interv. Author manuscript; available in PMC 2014 August 15. Published in final edited form as: Med Image Comput Comput Assist Interv.
More informationJournal of Statistical Software
JSS Journal of Statistical Software January 2014, Volume 56, Issue 9. http://www.jstatsoft.org/ POBE: A Computer Program for Optimal Design of Multi-Subject Blocked fmri Experiments Bärbel Maus University
More informationOptimally-Discriminative Voxel-Based Analysis
Optimally-Discriminative Voxel-Based Analysis Tianhao Zhang and Christos Davatzikos Section of Biomedical Image Analysis, Department of Radiology, University of Pennsylvania, Philadelphia, PA 19104, USA
More informationOnline Statistical Inference for Quantifying Mandible Growth in CT Images
Online Statistical Inference for Quantifying Mandible Growth in CT Images Moo K. Chung 1,2, Ying Ji Chuang 2, Houri K. Vorperian 2 1 Department of Biostatistics and Medical Informatics 2 Vocal Tract Development
More informationAn independent component analysis based tool for exploring functional connections in the brain
An independent component analysis based tool for exploring functional connections in the brain S. M. Rolfe a, L. Finney b, R. F. Tungaraza b, J. Guan b, L.G. Shapiro b, J. F. Brinkely b, A. Poliakov c,
More informationTutorial MICCAI'04. fmri data analysis: state of the art and future challenges. Jean-Baptiste Poline Philippe Ciuciu Alexis Roche
Tutorial MICCAI'04 fmri data analysis: state of the art and future challenges Jean-Baptiste Poline Philippe Ciuciu Alexis Roche CEA/SHFJ, Orsay (France) Goal of this tutorial and plan I. Talk 1: will orientate
More informationAdaptive statistical parametric mapping for fmri
Statistics and Its Interface Volume 3 (2010) 33 43 Adaptive statistical parametric mapping for fmri Ping Bai, Haipeng Shen, Jianhua Z Huang and Young K Truong, Brain activity is accompanied by changes
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationMonte Carlo for Spatial Models
Monte Carlo for Spatial Models Murali Haran Department of Statistics Penn State University Penn State Computational Science Lectures April 2007 Spatial Models Lots of scientific questions involve analyzing
More informationPattern Recognition for Neuroimaging Data
Pattern Recognition for Neuroimaging Data Edinburgh, SPM course April 2013 C. Phillips, Cyclotron Research Centre, ULg, Belgium http://www.cyclotron.ulg.ac.be Overview Introduction Univariate & multivariate
More informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationResidual Analysis for Detecting Mis-modeling in fmri. (Running title: Detecting Mis-modeling in fmri)
Residual Analysis for Detecting Mis-modeling in fmri (Running title: Detecting Mis-modeling in fmri) J. M. Loh 1, M. A. Lindquist and T. D. Wager 3 1 Dept of Statistics, Columbia U, New York; Email:meng@stat.columbia.edu,
More informationMultivariate pattern classification
Multivariate pattern classification Thomas Wolbers Space & Ageing Laboratory (www.sal.mvm.ed.ac.uk) Centre for Cognitive and Neural Systems & Centre for Cognitive Ageing and Cognitive Epidemiology Outline
More informationClassical and Bayesian Inference in fmri
1 Classical and Bayesian Inference in fmri William D. Penny and Karl J. Friston The Wellcome Department of Imaging Neuroscience, University College London Queen Square, London, UK WC1N 3BG Tel (44 020
More informationCorrection for multiple comparisons. Cyril Pernet, PhD SBIRC/SINAPSE University of Edinburgh
Correction for multiple comparisons Cyril Pernet, PhD SBIRC/SINAPSE University of Edinburgh Overview Multiple comparisons correction procedures Levels of inferences (set, cluster, voxel) Circularity issues
More informationAutomatic Generation of Training Data for Brain Tissue Classification from MRI
MICCAI-2002 1 Automatic Generation of Training Data for Brain Tissue Classification from MRI Chris A. Cocosco, Alex P. Zijdenbos, and Alan C. Evans McConnell Brain Imaging Centre, Montreal Neurological
More informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationBasic Introduction to Data Analysis. Block Design Demonstration. Robert Savoy
Basic Introduction to Data Analysis Block Design Demonstration Robert Savoy Sample Block Design Experiment Demonstration Use of Visual and Motor Task Separability of Responses Combined Visual and Motor
More information1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12
Contents 1 Introduction 10 1.1 Motivation and Aims....... 10 1.1.1 Functional Imaging.... 10 1.1.2 Computational Neuroanatomy... 12 1.2 Overview of Chapters... 14 2 Rigid Body Registration 18 2.1 Introduction.....
More informationAdaptive spatial smoothing of fmri images (Running title: Adaptive spatial smoothing)
Adaptive spatial smoothing of fmri images (Running title: Adaptive spatial smoothing) Yu (Ryan) Yue, Ji Meng Loh and M. A. Lindquist Deptartment of Statistics and CIS, Baruch College, New York, NY, 10010;
More informationTemporal and Cross-Subject Probabilistic Models for fmri Prediction Tasks
Temporal and Cross-Subject Probabilistic Models for fmri Prediction Tasks Alexis Battle Gal Chechik Daphne Koller Department of Computer Science Stanford University Stanford, CA 94305-9010 {ajbattle,gal,koller}@cs.stanford.edu
More informationArtifact detection and repair in fmri
Artifact detection and repair in fmri Paul K. Mazaika, Ph.D. Center for Interdisciplinary Brain Sciences Research (CIBSR) Division of Interdisciplinary Behavioral Sciences Stanford University School of
More informationMultivariate Pattern Classification. Thomas Wolbers Space and Aging Laboratory Centre for Cognitive and Neural Systems
Multivariate Pattern Classification Thomas Wolbers Space and Aging Laboratory Centre for Cognitive and Neural Systems Outline WHY PATTERN CLASSIFICATION? PROCESSING STREAM PREPROCESSING / FEATURE REDUCTION
More informationMedical Image Analysis
Medical Image Analysis Instructor: Moo K. Chung mchung@stat.wisc.edu Lecture 10. Multiple Comparisons March 06, 2007 This lecture will show you how to construct P-value maps fmri Multiple Comparisons 4-Dimensional
More informationSupplementary Figure 1
Supplementary Figure 1 BOLD and CBV functional maps showing EPI versus line-scanning FLASH fmri. A. Colored BOLD and CBV functional maps are shown in the highlighted window (green frame) of the raw EPI
More informationAn Intensity Consistent Approach to the Cross Sectional Analysis of Deformation Tensor Derived Maps of Brain Shape
An Intensity Consistent Approach to the Cross Sectional Analysis of Deformation Tensor Derived Maps of Brain Shape C. Studholme, V. Cardenas, A. Maudsley, and M. Weiner U.C.S.F., Dept of Radiology, VAMC
More informationarxiv: v2 [stat.ap] 11 Mar 2017
A Bayesian Heteroscedastic GLM with Application to fmr Data with Motion Spikes Anders Eklund a,b,c,, Martin A. Lindquist d, Mattias Villani* a a Division of Statistics & Machine Learning, Department of
More informationASAP_2.0 (Automatic Software for ASL Processing) USER S MANUAL
ASAP_2.0 (Automatic Software for ASL Processing) USER S MANUAL ASAP was developed as part of the COST Action "Arterial Spin Labelling Initiative in Dementia (AID)" by: Department of Neuroimaging, Institute
More informationStructural MRI analysis
Structural MRI analysis volumetry and voxel-based morphometry cortical thickness measurements structural covariance network mapping Boris Bernhardt, PhD Department of Social Neuroscience, MPI-CBS bernhardt@cbs.mpg.de
More informationNIH Public Access Author Manuscript Neuroimage. Author manuscript; available in PMC 2015 July 15.
NIH Public Access Author Manuscript Published in final edited form as: Neuroimage. 2014 July 15; 95: 162 175. doi:10.1016/j.neuroimage.2014.03.024. A Spatio-Temporal Nonparametric Bayesian Variable Selection
More informationIS MRI Based Structure a Mediator for Lead s Effect on Cognitive Function?
IS MRI Based Structure a Mediator for Lead s Effect on Cognitive Function? Brian Caffo, Sining Chen, Brian Schwartz Department of Biostatistics and Environmental Health Sciences Johns Hopkins University
More informationChapter 3 Set Redundancy in Magnetic Resonance Brain Images
16 Chapter 3 Set Redundancy in Magnetic Resonance Brain Images 3.1 MRI (magnetic resonance imaging) MRI is a technique of measuring physical structure within the human anatomy. Our proposed research focuses
More informationPBSI: A symmetric probabilistic extension of the Boundary Shift Integral
PBSI: A symmetric probabilistic extension of the Boundary Shift Integral Christian Ledig 1, Robin Wolz 1, Paul Aljabar 1,2, Jyrki Lötjönen 3, Daniel Rueckert 1 1 Department of Computing, Imperial College
More informationA TEMPORAL FREQUENCY DESCRIPTION OF THE SPATIAL CORRELATION BETWEEN VOXELS IN FMRI DUE TO SPATIAL PROCESSING. Mary C. Kociuba
A TEMPORAL FREQUENCY DESCRIPTION OF THE SPATIAL CORRELATION BETWEEN VOXELS IN FMRI DUE TO SPATIAL PROCESSING by Mary C. Kociuba A Thesis Submitted to the Faculty of the Graduate School, Marquette University,
More informationThe organization of the human cerebral cortex estimated by intrinsic functional connectivity
1 The organization of the human cerebral cortex estimated by intrinsic functional connectivity Journal: Journal of Neurophysiology Author: B. T. Thomas Yeo, et al Link: https://www.ncbi.nlm.nih.gov/pubmed/21653723
More informationImage analysis. Computer Vision and Classification Image Segmentation. 7 Image analysis
7 Computer Vision and Classification 413 / 458 Computer Vision and Classification The k-nearest-neighbor method The k-nearest-neighbor (knn) procedure has been used in data analysis and machine learning
More informationPOBE: A Computer Program for Optimal Design of Multi-Subject Blocked FMRI Experiments
POBE: A Computer Program for Optimal Design of Multi-Subject Blocked FMRI Experiments Bärbel Maus, University of Liège Gerard J.P. van Breukelen, Maastricht University March 20, 2013 Abstract For functional
More informationData mining for neuroimaging data. John Ashburner
Data mining for neuroimaging data John Ashburner MODELLING The Scientific Process MacKay, David JC. Bayesian interpolation. Neural computation 4, no. 3 (1992): 415-447. Model Selection Search for the best
More informationCOPULA MODELS FOR BIG DATA USING DATA SHUFFLING
COPULA MODELS FOR BIG DATA USING DATA SHUFFLING Krish Muralidhar, Rathindra Sarathy Department of Marketing & Supply Chain Management, Price College of Business, University of Oklahoma, Norman OK 73019
More informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
More informationMultiple comparisons problem and solutions
Multiple comparisons problem and solutions James M. Kilner http://sites.google.com/site/kilnerlab/home What is the multiple comparisons problem How can it be avoided Ways to correct for the multiple comparisons
More informationFeature Selection for fmri Classification
Feature Selection for fmri Classification Chuang Wu Program of Computational Biology Carnegie Mellon University Pittsburgh, PA 15213 chuangw@andrew.cmu.edu Abstract The functional Magnetic Resonance Imaging
More informationA Bayesian General Linear Modeling Approach to Cortical Surface fmri Data Analysis
A Bayesian General Linear Modeling Approach to Cortical Surface fmri Data Analysis Amanda F. Mejia a, Yu Ryan Yue b, David Bolin c, Finn Lindgren d and Martin A. Lindquist e arxiv:1706.00959v1 [stat.ap]
More informationMultiple Testing and Thresholding
Multiple Testing and Thresholding UCLA Advanced NeuroImaging Summer School, 2007 Thanks for the slides Tom Nichols! Overview Multiple Testing Problem Which of my 100,000 voxels are active? Two methods
More informationComputational Neuroanatomy
Computational Neuroanatomy John Ashburner john@fil.ion.ucl.ac.uk Smoothing Motion Correction Between Modality Co-registration Spatial Normalisation Segmentation Morphometry Overview fmri time-series kernel
More informationSimultaneous Estimation of Response Fields and Impulse Response Functions PhD Thesis Proposal
Simultaneous Estimation of Response Fields and Impulse Response Functions PhD Thesis Proposal Erich Huang Carnegie Mellon University - Department of Statistics February 11, 2009 Abstract Relevant issues
More informationHHS Public Access Author manuscript Med Image Comput Comput Assist Interv. Author manuscript; available in PMC 2018 March 20.
Online Statistical Inference for Large-Scale Binary Images Moo K. Chung 1,2, Ying Ji Chuang 2, and Houri K. Vorperian 2 1 Department of Biostatistics and Medical Informatics, University of Wisconsin, Madison,
More information