experimental design for functional mri david glahn updated by jll
TRANSCRIPT
![Page 1: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/1.jpg)
Experimental Designfor
Functional MRI
David GlahnUpdated by JLL
![Page 2: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/2.jpg)
General Experimental Design
- Neuropsychology -
• What is the question?• What are appropriate controls? • Which imaging modality?• Study style?
![Page 3: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/3.jpg)
Experimental Design: Terminology
• Variables– Independent vs. Dependent – Categorical vs. Continuous
• Contrasts– Experimental vs. Control– Parametric vs. subtractive
• Comparisons of subjects– Between- vs. Within-subjects
• Confounding factors • Randomization,
counterbalancing
From Scott Huettel, Duke
![Page 4: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/4.jpg)
Donder’s Method: Subtraction
• A random series of A’s and B’s presented and the subject must:– Task 1 - Respond whenever event A or B occurs (RT1)– Task 2 - Respond only to A not to B (RT2)– Task 3 - Respond X to A and Y to B (RT3)
RT = reaction time
• RT1 = T-detect + T-response• RT2 = T-detect + T-discrimination + T-response• RT3 = T-detect + T-discrimination + T-choice + T-response• T-discrimination = RT2 - RT1
• T-choice = RT3 - RT2
Example: How long does it take to choose between alternatives? (Mental Chronometry)
![Page 5: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/5.jpg)
Criticisms of Donder
• Assumes that adding components does not affect other components (i.e. assumption of pure insertion)
• One should pick tasks that differ along same dimension (time in our example)
• Although resting baseline is good to include, it may limit inference (e.g. Sternberg, 1964)
![Page 6: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/6.jpg)
What types of hypotheses are possible for fMRI data?
From Scott Huettel, Duke
![Page 7: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/7.jpg)
Experimental Design for fMRI
Must Account for Hemodynamic Response(HR)
Savoy et al., 1995
![Page 8: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/8.jpg)
Linear Systems Analysis Boynton et al. 1996
• The linear transform model of fMRI hypothesizes that responses are proportional to local average neural activity averaged over a period of time. – fMRI responses in human primary visual cortex (V1) depend on
both stimulus timing (8 Hz) and stimulus contrast (black/white). – Responses to long-duration stimuli can be predicted from a
hemodynamic response function (HRF) derived from shorter duration stimuli.
– The noise in the fMRI data is independent of stimulus contrast and stimulus temporal period.
• Because the linear transform model is consistent with our data, we proceeded to estimate the temporal fMRI response function and the underlying (presumably neural) contrast response function using HRF…
• Assumption is that HRF is linear and shift-invariant!
![Page 9: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/9.jpg)
Linearity of BOLD responseDale & Buckner, 1997
Sync each differential response to start of trial
Not quite linear but good enough for first order approximations
Reversing Checkerboard (8 Hz)
One-trial = 1 stimulus
Two-trial – 2 stimuli
Three-trial = 3 stimuli
Stim duration (SD) = 1 s
Inter-stim interval (ISI) = 2 s
![Page 10: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/10.jpg)
fMRI Design Types
1) Blocked Designs2) Event-Related Designs
a) Periodic Single Trial b) Jittered Single Trial
3) Mixed Designs- Combination blocked/event-
related
![Page 11: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/11.jpg)
Blocked Designs
![Page 12: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/12.jpg)
What are Blocked Designs?
• Blocked designs segregate different cognitive tasks into distinct time periods (blocks)
Task A Task B Task A Task B Task A Task B Task A Task B
Task A Task BREST REST Task A Task BREST REST
Paradigm – pattern or model; detailed plan for the experiment
fMRI brain images acquired continuously
![Page 13: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/13.jpg)
“Loose” vs. “Tight” Block Designs
• Loose: 1 Task, 1 contrast (with Baseline)
• Tight: more than 1 Task, multiple contrasts (including baseline)
![Page 14: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/14.jpg)
Types of Blocked Design
• Task A vs. Task B (… vs. Task C…)– Example: Squeezing Right Hand vs. Left Hand– Allows you to distinguish differential activation
between conditions– Does not allow identification of activity common to
both tasks• Can control for uninteresting activity
• Task A vs. No-task (… vs. Task C…)– Example: Squeezing Right Hand vs. Rest– Shows you activity associated with task– May introduce unwanted results if not matched
properly(example would be if rest acquired with eyes closed
but task had eyes open)
![Page 15: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/15.jpg)
Adapted from Gusnard & Raichle (2001)
(E - Bad Control Design)
![Page 16: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/16.jpg)
Adapted from Gusnard & Raichle (2001)
Oxygen Extractio
n Fraction
Cerebral Metabolic Rate of
O2
Cerebral Blood Flow
A True Baseline?
Depends on what is measured!
Different Areas may have different baselines
![Page 17: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/17.jpg)
Power in Blocked Designs
1. Summation of responses results in large signals then plateaus (~10 sec)
1. Response Duration does not plateau and onset does not change
Stimulus duration
and interval
compared with HRFISI = 1 sec
![Page 18: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/18.jpg)
Choosing Length of Blocks• Longer block lengths allow for stability of extended responses
– Hemodynamic response saturates following extended stimulation• After about 10s, activation reaches plateau
– Many tasks require extended intervals• Brain processing may differ throughout the task period
• Shorter block lengths move your signal to higher temporal frequencies– Away from low-frequency noise: scanner drift, etc.– Not possible in O-15 PET rCBF studies
• Periodic blocks may result in aliasing of other periodic signals in the data– Example: if the person breathes at a regular rate of 12 per min
and the blocks are 10s long (6 blocks/min)– Could be problem if the aliased signal falls within the range of
desired signals
From Scott Huettel, Duke
![Page 19: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/19.jpg)
What are the temporal limits?What is the shortest stimulus duration that fMRI can detect?
Blamire et al. (1992) – 2 secBandettini (1993): 0.5 secSavoy et al (1995): 34 msec
• With enough averaging, anything seems possible.
• Assume that the shape of the HRF is predictable.
• Event-related potentials (ERPs) are based on averaging small responses over many trials.
• Can we do the same thing with fMRI?
![Page 20: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/20.jpg)
Assumption of steady-state dynamics.
For block designs we assume that the BOLD effect remains constant across the epoch of interest.
For PET this assumption is valid given the half-life of the radiotracer used for CBF studies, task designs, and the time period for the image acquisition.
But the BOLD response is much more transient and more importantly may vary according to brain regions and stimulus durations and maybe even stimulus types.
Savoy et al., 1995
![Page 21: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/21.jpg)
Limitations of Blocked Designs
• Sensitive to signal drift or MR instability
• Poor choice of conditions/baseline may preclude meaningful conclusions
• Many tasks cannot be conducted well repeatedly
![Page 22: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/22.jpg)
Non-Task Brain Processing• In experiments activation can be greater in
baseline conditions than in task conditions!– Requires different processing for interpretation
• Suggests the idea of baseline/resting mental processes– Gathering/evaluation about the world around you– Awareness (of self)– Online monitoring of sensory information– Daydreaming– Neurons that are wired together fire together
• This collection of resting state brain processes is often called the “Default Mode Network” (DMN)
![Page 23: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/23.jpg)
Default Mode!
Damoiseaux 2006 analyzed separate 10-subject resting-state data sets, using
Independent Components analysis (ICA).
Vision.
Frontal Parietal
Resting State Networks (RSNs)Resting State Networks (RSNs)
![Page 24: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/24.jpg)
Event-Related Designs
![Page 25: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/25.jpg)
Buckner et al., 1998
Event RelatedEvent RelatedEvent RelatedEvent Related
![Page 26: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/26.jpg)
What are Event-Related Designs?
• Event-related designs associate brain processes with discrete events, which may occur at any point in the scanning session.
• Can detect transient BOLD responses• Supports adapting task to response such as changing
difficulty based on error rate
![Page 27: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/27.jpg)
Why use event-related designs?
• Some experimental tasks are naturally event-related (future stimuli based on response)
• Allows studying within-trial effects• Improves relation to behavioral
factors (behavior changes within blocks may be masked)
• Simple analyses– Selective averaging– General linear models (GLM)
![Page 28: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/28.jpg)
Same Event
Averaging
Sorting Into Common Groups
- Behavior
- Physiological Measure
- Outlier Rejection
- Transient vs. Task level Responses
![Page 29: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/29.jpg)
Periodic Single Trial Designs
• Stimulus events presented infrequently with long inter-stimulus intervals (ISIs)
500 ms 500 ms 500 ms 500 ms
18 s 18 s 18 s
![Page 30: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/30.jpg)
Trial Spacing Effects: Periodic Designs
ISI = 8sec (~12 trials) ISI = 4sec (~45 trials)
ISI = 20sec (9 trials) ISI = 12sec (15 trials)
A20
A4
A8
A12
Want to maximize amplitude times number of trials per study
![Page 31: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/31.jpg)
Bandettini & Cox, 2000 • The optimal inter-stimulus interval (ISI) for a stimulus duration (SD), was determined.
• Empirical Observation: For SD=2sec, ISI=12 to 14 sec.• Theory Predicts: For SD<=2 sec, the optimal repetition interval (RI=ISI+SD)• Theory Predicts: For SD>2sec, RI = 8+(2*SD).
• The statistical power of ER-fMRI relative to blocked-design was determined
• Empirical: For SD=2 sec, ER-fMRI was ~35% lower than that of blocked-design • Simulations that assumed a linear system demonstrated estimate ~65% reduction in power• Difference suggest that the ER-fMRI amplitude is greater than
that predicted by a linear shift-invariant system models.
![Page 32: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/32.jpg)
Jittered Single Trial Designs
• Varying the timing of trials within a run• Varying the timing of events within a trial
Trial 1 Trial 2 Trial 3 Trial 4
2 events 3 events 2 events1 event
![Page 33: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/33.jpg)
Effects of Jittering on Response
Stimulus
Response
Jittering allows us to sample BOLD response in more states
![Page 34: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/34.jpg)
Effects of ISI on Detectability
Birn et al, 2002
Jittered ISI
Constant ISI
Detectability
Estimated
Accuracy of
HRF
Max when ½ stims are task state and ½
stims are control state
![Page 35: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/35.jpg)
Dale and Buckner (1997)
Detecting Using Selective Averaging
Low Response
Fewer Samples
Mid Response
More Samples
Large Response
Most samples
Visual stim duration = 1 s; acquisition 240 sec
Trials subtracted then correlation analysis with predicted response
![Page 36: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/36.jpg)
Variability of HRF: EvidenceAguirre, Zarahn & D’Esposito, 1998• HRF shows considerable variability between subjects
• Within subjects, responses are more consistent, although there is still some variability between sessions
different subjects
same subject, same session same subject, different session
![Page 37: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/37.jpg)
Variability of HRF: ImplicationsAguirre, Zarahn & D’Esposito, 1998• Generic HRF models (gamma functions) account for 70% of variance• Subject-specific models account for 92% of the variance (22% more!)• Poor modeling reduces statistical power• Less of a problem for block designs than event-related (do you know why?)• Biggest problem with delay tasks where an inappropriate estimate of the initial and final components contaminates the delay component
• Possible solution: model the HRF individually for each subject
• Possible caveat: HRF may also vary between areas, not just subjects• Buckner et al., 1996:
• noted a delay of 0.5-1 sec between visual and prefrontal regions• vasculature difference?• processing latency?
• Bug or feature? • Menon & Kim – mental chronometry
![Page 38: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/38.jpg)
Post-Hoc Sorting of Trials
From Kim and Cabeza, 2007
Using information about fMRI activation at memory encoding to predict behavioral performance at
memory retrieval.
True Memory FormationTrue Memory Formation
vs.vs.
False Memory FormationFalse Memory Formation
![Page 39: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/39.jpg)
Limitations of Event-Related Designs
• Low power (maybe)– Collecting lots of data, many runs
• The key issues are:– Can my subjects perform the task as
designed?– Are the processes of interest independent
from each other (in time, amplitude, etc.)?
![Page 40: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/40.jpg)
Mixed Designs
![Page 41: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/41.jpg)
Mixed: Combination Blocked/Event
• Both blocked and event-related design aspects are used (for different purposes)– Blocked design: state-dependent effects – Event-related design: item-related effects
• Analyses can model these as separate phenomena, if cognitive processes are independent.– “Memory load effects” vs. “Item retrieval effects”
• Or, interactions can be modeled.– Effects of memory load on item retrieval activation.
![Page 42: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/42.jpg)
Blocked (solid)
Event-Related (dashed)
Event-related model reaches peak sooner…
… and returns to baseline more
slowly.
In this study, some language-related
regions were better modeled by event-
related.
From Mechelli, et al., 2003
You can model a block with events…
![Page 43: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/43.jpg)
Mixed Design
![Page 44: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/44.jpg)
Summary of Experiment Design
• Main Issues to Consider– What design constraints are induced by my task?– What am I trying to measure?– What sorts of non-task-related variability do I want to
avoid?
• Rules of thumb– Blocked Designs:
• Powerful for detecting activation• Useful for examining state changes
– Event-Related Designs: • Powerful for estimating time course of activity• Allows determination of baseline activity• Best for post hoc trial sorting
– Mixed Designs• Best combination of detection and estimation• Much more complicated analyses
![Page 45: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/45.jpg)
What is fMRI Experimental Design?
• Controlling the timing and quality of cognitive operations to influence brain activation
• What can we control?– Stimulus properties (what is presented?)– Stimulus timing (when is it presented?)– Subject instructions (what do subjects do with it?)
• What are the goals of experimental design?– To test specific hypotheses (i.e., hypothesis-driven)– To generate new hypotheses (i.e., data-driven)
![Page 46: Experimental Design for Functional MRI David Glahn Updated by JLL](https://reader036.vdocument.in/reader036/viewer/2022062423/56649ec45503460f94bce01f/html5/thumbnails/46.jpg)
Experimental Designfor
Functional MRI
David GlahnUpdated by JLL