the PyHRF package P. Ciuciu1,2 and T. Vincent1,2 Methods meeting at Neurospin 1: CEA/NeuroSpin/LNAO
|
|
- Hugh Clark
- 6 years ago
- Views:
Transcription
1 Joint detection-estimation of brain activity from fmri time series: the PyHRF package Methods meeting at Neurospin P. Ciuciu1,2 and T. Vincent1,2 1: CEA/NeuroSpin/LNAO 2: IFR49
2 Outline 2/61 I. Introduction II. A novel framework for intra-subject analysis in fmri III. Implementation : the PyHRF software IV. Conclusions and perspectives
3 3/61 Classical fmri analysis 1. Detect and localise brain activations Ex: In SPM [Friston et al, 1994], the BOLD response is modelled with: or Compute statistical activation Maps 2. Estimate the dynamics of activation [Goutte et al, IEEE TMI 2000; Marrelec et al, HBM 2003] HRF estimates Probe brain dynamics Time s
4 4/61 Classical fmri analysis GLM limitations: A single HRF shape is not physiologically appropriate Studies report HRF variability: within subject (between regions, sessions, conditions, trials) [Miezin et al., NIM 2000; Ciuciu et al., IEEE TMI 2003] [Neumann et al, NIM 2003, 2006; Smith et al, NIM 2005] between subjects [Handwerker et al., NIM 2004; Aguirre et al., NIM 1998] between groups (infants, patients,...)... [D'Esposito et al, NIM 1998, 2003; Richter and Richter, NIM 2003]
5 5/61 Motivations Two main issues in fmri Detection of brain activation and estimation of brain dynamics are addressed separately Any detection method supposes a given HRF shape Any estimation algorithm provides relevant results in activated voxels or regions only Address these two problems simultaneously in a joint detection-estimation (JDE) framework [Makni et al., IEEE SP 2005, ISBI'06, ICASSP'06, NIM, 2008; Vincent et al, ICASSP'07, EMBC'07, ISBI'08; Ciuciu et al, ISBI'08 ; Risser et al, sub. to NIPS'08; Vincent et al, sub. to IEEE TMI]
6 Outline 6/61 I. Introduction II. A novel framework for intra-subject analysis in fmri III. Implementation : the pyhrf software IV. Conclusions and perspectives
7 Forward BOLD model 7/61 [Makni, Ciuciu et al., IEEE SP 2005; Makni et al, Ciuciu, NeuroImage, 2008] Region definition: voxels Region of interest fmri time series in voxels of the ROI Experimental paradigm: Arrival times of M stimuli (onsets)
8 8/61 Forward BOLD model Temporal hypotheses: for standard ISIs Hyp. 1: linearity Hyp 2: stationarity + Convolution kernel Hyp 3: additivity Time in s TR 2TR3TR... Stimulus-varying NRLs:
9 Forward BOLD model 9/61 Spatial hypotheses: functionally homogeneous ROI Single HRF shape Voxel-dependent magnitudes of the BOLD response Neural Response Levels
10 Forward BOLD model 10/61 Spatial hypotheses: functionally homogeneous ROI Single HRF shape [Makni, Ciuciu et al., IEEE SP 2005; Makni et al, Ciuciu, NeuroImage, 2008] Voxel-dependent magnitudes of the BOLD response
11 11/61 Forward BOLD model Unknown parameters NRL in and for condition = = BOLD signal measured in voxel Arrival time of stimulus Noise statistics Drift coefficients in voxel HRF Orthonormal basis for low frequency drift modelling Known parameters
12 Bayes rule 12/61 likelihood How the data are generated from the parameters? Forward modeling
13 13/61 Likelihood definition Main hypothesis: noise decorrelated in space fmri time series are statistically independent in space: Temporal noise model: either white or serially correlated AR(1) [Makni et al, Ciuciu, NeuroImage 2008]
14 14/61 Bayes rule Prior What do we know about the parameters before the data are acquired? Prior modeling
15 HRF prior modeling 15/61 Information about the HRF shape nonparametric approaches Averaging: [Buckner et al, 1996] Selective averaging: [Dale et al, HBM 1997] FIR (Finite Impulse Response) Regularized FIR (smoothing prior): [Marrelec, Ciuciu et al, IPMI 03 Ciuciu et al., IEEE TMI 2003]
16 16/61 NRL prior modeling Spatial information [Vincent, Ciuciu et al, ICASSP'07] Independence between conditions: Spatial mixture model (SMM) for each Discrete Markov random Field (eg Potts/Ising): Partition function:
17 17/61 Prior model Two-class Gaussian mixture [Vincent et al., ICASSP' 07] Non-activating voxels: Activating voxels: Unknown hyper-parameters: Influence of
18 18/61 Bayes rule What do we know about the HRF, the NRLs and the hyper-parameters given the data? Keystone of learning scheme: simulating realizations of using Gibbs sampler
19 19/61 Real data sets Localizer fmri experiment: Experimental conditions under study: auditory and visual stimuli Event-related paradigm : - Short stimuli duration - Inter-stimulus interval : ~3s to 10s - Randomised sequence 125 scans with TR = 2.4s, scanning at 3T ROIs located in the primary auditory cortices
20 20/61 SPM vs. JDE Auditory Visual contrast: % Δ BOLD signal JDE SPM % Δ BOLD signal Time in s. Time in s. Bilateral activation detected along the gray matter from raw data sets (spatially unsmoothed)
21 21/61 Whole brain analysis First step: Segmentation of the gray/white matter interface from the T1 MRI Second step : brain parcellation based upon functional similarities and spatial connectivty [Flandin et al, ISBI'02; Thirion et al, HBM 2006] [Makni et al, Ciuciu, NIM 2008]
22 Sources of variability Subject, region 22/61 [Birn et al, NIM 2001; Marrelec et al, HBM 2003; Menon et al NIM 2003; Neumann et al NIM 2003; Handwerker et al, NIM 2004] [Makni et al, Ciuciu, NIM 2008]
23 Whole brain: supervised SMM =0.1 =0.3 23/61 =0.7
24 Outline 24/61 I. Introduction II. A novel framework for intra-subject analysis in fmri III. Implementation : the PyHRF software IV. Conclusions and perspectives
25 PyHRF package 25/61 Outline General framework Implementation and requirements Setup and use cases Viewer overview
26 fmri analysis with PyHRF 26/61 t t t PyHRF Parcellation = n-ary mask (0 label = background) from SPM beta maps from 4D fmri series
27 fmri analysis with PyHRF 27/61 t t t HRF PyHRF Parcellation = n-ary mask (0 label = background) from SPM beta maps from 4D fmri series contrasts
28 Package description 28/61 Mainly command line Setup with XML files Interactive dedicated viewer Python and C Time issues: parcellation : several minutes JDE : 1 to 2 hours (can be parallelized) Available at: $P4/neursopy-.../ $python setup.py install
29 29/61 Dependencies gsl scipy numpy fff PyHRF clustering viewer PyQt JDE PyXML pynifti SIP Qt nifticlib
30 30/61 Dependencies gsl scipy numpy fff PyHRF clustering viewer PyQt JDE PyXML pynifti SIP Qt nifticlib
31 Dependencies (versions) 31/61 fff (snapshot SVN may 2008) GSL 1.9 numpy 1.0 scipy PyXML nifticlib 1.0 pynifti PyQt Qt SIP 4.7
32 Use cases 32/61 Parcellation JDE Viewer
33 Parcellation 33/61 $pyhrf_clus_buildcfg clustering.xml $pyhrf_clus_run parcellation.nii
34 Use cases 34/61 Parcellation JDE Viewer
35 35/61 Joint Detection-Estimation $pyhrf_jde_buildcfg detectestim.xml $pyhrf_jde_estim -v1 session0_pm_hrf.nii session0_pmnrl_condition_audio.nii session0_pmnrl_condition_audio.nii...
36 Joint Detection-Estimation 36/61 $pyhrf_jde_buildcfg detectestim.xml Data definition Model parameters $pyhrf_jde_estim -v1 session0_pm_hrf.nii session0_pmnrl_condition_audio.nii session0_pmnrl_condition_audio.nii...
37 Joint Detection-Estimation 37/61 detectestim.xml Data definition - seconds - separator : space parcellation with positions excluded during analysis (output file) Input parcellation BOLD data - seconds - separator : space - if empty : spike stimuli
38 Joint Detection-Estimation 38/61 detectestim.xml Model definition (1) Drift base type ('polynomial' or 'cosine') Drift parameter: polynom order or cut-off frequency
39 Joint Detection-Estimation 39/61 detectestim.xml Définition du modèle (2) Amount of spatial correlation (between 0.3 et 0.7) Contrast definition Separator: ; If no contrast: ; Boolean flags for outputs Constraint of HRF=0 at first and last time position (boolean) Boolean flag for the HRF estimation. If 0 : use of a fixed canonical HRF
40 40/61 Joint Detection-Estimation Multi-session aspect $pyhrf_jde_buildcfg -s homoscedastic $pyhrf_jde_buildcfg -s independent Two different xml setup files
41 41/61 Joint Detection-Estimation $pyhrf_jde_buildcfg detectestim.xml $pyhrf_jde_estim -v1 session0_pm_hrf.nii session0_pmnrl_condition_audio.nii session0_pmnrl_condition_audio.nii...
42 Use cases 42/61 Parcellation JDE Viewer
43 Viewer 43/61 $pyhrf_view session0*
44 44/61 Viewer $pyhrf_view session0* Data browser (slice definition) current slice view
45 45/61 Viewer Data browser (slice definition) current slice view Slice dimension
46 46/61 Viewer Data browser (slice definition) current slice view Crop to mask
47 47/61 Viewer Data browser (slice definition) current slice view Slice axes selection
48 48/61 Viewer Data browser (slice definition) current slice view Mode 2D: image (matrix) Mode 1D: curves
49 PyHRF package in factory 49/61 Methodological developments Unsupervised spatial regularization Habituation modeling Model comparison Software developments Integration to the fmri Brainvisa toolbox User friendly parallel system Scriptable API Port of the matlab HRF toolbox: L. Risser's work.
50 Conclusions 50/61 The joint detection-estimation framework: directly accounts for different sources of variability provides both region-based HRF time courses and contrast maps No spatial filtering/smoothing of fmri datasets depends on an input parcellation: [Vincent et al, ISBI'08] gives improved RFX maps in comparison with SPM First release of Pyhrf package (v 1.0) available either in the BrainVisa software or as a standalone toolbox.
51 Extract temporal information 51/61 Multiple conditions Primary auditory cortex Time in s Planum temporale Time in s [Ciuciu et al, IEEE TMI 2003]
52 Sensitivity analysis 52/61 [Vincent, Ciuciu, ISBI 2008] Random parcellation JDE Single JDE run GMM parcellation JDE
53 Sensitivity analysis 53/61 [Vincent, Ciuciu, ISBI 2008] Random parcellations JDE JDE JDE JDE... Single JDE run GMM parcellation JDE
54 Whole brain: supervised SMM =0.1 =0.3 54/61 =0.7
55 55/61 Supervised vs unsupervised SMM unsup. =0.1 =0.7
56 Adaptive spatial regularization 56/61 RW Metropolis-Hastings step for sampling Preliminary estimation of the partition function Fondamental identity: Practical implementation: Tabulate over a fine grid For Generate of (SW scheme) Compute [Meng and Rubin BIOM 1998, Higdon et al, IEEE TMI 1997]
57 Why modeling the HRF? 57/61 Motivations Extract temporal information Account for variability sources Gain in sensitivity of detection [Marrelec et al, HBM 2003] Compare BOLD fmri sequences (EPI vs EVI) [Rabrait et al, JMRI 2008] Coherent approach for EEG/fMRI fusion [Daunizeau et al, IEEE SP 2005, NIM 2006/07]
58 58/61 Adaptive spatial regularization Remarks: Extrapolation method for ROI of variable size and shape Parcel-dependent regularization factor [Risser, Vincent and Ciuciu, sub. to NIPS 08]
59 59/61 Supervised vs unsupervised SMM unsup. =0.1 =0.7
60 60/61 fmri analysis with PyHRF t t t PyHRF SPM = Parcellation
61 61/61 fmri analysis with PyHRF t t t HRF PyHRF SPM = contrasts Parcellation
62 62/61 Implemented model Supervised spatial mixture modeling for neural responses White noise Polynomial/cosine low frequency drift (no output, integrated out of the problem) Gausian HRF which covariance is constraint by its second derivative (smooth constraint) Homoscedastic multisession modeling
Analysis of fmri data within Brainvisa Example with the Saccades database
Analysis of fmri data within Brainvisa Example with the Saccades database 18/11/2009 Note : All the sentences in italic correspond to informations relative to the specific dataset under study TP participants
More informationEvent-related design efficiency and How to plot fmri time series. Sepideh Sadaghiani NeuroSpin Methods Meeting 15. Sept. 2008
Event-related design efficiency and How to plot fmri time series Sepideh Sadaghiani NeuroSpin Methods Meeting 15. Sept. 2008 Event-related averaging or FIR? event-related fmri ability to average responses
More informationFirst-level fmri modeling
First-level fmri modeling Monday, Lecture 3 Jeanette Mumford University of Wisconsin - Madison What do we need to remember from the last lecture? What is the general structure of a t- statistic? How about
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 informationData Visualisation in SPM: An introduction
Data Visualisation in SPM: An introduction Alexa Morcom Edinburgh SPM course, April 2010 Centre for Cognitive & Neural Systems/ Department of Psychology University of Edinburgh Visualising results remembered
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 informationData Visualisation in SPM: An introduction
Data Visualisation in SPM: An introduction Alexa Morcom Edinburgh SPM course, April 2015 SPMmip [-30, 3, -9] 3 Visualising results remembered vs. fixation contrast(s) < < After the results table - what
More informationTutorial BOLD Module
m a k i n g f u n c t i o n a l M R I e a s y n o r d i c B r a i n E x Tutorial BOLD Module Please note that this tutorial is for the latest released nordicbrainex. If you are using an older version please
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 informationA fully Bayesian approach to the parcel-based detection-estimation of brain activity in fmri.
A fully Bayesian approach to the parcel-based detection-estimation of brain activity in fmri. Salima Makni, Jérôme Idier, Thomas Vincent, Bertrand Thirion, Ghislaine Dehaene-Lambertz, Philippe Ciuciu To
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 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 informationFunction-Structure Integration in FreeSurfer
Function-Structure Integration in FreeSurfer Outline Function-Structure Integration Function-Structure Registration in FreeSurfer fmri Analysis Preprocessing First-Level Analysis Higher-Level (Group) Analysis
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 informationSingle Subject Demo Data Instructions 1) click "New" and answer "No" to the "spatially preprocess" question.
(1) conn - Functional connectivity toolbox v1.0 Single Subject Demo Data Instructions 1) click "New" and answer "No" to the "spatially preprocess" question. 2) in "Basic" enter "1" subject, "6" seconds
More informationFunctional MRI data preprocessing. Cyril Pernet, PhD
Functional MRI data preprocessing Cyril Pernet, PhD Data have been acquired, what s s next? time No matter the design, multiple volumes (made from multiple slices) have been acquired in time. Before getting
More informationVersion. Getting Started: An fmri-cpca Tutorial
Version 11 Getting Started: An fmri-cpca Tutorial 2 Table of Contents Table of Contents... 2 Introduction... 3 Definition of fmri-cpca Data... 3 Purpose of this tutorial... 3 Prerequisites... 4 Used Terms
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 informationFMRI Pre-Processing and Model- Based Statistics
FMRI Pre-Processing and Model- Based Statistics Brief intro to FMRI experiments and analysis FMRI pre-stats image processing Simple Single-Subject Statistics Multi-Level FMRI Analysis Advanced FMRI Analysis
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 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 informationMULTIVARIATE ANALYSES WITH fmri DATA
MULTIVARIATE ANALYSES WITH fmri DATA Sudhir Shankar Raman Translational Neuromodeling Unit (TNU) Institute for Biomedical Engineering University of Zurich & ETH Zurich Motivation Modelling Concepts Learning
More informationExperimental Design for fmri OHBM Advanced fmri Educational Course Thomas Liu UCSD Center for Func6onal MRI
Experimental Design for fmri OHBM Advanced fmri Educational Course 2014 Thomas Liu UCSD Center for Func6onal MRI Experimental Design Design 1 Condition 1 Condition 2 Condition 3 Design 2 Why worry about
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 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 informationNA-MIC National Alliance for Medical Image Computing fmri Data Analysis
NA-MIC fmri Data Analysis Sonia Pujol, Ph.D. Wendy Plesniak, Ph.D. Randy Gollub, M.D., Ph.D. Acknowledgments NIH U54EB005149 Neuroimage Analysis Center NIH P41RR013218 FIRST Biomedical Informatics Research
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
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 information3488 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 53, NO. 9, SEPTEMBER 2005
3488 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 53, NO. 9, SEPTEMBER 2005 Joint Detection-Estimation of Brain Activity in Functional MRI: A Multichannel Deconvolution Solution Salima Makni, Philippe
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 informationSPM8 for Basic and Clinical Investigators. 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 informationMulti-voxel pattern analysis: Decoding Mental States from fmri Activity Patterns
Multi-voxel pattern analysis: Decoding Mental States from fmri Activity Patterns Artwork by Leon Zernitsky Jesse Rissman NITP Summer Program 2012 Part 1 of 2 Goals of Multi-voxel Pattern Analysis Decoding
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 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 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 informationAFNI Preprocessing: Outline, Recommendations, and New(ish) Stuff. Robert W Cox SSCC / NIMH & NINDS / NIH / DHHS / USA / EARTH
AFNI Preprocessing: Outline, Recommendations, and New(ish) Stuff Robert W Cox SSCC / NIMH & NINDS / NIH / DHHS / USA / EARTH HBM 2016 As a work of a US Government official, this presentation is not copyrighted
More informationJournal of Articles in Support of The Null Hypothesis
Data Preprocessing Martin M. Monti, PhD UCLA Psychology NITP 2016 Typical (task-based) fmri analysis sequence Image Pre-processing Single Subject Analysis Group Analysis Journal of Articles in Support
More informationGLM for fmri data analysis Lab Exercise 1
GLM for fmri data analysis Lab Exercise 1 March 15, 2013 Medical Image Processing Lab Medical Image Processing Lab GLM for fmri data analysis Outline 1 Getting Started 2 AUDIO 1 st level Preprocessing
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 informationNEURO M203 & BIOMED M263 WINTER 2014
NEURO M203 & BIOMED M263 WINTER 2014 MRI Lab 2: Neuroimaging Connectivity Lab In today s lab we will work with sample diffusion imaging data and the group averaged fmri data collected during your scanning
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 informationSPM Course! Single Subject Analysis
SPM Course! Single Subject Analysis Practical Session Dr. Jakob Heinzle & Dr. Frederike Petzschner & Dr. Lionel Rigoux Hands up: Who has programming experience with Matlab? Who has analyzed an fmri experiment
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 informationEfficiency and design optimization
Efficiency and design optimization Tuesday, Lecture 3 Jeanette Mumford University of Wisconsin - Madison Thanks to Tom Liu for letting me use some of his slides! What is the best way to increase your power?
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 informationSurface fmri data processing using BrainVISA (March 4th 2011)
Surface fmri data processing using BrainVISA (March 4th 2011) Index 1. Introduction and general goal 2. Running the right pipeline 3. Practical details for computing surface processing 4. Visualization
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 informationBayesian 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 informationDetecting Brain Activations in Functional Magnetic Resonance Imaging (fmri) Experiments with a Maximum Cross-Correlation Statistic
Journal of Data Science 10(2012), 403-418 Detecting Brain Activations in Functional Magnetic Resonance Imaging (fmri) Experiments with a Maximum Cross-Correlation Statistic Kinfemichael Gedif 1, William
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 informationIntroduction to fmri. Pre-processing
Introduction to fmri Pre-processing Tibor Auer Department of Psychology Research Fellow in MRI Data Types Anatomical data: T 1 -weighted, 3D, 1/subject or session - (ME)MPRAGE/FLASH sequence, undistorted
More informationSPM99 fmri Data Analysis Workbook
SPM99 fmri Data Analysis Workbook This file is a description of the steps needed to use SPM99 analyze a fmri data set from a single subject using a simple on/off activation paradigm. There are two parts
More informationImage Registration + Other Stuff
Image Registration + Other Stuff John Ashburner Pre-processing Overview fmri time-series Motion Correct Anatomical MRI Coregister m11 m 21 m 31 m12 m13 m14 m 22 m 23 m 24 m 32 m 33 m 34 1 Template Estimate
More informationfmri Image Preprocessing
fmri Image Preprocessing Rick Hoge, Ph.D. Laboratoire de neuroimagerie vasculaire (LINeV) Centre de recherche de l institut universitaire de gériatrie de Montréal, Université de Montréal Outline Motion
More informationAnalysis of group fmri datasets: Voxel-based techniques and beyond
Jirfni 2009, Marseille Analysis of group fmri datasets: Voxel-based techniques and beyond Bertrand Thirion INRIA Saclay Ile de France, PARIETAL team, Neurospin bertrand.thirion@inria.fr 1/80 Outline Background
More informationExploratory fmri analysis without spatial normalization
Exploratory fmri analysis without spatial normalization The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher
More informationCONN fmri functional connectivity toolbox
CONN fmri functional connectivity toolbox Gabrieli Lab. McGovern Institute for Brain Research Massachusetts Institute of Technology http://www.nitrc.org/projects/conn Susan Whitfield-Gabrieli Alfonso Nieto-Castanon
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 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 informationfmri Analysis Sackler Ins2tute 2011
fmri Analysis Sackler Ins2tute 2011 How do we get from this to this? How do we get from this to this? And what are those colored blobs we re all trying to see, anyway? Raw fmri data straight from the scanner
More informationSurface-based Analysis: Inter-subject Registration and Smoothing
Surface-based Analysis: Inter-subject Registration and Smoothing Outline Exploratory Spatial Analysis Coordinate Systems 3D (Volumetric) 2D (Surface-based) Inter-subject registration Volume-based Surface-based
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 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 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 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 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 informationfmri Basics: Single Subject Analysis
fmri Basics: Single Subject Analysis This session is intended to give an overview of the basic process of setting up a general linear model for a single subject. This stage of the analysis is also variously
More informationFMRI data: Independent Component Analysis (GIFT) & Connectivity Analysis (FNC)
FMRI data: Independent Component Analysis (GIFT) & Connectivity Analysis (FNC) Software: Matlab Toolbox: GIFT & FNC Yingying Wang, Ph.D. in Biomedical Engineering 10 16 th, 2014 PI: Dr. Nadine Gaab Outline
More informationEmpirical analyses of null-hypothesis perfusion FMRI data at 1.5 and 4 T
NeuroImage 19 (2003) 1449 1462 www.elsevier.com/locate/ynimg Empirical analyses of null-hypothesis perfusion FMRI data at 1.5 and 4 T Jiongjiong Wang, a,b, * Geoffrey K. Aguirre, a,c Daniel Y. Kimberg,
More informationNonparametric Mean Shift Functional Detection in the Functional Space for Task and Resting-state fmri
Nonparametric Mean Shift Functional Detection in the Functional Space for Task and Resting-state fmri Jian Cheng 1,2, Feng Shi 3, Kun Wang 1, Ming Song 1, Jiefeng Jiang 1, Lijuan Xu 1, Tianzi Jiang 1 1
More informationSelected topics on fmri image processing
fmri Symposium, 24 January, 2006, Ghent Selected topics on fmri image processing Jan Sijbers The fmri folks of the Vision Lab UZA Bio-Imaging Lab Vision Lab Univ. Maastricht T.U. Delft 1 Overview generalized
More informationBESA Research. CE certified software package for comprehensive, fast, and user-friendly analysis of EEG and MEG
BESA Research CE certified software package for comprehensive, fast, and user-friendly analysis of EEG and MEG BESA Research choose the best analysis tool for your EEG and MEG data BESA Research is the
More informationGetting Started Guide
Getting Started Guide Version 2.5 for BVQX 1.9 Rainer Goebel, Henk Jansma and Jochen Seitz Copyright 2007 Brain Innovation B.V. 2 Contents Preface...4 The Objects Tutorial...5 Scanning session information...5
More informationINDEPENDENT COMPONENT ANALYSIS APPLIED TO fmri DATA: A GENERATIVE MODEL FOR VALIDATING RESULTS
INDEPENDENT COMPONENT ANALYSIS APPLIED TO fmri DATA: A GENERATIVE MODEL FOR VALIDATING RESULTS V. Calhoun 1,2, T. Adali, 2 and G. Pearlson 1 1 Johns Hopkins University Division of Psychiatric Neuro-Imaging,
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 informationEstimating Noise and Dimensionality in BCI Data Sets: Towards Illiteracy Comprehension
Estimating Noise and Dimensionality in BCI Data Sets: Towards Illiteracy Comprehension Claudia Sannelli, Mikio Braun, Michael Tangermann, Klaus-Robert Müller, Machine Learning Laboratory, Dept. Computer
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 informationfmri pre-processing Juergen Dukart
fmri pre-processing Juergen Dukart Outline Why do we need pre-processing? fmri pre-processing Slice time correction Realignment Unwarping Coregistration Spatial normalisation Smoothing Overview fmri time-series
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 informationStatistical Methods in functional MRI. False Discovery Rate. Issues with FWER. Lecture 7.2: Multiple Comparisons ( ) 04/25/13
Statistical Methods in functional MRI Lecture 7.2: Multiple Comparisons 04/25/13 Martin Lindquist Department of iostatistics Johns Hopkins University Issues with FWER Methods that control the FWER (onferroni,
More informationOutline. Data integration/fusion: Previous Work. Data Fusion vs. Data Integration
Outline / data Three Approaches to integration/fusion Prediction Constraints 2 nd Level Fusion Conclusions Approaches to -MRI data integration Data Integration through: (i) Prediction some features of
More informationFSL Pre-Processing Pipeline
The Art and Pitfalls of fmri Preprocessing FSL Pre-Processing Pipeline Mark Jenkinson FMRIB Centre, University of Oxford FSL Pre-Processing Pipeline Standard pre-processing: Task fmri Resting-state fmri
More informationPixels to Voxels: Modeling Visual Representation in the Human Brain
Pixels to Voxels: Modeling Visual Representation in the Human Brain Authors: Pulkit Agrawal, Dustin Stansbury, Jitendra Malik, Jack L. Gallant Presenters: JunYoung Gwak, Kuan Fang Outlines Background Motivation
More informationBasic principles of MR image analysis. Basic principles of MR image analysis. Basic principles of MR image analysis
Basic principles of MR image analysis Basic principles of MR image analysis Julien Milles Leiden University Medical Center Terminology of fmri Brain extraction Registration Linear registration Non-linear
More informationDecoding the Human Motor Cortex
Computer Science 229 December 14, 2013 Primary authors: Paul I. Quigley 16, Jack L. Zhu 16 Comment to piq93@stanford.edu, jackzhu@stanford.edu Decoding the Human Motor Cortex Abstract: A human being s
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 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 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 informationLab 5: Multi Subject fmri
USA SPM Short Course April 6-8, 2005 Yale University Lab 5: Multi Subject fmri Rik Henson mailto:rik.henson@mrc-cbu.cam.ac.uk MRC Cognition & Brain Sciences Unit, Cambridge, UK CONTENTS Goals Introduction
More informationBrainVoyager TM. Getting Started Guide. Version 3.0. for BV 21
BrainVoyager TM Getting Started Guide Version 3.0 for BV 21 Rainer Goebel, Henk Jansma, Caroline Benjamins, Judith Eck, Hester Breman and Armin Heinecke Copyright 2018 Brain Innovation B.V. Contents About
More informationVisualization of group inference data in functional neuroimaging
Visualization of group inference data in functional neuroimaging Jan Gläscher Division of Humanities and Social Sciences California Institute of Technology 1200 E California Blvd, Pasadena, CA 91125 phone:
More informationMarkov chain Monte Carlo methods
Markov chain Monte Carlo methods (supplementary material) see also the applet http://www.lbreyer.com/classic.html February 9 6 Independent Hastings Metropolis Sampler Outline Independent Hastings Metropolis
More informationfmri Data Analysis Techniques and the Self-Organizing Maps Approach
fmri Data Analysis Techniques and the Self-Organizing Maps Approach Tiago José de Oliveira Jordão Instituto Superior Técnico Pólo do Taguspark, Porto Salvo, Portugal Functional Magnetic Resonance Imaging
More information- Graphical editing of user montages for convenient data review - Import of user-defined file formats using generic reader
Data review and processing Source montages and 3D whole-head mapping Onset of epileptic seizure with 3D whole-head maps and hemispheric comparison of density spectral arrays (DSA) Graphical display of
More informationSegmenting Lesions in Multiple Sclerosis Patients James Chen, Jason Su
Segmenting Lesions in Multiple Sclerosis Patients James Chen, Jason Su Radiologists and researchers spend countless hours tediously segmenting white matter lesions to diagnose and study brain diseases.
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 informationIndependent Component Analysis of fmri Data
Independent Component Analysis of fmri Data Denise Miller April 2005 Introduction Techniques employed to analyze functional magnetic resonance imaging (fmri) data typically use some form of univariate
More informationUnderstanding multivariate pattern analysis for neuroimaging applications
Understanding multivariate pattern analysis for neuroimaging applications Maria Giulia Preti Dimitri Van De Ville Medical Image Processing Lab, Institute of Bioengineering, Ecole Polytechnique Fédérale
More informationPattern Recognition for Neuroimaging Toolbox: PRoNTo
Click to edit Master title style Pattern Recognition for Neuroimaging Toolbox: PRoNTo Jessica Schrouff PRNI 2018 June 14 th NUS, Singapore Click Outline to edit Master title style PRoNTo s goals and history
More informationMR IMAGE SEGMENTATION
MR IMAGE SEGMENTATION Prepared by : Monil Shah What is Segmentation? Partitioning a region or regions of interest in images such that each region corresponds to one or more anatomic structures Classification
More information