- 87 Views
- Uploaded on
- Presentation posted in: General

Experimental design of fMRI studies

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Experimental design of fMRI studies

Klaas Enno Stephan Laboratory for Social and Neural Systems Research

Institute for Empirical Research in Economics

University of Zurich

Functional Imaging Laboratory (FIL)

Wellcome Trust Centre for Neuroimaging

University College London

With many thanks for slides & images to:

FIL Methods group,

particularly Christian Ruff

Methods & models for fMRI data analysis12 November 2008

Overview of SPM

Statistical parametric map (SPM)

Design matrix

Image time-series

Kernel

Realignment

Smoothing

General linear model

Gaussian

field theory

Statistical

inference

Normalisation

p <0.05

Template

Parameter estimates

- Categorical designs
- Subtraction - Pure insertion, evoked / differential responses
- Conjunction - Testing multiple hypotheses
- Parametric designs
- Linear - Adaptation, cognitive dimensions
- Nonlinear- Polynomial expansions, neurometric functions
- Factorial designs
- Categorical- Interactions and pure insertion
- Parametric- Linear and nonlinear interactions
- - Psychophysiological Interactions

- Example:

- Neuronal structures computing face recognition?

- Aim:
- Neuronal structures underlying a single process P?

- Procedure:
- Contrast: [Task with P] – [control task without P ] = P
- the critical assumption of „pure insertion“

- „Related“ stimuli

- P implicit in control task ?

„Queen!“ „Aunt Jenny?“

- Same stimuli, different task

- Interaction of process and task?

Name Person! Name Gender!

- „Distant“ stimuli

- Several components differ!

SPM{F} testing for evoked responses

Parahippocampal responses to words

- “Baseline” here corresponds to session mean
- Estimation depends crucially on inclusion of null events or long SOAs
- “Cognitive” interpretation difficult, but useful to define regions generally involved in the task

BOLD EPI fMRI at 2T, TR 3.2sec. Words presented every 16 secs; (i) studied words or (ii) new words

SPM{F} testing for evoked responses

Parahippocampal responses to words

SPM{F} testing for differences

studied words

new words

BOLD EPI fMRI at 2T, TR 3.2sec. Words presented every 16 secs; (i) studied words or (ii) new words

Peri-stimulus time {secs}

Experimental design

Word generationG

Word repetitionR

R G R G R G R G R G R G

G - R = Intrinsic word generation

…under assumption of pure insertion

- Categorical designs
- Subtraction - Pure insertion, evoked / differential responses
- Conjunction - Testing multiple hypotheses
- Parametric designs
- Linear - Adaptation, cognitive dimensions
- Nonlinear- Polynomial expansions, neurometric functions
- Factorial designs
- Categorical- Interactions and pure insertion
- Parametric- Linear and nonlinear interactions
- - Psychophysiological Interactions

- One way to minimise the baseline/pure insertion problem is to isolate the same process by two or more separate comparisons, and inspect the resulting simple effects for commonalities
- A test for such activation common to several independent contrasts is called “conjunction”
- Conjunctions can be conducted across a whole variety of different contexts:
- tasks
- stimuli
- senses (vision, audition)
- etc.

- Note: the contrasts entering a conjunction must be orthogonal !

Task (1/2)

Viewing Naming

A1

A2

Objects Colours

Stimuli (A/B)

B1

B2

Common object

recognition response (R)

Price et al. 1997

Example: Which neural structures support object recognition, independent of task (naming vs viewing)?

Visual ProcessingV

Object Recognition R

Phonological RetrievalP

Object viewingV,R

Colour viewingV

Object namingP,V,R

Colour namingP,V

(Object - Colour viewing) [1 -1 0 0]

&

(Object - Colour naming) [0 0 1 -1]

[ V,R - V ] & [ P,V,R - P,V ] = R & R = R

B1-B2

A1-A2

Two types of conjunctions

- Test of global null hypothesis: Significant set of consistent effects
- “Which voxels show effects of similar direction (but not necessarily individual significance) across contrasts?”
- Null hypothesis: No contrast is significant: k = 0
- does not correspond to a logical AND !
- Test of conjunction null hypothesis: Set of consistently significant effects
- “Which voxels show, for each specified contrast, significant effects?”
- Null hypothesis: Not all contrasts are significant: k < n
- corresponds to a logical AND

p(A1-A2) <

+

+

p(B1-B2) <

Friston et al.(2005). Neuroimage, 25:661-667.

Nichols et al. (2005). Neuroimage, 25:653-660.

specificity

specificity

Global null:

k = 0

Conjunction null:

k < n

Friston et al. 2005, Neuroimage, 25:661-667.

Friston et al. 2005, Neuroimage, 25:661-667.

