Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco.

Size: px
Start display at page:

Download "Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco."

Transcription

1 Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu Jenni Pacheco

2 Outline Processing Stages Command line Stream Assemble Data (mris_preproc, mri_surf2surf) Design/Contrast (GLM Theory) (design matrix) Analyze (mri_glmfit) Visualize (tksurfer, QDEC) Interac)ve/Automated GUI (QDEC) Correc)on for mul)ple comparisons (mri_surfcluster) 2

3 Surface-based Study (Thickness) 3

4 Example: Thickness Study 1. $SUBJECTS_DIR/bert/surf/lh.thickness 2. $SUBJECTS_DIR/fred/surf/lh.thickness 3. $SUBJECTS_DIR/sara/surf/lh.thickness 4. $SUBJECTS_DIR/margaret/surf/lh.thickness 5. lh.thickness is in the subject s surface space, Ie, not a common surface space.

5 Processing Stages Specify Subjects and Surface measures Assemble Data: Resample into common space Smooth Concatenate into one file Model and Contrasts (GLM) Fit Model (Es)mate) Correct for mul)ple comparisons Visualize 5

6 Subject 1 Native Inter Subject Averaging Spherical GLM Spherical Subject 2 Demographics Surface-to- Surface Surface-to- Surface

7 GLM Theory y is surface-based measure X is design matrix Beta are regression coefficients C is contrast matrix GLM: fit for beta Null Hypotheses: gamma = 0 Performed at each surface vertex Need: y X C Conceptually no different than volumebased Output: beta gamma (cope) t/f/p

8 Surface based Sta)s)cal Analysis Making an average subject from your set of subjects (or use fsaverage) Construc)ng a FreeSurfer Group Descriptor File (FSGD) Preprocessing the group data Construc)ng the design matrix Construc)ng contrast matrices to test hypotheses Correc)ng for mul)ple comparisons

9 Average(Target) Subject Use the pre made average subject, fsaverage Average of 40 subjects Create your own average subjects from your data set Use one individual subject from your group Something else

10 Surface based Sta)s)cal Analysis Making an average subject from your set of subjects (or use fsaverage) Construc)ng a FreeSurfer Group Descriptor File (FSGD) Preprocessing the group data Construc)ng the design matrix Construc)ng contrast matrices to test hypotheses Correc)ng for mul)ple comparisons

11 FSGD File FSGDF = FreeSurfer Group Descriptor File, Supplies (most of) X and y. GroupDescriptorFile 1 Class Male red square Class Female blue circle Variables Age Weight IQ Input bert Male Input fred Male Input sara Female Input margaret Female One Discrete Factor (Gender) with Two Levels (M&F) Three Continuous Variables: Age, Weight, IQ

12 Surface based Sta)s)cal Analysis Making an average subject from your set of subjects (or use fsaverage) Construc)ng a FreeSurfer Group Descriptor File (FSGD) Preprocessing the group data Construc)ng the design matrix Construc)ng contrast matrices to test hypotheses Correc)ng for mul)ple comparisons

13 mris_preproc Assembles your subjects into a common space (spherical) and gathers the information from the measurement you are using (thickness) into one file. Will need to specify: fsgd file Hemisphere Measure (thickness, curv, sulc, functional values, etc ) Target subject Input: y X C

14 mris_preproc Command mris_preproc fsgd gender_age.txt target average hemi lh meas thickness out lh.gender_age.thickness.mgh

15 Surface Smoothing Uses output from mris_preproc lh.thickness.mgh 2D surface based smoothing (eg, fwhm = 10 mm) Saves lh.thickness.sm10.mgh Why should you smooth? Improve CNR Improve inter-subject registration How much smoothing? Blob-size Typically mm FWHM More forgiving than volume-based Input: y X C

16 mri_surf2surf Command mri_surf2surf hemi lh s average sval lh.gender_age.thickness.mgh fwhm 10 tval lh.gender_age.thickness.10.mgh

