Linear Models in Medical Imaging. John Kornak MI square February 21, 2012
|
|
- Elwin Arnold
- 6 years ago
- Views:
Transcription
1 Linear Models in Medical Imaging John Kornak MI square February 21, 2012
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 Imaging data statistical methods to look for regional effects Tissue differences between groups or over time VBM, TBM (voxel/tensor-based morphometry) PET (positron emmission tomography), fmri (functional MRI) determine activation in the brain due to thought, stimulus or task Diffusion (DWI, DTI, tractography), Bone mineral density etc. etc.
5 FMRI Data: Set of Volumes (over time) or Set of Time-Series (over space) Serial Snapshots of Volunteers brain Time Time Active Passive Baseline
6 Software etc. SPM PET, fmri, VBM and TBM, EEG/MEG ( needs Matlab) FSL fmri primarily + DTI ( R AnalyzeFMRI package + linear models in general ( and then go to your nearest CRAN mirror) Also, check Venables and Ripley Splus book + many R books (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 ) = (variables 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. age, blood pressure etc., (random or controlled) Factors: e.g. diagnosis, gender, drinking level (low, medium, high) etc.
12 The (General) Linear Model A simple linear model might take the form: e.g. y i = β 1 + x i2 β 2 + x i3 β x im β m + ε i i.i.d. = independently and identically distributed
13 The (General) Linear Model For univariate data: y i = x i T β + ε i, i = 1,...,n β = (β 1,...,β m ) T or in matrix notation y = X T β + ε is a set of unknown parameters 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)
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 y i = x age exp(β 1 + x height β 2 ) i = 1,...,n i = 1,...,n
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 i.e. y i = β 1 + β 2 x i + ε i β 1 (0,0,1) β 2 β 3 = 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 Are 2 different conditions equivalent? E.g. is the activation effect: reading a word vs imagining the object equal?
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 of 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 Time Active Passive Baseline
32 image data design matrix kernel 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 Active Passive Baseline β 1 + β 2 + β 3
34 Parameter Estimates Same model for all voxels Different parameters for each voxel beta_0001.img Time-series ˆβ = ˆβ = 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) 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 design matrix kernel 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 significant? α =0.05 5,000 false positive voxels
44 How can we determine a sensible threshold level? 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 Poor Power (risk of false negatives) Poor Specificity (risk of false positives) Good Power
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 #blobs - #holes At high thresholds, just counts blobs Random Field Threshold FWER = Pr(Max voxel u H o ) No holes = Pr(One or more blobs H o ) Never more than 1 blob Pr(χ u 1 H o ) E(χ u H o ) 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
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% 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 Journal of the Royal Statistical Society Series B (1995) 57: 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 NB, no spatial consideration i/v p (i) i/v q/c(v)
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 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 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 u 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 x Model Illustration 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
66 Plug for: February 29 th SFASA Seminar Speaker: Time: Yoav Benjamini, PhD., Professor of Statistics Wednesday, 5pm - 6pm (4:30-5pm pre-seminar social) Location: UCSF China Basin Landing, Room 6702 (specific directions to classrooms: Title: Hierarchical Testing of Families of Hypotheses Abstract: As the size of large testing problems encountered in genomic research keeps increasing, more of these problems have further structure where the set of hypotheses can be partitioned into families of the hypotheses, and the true state of the tested signals tends to be more similar within these subsets than across the subsets. Moreover, interest may lie with a discovery of a family with some signal in it, on top of the discovery of a signal in each of the many hypotheses on its own. The challenges in the analysis of such multiple testing problems will be discussed. We then present the concept of the control on the average over the selected families of a desired error-rate, be it the family-wise error rate, the False Discovery Rate, or their generalizations. We discuss the various considerations involved using the genomic part of a Norwegian epidemiological study of breast cancer, and a study involving genomics and brain imaging. Transportation/Parking: 1) Right at the Caltrain station. 2) If driving, 2-hour parking in the neighborhood as well as reasonably priced parking in lot A for the Giants. The building parking is expensive. 3) If by BART, transfer at Embarcadero to the N-Judah or T-train (both inbound) 66
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationRecap: Gaussian (or Normal) Distribution. Recap: Minimizing the Expected Loss. Topics of This Lecture. Recap: Maximum Likelihood Approach
Truth Course Outline Machine Learning Lecture 3 Fundamentals (2 weeks) Bayes Decision Theory Probability Density Estimation Probability Density Estimation II 2.04.205 Discriminative Approaches (5 weeks)
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 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 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 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 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 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 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 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 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 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 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 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 informationClassification. Vladimir Curic. Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University
Classification Vladimir Curic Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University Outline An overview on classification Basics of classification How to choose appropriate
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationSolution Sketches Midterm Exam COSC 6342 Machine Learning March 20, 2013
Your Name: Your student id: Solution Sketches Midterm Exam COSC 6342 Machine Learning March 20, 2013 Problem 1 [5+?]: Hypothesis Classes Problem 2 [8]: Losses and Risks Problem 3 [11]: Model Generation
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 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 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 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 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 informationLast time... Coryn Bailer-Jones. check and if appropriate remove outliers, errors etc. linear regression
Machine learning, pattern recognition and statistical data modelling Lecture 3. Linear Methods (part 1) Coryn Bailer-Jones Last time... curse of dimensionality local methods quickly become nonlocal as
More informationThe Anatomical Equivalence Class Formulation and its Application to Shape-based Computational Neuroanatomy
The Anatomical Equivalence Class Formulation and its Application to Shape-based Computational Neuroanatomy Sokratis K. Makrogiannis, PhD From post-doctoral research at SBIA lab, Department of Radiology,
More informationINF 4300 Classification III Anne Solberg The agenda today:
INF 4300 Classification III Anne Solberg 28.10.15 The agenda today: More on estimating classifier accuracy Curse of dimensionality and simple feature selection knn-classification K-means clustering 28.10.15
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 informationInterpreting predictive models in terms of anatomically labelled regions
Interpreting predictive models in terms of anatomically labelled regions Data Feature extraction/selection Model specification Accuracy and significance Model interpretation 2 Data Feature extraction/selection
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 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 informationEvaluating Machine Learning Methods: Part 1
Evaluating Machine Learning Methods: Part 1 CS 760@UW-Madison Goals for the lecture you should understand the following concepts bias of an estimator learning curves stratified sampling cross validation
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 informationData Mining. ❷Chapter 2 Basic Statistics. Asso.Prof.Dr. Xiao-dong Zhu. Business School, University of Shanghai for Science & Technology
❷Chapter 2 Basic Statistics Business School, University of Shanghai for Science & Technology 2016-2017 2nd Semester, Spring2017 Contents of chapter 1 1 recording data using computers 2 3 4 5 6 some famous
More informationAn Introduction To Automatic Tissue Classification Of Brain MRI. Colm Elliott Mar 2014
An Introduction To Automatic Tissue Classification Of Brain MRI Colm Elliott Mar 2014 Tissue Classification Tissue classification is part of many processing pipelines. We often want to classify each voxel
More 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 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 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 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 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 informationModern Medical Image Analysis 8DC00 Exam
Parts of answers are inside square brackets [... ]. These parts are optional. Answers can be written in Dutch or in English, as you prefer. You can use drawings and diagrams to support your textual answers.
More informationEvaluating Machine-Learning Methods. Goals for the lecture
Evaluating Machine-Learning Methods Mark Craven and David Page Computer Sciences 760 Spring 2018 www.biostat.wisc.edu/~craven/cs760/ Some of the slides in these lectures have been adapted/borrowed from
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 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 informationPATTERN CLASSIFICATION AND SCENE ANALYSIS
PATTERN CLASSIFICATION AND SCENE ANALYSIS RICHARD O. DUDA PETER E. HART Stanford Research Institute, Menlo Park, California A WILEY-INTERSCIENCE PUBLICATION JOHN WILEY & SONS New York Chichester Brisbane
More information