Statistical Methods in functional MRI. False Discovery Rate. Issues with FWER. Lecture 7.2: Multiple Comparisons ( ) 04/25/13
|
|
- Sara Parsons
- 5 years ago
- Views:
Transcription
1 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, RF, Permutation ests) provide a strong control over the number of false positives. While this is appealing, the resulting thresholds often lead to tests that suffer from low power. Power is critical in fmri applications because the most interesting effects are usually at the edge of detection. False Discovery Rate he false discovery rate (FDR) is a recent development in multiple comparison problems due to enjamini and Hochberg (1995). While the FWER controls the probability of any false positives, the FDR controls the proportion of false positives among all rejected tests. Notation Suppose we perform tests on m voxels. Declared Inactive Declared ctive ruly inactive U V m 0 ruly active S m-m 0 m-r R m In this notation: Definitions ( ) FWER = P V 1 False discovery rate: V FDR = E R U, V, and S are unobservable random variables. R is an observable random variable. he FDR is defined to be 0 if R=0. 1
2 Properties procedure controlling the FDR ensures that on average the FDR is no bigger than a prespecified rate q which lies between 0 and 1. H Procedure 1. Select desired limit q on FDR (e.g., 0.05) 2. Rank p-values, p (1) p (2)... p (m) However, for any given data set the FDR need not be below the bound. n FDR-controlling technique guarantee controls of the FDR in the sense that FDR q. 3. Let r be largest i such that p (i) i/m q 4. Reject all hypotheses corresponding to p (1),..., p (r). p-value 0 1 p (i) i/m q 0 1 Comments he H procedure is adaptive in the sense that the larger the signal, the lower the threshold. If all null hypothesis are true, the FDR is equivalent to the FWER. Low signal High signal ny procedure that controls the FWER also controls the FDR. procedure that controls the FDR only can be less stringent and lead to a gain in power. q m q Since FDR controlling procedures work only on the p-values and not on the actual test statistics, it can be applied to any valid statistical test. Example Signal + Noise = Signal + Noise α=0.10, No correction Percentage of false positives FWER control at 10% FDR control at 10% FWER Occurrence of false positive Percentage of active voxels that are false positives 2
3 Uncorrected hresholds Most published PE and fmri studies use arbitrary uncorrected thresholds (e.g., p<0.001). With available sample sizes, corrected thresholds are so stringent that power is extremely low. Using uncorrected thresholds is problematic when interpreting conclusions from individual studies, as many activated regions may be false positives. Null findings are hard to disseminate, hence it is difficult to refute false positives established in the literature. Extent hreshold Sometimes an arbitrary extent threshold is used when reporting results. Here a voxel is only deemed truly active if it belongs to a cluster of k contiguous active voxels (e.g., p<0.001, 10 contingent voxels). Unfortunately, this does not necessarily correct the problem because imaging data are spatially smooth and therefore false positives may appear in clusters. Example ctivation maps with spatially correlated noise thresholded at three different significance levels. Due to the smoothness, the false-positive activation are contiguous regions of multiple voxels. α=0.10 α=0.01 α=0.001 Example Similar activation maps using null data. α=0.10 α=0.01 α=0.001 Note: ll images smoothed with FWHM=12mm Note: ll images smoothed with FWHM=12mm Lecture 8: Functional Connectivity 04/25/13 3
4 Data Processing Pipeline rain Networks Experimental Design It has become common practice to talk about brain networks, i.e. sets of interconnected brain regions with information transfer among regions. Data cquisition Reconstruction Preprocessing Slice-time Correction Motion Correction, Co-registration & Normalization Data nalysis Localizing rain ctivity Connectivity o construct a network: Define a set of nodes (e.g., ROIs) Estimate the set of connections, or edges, between the nodes. Spatial Smoothing Prediction C C C Network Methods number of methods have been suggested in the neuroimaging literature to quantify the relationship between nodes/regions. heir appropriateness depend upon: what type of conclusions one is interested in making; what type of assumptions one is willing to make; and the level of the analysis and modality. rain Connectivity Functional Connectivity Undirected association between two or more fmri time series. Makes statements about the structure of relationships among brain regions. VMPFC DLPFC dcc MG rain Connectivity Effective Connectivity Directed influence of one brain region on the physiological activity recorded in other brain regions. Makes statements about causal effects among tasks and regions. V5 Methods: Functional Connectivity Seed analysis Inverse covariance methods Multivariate decomposition methods Principle Components nalysis Independent Components nalysis Partial Least Squares V1 PPC Mediation analysis Psychophysiological interaction (PPI) analysis 4
5 Methods: Effective Connectivity Structural Equation Modeling Granger Causality Dynamic Causal Modeling ayes Net Mediation analysis Psychophysiological interaction (PPI) analysis Levels of nalysis Functional connectivity can be applied at different levels of analysis, with different interpretations at each. Connectivity across time can reveal networks that are dynamically activated across time. Connectivity across trials can identify coherent networks of task related activations. Levels of nalysis ivariate Connectivity Connectivity across subjects can reveal patterns of coherent individual differences. Connectivity across studies can reveal tendencies for studies to co-activate within sets of regions. Simple functional connectivity Region is correlated with Region. Provides information about relationships among regions. Can be performed on time series data within a subject, or individual differences (contrast maps, one per subject). ime Series Connectivity Calculate the cross-correlation between time series from two separate brain regions. Region 1 Region 2 Seed nalysis In seed analysis the cross-correlation is computed between the time course from a predetermined region (seed region) and all other voxels. Subject 1 Subject 2 Subject n Group nalysis his allows researchers to find regions correlated with the activity in the seed region. he seed time course can also be a performance or physiological variable 5
6 Correlations between brain activity and heart-rate VMPFC Issues One of the main problems with time series connectivity is the fact that there may be different hemodynamic lags in different regions: ime series from different regions may not match up, even if neural activity patterns match up. ime (Rs, 2 s) hreshold: p <.005 If lags are estimated from data, temporal order may be caused by vascular (uninteresting) or neural (interesting) response. verage within-subject correlation (r) eta Series eta Series he beta series approach can be used to minimize issues of inter-region neurovascular coupling. Region 1 Region 2 Procedure: Fit a GLM to obtain separate parameter estimates for each individual trial. Compute the correlation between these estimates across voxels. Subject 1 Subject 2 Subject n Group nalysis Individual Differences Partial Correlation Subject Contrast Image Seed Value Group Results 1 x 1 C 2.. x 2 Partial Correlation Correlation between two regions, after the effect of all other regions have been removed. Helps protect against illusory correlations between regions (e.g., and C uncorrelated after controlling for ). N x N 6
7 Inverse Covariance Methods For multivariate normal data there exists a duality between the inverse covariance (precision) matrix and the graph representing relationships between regions. Conditional independence between variables (regions) corresponds to zero entries in the precision matrix. Graphical lasso (GLSSO) can be used to estimate sparse precision matrices and graphs. 0 Σ 1 = C 0 Mediation Mediation he relationship between regions and is mediated by M Can identify functional pathways spanning > 2 regions Can be performed on time series data within a subject, or individual differences (contrast maps, one per subject) lso: est of whether task-related activations in are mediated, or explained, by M. ask M M Demonstrating Mediation Decomposition of Effects Full model, with mediator x a m c b m = i m + ax + e m y = i y + bm + c'x + e y y Reduced model, without mediator x y c y = i y + cx + e y he mediation framework allows us to decompose the total effect of x on y as follows: c = c' + ab otal effect = Direct effect + Mediated effect Does m explain some of the x-y relationship? est c c, which is equivalent to significance of the ab product. Sobel test or bootstrap test. Functional Mediation Pain Data X ( m t), m ( t), m ( )) 1( 2 n t M α (t) β (t) otal: γ Direct: γ ' Y ( x, x, 1 2 x n ) ( y, y, 1 2 y n ) y = γ x + ε ( ) i m ( t) = α ( t) x + ε, ( t) i i i m yi = α ( s) mi ( s) ds + γ ' xi + εi, y ( t) 0 i i, x t γ = γ ' + α( s) β ( s) ds 0 ctivation in the right anterior insula mediates the relationship between temperature and pain rating. α pathway function he key time interval driving the mediation is between seconds following activation. emp rain Response αβ pathway function β pathway function Rating 7
8 Moderation PPI M Moderation he relationship between regions and is moderated by M Connectivity between and depends on state (level) of M Can be performed on time series data within a subject, or individual differences (contrast maps, one per subject) M can be task state or other variable In the psychophysiological interaction (PPI) approach, the standard GLM model is supplemented with additional regressors that model the interaction between the task and the time course in a seed region. Y = β + β X + β R + β R* X + ε In SPM, on time series data: Psychophysiological interaction (PPI). ask ime course from seed region Interaction term PPI PPI can be used to determine whether the correlation between two brain areas is altered by different psychological contexts. he interaction term reflects the modulation of the slope of the linear relationship with the seed voxel depending on the variable used to create the interaction. Decomposition Methods We often use multivariate decomposition methods to study functional connectivity. Provides a decomposition of the data into separate components. Can be used to find coherent brain networks. Provides information on how different brain regions interact with one another. he most common decomposition methods are principal components analysis and independent components analysis. Data Organization Principal Components nalysis hroughout we organize the fmri data in a N matrix X. he row dimension is the number of time points and the column dimension the number of voxels. ime Voxels X Principal components analysis involves finding spatial modes, or eigenimages, in the data. hese are the patterns that account for most of the variance-covariance structure in the data. hey are ranked in order of the amount of variation they explain. he eigenimages can be obtained using singular value decomposition (SVD), which decomposes the data into two sets of orthogonal vectors that correspond to patterns in space and time. 8
9 Using SVD, we can decompose the matrix X as: X = USV where U and V are unitary orthogonal matrices and S is a diagonal matrix consisting of ranked singular values. ime Voxels = ime courses Eigenimages Each column of V defines a distributed brain region that can be displayed as an image (eigenimages). X = USV Each column of U correspond to the time-dependent profiles associated with each eigenimage. X = s1u 1v1 + s2u2v2 + + s N u N v N v 1 v 2 PPROX. = s 1 u 1 OF Y + s PPROX. 2 u OF Y Worsley Independent Components nalysis Independent Components nalysis (IC) is a family of techniques used to extract independent signals from some source signal. IC provides a method to blindly separate the data into spatially independent components. Cocktail Party Problem wo people are talking simultaneously in a room with two microphones. Speakers: s 1 (t) and s 2 (t). Microphones: x 1 (t) and x 2 (t) he key assumption is that the data set consists of p spatially independent components, which are linearly mixed and spatially fixed. x ( t) = a 1 x ( t) = a 2 s ( t) + a 11 1 s ( t) + a 21 1 s ( t) 12 2 s ( t) 22 2 X = S Mixing matrix Source matrix 9
10 ssumptions If the mixing matrix is known, the problem is straight forward. However, IC solves this problem without knowing the mixing parameters. IC Estimation We can find the independent components using a variety of different approaches. Maximizing non-gaussianity Minimizing the mutual information Maximum likelihood estimation Projection pursuit Instead it exploits some key assumptions: Linear mixing of sources. he components s i are statistically independent. he components s i are non-gaussian. IC for fmri IC for fmri It is assumed that the fmri data can be modeled by identifying sets of voxels whose activity both vary together over time and are different from the activity in other sets. Decompose the data set into a set of spatially independent component maps with a set of corresponding time-courses. fmri data Source 1 s = [ s1 s2] + ime course 1 fmri data is assumed to be a linear mixture of statistically independent sources, s. Source 2 ime course 2 Vince Calhoun IC for fmri Overview We seek to decompose X as follows: Voxels ime Courses Spatially independent Components X = S ime = where the matrix S contains statistically independent maps in its rows each with an internally consistent time-course contained in the associated column of the mixing matrix. Data X = S Mixing Matrix Components Use an IC algorithm to find and S. 10
11 Comments Unlike PC which assumes an orthonormality constraint, IC assumes statistical independence among a collection of spatial patterns. Independence is a stronger requirement than orthonormality. However, in IC the spatially independent components are not ranked in order of importance as they are when performing PC. ypes of IC n IC that decomposes the original data into spatially statistically independent components is called spatial IC (sic). It is possible to switch the order and make the temporal dynamics independent. his is called temporal IC (tic). Spatial IC is more common in fmri data analysis. Multi-subject nalysis Using IC to analyze fmri data from multiple subjects raises several questions. How should components be combined across subjects? How should the final results be thresholded and/or presented? here are several approaches: Stack time courses (forces time courses to be the same) Stack images and back-reconstruct (allows time courses to vary, allows some flexibility in images) Stack into a cube (forces images and time courses to be the same) McKeown, et. al. Group IC Group IC Group IC is based on temporal concatenation. It decomposes the group matrix, and estimates through back-reconstruction the spatial weights for each subject for a component of interest. For each subject the spatial weights at each voxel are treated as random variables, and a t- test is used to test whether that voxel loaded significantly on that component in the group. Data Subject 1 X Subject N IC = 1 S_agg N ack-reconstruction 1-1 Subject i i = S i 11
Statistical 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 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 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 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 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 informationControlling for multiple comparisons in imaging analysis. Wednesday, Lecture 2 Jeanette Mumford University of Wisconsin - Madison
Controlling for multiple comparisons in imaging analysis Wednesday, Lecture 2 Jeanette Mumford University of Wisconsin - Madison Motivation Run 100 hypothesis tests on null data using p
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 informationLinear Methods for Regression and Shrinkage Methods
Linear Methods for Regression and Shrinkage Methods Reference: The Elements of Statistical Learning, by T. Hastie, R. Tibshirani, J. Friedman, Springer 1 Linear Regression Models Least Squares Input vectors
More informationUnsupervised Learning
Unsupervised Learning Learning without Class Labels (or correct outputs) Density Estimation Learn P(X) given training data for X Clustering Partition data into clusters Dimensionality Reduction Discover
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 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 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 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 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 informationMultiple Testing and Thresholding
Multiple Testing and Thresholding NITP, 2010 Thanks for the slides Tom Nichols! Overview Multiple Testing Problem Which of my 100,000 voxels are active? Two methods for controlling false positives Familywise
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 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 informationStatistical Methods in functional MRI. Standard Analysis. Data Processing Pipeline. Multiple Comparisons Problem. Multiple Comparisons Problem
Statistical Methods in fnctional MRI Lectre 7: Mltiple Comparisons 04/3/13 Martin Lindqist Department of Biostatistics Johns Hopkins University Data Processing Pipeline Standard Analysis Data Acqisition
More informationMultiple Testing and Thresholding
Multiple Testing and Thresholding NITP, 2009 Thanks for the slides Tom Nichols! Overview Multiple Testing Problem Which of my 100,000 voxels are active? Two methods for controlling false positives Familywise
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationFSL Workshop Session 3 David Smith & John Clithero
FSL Workshop 12.09.08 Session 3 David Smith & John Clithero What is MELODIC? Probabilistic ICA Improves upon standard ICA Allows for inference Avoids over-fitting Three stage process ( ppca ) 1.) Dimension
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 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 informationLecture 3: Camera Calibration, DLT, SVD
Computer Vision Lecture 3 23--28 Lecture 3: Camera Calibration, DL, SVD he Inner Parameters In this section we will introduce the inner parameters of the cameras Recall from the camera equations λx = P
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 information3. Data Structures for Image Analysis L AK S H M O U. E D U
3. Data Structures for Image Analysis L AK S H M AN @ O U. E D U Different formulations Can be advantageous to treat a spatial grid as a: Levelset Matrix Markov chain Topographic map Relational structure
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 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 informationI How does the formulation (5) serve the purpose of the composite parameterization
Supplemental Material to Identifying Alzheimer s Disease-Related Brain Regions from Multi-Modality Neuroimaging Data using Sparse Composite Linear Discrimination Analysis I How does the formulation (5)
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 informationNew Approaches for EEG Source Localization and Dipole Moment Estimation. Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine
New Approaches for EEG Source Localization and Dipole Moment Estimation Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine Outline Motivation why EEG? Mathematical Model equivalent
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 informationDetecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference
Detecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference Minh Dao 1, Xiang Xiang 1, Bulent Ayhan 2, Chiman Kwan 2, Trac D. Tran 1 Johns Hopkins Univeristy, 3400
More informationEE613 Machine Learning for Engineers LINEAR REGRESSION. Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov.
EE613 Machine Learning for Engineers LINEAR REGRESSION Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov. 4, 2015 1 Outline Multivariate ordinary least squares Singular value
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationUnsupervised learning in Vision
Chapter 7 Unsupervised learning in Vision The fields of Computer Vision and Machine Learning complement each other in a very natural way: the aim of the former is to extract useful information from visual
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 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 informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
More informationHistograms. h(r k ) = n k. p(r k )= n k /NM. Histogram: number of times intensity level rk appears in the image
Histograms h(r k ) = n k Histogram: number of times intensity level rk appears in the image p(r k )= n k /NM normalized histogram also a probability of occurence 1 Histogram of Image Intensities Create
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 informationData Mining Chapter 3: Visualizing and Exploring Data Fall 2011 Ming Li Department of Computer Science and Technology Nanjing University
Data Mining Chapter 3: Visualizing and Exploring Data Fall 2011 Ming Li Department of Computer Science and Technology Nanjing University Exploratory data analysis tasks Examine the data, in search of structures
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 informationEE613 Machine Learning for Engineers LINEAR REGRESSION. Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov.
EE613 Machine Learning for Engineers LINEAR REGRESSION Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov. 9, 2017 1 Outline Multivariate ordinary least squares Matlab code:
More informationDimension Reduction CS534
Dimension Reduction CS534 Why dimension reduction? High dimensionality large number of features E.g., documents represented by thousands of words, millions of bigrams Images represented by thousands of
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 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 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 informationNetwork statistics and thresholding
Network statistics and thresholding Andrew Zalesky azalesky@unimelb.edu.au HBM Educational Course June 25, 2017 Network thresholding Unthresholded Moderate thresholding Severe thresholding Strong link
More informationCSC 411 Lecture 18: Matrix Factorizations
CSC 411 Lecture 18: Matrix Factorizations Roger Grosse, Amir-massoud Farahmand, and Juan Carrasquilla University of Toronto UofT CSC 411: 18-Matrix Factorizations 1 / 27 Overview Recall PCA: project data
More informationBayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis
Bayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis Xavier Le Faucheur a, Brani Vidakovic b and Allen Tannenbaum a a School of Electrical and Computer Engineering, b Department of Biomedical
More informationStatistical inference on images
7 Statistical inference on images The goal of statistical inference is to make decisions based on our data, while accounting for uncertainty due to noise in the data. From a broad perspective, statistical
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 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 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 informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationControlling for mul2ple comparisons in imaging analysis. Where we re going. Where we re going 8/15/16
Controlling for mul2ple comparisons in imaging analysis Wednesday, Lecture 2 Jeane?e Mumford University of Wisconsin - Madison Where we re going Review of hypothesis tes2ng introduce mul2ple tes2ng problem
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 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 informationEstimating image bases for visual image reconstruction from human brain activity
Estimating image bases for visual image reconstruction from human brain activity Yusuke Fujiwara 1 Yoichi Miyawaki 2,1 Yukiyasu Kamitani 1 1 ATR Computational Neuroscience Laboratories 2 National Institute
More informationSparse & Functional Principal Components Analysis
Sparse & Functional Principal Components Analysis Genevera I. Allen Department of Statistics and Electrical and Computer Engineering, Rice University, Department of Pediatrics-Neurology, Baylor College
More informationPreface to the Second Edition. Preface to the First Edition. 1 Introduction 1
Preface to the Second Edition Preface to the First Edition vii xi 1 Introduction 1 2 Overview of Supervised Learning 9 2.1 Introduction... 9 2.2 Variable Types and Terminology... 9 2.3 Two Simple Approaches
More informationWeek 7 Picturing Network. Vahe and Bethany
Week 7 Picturing Network Vahe and Bethany Freeman (2005) - Graphic Techniques for Exploring Social Network Data The two main goals of analyzing social network data are identification of cohesive groups
More informationCLASSIFICATION AND CHANGE DETECTION
IMAGE ANALYSIS, CLASSIFICATION AND CHANGE DETECTION IN REMOTE SENSING With Algorithms for ENVI/IDL and Python THIRD EDITION Morton J. Canty CRC Press Taylor & Francis Group Boca Raton London NewYork CRC
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 informationMotion Interpretation and Synthesis by ICA
Motion Interpretation and Synthesis by ICA Renqiang Min Department of Computer Science, University of Toronto, 1 King s College Road, Toronto, ON M5S3G4, Canada Abstract. It is known that high-dimensional
More informationShort on camera geometry and camera calibration
Short on camera geometry and camera calibration Maria Magnusson, maria.magnusson@liu.se Computer Vision Laboratory, Department of Electrical Engineering, Linköping University, Sweden Report No: LiTH-ISY-R-3070
More informationMapping of Hierarchical Activation in the Visual Cortex Suman Chakravartula, Denise Jones, Guillaume Leseur CS229 Final Project Report. Autumn 2008.
Mapping of Hierarchical Activation in the Visual Cortex Suman Chakravartula, Denise Jones, Guillaume Leseur CS229 Final Project Report. Autumn 2008. Introduction There is much that is unknown regarding
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 informationSnPM is an SPM toolbox developed by Andrew Holmes & Tom Nichols
1 of 14 3/30/2005 9:24 PM SnPM A Worked fmri Example SnPM is an SPM toolbox developed by Andrew Holmes & Tom Nichols This page... introduction example data background design setup computation viewing results
More informationControlling for mul-ple comparisons in imaging analysis. Wednesday, Lecture 2 Jeane:e Mumford University of Wisconsin - Madison
Controlling for mul-ple comparisons in imaging analysis Wednesday, Lecture 2 Jeane:e Mumford University of Wisconsin - Madison Where we re going Review of hypothesis tes-ng introduce mul-ple tes-ng problem
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 informationContents. comparison. Multiple comparison problem. Recap & Introduction Inference & multiple. «Take home» message. DISCOS SPM course, CRC, Liège, 2009
DISCOS SPM course, CRC, Liège, 2009 Contents Multiple comparison problem Recap & Introduction Inference & multiple comparison «Take home» message C. Phillips, Centre de Recherches du Cyclotron, ULg, Belgium
More informationLRLW-LSI: An Improved Latent Semantic Indexing (LSI) Text Classifier
LRLW-LSI: An Improved Latent Semantic Indexing (LSI) Text Classifier Wang Ding, Songnian Yu, Shanqing Yu, Wei Wei, and Qianfeng Wang School of Computer Engineering and Science, Shanghai University, 200072
More informationD-Optimal Designs. Chapter 888. Introduction. D-Optimal Design Overview
Chapter 888 Introduction This procedure generates D-optimal designs for multi-factor experiments with both quantitative and qualitative factors. The factors can have a mixed number of levels. For example,
More informationFault detection with principal component pursuit method
Journal of Physics: Conference Series PAPER OPEN ACCESS Fault detection with principal component pursuit method Recent citations - Robust Principal Component Pursuit for Fault Detection in a Blast Furnace
More informationFNC Toolbox Walk Through
FNC Toolbox Walk Through Date: August 10 th, 2009 By: Nathan Swanson, Vince Calhoun at The Mind Research Network Email: nswanson@mrn.org Introduction The FNC Toolbox is an extension of the GIFT toolbox
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 informationChapter 6: Linear Model Selection and Regularization
Chapter 6: Linear Model Selection and Regularization As p (the number of predictors) comes close to or exceeds n (the sample size) standard linear regression is faced with problems. The variance of the
More informationCSC321: Neural Networks. Lecture 13: Learning without a teacher: Autoencoders and Principal Components Analysis. Geoffrey Hinton
CSC321: Neural Networks Lecture 13: Learning without a teacher: Autoencoders and Principal Components Analysis Geoffrey Hinton Three problems with backpropagation Where does the supervision come from?
More informationSegmentation of MR Images of a Beating Heart
Segmentation of MR Images of a Beating Heart Avinash Ravichandran Abstract Heart Arrhythmia is currently treated using invasive procedures. In order use to non invasive procedures accurate imaging modalities
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 information3D Models and Matching
3D Models and Matching representations for 3D object models particular matching techniques alignment-based systems appearance-based systems GC model of a screwdriver 1 3D Models Many different representations
More informationThe Singular Value Decomposition: Let A be any m n matrix. orthogonal matrices U, V and a diagonal matrix Σ such that A = UΣV T.
Section 7.4 Notes (The SVD) The Singular Value Decomposition: Let A be any m n matrix. orthogonal matrices U, V and a diagonal matrix Σ such that Then there are A = UΣV T Specifically: The ordering of
More information