Friston et al. 2005, Neuroimage, 25:661-667.

- Categorical designs
- Subtraction - Pure insertion, evoked / differential responses
- Conjunction - Testing multiple hypotheses
- Parametric designs
- Linear - Adaptation, cognitive dimensions
- Nonlinear- Polynomial expansions, neurometric functions
- Factorial designs
- Categorical- Interactions and pure insertion
- Parametric- Linear and nonlinear interactions
- - Psychophysiological Interactions

- Parametric designs approach the baseline problem by:
- Varying the stimulus-parameter of interest on a continuum, in multiple (n>2) steps...
- ... and relating measured BOLD signal to this parameter

- Possible tests for such relations are manifold:
- Linear
- Nonlinear: Quadratic/cubic/etc. (polynomial expansion)
- Model-based (e.g. predictions from learning models)

Quadratic

Linear

SPM{F}

F-contrast [0 1 0] on

quadratic parameter

→ inverted ‘U’ response to

increasing word presentation

rate in the DLPFC

Polynomial expansion:

f(x) ~ b1 x + b2 x2 + b3 x3 ...

SPM5 offers polynomial expansion as option for parametric modulation of regressors

Büchel et al. 1996, NeuroImage 4:60-66

Büchel et al. 1998, NeuroImage 8:140-148

Polynomial expansion

&

orthogonalisation

Büchel et al. 1998, NeuroImage 8:140-148

Stimulus

awareness

Stimulus

intensity

Pain

intensity

Pain threshold: 410 mJ

- P0-P4: Variation of intensity of a laser stimulus applied to the right hand (0, 300, 400, 500, and 600 mJ)

P3

P2

P1

P4

Büchel et al. 2002, J. Neurosci. 22:970-976

Stimulus intensity

Stimulus presence

Pain intensity

Büchel et al. 2002, J. Neurosci. 22:970-976

- general idea:generate predictions from a computational model, e.g. of learning or decision-making
- Commonly used models:
- Rescorla-Wagner learning model
- temporal difference (TD) learning model

- use these predictions to define regressors
- include these regressors in a GLM and test for significant correlations with voxel-wise BOLD responses

TD model of reinforcement learning

- appetitive conditing with pleasant taste reward
- activity in ventral striatum and OFC correlated with TD prediction error at the time of the CS

O‘Doherty et al. 2003, Neuron

- Categorical designs
- Subtraction - Pure insertion, evoked / differential responses
- Conjunction - Testing multiple hypotheses
- Parametric designs
- Linear - Adaptation, cognitive dimensions
- Nonlinear- Polynomial expansions, neurometric functions
- Factorial designs
- Categorical- Interactions and pure insertion
- Parametric- Linear and nonlinear interactions
- - Psychophysiological Interactions

Task (1/2)

Viewing Naming

A1

A2

Objects Colours

Stimuli (A/B)

B1

B2

- Main effect of task:(A1 + B1) – (A2 + B2)
- Main effect of stimuli: (A1 + A2) – (B1 + B2)
- Interaction of task and stimuli:Can show a failure of pure insertion
- (A1 – B1) – (A2 – B2)

interaction effect

(Stimuli x Task)

Colours Objects Colours Objects

Viewing

Naming

V+

V-

V+

V-

V+

V-

V+

V-

A+

A-

A+

A-

Rescorla-Wagner learning

curves (=0.075):

Primary visual cortex

Significant four-way interaction in V1 and putamen:

V1

- during learning, predictive tones
- increasingly activate V1 & PUT in the absence of predicted visual stimuli,
- increasingly deactivate V1 & PUT in the presence of predicted vis. stimuli

PUT

p < 0.05, corrected

random effects, n=16

den Ouden et al. 2009, Cereb. Cortex

Task factor

Task B

Task A

TA/S1

TB/S1

Stim 1

Stimulus factor

Stim 2

TB/S2

TA/S2

Psycho-physiological interactions (PPI)

GLM of a 2x2 factorial design:

main effect

of task

main effect

of stim. type

interaction

main effect

of task

We can replace one main effect in the GLM by the time series of an area that shows this main effect.

E.g. let's replace the main effect of stimulus type by the time series of area V1:

V1 time series main effect

of stim. type

psycho-

physiological

interaction

SPM{Z}

V5 activity

time

V1

V5

V5

attention

V5 activity

no attention

V1 activity

Attention

=

V1 x Att.

Friston et al. 1997, NeuroImage 6:218-229

Büchel & Friston 1997, Cereb. Cortex 7:768-778

V1

V5

V5

V1

attention

attention

Two possible interpretations of the PPI term:

V1

V1

Modulation of V1V5 by attention

Modulation of the impact of attention on V5 by V1.