17 Surface based Sta)s)cal Analysis Making an average subject from your set of subjects (or use fsaverage) Construc)ng a FreeSurfer Group Descriptor File (FSGD) Preprocessing the group data Construc)ng the design matrix Construc)ng contrast matrices to test hypotheses Correc)ng for mul)ple comparisons

18 DOSS vs DODS Different Offset, Same Slope Different Offset, Different Slope

19 Design matrix and contrasts DOSS Different Offset, Same Slope Female Class Male Class X = Age for Males and Females C = [ ] } Tests for the difference in offset between groups Input: y X C #Regressors = Nv+Nc = 3+2=5 Fewer regressors than DODS DOF = #Rows - #Regressors

20 Design matrix and contrasts DODS Different Offset, Different Slope Female Group Male Group X = Male Age Female Age C = [ ] } Same test, different vector #Regressors = (Nv+1)*Nc = (3+1)*2=8 DOF = #Rows - #Regressors Input: y X C

21 mri_glmfit Reads in FSGD File and constructs X Reads in your contrasts (C1, C2, etc) Uses output from mris_surf2surf lh.thickness.sm10.mgh Fits GLM Computes contrasts (gamma) t or F ra)os and significances

22 mri_glmfit Command mri_glmfit y lh.thickness.sm10.mgh fsgd gender_age.txt doss surf fsaverage lh glmdir lh.gender_age.glmdir C age.mat C gender.mat

23 mri_glmfit Output Creates: lh.gender_age.glmdir/ beta.mgh parameter estimates eres.mgh residual error rvar.mgh residual error variance y.fsgd fsgd file used for scatter plot etc age/ sig.mgh -log10(p) gamma.mgh, F.mgh gender/ sig.mgh -log10(p) gamma.mgh, F.mgh

24 Visualization with tksurfer Threshold: -log10(p), Eg, 2=.01 Saturation: -log10(p), Eg, 5= False Discovery Rate Eg,.01 View->Configure->Overlay File->LoadOverlay

25 Visualization with tksurfer File-> Load Group Descriptor File

26 Problem of Mul)ple Comparisons p < 0.10 p < 0.01 p <

27 Correc)on for Mul)ple Comparisons Cluster based Monte Carlo simula)on Permuta)on Tests Surface Gaussian Random Fields (GRF) There but not fully tested False Discovery Rate (FDR) built into tksurfer and QDEC. (Genovese, et al, NI 2002) 27

28 Clustering 1. Choose a vertex wise threshold Eg, 2 (p<.01), or 3 (p<.001) Sign (pos, neg, abs) 2. A cluster is a group of connected (neighboring) ver)ces above threshold 3. Cluster has a size (area in mm 2 ) p<.01 (-log10(p)=2) Negative p<.0001 (-log10(p)=4) Negative 28

29 Cluster based Correc)on for Mul)ple Comparisons 1. Simulate data under Null Hypothesis: Synthesize Gaussian noise and then smooth (Monte Carlo) Permute rows of design matrix (Permuta)on, orthog) 2. Analyze, threshold, cluster, max cluster size 3. Repeat 10,000 )mes 4. Analyze real data, get cluster sizes 5. P(cluster) = #MaxClusterSize > ClusterSize/ mri_glmfit-sim 29

Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco.

Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco. Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu Jenni Pacheco jpacheco@mail.utexas.edu fmri Integra)on Visualize individual fmri results on surface volume ROI Volume Study: Count number of voxels

More information

Surface-based Analysis: Inter-subject Registration and Smoothing

Surface-based Analysis: Inter-subject Registration and Smoothing Surface-based Analysis: Inter-subject Registration and Smoothing Outline Exploratory Spatial Analysis Coordinate Systems 3D (Volumetric) 2D (Surface-based) Inter-subject registration Volume-based Surface-based

More information

Function-Structure Integration in FreeSurfer

Function-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 information

