The General Linear Model GLM Klaas Enno Stephan

  • Slides: 28
Download presentation
The General Linear Model (GLM) Klaas Enno Stephan Branco Weiss Laboratory (BWL) Institute for

The General Linear Model (GLM) Klaas Enno Stephan Branco Weiss Laboratory (BWL) Institute for Empirical Research in Economics University of Zurich Functional Imaging Laboratory (FIL) Wellcome Trust Centre for Neuroimaging University College London SPM Short Course, Wellcome Trust Centre for Neuroimaging October 2008

Overview of SPM Image time-series Realignment Kernel Smoothing Design matrix Statistical parametric map (SPM)

Overview of SPM Image time-series Realignment Kernel Smoothing Design matrix Statistical parametric map (SPM) General linear model Statistical inference Normalisation Gaussian field theory p <0. 05 Template Parameter estimates

A very simple f. MRI experiment One session Passive word listening versus rest 7

A very simple f. MRI experiment One session Passive word listening versus rest 7 cycles of rest and listening Blocks of 6 scans with 7 sec TR Stimulus function Question: Is there a change in the BOLD response between listening and rest?

Modelling the measured data Why? Make inferences about effects of interest How? 1. Decompose

Modelling the measured data Why? Make inferences about effects of interest How? 1. Decompose data into effects and error 2. Form statistic using estimates of effects and error stimulus function data linear model effects estimate error estimate statistic

Voxel-wise time series analysis model specification Time parameter estimation hypothesis statistic e m Ti

Voxel-wise time series analysis model specification Time parameter estimation hypothesis statistic e m Ti BOLD signal single voxel time series SPM

Time = 1 BOLD signal + 2 x 1 + x 2 error Single

Time = 1 BOLD signal + 2 x 1 + x 2 error Single voxel regression model e

Mass-univariate analysis: voxel-wise GLM y = X + Model is specified by 1. Design

Mass-univariate analysis: voxel-wise GLM y = X + Model is specified by 1. Design matrix X 2. Assumptions about e N: number of scans p: number of regressors The design matrix embodies all available knowledge about experimentally controlled factors and potential confounds.

GLM: mass-univariate parametric analysis • parametric • • • one sample t-test two sample

GLM: mass-univariate parametric analysis • parametric • • • one sample t-test two sample t-test paired t-test Anova An. Cova correlation linear regression multiple regression F-tests etc… all cases of the General Linear Model • non-parametric? Sn. PM

GLM assumes Gaussian “spherical” (i. i. d. ) errors sphericity = iid: error covariance

GLM assumes Gaussian “spherical” (i. i. d. ) errors sphericity = iid: error covariance is scalar multiple of identity matrix: Cov(e) = 2 I Examples for non-sphericity: non-identity non-independence

Parameter estimation Objective: estimate parameters to minimize = y + X Ordinary least squares

Parameter estimation Objective: estimate parameters to minimize = y + X Ordinary least squares estimation (OLS) (assuming i. i. d. error):

A geometric perspective on the GLM OLS estimates Residual forming matrix R y e

A geometric perspective on the GLM OLS estimates Residual forming matrix R y e x 2 x 1 Design space defined by X Projection matrix P

Correlated and orthogonal regressors y x 2* Correlated regressors = explained variance is shared

Correlated and orthogonal regressors y x 2* Correlated regressors = explained variance is shared between regressors x 2 x 1 When x 2 is orthogonalized with regard to x 1, only the parameter estimate for x 1 changes, not that for x 2!

What are the problems of this model? 1. BOLD responses have a delayed and

What are the problems of this model? 1. BOLD responses have a delayed and dispersed form. 2. The BOLD signal includes substantial amounts of low-frequency noise. 3. The data are serially correlated (temporally autocorrelated) this violates the assumptions of the noise model in the GLM HRF

Problem 1: Shape of BOLD response Solution: Convolution model hemodynamic response function (HRF) The

Problem 1: Shape of BOLD response Solution: Convolution model hemodynamic response function (HRF) The response of a linear time-invariant (LTI) system is the convolution of the input with the system's response to an impulse (delta function). expected BOLD response = input function impulse response function (HRF)

Convolution model of the BOLD response Convolve stimulus function with a canonical hemodynamic response

Convolution model of the BOLD response Convolve stimulus function with a canonical hemodynamic response function (HRF): HRF

Problem 2: Low-frequency noise Solution: High pass filtering S = residual forming matrix of

Problem 2: Low-frequency noise Solution: High pass filtering S = residual forming matrix of DCT set discrete cosine transform (DCT) set

High pass filtering: example blue = data black = mean + low-frequency drift green

High pass filtering: example blue = data black = mean + low-frequency drift green = predicted response, taking into account low-frequency drift red = predicted response, NOT taking into account low-frequency drift

Problem 3: Serial correlations with 1 st order autoregressive process: AR(1) autocovariance function

Problem 3: Serial correlations with 1 st order autoregressive process: AR(1) autocovariance function

Dealing with serial correlations • Pre-colouring: impose some known autocorrelation structure on the data

Dealing with serial correlations • Pre-colouring: impose some known autocorrelation structure on the data (filtering with matrix W) and use Satterthwaite correction for df’s. • Pre-whitening: 1. Use an enhanced noise model with multiple error covariance components. 2. Use estimated autocorrelation to specify filter matrix W for whitening the data.

How do we define V ? • Enhanced noise model • Remember linear transform

How do we define V ? • Enhanced noise model • Remember linear transform for Gaussians • Choose W such that error covariance becomes spherical • Conclusion: W is a function of V so how do we estimate V ?

Multiple covariance components enhanced noise model V = 1 error covariance components Q and

Multiple covariance components enhanced noise model V = 1 error covariance components Q and hyperparameters Q 1 + 2 Q 2 Estimation of hyperparameters with Re. ML (restricted maximum likelihood).

Contrasts & statistical parametric maps c=100000 Q: activation during listening ? Null hypothesis:

Contrasts & statistical parametric maps c=100000 Q: activation during listening ? Null hypothesis:

t-statistic based on ML estimates c=100000 For brevity: Re. MLestimates

t-statistic based on ML estimates c=100000 For brevity: Re. MLestimates

Physiological confounds • head movements • arterial pulsations • breathing • eye blinks •

Physiological confounds • head movements • arterial pulsations • breathing • eye blinks • adaptation affects, fatigue, fluctuations in concentration, etc.

Outlook: further challenges • correction for multiple comparisons • variability in the HRF across

Outlook: further challenges • correction for multiple comparisons • variability in the HRF across voxels • slice timing • limitations of frequentist statistics Bayesian analyses • GLM ignores interactions among voxels models of effective connectivity These issues are discussed in future lectures.

Correction for multiple comparisons • Mass-univariate approach: We apply the GLM to each of

Correction for multiple comparisons • Mass-univariate approach: We apply the GLM to each of a huge number of voxels (usually > 100, 000). • Threshold of p<0. 05 more than 5000 voxels significant by chance! • Massive problem with multiple comparisons! • Solution: Gaussian random field theory

Variability in the HRF • HRF varies substantially across voxels and subjects • For

Variability in the HRF • HRF varies substantially across voxels and subjects • For example, latency can differ by ± 1 second • Solution: use multiple basis functions • See talk on event-related f. MRI

Summary • Mass-univariate approach: same GLM for each voxel • GLM includes all known

Summary • Mass-univariate approach: same GLM for each voxel • GLM includes all known experimental effects and confounds • Convolution with a canonical HRF • High-pass filtering to account for low-frequency drifts • Estimation of multiple variance components (e. g. to account for serial correlations) • Parametric statistics