Linear Models in Medical Imaging. John Kornak MI square February 19, 2013
|
|
- Estella Baldwin
- 5 years ago
- Views:
Transcription
1 Linear Models in Medical Imaging John Kornak MI square February 19, 2013
2 Acknowledgement / Disclaimer Many of the slides in this lecture have been adapted from slides available in talks available on the SPM web site.
3 Overview Motivation Linear model formulation Region of interest analyses Pixel/voxel based analyses Multiple comparisons for images Bayesian image analysis methods
4 Motivation Linear Models = cornerstone of statistical methods - linear regression to complex models Imaging data statistical methods to look for regional effects or at the voxel level Across subjects: Tissue differences between groups (e.g., voxel/tensor-based morphometry) Within subjects: PET (positron emission tomography), fmri (functional MRI) activation in brain due to thought, stimulus or task Non-brain: Diffusion (DWI, DTI, tractography), Bone mineral density etc., etc.
5 FMRI Data Serial Snapshots of Volunteers brain Time Active Passive Baseline
6 Software SPM PET, fmri, VBM and TBM, EEG/MEG ( needs Matlab) FSL fmri primarily + DTI ( R AnalyzeFMRI package + see medical imaging view + linear models in general ( Also, many R books and free ones (see R web site) + online tutorials
7 Challenges Generating suitable statistical imaging models Dealing with highly multivariate responses (curse of dimensionality) Defining imaging hypotheses Creating computationally efficient analysis procedures
8 Aims of Statistical Modeling Summarize data Estimation: point and interval estimates Inference: hypotheses / relationships Prediction
9 Aims of Statistical Modeling Summarize data Estimation: point and interval estimates Inference: hypotheses / relationships Prediction
10 Statistical Modeling Strategy Propose a model for the data Fit the model Assess the model s adequacy Fit other plausible models Compare all fitted models Interpret the best model
11 Statistical Models: Definitions Univariate response variable y i (for exp. unit i) Covariates (x i1, x i2,..., x ik ) = (predictors of interest and nuisance variables) x i T Data is: { y i,x T i ;i = 1,...,n}, n experimental units Continuous covariates: e.g. blood pressure, druglevel, age (random or controlled) Factors: e.g. diagnosis, gender, drinking level (low, medium, high)
12 The (General) Linear Model A simple linear model might take the form: e.g. y i = β i + x i2 β 2 + x i 3 β x im β m + ε i y i = β mean + x i,age β age + x i,gender β gender + x i,diagnosis β diagnosis + ε i ε i ~ N(0,σ 2 ), i.i.d., i =1,.,n i.i.d. = independently and identically distributed
13 The (General) Linear Model For univariate data: y i = x i T β + ε i, β = (β 1,...,β m ) T i = 1,...,n is a set of unknown parameters or in matrix notation y = X T β + ε This can be extended to a multivariate response Y = X T B + E
14 Ex. Hippocampal Volume HCV ~ Age + Diagnosis (Wilkinson notation) Diagnosis can be normal control (NC) or Alzheimer s disease (AD)
15 Ex. Hippocampal Volume HCV ~ Age + Diagnosis + Age*Diagnosis (Wilkinson notation) Diagnosis can be normal control (NC) or Alzheimer s disease (AD)
16 Structural T1 weighted MRI s Hippocampal volumes manually traced Volume measure = response for each subject Disease status encoded 1 for AD and 0 for NC (the term) x diagnosis = x diag.
17 y i = β 1 + x i,age β age + x i,diag. β diag. + x i.age x i,diag. β inter + ε i Case 1 HCV age
18 y i = β 1 + x i,age β age + x i,diag. β diag. + x i.age x i,diag. β inter + ε i Case 2 HCV age
19 Case 3 HCV NC AD age
20 Case 4 HCV NC AD age
21 y i = β 1 + x i,age β age + x i,diag. β diag. + x i.age x i,diag. β inter + ε i Case 4 HCV NC AD age
22 Linear models can be more general - only needs to be linear in the parameters: We can have: y = x β + x 2 β + exp(x )β + x π x β + ε i age 1 age 2 height 3 age height 4 i But not i = 1,...,n y i = x age β 1 + exp( x height β 2 )
23 Estimation Minimize squared error (Least Squares Error) = Maximum Likelihood Estimation for linear model ˆβ = (X T X) 1 X T y E( ˆβ) = β V ( ˆβ) = σ 2 (X T X) 1 Estimate by sum of squares error ˆσ 2 = n or divide by n-1 for unbiased estimate
24 Inference Model Comparison Take linear model y = X T β + ε And add constraint this defines a new model that is a simplification of the previous one
25 Inference Model Comparison E.g., cf. model y i = β 1 + β 2 x i1 + β 3 x i2 + ε i to simplification with β 3 = 0 (0,0,1) β 1 β 2 β 3 i.e. y i = β 1 + β 2 x i + ε i = 0 i.e. Aβ = c
26 What about &? β 2 = 0 β 3 = 0 Aβ = c β 1 β 2 β 3 = 0 0
27 And what about? β 2 = β 3 β 1 ( ) Aβ = c β 2 β 3 = 0 Is there a difference between 2 conditions? E.g. is the activation effect of reading a word vs. imagining the object different?
28 Definition: Linear model nested in another if 1 st model can be obtained by linear constraint on the 2 nd Nesting tree: age + gender age gender null
29 F-test for General Linear Hypothesis ( ) y = X T β + ε ε N n 0,σ 2 I n Consider This is the General Linear Hypothesis
30 Under, i.e., Aβ = c F = (SSE nested SSE larger ) / ( p larger p nested ) (SSE larger ) / (n p larger ) F plarger p nested,n p larger p denotes the number of model parameters n denotes the number of data points SSE = Deviance = sum of squares error /residuals Tests whether ratio of variances = 1
31 FMRI Data: Set of Volumes (over time) or Set of Time-Series (over space) Serial Snapshots of Volunteers brain Time Active Passive Baseline
32 image data kernel design matrix parameter estimates realignment & motion correction smoothing Linear Model model fitting statistic image Thresholding & Random Field Theory normalisation anatomical reference Statistical Parametric Map (test statistics) Corrected thresholds & p-values
33 Estimation The estimation entails finding the parameter values such that the linear combination best fits the data at each voxel Active Passive Baseline β 1 + β 2 + β 3
34 Parameter Estimates Same model for all voxels Different parameters for each voxel 0.83 ˆβ = ˆβ = beta_0001.img beta_0002.img ˆβ = beta_0003.img
35 y X T β SPM View β 1 β 2 β 1 +β 2 +β 3 β 1 +β 2 +β 3 β 3 Note: We trust: Long series with large effects and small error
36 Spatial Modeling
37 Spatial Hypotheses Question - how do we extend from standard univariate hypotheses to answering spatially motivated questions? Not easy - curse of dimensionality (millions of voxels) A B A 2D 1D B e.g. in 1D it makes sense to infer A is less than B, but what is the equivalent in 2D?
38 Spatial Testing Solutions Summarize the image into one dimensional quantities for testing (e.g. region of interest analysis) Consider the overall test as a combination of individual voxel tests (voxel based analysis) Perform shape/object analysis on objects defined via landmarks Build Bayesian image analysis models
39 Spatial Testing Solutions Summarize the image into one dimensional quantities for testing (e.g. region of interest analysis) Consider the overall test as a combination of individual voxel tests (voxel based analysis) Perform shape/object analysis on objects defined via landmarks Build Bayesian image analysis models
40 Voxel based analysis Each voxel obtains a test statistic from the linear model, e.g. t or F Forms statistical maps of the statistics
41 image data kernel design matrix parameter estimates realignment & motion correction smoothing Linear Model model fitting statistic image Thresholding & Random Field Theory normalisation anatomical reference Statistical Parametric Map (test statistics) Corrected thresholds & p-values
42 Hypothesis Testing Null Hypothesis H 0 Test statistic T t observed realization of T α-level Acceptable false positive risk Level α = Pr( T>u α H 0 ) Threshold u α controls false positive risk at level α u α Null Distribution of T α
43 Multiple Comparisons Problem Which of 100,000 voxels are statistically significant? α =0.05 5,000 false positive voxels
44 Assessing Statistic Images Where s the signal or change? High Threshold Med. Threshold Low Threshold t > 5.5 t > 3.5 t > 0.5 Good Specificity Good Sensitivity Poor Sensitivity (risk of false negatives) Poor Specificity (risk of false positives) How can we determine a sensible threshold level?
45 Multiple Comparison Solutions : Measuring False Positives Familywise Error Rate (FWER) Familywise Error Existence of one or more false positives False Discovery Rate (FDR) FDR = E[FP/(TP+FP)] TP+FP voxels declared active, FP falsely so Realized false discovery rate: FP/(TP+FP)
46 Bonferroni Correction FWE, α, for N independent voxels is approximately α = Nv (v = voxel-wise error rate) To control FWE set v = α / N Independent Voxels Spatially Correlated Voxels Bonferroni is too conservative for brain images
47 FWER MCP Solutions: Random Field Theory Euler Characteristic χ u Topological Measure No holes Never more than 1 blob #blobs - #holes At high thresholds, just counts blobs Random Field FWER = Pr(Max voxel u H o ) = Pr(One or more blobs H o ) Pr(χ u 1 H o ) E(χ u H o ) Threshold See description at Suprathreshold Sets
48 Random Field Theory Limitations Multivariate normality (Gaussianity) Virtually impossible to check Sufficient smoothness FWHM smoothness 3-4 voxel size Smoothness estimation Estimate is biased when images not sufficiently smooth Several layers of approximations Lattice Image Data Continuous Random Field
49 Multiple Comparison Solutions: Measuring False Positives Familywise Error Rate (FWER) Familywise Error Existence of one or more false positives False Discovery Rate (FDR) FDR = E[FP/(TP+FP)] TP+FP voxels declared active, FP falsely so Realized false discovery rate: FP/(TP+FP)
50 False Discovery Rate For any threshold, all voxels can be crossclassified: Accept Null Reject Null Null True TN FP Null False FN TP N A N R Realized FDR rfdr = FP /(TP+FP) = FP /N R Special case: if N R = 0, rfdr = 0 But only can observe N R, don t know TP & FP We therefore control the expected rfdr FDR = E(rFDR)
51 False Discovery Rate Illustration: Noise Signal Signal+Noise
52 Control of Per Comparison Rate at 10% 11.3% 11.3% 12.5% 10.8% 11.5% 10.0% 10.7% 11.2% 10.2% 9.5% Percentage of Null Pixels that are False Positives Control of Familywise Error Rate at 10% Occurrence of Familywise Error FWE Control of False Discovery Rate at 10% 6.7% 10.4% 14.9% 9.3% 16.2% 13.8% 14.0% 10.5% 12.2% 8.7% Percentage of Observed Above Threshold Pixels that are False Positives
53 Benjamini & Hochberg Procedure Select desired limit q on FDR Order p-values, p (1) p (2)... p (V) Let r be largest i such that p (i) i/v q/c(v) Reject all hypotheses corresponding to p (1),..., p (r) p-value NB, no spatial consideration 0 1 Journal of the Royal Statistical Society Series B (1995) 57: i/v p (i) i/v q/c(v) 0 1
54 Also, Non-Parametric Testing If H 0 is true then time order irrelevant (if noise really white) Therefore permute the timepoints and obtain test statistics If true test statistic is extreme compared to others then reject H 0
55 Types of Spatial Inference Individual voxel level Cluster level Set level Bayesian model based
56 Voxel-level Inference Retain voxels above α-level threshold u α Gives best spatial specificity H 0 at a single voxel can be rejected u α space Significant Voxels No significant Voxels
57 Cluster-level Inference Two step-process Define clusters by arbitrary threshold u clus Retain clusters larger than α-level threshold k α u clus space Cluster not significant k α k α Cluster significant
58 Cluster-level Inference Typically better sensitivity Worse spatial specificity The null hyp. of entire cluster is rejected Only means that one or more of voxels in cluster active u clus space Cluster not significant k α k α Cluster significant
59 Set-level Inference Count number of blobs c uses minimum blob size k to count significant activity if number of blobs > n(k,u ) Worst spatial specificity Only can reject global null hypothesis clus u clus k k space Here c = 1; only 1 cluster larger than k
60 A flexible Bayesian Approach Model the form of activity Provides an adaptive thresholding approach space Active voxels
61 Bayesian Model y = data, parameter estimates of statistics z = binary activation map modeled as a MRF x = activation level field modeled as a MRF = residual error MRF = Markov Random Field (similar random field but defined on a lattice)
62 Model Illustration x zx + z = 0 z = 1 z = 0
63 Model Illustration x zx fit z = 0 z = 1 z = 0
64 y x z zx + ε
65 Other Topics and Omissions Hemodynamic response function Multiple subjects (random and mixed effects models) PCA, ICA Multivariate analysis with variogram modeling Space-time modeling
Linear 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationBayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri
Bayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri Galin L. Jones 1 School of Statistics University of Minnesota March 2015 1 Joint with Martin Bezener and John Hughes Experiment
More 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 informationStatistical Methods in functional MRI. False Discovery Rate. Issues with FWER. Lecture 7.2: Multiple Comparisons ( ) 04/25/13
Statistical Methods in functional MRI Lecture 7.2: Multiple Comparisons 04/25/13 Martin Lindquist Department of iostatistics Johns Hopkins University Issues with FWER Methods that control the FWER (onferroni,
More informationZurich SPM Course Voxel-Based Morphometry. Ged Ridgway (Oxford & UCL) With thanks to John Ashburner and the FIL Methods Group
Zurich SPM Course 2015 Voxel-Based Morphometry Ged Ridgway (Oxford & UCL) With thanks to John Ashburner and the FIL Methods Group Examples applications of VBM Many scientifically or clinically interesting
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 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 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 informationModelling Neuroimaging Data Using the General Linear Model (GLM Karl) Jesper Andersson KI, Stockholm & BRU, Helsinki
Modelling Neuroimaging Data Using the General Linear Model (GLM Karl) Jesper Andersson KI, Stockholm & BRU, Helsinki No escaping this one fmri time-series Design matrix Parameter Estimates Motion Correction
More informationStatistical Analysis of Neuroimaging Data. Phebe Kemmer BIOS 516 Sept 24, 2015
Statistical Analysis of Neuroimaging Data Phebe Kemmer BIOS 516 Sept 24, 2015 Review from last time Structural Imaging modalities MRI, CAT, DTI (diffusion tensor imaging) Functional Imaging modalities
More informationA Non-Parametric Approach
Andrew P. Holmes. Ph.D., 1994. Chapter Six A Non-Parametric Approach In this chapter, a non-parametric approach to assessing functional mapping experiments is presented. A multiple comparisons randomisation
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 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 information1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12
Contents 1 Introduction 10 1.1 Motivation and Aims....... 10 1.1.1 Functional Imaging.... 10 1.1.2 Computational Neuroanatomy... 12 1.2 Overview of Chapters... 14 2 Rigid Body Registration 18 2.1 Introduction.....
More informationStatistical Analysis of Metabolomics Data. Xiuxia Du Department of Bioinformatics & Genomics University of North Carolina at Charlotte
Statistical Analysis of Metabolomics Data Xiuxia Du Department of Bioinformatics & Genomics University of North Carolina at Charlotte Outline Introduction Data pre-treatment 1. Normalization 2. Centering,
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 informationAdvances in FDR for fmri -p.1
Advances in FDR for fmri Ruth Heller Department of Statistics, University of Pennsylvania Joint work with: Yoav Benjamini, Nava Rubin, Damian Stanley, Daniel Yekutieli, Yulia Golland, Rafael Malach Advances
More informationStatistical Methods in functional MRI. 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 informationPart I. Hierarchical clustering. Hierarchical Clustering. Hierarchical clustering. Produces a set of nested clusters organized as a
Week 9 Based in part on slides from textbook, slides of Susan Holmes Part I December 2, 2012 Hierarchical Clustering 1 / 1 Produces a set of nested clusters organized as a Hierarchical hierarchical clustering
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 informationBayesian Methods in Functional Magnetic Resonance Imaging
Bayesian Methods in Functional Magnetic Resonance Imaging Galin L. Jones Kuo-Jung Lee Brian S. Caffo Susan Spear Bassett Abstract: One of the major objectives of functional magnetic resonance imaging studies
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 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 informationPower analysis. Wednesday, Lecture 3 Jeanette Mumford University of Wisconsin - Madison
Power analysis Wednesday, Lecture 3 Jeanette Mumford University of Wisconsin - Madison Power Analysis-Why? To answer the question. How many subjects do I need for my study? How many runs per subject should
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 informationStatistics 202: Data Mining. c Jonathan Taylor. Outliers Based in part on slides from textbook, slides of Susan Holmes.
Outliers Based in part on slides from textbook, slides of Susan Holmes December 2, 2012 1 / 1 Concepts What is an outlier? The set of data points that are considerably different than the remainder of the
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 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 informationSplines and penalized regression
Splines and penalized regression November 23 Introduction We are discussing ways to estimate the regression function f, where E(y x) = f(x) One approach is of course to assume that f has a certain shape,
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 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 informationVoxel-Based Morphometry & DARTEL. Ged Ridgway, London With thanks to John Ashburner and the FIL Methods Group
Zurich SPM Course 2012 Voxel-Based Morphometry & DARTEL Ged Ridgway, London With thanks to John Ashburner and the FIL Methods Group Aims of computational neuroanatomy * Many interesting and clinically
More informationThe Optimal Discovery Procedure: A New Approach to Simultaneous Significance Testing
UW Biostatistics Working Paper Series 9-6-2005 The Optimal Discovery Procedure: A New Approach to Simultaneous Significance Testing John D. Storey University of Washington, jstorey@u.washington.edu Suggested
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 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 informationMachine Learning. Topic 4: Linear Regression Models
Machine Learning Topic 4: Linear Regression Models (contains ideas and a few images from wikipedia and books by Alpaydin, Duda/Hart/ Stork, and Bishop. Updated Fall 205) Regression Learning Task There
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 informationCPSC 340: Machine Learning and Data Mining. Principal Component Analysis Fall 2017
CPSC 340: Machine Learning and Data Mining Principal Component Analysis Fall 2017 Assignment 3: 2 late days to hand in tonight. Admin Assignment 4: Due Friday of next week. Last Time: MAP Estimation MAP
More information[Programming Assignment] (1)
http://crcv.ucf.edu/people/faculty/bagci/ [Programming Assignment] (1) Computer Vision Dr. Ulas Bagci (Fall) 2015 University of Central Florida (UCF) Coding Standard and General Requirements Code for all
More informationDATA MINING AND MACHINE LEARNING. Lecture 6: Data preprocessing and model selection Lecturer: Simone Scardapane
DATA MINING AND MACHINE LEARNING Lecture 6: Data preprocessing and model selection Lecturer: Simone Scardapane Academic Year 2016/2017 Table of contents Data preprocessing Feature normalization Missing
More informationDecoding analyses. Neuroimaging: Pattern Analysis 2017
Decoding analyses Neuroimaging: Pattern Analysis 2017 Scope Conceptual explanation of decoding/machine learning (ML) and related concepts; Not about mathematical basis of ML algorithms! Learning goals
More informationLecture 7: Linear Regression (continued)
Lecture 7: Linear Regression (continued) Reading: Chapter 3 STATS 2: Data mining and analysis Jonathan Taylor, 10/8 Slide credits: Sergio Bacallado 1 / 14 Potential issues in linear regression 1. Interactions
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 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 informationNorbert Schuff VA Medical Center and UCSF
Norbert Schuff Medical Center and UCSF Norbert.schuff@ucsf.edu Medical Imaging Informatics N.Schuff Course # 170.03 Slide 1/67 Objective Learn the principle segmentation techniques Understand the role
More informationBig Data Methods. Chapter 5: Machine learning. Big Data Methods, Chapter 5, Slide 1
Big Data Methods Chapter 5: Machine learning Big Data Methods, Chapter 5, Slide 1 5.1 Introduction to machine learning What is machine learning? Concerned with the study and development of algorithms that
More information3 Feature Selection & Feature Extraction
3 Feature Selection & Feature Extraction Overview: 3.1 Introduction 3.2 Feature Extraction 3.3 Feature Selection 3.3.1 Max-Dependency, Max-Relevance, Min-Redundancy 3.3.2 Relevance Filter 3.3.3 Redundancy
More informationIntroduction to machine learning, pattern recognition and statistical data modelling Coryn Bailer-Jones
Introduction to machine learning, pattern recognition and statistical data modelling Coryn Bailer-Jones What is machine learning? Data interpretation describing relationship between predictors and responses
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 informationWeek 5: Multiple Linear Regression II
Week 5: Multiple Linear Regression II Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Adjusted R
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 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 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 informationLecture 16: High-dimensional regression, non-linear regression
Lecture 16: High-dimensional regression, non-linear regression Reading: Sections 6.4, 7.1 STATS 202: Data mining and analysis November 3, 2017 1 / 17 High-dimensional regression Most of the methods we
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 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 informationDimension Reduction for Big Data Analysis. Dan Shen. Department of Mathematics & Statistics University of South Florida.
Dimension Reduction for Big Data Analysis Dan Shen Department of Mathematics & Statistics University of South Florida danshen@usf.edu October 24, 2014 1 Outline Multiscale weighted PCA for Image Analysis
More informationTopics in Machine Learning
Topics in Machine Learning Gilad Lerman School of Mathematics University of Minnesota Text/slides stolen from G. James, D. Witten, T. Hastie, R. Tibshirani and A. Ng Machine Learning - Motivation Arthur
More informationGeneralized additive models I
I Patrick Breheny October 6 Patrick Breheny BST 764: Applied Statistical Modeling 1/18 Introduction Thus far, we have discussed nonparametric regression involving a single covariate In practice, we often
More informationCPSC 340: Machine Learning and Data Mining. Principal Component Analysis Fall 2016
CPSC 340: Machine Learning and Data Mining Principal Component Analysis Fall 2016 A2/Midterm: Admin Grades/solutions will be posted after class. Assignment 4: Posted, due November 14. Extra office hours:
More informationDetecting Changes In Non-Isotropic Images
Detecting Changes In Non-Isotropic Images K.J. Worsley 1, M. Andermann 1, T. Koulis 1, D. MacDonald, 2 and A.C. Evans 2 August 4, 1999 1 Department of Mathematics and Statistics, 2 Montreal Neurological
More 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 information08 An Introduction to Dense Continuous Robotic Mapping
NAVARCH/EECS 568, ROB 530 - Winter 2018 08 An Introduction to Dense Continuous Robotic Mapping Maani Ghaffari March 14, 2018 Previously: Occupancy Grid Maps Pose SLAM graph and its associated dense occupancy
More informationDATA ANALYSIS USING HIERARCHICAL GENERALIZED LINEAR MODELS WITH R
DATA ANALYSIS USING HIERARCHICAL GENERALIZED LINEAR MODELS WITH R Lee, Rönnegård & Noh LRN@du.se Lee, Rönnegård & Noh HGLM book 1 / 24 Overview 1 Background to the book 2 Crack growth example 3 Contents
More informationCLASSIFICATION WITH RADIAL BASIS AND PROBABILISTIC NEURAL NETWORKS
CLASSIFICATION WITH RADIAL BASIS AND PROBABILISTIC NEURAL NETWORKS CHAPTER 4 CLASSIFICATION WITH RADIAL BASIS AND PROBABILISTIC NEURAL NETWORKS 4.1 Introduction Optical character recognition is one of
More informationWildfire Chances and Probabilistic Risk Assessment
Wildfire Chances and Probabilistic Risk Assessment D R Brillinger, H K Preisler and H M Naderi Statistics Department University of California Berkeley, CA, 97-386 Pacific Southwest Research Station USDA
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 10 130221 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Canny Edge Detector Hough Transform Feature-Based
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 informationOnline Statistical Inference for Quantifying Mandible Growth in CT Images
Online Statistical Inference for Quantifying Mandible Growth in CT Images Moo K. Chung 1,2, Ying Ji Chuang 2, Houri K. Vorperian 2 1 Department of Biostatistics and Medical Informatics 2 Vocal Tract Development
More informationfmri Basics: Single Subject Analysis
fmri Basics: Single Subject Analysis This session is intended to give an overview of the basic process of setting up a general linear model for a single subject. This stage of the analysis is also variously
More informationSelected topics on fmri image processing
fmri Symposium, 24 January, 2006, Ghent Selected topics on fmri image processing Jan Sijbers The fmri folks of the Vision Lab UZA Bio-Imaging Lab Vision Lab Univ. Maastricht T.U. Delft 1 Overview generalized
More informationMachine Learning / Jan 27, 2010
Revisiting Logistic Regression & Naïve Bayes Aarti Singh Machine Learning 10-701/15-781 Jan 27, 2010 Generative and Discriminative Classifiers Training classifiers involves learning a mapping f: X -> Y,
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 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 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 informationStatistical Modeling of Neuroimaging Data: Targeting Activation, Task-Related Connectivity, and Prediction
Statistical Modeling of Neuroimaging Data: Targeting Activation, Task-Related Connectivity, and Prediction F. DuBois Bowman Department of Biostatistics and Bioinformatics Emory University, Atlanta, GA,
More informationFast or furious? - User analysis of SF Express Inc
CS 229 PROJECT, DEC. 2017 1 Fast or furious? - User analysis of SF Express Inc Gege Wen@gegewen, Yiyuan Zhang@yiyuan12, Kezhen Zhao@zkz I. MOTIVATION The motivation of this project is to predict the likelihood
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 informationCSE 158. Web Mining and Recommender Systems. Midterm recap
CSE 158 Web Mining and Recommender Systems Midterm recap Midterm on Wednesday! 5:10 pm 6:10 pm Closed book but I ll provide a similar level of basic info as in the last page of previous midterms CSE 158
More informationTutorial MICCAI'04. fmri data analysis: state of the art and future challenges. Jean-Baptiste Poline Philippe Ciuciu Alexis Roche
Tutorial MICCAI'04 fmri data analysis: state of the art and future challenges Jean-Baptiste Poline Philippe Ciuciu Alexis Roche CEA/SHFJ, Orsay (France) Goal of this tutorial and plan I. Talk 1: will orientate
More 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 information