Introduc)on to fmri. Natalia Zaretskaya

Introduc)on to fmri. Natalia Zaretskaya Introduc)on to fmri Natalia Zaretskaya Content fmri signal fmri versus neural ac)vity A classical experiment: flickering checkerboard Preprocessing Univariate analysis Single- subject analysis Group analysis

More information

Extending the GLM. Outline. Mixed effects motivation Evaluating mixed effects methods Three methods. Conclusions. Overview

Extending 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 information

New and best-practice approaches to thresholding

New 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 information

Multiple Testing and Thresholding

Multiple 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 information

Multiple Testing and Thresholding

Multiple 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 information

Introductory Concepts for Voxel-Based Statistical Analysis

Introductory 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 information

Multiple comparisons problem and solutions

Multiple 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 information

Controlling for mul2ple comparisons in imaging analysis. Where we re going. Where we re going 8/15/16

Controlling 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 information

Journal of Articles in Support of The Null Hypothesis

Journal 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 information

Introduction to Statistical Analyses in SAS

Introduction to Statistical Analyses in SAS Introduction to Statistical Analyses in SAS Programming Workshop Presented by the Applied Statistics Lab Sarah Janse April 5, 2017 1 Introduction Today we will go over some basic statistical analyses in

More information

Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco.

Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco. Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu Jenni Pacheco jpacheco@mail.utexas.edu Outline FreeSurfer ROI Terminology ROI Sta)s)cs Files ROI Studies Volumetric/Area Intensity FreeSurfer

More information

Statistical Methods in functional MRI. False Discovery Rate. Issues with FWER. Lecture 7.2: Multiple Comparisons ( ) 04/25/13

Statistical 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 information

seam Documentation Release 0.0 Scott Burns

seam Documentation Release 0.0 Scott Burns seam Documentation Release 0.0 Scott Burns April 08, 2014 Contents 1 Contents 3 1.1 Installation................................................ 3 1.2 Philosophy................................................

More information

Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco.

Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu. Jenni Pacheco. Introduc)on to FreeSurfer h0p://surfer.nmr.mgh.harvard.edu Jenni Pacheco jpacheco@mail.utexas.edu Analyzing the individual subject What happens? How do I do that? Now what? Tools you will use recon all

More information

Group Sta*s*cs in MEG/EEG

Group Sta*s*cs in MEG/EEG Group Sta*s*cs in MEG/EEG Will Woods NIF Fellow Brain and Psychological Sciences Research Centre Swinburne University of Technology A Cau*onary tale. A Cau*onary tale. A Cau*onary tale. Overview Introduc*on

More information

Medical Image Analysis

Medical 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 information

Brief Guide on Using SPSS 10.0

Brief Guide on Using SPSS 10.0 Brief Guide on Using SPSS 10.0 (Use student data, 22 cases, studentp.dat in Dr. Chang s Data Directory Page) (Page address: http://www.cis.ysu.edu/~chang/stat/) I. Processing File and Data To open a new

More information

Linear Models in Medical Imaging. John Kornak MI square February 19, 2013

Linear 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 information

Linear Models in Medical Imaging. John Kornak MI square February 22, 2011

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 information

CHAPTER 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 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 information

Multiple Testing and Thresholding

Multiple 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 information

fmri Analysis Sackler Ins2tute 2011

fmri Analysis Sackler Ins2tute 2011 fmri Analysis Sackler Ins2tute 2011 How do we get from this to this? How do we get from this to this? And what are those colored blobs we re all trying to see, anyway? Raw fmri data straight from the scanner

More information

Sta$s$cs & Experimental Design with R. Barbara Kitchenham Keele University

Sta$s$cs & Experimental Design with R. Barbara Kitchenham Keele University Sta$s$cs & Experimental Design with R Barbara Kitchenham Keele University 1 Analysis of Variance Mul$ple groups with Normally distributed data 2 Experimental Design LIST Factors you may be able to control

More information

fmri Basics: Single Subject Analysis

fmri 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 information

5. Key ingredients for programming a random walk

5. Key ingredients for programming a random walk 1. Random Walk = Brownian Motion = Markov Chain (with caveat that each of these is a class of different models, not all of which are exactly equal) 2. A random walk is a specific kind of time series process

More information

Correction for multiple comparisons. Cyril Pernet, PhD SBIRC/SINAPSE University of Edinburgh

Correction 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 information

Cluster failure: Why fmri inferences for spatial extent have inflated false positive rates

Cluster failure: Why fmri inferences for spatial extent have inflated false positive rates Supporting Information Appendix Cluster failure: Why fmri inferences for spatial extent have inflated false positive rates Anders Eklund, Thomas Nichols, Hans Knutsson Methods Resting state fmri data Resting

More information

Single Subject Demo Data Instructions 1) click "New" and answer "No" to the "spatially preprocess" question.

