Multivariate Analysis of fmri Group Data Using Independent Vector Analysis
|
|
- Dustin Goodwin
- 5 years ago
- Views:
Transcription
1 Multivariate Analysis of fmri Group Data Using Independent Vector Analysis Jong-Hwan Lee 1, Te-Won Lee 2, Ferenc A. Jolesz 1, and Seung-Schik Yoo 1,3 1 Department of Radiology, Brigham and Women s Hospital, Harvard Medical School, MA, USA jhlee@bwh.harvard.edu 2 Institute of Neural Computation, University of California at San Diego, La Jolla, CA, USA 3 Dept. of BioSystems, Korea Advanced Insitute of Science and Techonology, Daejeon, Korea Abstract. A multivariate non-parametric approach for the processing of fmri group data is important to address variability of hemodynamic responses across subjects, sessions, and brain regions. Independent component analysis (ICA) has a limitation during the inference of group effects due to a permutation problem of independent components. In order to address this limitation, we present an independent vector analysis (IVA) for the processing of fmri group data. Compared to the ICA, the IVA offers an extra dimension for the dependent parameters, which can be assigned for the automated grouping of dependent activation patterns across subjects. The IVA was applied to the fmri data obtained from 12 subjects performing a left-hand motor task. In comparison with conventional univariate methods, IVA successfully characterized the group-representative activation time courses (as component vectors) without extra data processing schemes to circumvent the permutation problem, while effectively detecting the areas with hemodynamic responses deviating from canonical, model-driven ones. Keywords: Independent Vector Analysis, fmri, Neuroimaging, Group Study, Multivariate Analysis. 1 Introduction Functional MRI (fmri) measures the blood-oxygenation-level-dependent (BOLD) signal changes associated with neural activity. Thus, temporal dynamics of the BOLD signal (time series; TS), called as hemodynamic response function (HRF), is the key element in analyzing fmri data. Typically, univariate approaches such as the generalized linear model (GLM) or regression analysis, are performed to estimate the conformity of a measured voxel-wise BOLD TS to the fixed, canonical HRF [1]. However, the BOLD TS may not be fully appreciated from the univariate methods due to the variations across subjects, scans, and brain regions [2]. ICA [3], as one of the multivariate approaches, has provided flexibility in data processing compared to the hypothesis-driven, univariate methods. Such flexibility applies especially when observed hemodynamic responses deviated from expected (hypothesized) HRF [4]. Task-related activation components, often similar in their M.E. Davies et al. (Eds.): ICA 2007, LNCS 4666, pp , Springer-Verlag Berlin Heidelberg 2007
2 634 J.-H. Lee et al. spatial/temporal patterns across subjects/sessions, may not be inherently generalized from individual to group level analysis since the ICA algorithm permutes the order of output components. Therefore, the task-related components across subjects/sessions were manually inspected and grouped in the previous study [5]. This method may require careful selection of the component-of-interest from the large number of subjects/sessions. In the present study, we propose a novel fmri analysis method for group processing based on independent vector analysis (IVA) [6]. IVA was originally proposed to address the limitation of the conventional ICA approach during the blind source/signal separation (BSS) in the frequency domain (i.e. permutation of extracted independent components across frequency bins). IVA correctly indexed the independent components (ICs) that were identified during the BSS in the frequency domain by utilizing the mutual dependency among the extracted ICs across frequency bins. Intuitively, the IVA model offers an extra dimension for processing dependent components, compared to the ICA model. In the fmri study, this extra dimension can be assigned for automated grouping of similar IC maps across subjects. Fundamentally, IVA is an extension of ICA, whereby the component of an input and an output stage forms a vector (instead of a scalar value as in the case of ICA). IVA assumes and, therefore, attempts to increase independency across output vector components, while maintaining dependency among scalar elements within each vector. The dependency in the fmri study is analogous to mutual activation patterns across the subjects, comparable to the group trend in similar spatial activation patterns. Using the IVA algorithm, the spatially-similar trend in activation maps across subjects (dependent, thus representing group-trend) can be derived as the output vector components. As a result, one can avoid the complications of manual selection of task-related IC maps/time-courses (TCs) across subjects, thus rendering the whole process user-independent. In order to show the utility of the proposed method, we implemented and applied the IVA algorithm to analyze fmri data of a left-hand motor clenching task using a short-time trial-based paradigm design. The obtained result was compared with the result from the generalized linear model (GLM) in SPM2 (Wellcome Department of Imaging Neuroscience, University College London, London, UK; 2 Methods 2.1 IVA Model and Learning Algorithm Figure 1 shows a schematic diagram of concurrent synthesis (generative) and analysis models for group data processing using IVA. In the synthesis model, the weighting values at the v th voxels associated with IC maps (assuming mutual dependence across subjects) were grouped into a single vector array assuming the spatial similarity across the subjects. Through the mixing matrix A (TCs represented by each column), the arrays of these vectors are linearly combined to form sets of other vector arrays (equal to the number of temporal points of the BOLD TS), each containing the measured BOLD signal from a specific (v th ) voxel location across all of the subjects. In the analysis model, the unknown weighting values (also associated with IC map) can be estimated via a corresponding unmixing matrix W. The matrix W can be
3 Multivariate Analysis of fmri Group Data Using Independent Vector Analysis 635 obtained by applying a learning rule to increase independence among output vector arrays (thus, sorting out the activation patterns). This is accompanied by maintaining dependence of weighting values within each output vector array, thus deriving mutual activations across the subjects for group inference. Fig. 1. The schematic diagrams of synthesis and analysis models for fmri group data pro- (1) ( M ) (1) ( M ) = a, L, a w = w, L, w. cessing using independent vector analysis. a [ ] and [ ] ji From the synthesis model in Fig. 1, the measured BOLD TS at the v th voxel of subject m can be represented as ji ji x = A c. (1) Each subject has its own mixing matrix (i.e. TCs) with independent activations (i.e. weighting values associated with IC maps) for the measured BOLD TS and does not share a mixing matrix with other subjects. In addition, the dependence among weighting values across subjects, within each unknown vector component, is assumed by the multivariate probability density function (p.d.f.) (1) ( M ) p( c i ) ( = p( ci, L, ci ) ). The boldface and lightface represent a vector/matrix and scalar value, respectively. The superscript and subscript denote the indices of the subjects and the IC map/fmri volumes, respectively. M is the number of subjects ( m { 1, L, M }), N is the number of fmri volume acquisitions corresponding to the unknown IC maps ( i, j, k { 1, L, N} ), and V is the number of voxels within a brain region ( v { 1, L,V} ). The number of unknown IC maps was assumed to be the same as the number of volume acquisitions in Fig. 1. The assumed number of IC maps can be further reduced using a dimension reduction scheme such as principle component analysis (PCA) [3,5]. Then, by applying an unmixing matrix in the analysis model (i.e. inverse of TCs), the weighting values at the v th voxel (across IC maps) of subject m can be estimated as ˆ ) ( m c = W x. (2) kj kj kj
4 636 J.-H. Lee et al. In order to derive learning algorithm, as presented in [6], the Kullback-Leibler (KL) divergence was adopted as a measure of independence of output vector components. Also, variance-dependent multivariate p.d.f. was utilized as a measure of dependence among elements within output vector component. By following the procedure in [6], (m) the algorithm for calculating an update term, Δ W corresponding to the unmixing matrix of the subject m, can be derived as T ( ) [ I ϕ( cˆ )( cˆ ) ] W, m Δ W (3) [ ] T where I is an identity matrix ( N N ), (m) ϕ( c ˆ ) = ϕ( cˆ 1 ) L ϕ( cˆ N ), and M ( l) 2 ϕ ( cˆ ( )) = ˆ ( ) = ˆ k v ck v c ( ). l 1 k v Applying Eq. (3), the unmixing matrix can be iteratively updated for the data from all voxels ( v = 1, L, V ). The only difference compared to the Infomax-based ICA (for the processing of a single subject data [3]) is the application of a nonlinear function ϕ ( cˆ ) (c.f. score function in ICA), which is dependent across subjects in IVA. 2.2 Application to Trial-Based fmri Data This study was approved by the local Institutional Review Board. Twelve righthanded subjects (aged , 5 females) performed one session of a left hand (LH) clenching (2 times/sec) task based on a short-time trial-based paradigm design (65-sec duration excluding 10-sec of dummy scans; task onset occurred at 15-sec followed by a 3-sec task-period). For the start/end of the task, a pre-recorded sound cue was played to the subject in the MRI system via an auditory headset (Avotec, Jensen Beach, FL). The fmri data was obtained in a 3-Tesla clinical scanner (Signa VH, GE Medical Systems) using a single channel, standard birdcage, head coil. To obtain functional data, an EPI sequence was applied to image most of the brain volume (13 axial slices, flip angle=80, TE/TR=40/1000msec, 64 frequency and phase encoding: in-plane voxels, 5mm thickness with a 1mm gap, 240mm square field-of-view) for detection of the BOLD TS associated with neural activity. Prior to group processing, individual EPI data was standardized to the MNI (Montreal Neurological Institute) space by following preprocessing steps in SPM2 (i.e. in order: slice timing correction, realignment, normalization, and spatial smoothing using an 8mm full-width-at-half-maximum 3-D Gaussian kernel). Before processing using IVA algorithm, a PCA-based dimension reduction scheme [3,5] was applied to reduce the number of IC maps/tcs to 50. The sum of corresponding 50 eigenvalues was more than 99% of a sum of total 65 eigenvalues for each subject. Using the IVA algorithm in Eq. (3) to a normalized set of dimension-reduced fmri group data, a semi-batch learning scheme [3] was applied to update the unmixing matrix of each subject. A batch of a mm 3 isotropic cluster (5 5 5=125 voxels due to the 2 2 2mm 3 isotropic voxel) was used, assuming dependencies of neural activations within this cluster across subjects. The learning rate (η ) was set to 10-3 throughout iterations. The iteration was stopped when a ratio of weight change ( ) ( ( ) m η ΔW m W ) was stabilized ( ~ ). After the algorithm converged, the resulting IC map (i.e. weighting values of TC) was transformed into a z-scored map by subtracting a mean value and dividing by the standard deviation [3].
5 Multivariate Analysis of fmri Group Data Using Independent Vector Analysis 637 Sign ambiguity of IC maps/tcs that typically arise from ICA-based methods [3, 4, 5, 7] also applied to the case of IVA. In IVA, the voxel-wise correlation coefficients between (1) the IC z-map within the activated regions ( z >threshold) and (2) the original fmri volumes were utilized. If the averaged (across time points) value of the correlation coefficients was negative, the sign of the IC z-map (& corresponding TC) was inverted. In order to find task-related components, a spatial sorting scheme was employed, whereby a voxel-wise correlation value between the sign-corrected IC z-map and cross-correlation (CC) map (obtained from the temporal correlation coefficients using the task-related canonical HRF) was regarded as the degree of task relevance of each IC z-map for each subject. This similarity measure (i.e. the voxel-wise correlation value) was then averaged across all subjects for each output vector array of IVA and subsequently sorted in descending order. The resulting five most task-related IC z-maps across subjects were processed using one-sample t-test implemented in SPM2 by considering a random effect (RFX) model [8]. These resulting task-related group activation maps of IVA were compared to the group activation map of GLM obtained from SPM2. After obtaining a group inference, the group activation maps were qualitatively compared in terms of location of activation. The areas with an activation volume greater than 5 5 5mm 3 (p<10-3 ) were identified and labeled from the Brodmann s Area (BA) and Automated Anatomical Labeling (AAL) templates (provided by MRIcro; 3 Results Figure 2 shows the task-related group activation maps obtained from GLM and IVA. In the results of GLM (Fig.2A), which is a univariate approach, a single task-related group activation map was acquired. In the results of IVA, the 5 most task-related components were selected after the reordering of output components based on the spatial sorting strategy explained in Section 2.2. The group activations (p<10-2 ~ p<10-7 ) were coded with a color gradient. The labeled anatomical areas of the group activations were listed in Table 1. From the analysis by GLM, the directly task-related areas (e.g. right primary motor area: M1, supplementary motor area: SMA, primary sensory area: S1, and cingulate gyrus) and paradigm-related areas (auditory; superior temporal) showed more significant activations compared to basal ganglia (caudate/putamen/pellidum) and thalamus (see Fig. 2A). In the results of IVA, the task-related activations (e.g. the right M1, SMA, S1, cingulate gyrus, and sup./mid. temporal gyrus) were extracted as the 1 st task-related group activation map ĉ 1. The group activations in the remaining components were dominant in the primary auditory area ( ĉ ), basal ganglia & thalamus ( 2 ĉ 3 ), inf. frontal & insular cortex ( ĉ ), and mid. 4 frontal & cingulate gyrus ( ĉ 5 ). From the comparison of group activation maps between GLM and IVA (Fig. 2), some of the activations revealed by IVA were underestimated by GLM. For example, the size of activated regions obtained in the auditory area by GLM (examples were shown with green circles in Fig. 2A) was reduced compared to the area detected by IVA (marked as green in Fig. 2B). Also, the activations in the basal ganglia/thalamus
6 638 J.-H. Lee et al. identified from IVA (blue area in Fig. 2B) were not detected when processed with GLM (examples were shown with blue circles in Fig. 2A). Fig. 2. Group activation maps obtained by (A) GLM and (B) IVA methods. For IVA, five most task-related component maps ( cˆ 1 ~ cˆ 5 ) were color-coded. Table 1. The strongly activated cortical areas inside group activation map. If any activated cluster (p<10-3) was bigger than 5 5 5mm3, the center of this cluster was registered using BA and AAL indices (provided by MRIcro; GLM Left-Hemisphere Inf. frontal / SMA / Insula / Mid. Cingulate / S1 / SupraMarginal / Sup. Temp. *BA: 6, 23, 24, 32, 42, 47, 48 Right-Hemisphere M1 / Inf. Frontal / SMA / Insula / Mid. Cingulate / S1 / Parietal (Sup., Inf) / SupraMarginal / Sup. Temp. *BA: 2, 3, 4, 6, 8, 23, 24, 32, 38, 40, 42, 44, 47, 48 c 1 : SMA / Mid. Cingulate / S1 / SupraMarginal / Temp. (Sup., Mid.) *BA: 3, 6, 8, 22, 23, 24, 32, 41, 42, 48 c 1 : M1 / Frontal (Sup., Mid.) / SMA / Mid. Cingulate / S1 / Parietal (Sup., Inf) / SupraMarginal / Sup. Temp. *BA: 2, 3, 4, 6, 23, 24, 32, 40, 48 c 2 : Insula / SupraMarginal / Heschl / Temporal (Sup., Mid.) *BA: 22, 42, 48 c 3 IVA : Frontal (Sup., Inf.) / Insula / Hippocampus / Amygdala / Putamen / Pallidum / Thalamus / Heschl / Sup.Temp. *BA: 11, 20, 25, 34, 38, 47, 48 c 2 : Heschl / Sup. Temp. *BA: 22, 48 c 3 : Insula / Hippocampus / Amygdala / Caudate / Putamen / Pallidum Thalamus *BA: 11, 20, 27, 34, 37, 38, 48 / c 4 / *BA: 6, 48 : Inf. Frontal / Insula SupraMarginal / Putamen / Heschl *BA: 45, 47, 48 c 5 : Mid. Frontal / Cingulate (Ant., Mid.) c 5 : / Sup. Temp. Sup. Temp. *BA: 24, 32, 48 c 4 : Inf. Frontal / Insula / Putamen *BA: 24, 32, 38, 45, 46 Insula / Cingulate (Ant., Mid.) /
7 Multivariate Analysis of fmri Group Data Using Independent Vector Analysis 639 We conjectured that the differences in activation maps between GLM and IVA were caused by the individual differences in temporal patterns of hemodynamic responses from the areas. Therefore, individual differences in the obtained TCs were further compared across subjects (Note that the TC represents a dominant feature of the BOLD TS corresponding to the activated voxels in the IC map). Figure 3 shows the individual TCs corresponding to three highly task-related group activation maps ĉ in Fig. 2B) by IVA. Here, we adopted the convention introduced by Duann et ( 1 ĉ ~ 3 al., [8], whereby the normalized TC (0~1) was coded in gray scale (black: 0 & white: 1) so that relative amplitudes can be readily discriminated within/across subjects. First of all, the TCs (Fig. 3A) corresponding to the 1 st (highly) task-related group activation map ĉ was in good agreement with the hypothesized HRF in GLM (yellow box: 1 task-related period; correlation coefficient between the averaged TC across subjects and the hypothesized HRF: 0.88). On the other hand, the TCs (Fig. 3B&C) corresponding to ĉ 2 & ĉ 3 showed some degree of variations in peak position within the task-related period. It is also notable that additional peaks during rest-periods were observed (examples are shown with arrows) with reduced correlation coefficients of 0.56 & 0.63, compared to the ĉ. Because of these large variations between the actual 1 hemodynamic responses (analogous to TCs) and the hypothesized HRF across subjects, the activations in the corresponding areas may not be detected by GLM. Fig. 3. Image plots of TCs across subjects corresponding to 3 highly task-related group activation maps ( ĉ ~ 1 ĉ 3 ). Each TC was normalized between 0 (black) and 1 (white). A yellow box indicates the period of task-related response. Green arrows indicate examples of the peaks during the rest-period. The plots in the bottom are the hypothesized HRF (green line), along with the averaged TC (blue line), and standard deviations (red bars) across all subjects. A correlation coefficient between the averaged TC and the hypothesized HRF is shown in the topright corner of each averaged time plot. A task-period (3-sec) is marked with thick black bar. 4 Conclusion In this study, we have proposed the use of IVA to infer the group-activation pattern from fmri data. The IVA algorithm was applied for multiple subjects BOLD signals
8 640 J.-H. Lee et al. and the spatially-similar trend in activation maps across subjects (dependent, thus representing group-trend) were derived as the output vector components. From the application of the proposed method to fmri data, the resulting IC maps/tcs of individual subjects provided reliable task-related information for a further group level inference. In addition, IVA provided more robust activation patterns than GLM (based on the hypothesized univariate HRF), especially when the HRF deviated from the hypothesized HRF (e.g. from the basal ganglia/thalamus). These results show the feasibility of IVA to be used in fmri group studies as potential alternative to conventional univariate approach. The proposed model may also be adopted to find out multivariate common activation patterns across multiple trials/sessions from a single subject data by substituting the index of a subject for the index of a trial/session. The IVA approach demands computational load since an individual unmixing matrix is iteratively trained using the results from other subjects (represented as nonlinear function, ϕ ( cˆ )), and thus, all of the unmixing matrices should be parallely updated. This increased computational demand can be alleviated by the increasing hardware memory. In order to achieve the fully-automated group processing using IVA, elaborate sets of optimization in terms of learning parameters and sorting schemes (for the selection of task-related features) are needed. Acknowledgments. This work was partially supported in part by grants from NIH (R01-NS to Yoo, SS and NIH U41RR to Jolesz FA). References 1. Worsley, K.J., Friston, K.J.: Analysis of fmri time-series revisited again. Neuroimage. 2, (1995) 2. Aguirre, G.K., Zarahn, E., D Esposito, M.: The variability of human, BOLD hemodynamic responses. NeuroImage 8, (1998) 3. McKeown, M.J., Makeig, S., Brown, G.G., Jung, T.P., Kindermann, S.S., Bell, A.J., Sejnowski, T.J.: Analysis of fmri data by blind separation into independent spatial components. Hum. Brain Mapp. 6, (1998) 4. Quigley, M.A., Haughton, V.M., Carew, J., Cordes, D., Moritz, C.H., Meyerand, M.E.: Comparison of independent component analysis and conventional hypothesis-driven analysis for clinical functional MR image processing. AJNR Am. J. Neuroradiol. 23, (2002) 5. Calhoun, V.D., Adali, T., McGinty, V.B., Pekar, J.J., Watson, T.D., Pearlson, G.D.: fmri activation in a visual-perception task: network of areas detected using the generalized linear model and independent component analysis. Neuroimage 14, (2001) 6. Kim, T.S., Attias, H.T., Lee, S.Y., Lee, T.W.: Blind source separation exploiting higherorder frequency dependencies. IEEE Trans. Audio, Speech and Language Process. 15, (2007) 7. Duann, J.R., Jung, T.P., Kuo, W.J., Yeh, T.C., Makeig, S., Hsieh, J.C., Sejnowski, T.J.: Single-trial variability in event-related BOLD signals. Neuroimage 15, (2002) 8. Friston, K.J., Holmes, A.P., Worsley, K.J.: How many subjects constitute a study? Neuroimage 10, 1 5 (1999)
Independent vector analysis (IVA): Multivariate approach for fmri group study
www.elsevier.com/locate/ynimg NeuroImage 40 (2008) 86 109 Independent vector analysis (IVA): Multivariate approach for fmri group study Jong-Hwan Lee, a Te-Won Lee, b Ferenc A. Jolesz, a and Seung-Schik
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 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 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 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 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 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 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 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 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 informationThis Week. Last Week. Independent Component Analysis. ICA vs. PCA. ICA vs PCA. Blind Source Separation: The Cocktail Party Problem
1 Last Week Functional Connectivity This Week (Brief) ICA Intro ICA of fmri Sorting/Calibration Validation Demo single-subject subject ICA 2 Blind Source Separation: The Cocktail Party Problem Observations
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 informationAttention modulates spatial priority maps in human occipital, parietal, and frontal cortex
Attention modulates spatial priority maps in human occipital, parietal, and frontal cortex Thomas C. Sprague 1 and John T. Serences 1,2 1 Neuroscience Graduate Program, University of California San Diego
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 informationINDEPENDENT COMPONENT ANALYSIS WITH FEATURE SELECTIVE FILTERING
INDEPENDENT COMPONENT ANALYSIS WITH FEATURE SELECTIVE FILTERING Yi-Ou Li 1, Tülay Adalı 1, and Vince D. Calhoun 2,3 1 Department of Computer Science and Electrical Engineering University of Maryland Baltimore
More informationFigure 1. Comparison of the frequency of centrality values for central London and our region of Soho studied. The comparison shows that Soho falls
A B C Figure 1. Comparison of the frequency of centrality values for central London and our region of Soho studied. The comparison shows that Soho falls well within the distribution for London s streets.
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 informationA Modified Infomax ICA Algorithm for fmri Data Source Separation
Research Journal of Applied Sciences, Engineering and echnology 5(20): 4862-4868, 2013 ISSN: 2040-7459; e-issn: 2040-7467 Maxwell Scientific Organization, 2013 Submitted: December 17, 2012 Accepted: January
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 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 informationX A S. Overview. Modeling the Brain? ICA vs PCA. Independent Component Analysis
fmri Course Lecture 9: Introduction to Independent Component Analysis Vince D. Calhoun, Ph.D. Director, Image Analysis & MR Research The Mind Research Network Associate Professor, Electrical and Computer
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 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 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 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 informationI.e. Sex differences in child appetitive traits and Eating in the Absence of Hunger:
Supplementary Materials I. Evidence of sex differences on eating behavior in children I.e. Sex differences in child appetitive traits and Eating in the Absence of Hunger: Table 2. Parent Report for Child
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 informationThis exercise uses one anatomical data set (ANAT1) and two functional data sets (FUNC1 and FUNC2).
Exploring Brain Anatomy This week s exercises will let you explore the anatomical organization of the brain to learn some of its basic properties, as well as the location of different structures. The human
More informationFmri Spatial Processing
Educational Course: Fmri Spatial Processing Ray Razlighi Jun. 8, 2014 Spatial Processing Spatial Re-alignment Geometric distortion correction Spatial Normalization Smoothing Why, When, How, Which Why is
More informationSpatial Distribution of Temporal Lobe Subregions in First Episode Schizophrenia: Statistical Probability Atlas Approach
Spatial Distribution of Temporal Lobe Subregions in First Episode Schizophrenia: Statistical Probability Atlas Approach Hae-Jeong Park, Martha E. Shenton, James Levitt, Dean Salisbury, Marek Kubicki, Ferenc
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 informationA Data-Driven fmri Neuronal Activation Analysis Method Using Temporal Clustering Technique and an Adaptive Voxel Selection Criterion
A Data-Driven fmri Neuronal Activation Analysis Method Using Temporal Clustering Technique and an Adaptive Voxel Selection Criterion Sarah Lee, Fernando Zelaya, Stephanie A. Amiel and Michael J. Brammer
More informationFunctional Connectivity in the Thalamus and Hippocampus Studied with Functional MR Imaging
AJNR Am J Neuroradiol 21:1397 1401, September 2000 Functional Connectivity in the Thalamus and Hippocampus Studied with Functional MR Imaging Thor Stein, Chad Moritz, Michelle Quigley, Dietmar Cordes,
More informationAppendix E1. Supplementary Methods. MR Image Acquisition. MR Image Analysis
RSNA, 2015 10.1148/radiol.2015150532 Appendix E1 Supplementary Methods MR Image Acquisition By using a 1.5-T system (Avanto, Siemens Medical, Erlangen, Germany) under a program of regular maintenance (no
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 informationGiven a network composed of a set of nodes and edges, shortest path length, dij, is the
SUPPLEMENTARY METHODS Efficiency Given a network composed of a set of nodes and edges, shortest path length, dij, is the minimum number of edges between two nodes i and j, and the efficiency, eij, between
More informationSupplementary methods
Supplementary methods This section provides additional technical details on the sample, the applied imaging and analysis steps and methods. Structural imaging Trained radiographers placed all participants
More informationIndependent Component Analysis of Functional Magnetic Resonance Imaging (fmri) Data: A simple Approach
Research Journal of Applied Sciences, Engineering and Technology 5(24): 5494-552, 213 ISSN: 24-7459; e-issn: 24-7467 Maxwell Scientific Organization, 213 Submitted: August 17, 212 Accepted: September 8,
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 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 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 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 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 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 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 informationIndependent components of face images: A represent at ion for face recognit ion
4th Joint Symposium on Neural Computation Proceedings 1997 Independent components of face images: A represent at ion for face recognit ion Marian Stewart Bartlett Terrence J. Sejnowski University of California
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 informationSuper-resolution Reconstruction of Fetal Brain MRI
Super-resolution Reconstruction of Fetal Brain MRI Ali Gholipour and Simon K. Warfield Computational Radiology Laboratory Children s Hospital Boston, Harvard Medical School Worshop on Image Analysis for
More informationCocozza S., et al. : ALTERATIONS OF FUNCTIONAL CONNECTIVITY OF THE MOTOR CORTEX IN FABRY'S DISEASE: AN RS-FMRI STUDY
ALTERATIONS OF FUNCTIONAL CONNECTIVITY OF THE MOTOR CORTEX IN FABRY'S DISEASE: AN RS-FMRI STUDY SUPPLEMENTARY MATERIALS Sirio Cocozza, MD 1*, Antonio Pisani, MD, PhD 2, Gaia Olivo, MD 1, Francesco Saccà,
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 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 informationFrequencies Contributing to Functional Connectivity in the Cerebral Cortex in Resting-state Data
AJNR Am J Neuroradiol 22:1326 1333, August 2001 Frequencies Contributing to Functional Connectivity in the Cerebral Cortex in Resting-state Data Dietmar Cordes, Victor M. Haughton, Konstantinos Arfanakis,
More informationAdvances in Neural Information Processing Systems, 1999, In press. Unsupervised Classication with Non-Gaussian Mixture Models using ICA Te-Won Lee, Mi
Advances in Neural Information Processing Systems, 1999, In press. Unsupervised Classication with Non-Gaussian Mixture Models using ICA Te-Won Lee, Michael S. Lewicki and Terrence Sejnowski Howard Hughes
More informationA GPU Accelerated Interactive Interface for Exploratory Functional Connectivity Analysis of fmri Data
A GPU Accelerated Interactive Interface for Exploratory Functional Connectivity Analysis of fmri Data Anders Eklund, Ola Friman, Mats Andersson and Hans Knutsson Linköping University Post Print N.B.: When
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 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 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 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 informationModeling the Other Race Effect with ICA
Modeling the Other Race Effect with ICA Marissa Grigonis Department of Cognitive Science University of California, San Diego La Jolla, CA 92093 mgrigoni@cogsci.uscd.edu Abstract Principal component analysis
More informationSupplementary Data. in residuals voxel time-series exhibiting high variance, for example, large sinuses.
Supplementary Data Supplementary Materials and Methods Step-by-step description of principal component-orthogonalization technique Below is a step-by-step description of the principal component (PC)-orthogonalization
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 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 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 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 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 informationAnalysis 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 informationICA mixture models for image processing
I999 6th Joint Sy~nposiurn orz Neural Computation Proceedings ICA mixture models for image processing Te-Won Lee Michael S. Lewicki The Salk Institute, CNL Carnegie Mellon University, CS & CNBC 10010 N.
More informationGroup ICA of FMRI: Introduction and Review of Current Work Vince D. Calhoun, Ph.D.
Group ICA of FMRI: Introduction and Review of Current Work Vince D. Calhoun, Ph.D. Director, Image Analysis & MR Research The Mind Research Network Associate Professor, Electrical and Computer Engineering
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 informationEvaluation of multiple voxel-based morphometry approaches and applications in the analysis of white matter changes in temporal lobe epilepsy
Evaluation of multiple voxel-based morphometry approaches and applications in the analysis of white matter changes in temporal lobe epilepsy Wenjing Li a, Huiguang He a, Jingjing Lu b, Bin Lv a, Meng Li
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 informationIndependent Component Analysis of fmri Data in the Complex Domain
Magnetic Resonance in Medicine 48:180 192 (2002) Independent Component Analysis of fmri Data in the Complex Domain V.D. Calhoun, 1,4 * T. Adalı, 4 G.D. Pearlson, 1 P.C.M. van Zijl, 2,3 and J.J. Pekar 2,3
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 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 informationNeuroimage Processing
Neuroimage Processing Instructor: Moo K. Chung mkchung@wisc.edu Lecture 2-3. General Linear Models (GLM) Voxel-based Morphometry (VBM) September 11, 2009 What is GLM The general linear model (GLM) is a
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 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 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 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 informationTopographic Mapping with fmri
Topographic Mapping with fmri Retinotopy in visual cortex Tonotopy in auditory cortex signal processing + neuroimaging = beauty! Topographic Mapping with fmri Retinotopy in visual cortex Tonotopy in auditory
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 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 informationEffect of Focal and Nonfocal Cerebral Lesions on Functional Connectivity Studied with MR Imaging
AJNR Am J Neuroradiol 22:294 300, February 2001 Effect of Focal and Nonfocal Cerebral Lesions on Functional Connectivity Studied with MR Imaging Michelle Quigley, Dietmar Cordes, Gary Wendt, Pat Turski,
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 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 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 informationAutomatic Optimization of Segmentation Algorithms Through Simultaneous Truth and Performance Level Estimation (STAPLE)
Automatic Optimization of Segmentation Algorithms Through Simultaneous Truth and Performance Level Estimation (STAPLE) Mahnaz Maddah, Kelly H. Zou, William M. Wells, Ron Kikinis, and Simon K. Warfield
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 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 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 informationEVALUATION OF WEIGHT SPARSITY REGULARIZION SCHEMES OF DEEP NEURAL NETWORKS APPLIED TO FUNCTIONAL NEUROIMAGING DATA. Hyun-Chul Kim, Jong-Hwan Lee
EVALUATION OF WEIGHT SPARSITY REGULARIZION SCHEMES OF DEEP NEURAL NETWORKS APPLIED TO FUNCTIONAL NEUROIMAGING DATA Hyun-Chul Kim, Jong-Hwan Lee Department of Brain and Cognitive Engineering, Korea University,
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 informationA Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction
Tina Memo No. 2007-003 Published in Proc. MIUA 2007 A Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction P. A. Bromiley and N.A. Thacker Last updated 13 / 4 / 2007 Imaging Science and
More informationNIH Public Access Author Manuscript Conf Proc IEEE Eng Med Biol Soc. Author manuscript; available in PMC 2009 December 11.
NIH Public Access Author Manuscript Published in final edited form as: Conf Proc IEEE Eng Med Biol Soc. 2004 ; 2: 1521 1524. doi:10.1109/iembs.2004.1403466. An Elliptic PDE Approach for Shape Characterization
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 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 informationVoxel Selection in fmri Data Analysis Based on Sparse Representation
IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 56, NO. 10, OCTOBER 2009 2439 Voxel Selection in fmri Data Analysis Based on Sparse Representation Yuanqing Li, Member, IEEE, Praneeth Namburi, Zhuliang
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationSupplementary Information. Task-induced brain state manipulation improves prediction of individual traits. Greene et al.
Supplementary Information Task-induced brain state manipulation improves prediction of individual traits Greene et al. Supplementary Note 1 Analyses of effects of gf measurement technique PNC CPM results
More informationCHAPTER 9: Magnetic Susceptibility Effects in High Field MRI
Figure 1. In the brain, the gray matter has substantially more blood vessels and capillaries than white matter. The magnified image on the right displays the rich vasculature in gray matter forming porous,
More information