Single 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 information

FMRI Pre-Processing and Model- Based Statistics

FMRI 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

Decision Trees: Representa:on

Decision Trees: Representa:on Decision Trees: Representa:on Machine Learning Fall 2017 Some slides from Tom Mitchell, Dan Roth and others 1 Key issues in machine learning Modeling How to formulate your problem as a machine learning

More information

= secant lines of the n+1 unit circle divisions (d)

= secant lines of the n+1 unit circle divisions (d) The logarithm of is a complex number such that. In my theory, I represent as the complex number, which is the same as Euler s, where Taking it a step further with, you get a complex number and its conjugate

More information

Basic Introduction to Data Analysis. Block Design Demonstration. Robert Savoy

Basic 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 information

Controlling 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 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 information

HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006

HST.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 information

Lecture 6: Edge Detection

Lecture 6: Edge Detection #1 Lecture 6: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Options for Image Representation Introduced the concept of different representation or transformation Fourier Transform

More information

First-level fmri modeling

First-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 information

Package gee. June 29, 2015

Package gee. June 29, 2015 Title Generalized Estimation Equation Solver Version 4.13-19 Depends stats Suggests MASS Date 2015-06-29 DateNote Gee version 1998-01-27 Package gee June 29, 2015 Author Vincent J Carey. Ported to R by

More information

Weka ( )

Weka (  ) Weka ( http://www.cs.waikato.ac.nz/ml/weka/ ) The phases in which classifier s design can be divided are reflected in WEKA s Explorer structure: Data pre-processing (filtering) and representation Supervised

More information

Applied Regression Modeling: A Business Approach

Applied Regression Modeling: A Business Approach i Applied Regression Modeling: A Business Approach Computer software help: SAS SAS (originally Statistical Analysis Software ) is a commercial statistical software package based on a powerful programming

More information

Linear Models in Medical Imaging. John Kornak MI square February 23, 2010

Linear 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 information

FreeSurfer Register Download Platforms: Linux, Mac, Windows.

FreeSurfer  Register Download Platforms: Linux, Mac, Windows. FreeSurfer http://surfer.nmr.mgh.harvard.edu Register Download Platforms: Linux, Mac, Windows freesurfer@nmr.mgh.harvard.edu Big Picture Input: T1-weighted (MPRAGE,SPGR) 1mm 3 resolution (.dcm) Output:

More information

FSL Pre-Processing Pipeline

FSL Pre-Processing Pipeline The Art and Pitfalls of fmri Preprocessing FSL Pre-Processing Pipeline Mark Jenkinson FMRIB Centre, University of Oxford FSL Pre-Processing Pipeline Standard pre-processing: Task fmri Resting-state fmri

More information

7/15/2016 ARE YOUR ANALYSES TOO WHY IS YOUR ANALYSIS PARAMETRIC? PARAMETRIC? That s not Normal!

7/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 information

SYS 6021 Linear Statistical Models

SYS 6021 Linear Statistical Models SYS 6021 Linear Statistical Models Project 2 Spam Filters Jinghe Zhang Summary The spambase data and time indexed counts of spams and hams are studied to develop accurate spam filters. Static models are

More information

GLM II. Basic Modeling Strategy CAS Ratemaking and Product Management Seminar by Paul Bailey. March 10, 2015

GLM II. Basic Modeling Strategy CAS Ratemaking and Product Management Seminar by Paul Bailey. March 10, 2015 GLM II Basic Modeling Strategy 2015 CAS Ratemaking and Product Management Seminar by Paul Bailey March 10, 2015 Building predictive models is a multi-step process Set project goals and review background

More information

Bayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis

Bayesian 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 information

SYDE Winter 2011 Introduction to Pattern Recognition. Clustering

SYDE Winter 2011 Introduction to Pattern Recognition. Clustering SYDE 372 - Winter 2011 Introduction to Pattern Recognition Clustering Alexander Wong Department of Systems Design Engineering University of Waterloo Outline 1 2 3 4 5 All the approaches we have learned

More information

UNC 4D Infant Cortical Surface Atlases, from Neonate to 6 Years of Age

UNC 4D Infant Cortical Surface Atlases, from Neonate to 6 Years of Age UNC 4D Infant Cortical Surface Atlases, from Neonate to 6 Years of Age Version 1.0 UNC 4D infant cortical surface atlases from neonate to 6 years of age contain 11 time points, including 1, 3, 6, 9, 12,

More information

Design of Experiments

Design of Experiments Seite 1 von 1 Design of Experiments Module Overview In this module, you learn how to create design matrices, screen factors, and perform regression analysis and Monte Carlo simulation using Mathcad. Objectives

More information

Minitab 17 commands Prepared by Jeffrey S. Simonoff

Minitab 17 commands Prepared by Jeffrey S. Simonoff Minitab 17 commands Prepared by Jeffrey S. Simonoff Data entry and manipulation To enter data by hand, click on the Worksheet window, and enter the values in as you would in any spreadsheet. To then save

More information

COSC160: Detection and Classification. Jeremy Bolton, PhD Assistant Teaching Professor

COSC160: Detection and Classification. Jeremy Bolton, PhD Assistant Teaching Professor COSC160: Detection and Classification Jeremy Bolton, PhD Assistant Teaching Professor Outline I. Problem I. Strategies II. Features for training III. Using spatial information? IV. Reducing dimensionality

More information

Introduction. About this Document. What is SPSS. ohow to get SPSS. oopening Data

Introduction. About this Document. What is SPSS. ohow to get SPSS. oopening Data Introduction About this Document This manual was written by members of the Statistical Consulting Program as an introduction to SPSS 12.0. It is designed to assist new users in familiarizing themselves

More information

Robust Linear Regression (Passing- Bablok Median-Slope)

Robust Linear Regression (Passing- Bablok Median-Slope) Chapter 314 Robust Linear Regression (Passing- Bablok Median-Slope) Introduction This procedure performs robust linear regression estimation using the Passing-Bablok (1988) median-slope algorithm. Their

More information

Functional MRI in Clinical Research and Practice Preprocessing

Functional 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 information

Scott Smith Advanced Image Processing March 15, Speeded-Up Robust Features SURF

Scott Smith Advanced Image Processing March 15, Speeded-Up Robust Features SURF Scott Smith Advanced Image Processing March 15, 2011 Speeded-Up Robust Features SURF Overview Why SURF? How SURF works Feature detection Scale Space Rotational invariance Feature vectors SURF vs Sift Assumptions

More information

Statistical 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 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 information

Excel 2010 with XLSTAT

Excel 2010 with XLSTAT Excel 2010 with XLSTAT J E N N I F E R LE W I S PR I E S T L E Y, PH.D. Introduction to Excel 2010 with XLSTAT The layout for Excel 2010 is slightly different from the layout for Excel 2007. However, with

More information

Data Mining: Concepts and Techniques. (3 rd ed.) Chapter 3. Chapter 3: Data Preprocessing. Major Tasks in Data Preprocessing

Data Mining: Concepts and Techniques. (3 rd ed.) Chapter 3. Chapter 3: Data Preprocessing. Major Tasks in Data Preprocessing Data Mining: Concepts and Techniques (3 rd ed.) Chapter 3 1 Chapter 3: Data Preprocessing Data Preprocessing: An Overview Data Quality Major Tasks in Data Preprocessing Data Cleaning Data Integration Data

More information

Decision Trees: Representation

Decision Trees: Representation Decision Trees: Representation Machine Learning Fall 2018 Some slides from Tom Mitchell, Dan Roth and others 1 Key issues in machine learning Modeling How to formulate your problem as a machine learning

More information

EPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing

EPI 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 information

Linear Methods for Regression and Shrinkage Methods

Linear 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 information

SPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing

SPM8 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 information

Package PSTR. September 25, 2017

Package PSTR. September 25, 2017 Type Package Package PSTR September 25, 2017 Title Panel Smooth Transition Regression Modelling Version 1.1.0 Provides the Panel Smooth Transition Regression (PSTR) modelling. The modelling procedure consists

More information

MULTI-RESOLUTION STATISTICAL ANALYSIS ON GRAPH STRUCTURED DATA IN NEUROIMAGING

MULTI-RESOLUTION STATISTICAL ANALYSIS ON GRAPH STRUCTURED DATA IN NEUROIMAGING MULTI-RESOLUTION STATISTICAL ANALYSIS ON GRAPH STRUCTURED DATA IN NEUROIMAGING, Vikas Singh, Moo Chung, Nagesh Adluru, Barbara B. Bendlin, Sterling C. Johnson University of Wisconsin Madison Apr. 19, 2015

More information

WEEK 4 REVIEW. Graphing Systems of Linear Inequalities (3.1)

WEEK 4 REVIEW. Graphing Systems of Linear Inequalities (3.1) WEEK 4 REVIEW Graphing Systems of Linear Inequalities (3.1) Linear Programming Problems (3.2) Checklist for Exam 1 Review Sample Exam 1 Graphing Linear Inequalities Graph the following system of inequalities.

More information

BrainVoyager TM. Getting Started Guide. Version 3.0. for BV 21

BrainVoyager TM. Getting Started Guide. Version 3.0. for BV 21 BrainVoyager TM Getting Started Guide Version 3.0 for BV 21 Rainer Goebel, Henk Jansma, Caroline Benjamins, Judith Eck, Hester Breman and Armin Heinecke Copyright 2018 Brain Innovation B.V. Contents About

More information

Application of Hierarchical Clustering to Find Expression Modules in Cancer

Application of Hierarchical Clustering to Find Expression Modules in Cancer Application of Hierarchical Clustering to Find Expression Modules in Cancer T. M. Murali August 18, 2008 Innovative Application of Hierarchical Clustering A module map showing conditional activity of expression

More information

Time Series Analysis by State Space Methods

Time Series Analysis by State Space Methods Time Series Analysis by State Space Methods Second Edition J. Durbin London School of Economics and Political Science and University College London S. J. Koopman Vrije Universiteit Amsterdam OXFORD UNIVERSITY

More information

Basic fmri Design and Analysis. Preprocessing

Basic 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 information

Training-Free, Generic Object Detection Using Locally Adaptive Regression Kernels

Training-Free, Generic Object Detection Using Locally Adaptive Regression Kernels Training-Free, Generic Object Detection Using Locally Adaptive Regression Kernels IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIENCE, VOL.32, NO.9, SEPTEMBER 2010 Hae Jong Seo, Student Member,

More information

Digital Image Processing COSC 6380/4393

Digital Image Processing COSC 6380/4393 Digital Image Processing COSC 6380/4393 Lecture 21 Nov 16 th, 2017 Pranav Mantini Ack: Shah. M Image Processing Geometric Transformation Point Operations Filtering (spatial, Frequency) Input Restoration/

More information

COMPUTER AND ROBOT VISION

COMPUTER AND ROBOT VISION VOLUME COMPUTER AND ROBOT VISION Robert M. Haralick University of Washington Linda G. Shapiro University of Washington T V ADDISON-WESLEY PUBLISHING COMPANY Reading, Massachusetts Menlo Park, California

More information

Introduction to Neuroimaging Janaina Mourao-Miranda

Introduction 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 information

Smoothing Dissimilarities for Cluster Analysis: Binary Data and Functional Data

Smoothing Dissimilarities for Cluster Analysis: Binary Data and Functional Data Smoothing Dissimilarities for Cluster Analysis: Binary Data and unctional Data David B. University of South Carolina Department of Statistics Joint work with Zhimin Chen University of South Carolina Current

More information

Cross- Valida+on & ROC curve. Anna Helena Reali Costa PCS 5024

Cross- Valida+on & ROC curve. Anna Helena Reali Costa PCS 5024 Cross- Valida+on & ROC curve Anna Helena Reali Costa PCS 5024 Resampling Methods Involve repeatedly drawing samples from a training set and refibng a model on each sample. Used in model assessment (evalua+ng

More information

FreeSurfer Tutorial :02

FreeSurfer Tutorial :02 FreeSurfer Tutorial FreeSurfer is brought to you by the Martinos Center for Biomedical Imaging (supported by NCRR/P41RR14075), Massachusetts General Hospital, Boston, MA USA 2009-06-01 12:02 1 FSL/FreeSurfer

More information

Multiple Comparisons of Treatments vs. a Control (Simulation)

Multiple Comparisons of Treatments vs. a Control (Simulation) Chapter 585 Multiple Comparisons of Treatments vs. a Control (Simulation) Introduction This procedure uses simulation to analyze the power and significance level of two multiple-comparison procedures that

More information

Dynamic Geometry Processing

Dynamic Geometry Processing Dynamic Geometry Processing EG 2012 Tutorial Will Chang, Hao Li, Niloy Mitra, Mark Pauly, Michael Wand Tutorial: Dynamic Geometry Processing 1 Articulated Global Registration Introduction and Overview

More information

Advances in FDR for fmri -p.1

Advances 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 information

2. Linear Regression and Gradient Descent

2. Linear Regression and Gradient Descent Pattern Recognition And Machine Learning - EPFL - Fall 2015 Emtiyaz Khan, Timur Bagautdinov, Carlos Becker, Ilija Bogunovic & Ksenia Konyushkova 2. Linear Regression and Gradient Descent 2.1 Goals The

More information

Chapter 6: DESCRIPTIVE STATISTICS

Chapter 6: DESCRIPTIVE STATISTICS Chapter 6: DESCRIPTIVE STATISTICS Random Sampling Numerical Summaries Stem-n-Leaf plots Histograms, and Box plots Time Sequence Plots Normal Probability Plots Sections 6-1 to 6-5, and 6-7 Random Sampling

More information

Overview and Practical Application of Machine Learning in Pricing

Overview and Practical Application of Machine Learning in Pricing Overview and Practical Application of Machine Learning in Pricing 2017 CAS Spring Meeting May 23, 2017 Duncan Anderson and Claudine Modlin (Willis Towers Watson) Mark Richards (Allstate Insurance Company)

More information

CSE 158 Lecture 2. Web Mining and Recommender Systems. Supervised learning Regression

CSE 158 Lecture 2. Web Mining and Recommender Systems. Supervised learning Regression CSE 158 Lecture 2 Web Mining and Recommender Systems Supervised learning Regression Supervised versus unsupervised learning Learning approaches attempt to model data in order to solve a problem Unsupervised

More information

EE795: Computer Vision and Intelligent Systems

EE795: 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 information

Getting Started Guide

Getting Started Guide Getting Started Guide Version 2.5 for BVQX 1.9 Rainer Goebel, Henk Jansma and Jochen Seitz Copyright 2007 Brain Innovation B.V. 2 Contents Preface...4 The Objects Tutorial...5 Scanning session information...5

More information

Edge detection. Gradient-based edge operators

Edge detection. Gradient-based edge operators Edge detection Gradient-based edge operators Prewitt Sobel Roberts Laplacian zero-crossings Canny edge detector Hough transform for detection of straight lines Circle Hough Transform Digital Image Processing:

More information

Types of Edges. Why Edge Detection? Types of Edges. Edge Detection. Gradient. Edge Detection

Types of Edges. Why Edge Detection? Types of Edges. Edge Detection. Gradient. Edge Detection Why Edge Detection? How can an algorithm extract relevant information from an image that is enables the algorithm to recognize objects? The most important information for the interpretation of an image

More information

Fathom Dynamic Data TM Version 2 Specifications

Fathom Dynamic Data TM Version 2 Specifications Data Sources Fathom Dynamic Data TM Version 2 Specifications Use data from one of the many sample documents that come with Fathom. Enter your own data by typing into a case table. Paste data from other

More information

Blood vessel tracking in retinal images

Blood vessel tracking in retinal images Y. Jiang, A. Bainbridge-Smith, A. B. Morris, Blood Vessel Tracking in Retinal Images, Proceedings of Image and Vision Computing New Zealand 2007, pp. 126 131, Hamilton, New Zealand, December 2007. Blood

More information

Multivariate Capability Analysis

Multivariate Capability Analysis Multivariate Capability Analysis Summary... 1 Data Input... 3 Analysis Summary... 4 Capability Plot... 5 Capability Indices... 6 Capability Ellipse... 7 Correlation Matrix... 8 Tests for Normality... 8

More information

Data Management - 50%

Data Management - 50% Exam 1: SAS Big Data Preparation, Statistics, and Visual Exploration Data Management - 50% Navigate within the Data Management Studio Interface Register a new QKB Create and connect to a repository Define

More information

IMP:DisparityBox6. H. David Sheets, Dept. of Physics, Canisius College, Buffalo, NY 14208,

IMP:DisparityBox6. H. David Sheets, Dept. of Physics, Canisius College, Buffalo, NY 14208, IMP:DisparityBox6 H. David Sheets, Dept. of Physics, Canisius College, Buffalo, NY 14208, nual, resampling statistics) Miriam L. Zelditch, Museum of Paleontology, University of Michigan, Ann Arbor, Michigan

More information

SPM8 for Basic and Clinical Investigators. Preprocessing

SPM8 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 information

Detection Theory for Graphs

Detection Theory for Graphs Detection Theory for Graphs Benjamin A. Miller, Nadya T. Bliss, Patrick J. Wolfe, and Michelle S. Beard Graphs are fast emerging as a common data structure used in many scientific and engineering fields.

More information

Package EBglmnet. January 30, 2016

Package EBglmnet. January 30, 2016 Type Package Package EBglmnet January 30, 2016 Title Empirical Bayesian Lasso and Elastic Net Methods for Generalized Linear Models Version 4.1 Date 2016-01-15 Author Anhui Huang, Dianting Liu Maintainer

More information

OVERVIEW OF ESTIMATION FRAMEWORKS AND ESTIMATORS

OVERVIEW OF ESTIMATION FRAMEWORKS AND ESTIMATORS OVERVIEW OF ESTIMATION FRAMEWORKS AND ESTIMATORS Set basic R-options upfront and load all required R packages: > options(prompt = " ", digits = 4) setwd('c:/klaus/aaec5126/test/')#r Sweaves to the default

More information

SnPM is an SPM toolbox developed by Andrew Holmes & Tom Nichols

SnPM 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 information