Fmri Task Design Guide

Research Planning Protocol

Before executing the domain-specific steps below, you MUST:

1. State the research question — What specific question is this analysis/paradigm addressing?
2. Justify the method choice — Why is this approach appropriate? What alternatives were considered?
3. Declare expected outcomes — What results would support vs. refute the hypothesis?
4. Note assumptions and limitations — What does this method assume? Where could it mislead?
5. Present the plan to the user and WAIT for confirmation before proceeding.

For detailed methodology guidance, see the research-literacy skill.

⚠️ Verification Notice

This skill was generated by AI from academic literature. All parameters, thresholds, and citations require independent verification before use in research. If you find errors, please open an issue.

Design Type Selection

Comparison of Design Types

Decision Tree

What is the primary goal?
 |
 +-- Detect presence/absence of activation (localization)
 | |
 | +-- Is HRF shape estimation needed?
 | |
 | +-- NO --> Block design (maximum detection power)
 | |
 | +-- YES --> Mixed design (blocks + events within blocks)
 |
 +-- Estimate trial-by-trial neural responses
 | |
 | +-- Are there enough trials (>40 per condition)?
 | |
 | +-- YES --> Rapid event-related design (jittered ISI)
 | |
 | +-- NO --> Slow event-related design (ISI > 12 s)
 |
 +-- Separate sustained state vs. transient item effects
 --> Mixed design (Petersen & Dubis, 2012)

Block Design Parameters

• Optimal block duration: 15-20 seconds for maximum detection power (Maus et al., 2010; Bandettini et al., 1993). Shorter blocks (< 12 s) reduce sensitivity because the BOLD response does not reach steady state. Longer blocks (> 30 s) increase habituation and strategy effects (Poldrack et al., 2011, Ch. 3)
• Minimum block duration: 12 seconds to allow the BOLD signal to reach near-plateau (Bandettini et al., 1993)
• Number of blocks per condition: At least 4-6 blocks per condition per run for stable estimates (Poldrack et al., 2011, Ch. 3)
• Condition alternation: Alternate conditions (ABAB or ABCABC) rather than grouping (AAABBB), which confounds condition with time (Poldrack et al., 2011, Ch. 3)
• Rest blocks: Include rest/fixation blocks of at least 12-16 seconds between active blocks to allow BOLD signal return to baseline (Glover, 1999)

Event-Related Design Parameters

Inter-Stimulus Interval (ISI) and Jittering

The ISI between events is critical for statistical efficiency and BOLD signal separability.

Jittering strategies (from most to least recommended):

1. Optimized sequences: Use design optimization tools (optseq2, NeuroDesign) to maximize efficiency for specific contrasts (Dale, 1999; Durnez et al., 2017)
2. Truncated exponential distribution: More short ISIs, fewer long ISIs; near-optimal efficiency (Hagberg et al., 2001)
3. Uniform random: Equal probability across ISI range; acceptable but suboptimal
4. Fixed ISI: Avoid for rapid event-related designs; severely reduces design efficiency

Domain warning: Jittered designs can be over 10x more efficient than fixed-ISI designs with the same mean interval (Dale, 1999). Always jitter for event-related fMRI.

HRF Timing Constraints

The BOLD hemodynamic response imposes hard constraints on fMRI design timing:

• HRF peak: 4-6 seconds after neural event onset (Glover, 1999)
• Return to baseline: 12-16 seconds after a brief event (Glover, 1999)
• BOLD nonlinearity: Responses to stimuli separated by < 2 seconds sum nonlinearly (reduced amplitude), making them harder to separate (Glover, 1999; Wager & Nichols, 2003)

Trial Count Requirements

Domain insight: These are per-condition minimums. If comparing conditions (A vs. B), each condition needs this many trials. More conditions require longer scan sessions or fewer trials per condition, creating a power trade-off.

Design Efficiency

Efficiency Calculation

Design efficiency quantifies how well a given design matrix allows detection of specific contrasts:

Detection efficiency = 1 / trace(c' * inv(X'X) * c)

where c is the contrast vector and X is the design matrix (Dale, 1999; Liu et al., 2001).

Detection vs. Estimation Trade-off

• Detection power: Ability to detect whether an effect exists. Maximized by block designs and rapid event-related designs with high event density (Liu et al., 2001)
• Estimation efficiency: Ability to accurately characterize the HRF shape. Maximized by jittered designs with sufficient ISI variability (Liu et al., 2001)
• These are inherently in tension: block designs maximize detection but cannot estimate HRF shape

Design Optimization Tools

• optseq2 (FreeSurfer): Optimizes event ordering and null events for maximum efficiency (Dale, 1999)
• NeuroDesign (Python): Genetic algorithm-based optimization (Durnez et al., 2017)
• fMRIpower: Power calculations accounting for design and temporal autocorrelation (Mumford & Nichols, 2008)

Mixed Design Specification

Mixed designs combine sustained (block-level) and transient (event-level) components (Petersen & Dubis, 2012):

1. Sustained regressor: Models the tonic task state (e.g., "task block on" vs. "rest"), boxcar convolved with HRF
2. Transient regressors: Model individual trial onsets within each block, convolved with HRF
3. These regressors are separable because they operate at different temporal frequencies

Key parameters:

• Block duration: 20-40 seconds to provide enough events within each block (Petersen & Dubis, 2012)
• Events per block: At least 4-6 for stable transient estimates
• Inter-block rest: 16-20 seconds minimum for BOLD signal recovery (Glover, 1999)

Scanner-Related Timing Constraints

Rest/Fixation Baseline

• Allocate 25-30% of total scan time to rest/fixation baseline (Friston et al., 1999)
• Rest periods serve dual purpose: allow BOLD signal recovery and provide baseline estimate
• For block designs, rest blocks of 12-16 seconds between active blocks (Glover, 1999)
• For event-related designs, null trials (fixation) distributed throughout the sequence (Friston et al., 1999)

Condition Ordering

• m-sequences: Pseudo-random sequences with optimal counterbalancing properties (Buracas & Boynton, 2002)
• Optimized pseudo-random: Generated by optimization algorithms (optseq2, NeuroDesign) to maximize design efficiency while controlling trial-order effects (Dale, 1999)
• Avoid: Simple alternation (ABABAB) which confounds condition with time, or purely random sequences which may produce long runs of the same condition

Behavioral Task Constraints

Common Pitfalls

1. Fixed ISI in event-related designs: Dramatically reduces design efficiency compared to jittered designs. Always jitter ISI for event-related fMRI (Dale, 1999)
2. Too few trials per condition: Fewer than 20 events per condition yields unreliable single-subject estimates (Desmond & Glover, 2002). Plan for at least 30 per condition
3. Ignoring HRF recovery time: Events separated by < 2 seconds produce nonlinear BOLD summation, making responses difficult to separate (Glover, 1999)
4. No baseline/rest periods: Without rest periods, the model cannot estimate absolute activation levels and efficiency drops substantially (Friston et al., 1999)
5. Confounding condition with time: Presenting all trials of one condition before another confounds the effect with scanner drift and fatigue
6. Not counterbalancing response mappings: Lateralized motor responses (left vs. right hand) produce motor cortex activation that confounds task effects
7. Ceiling/floor performance: If accuracy is near 100% or chance, there is no behavioral variance to correlate with brain activity
8. Not optimizing the design matrix: Using arbitrary event timing instead of optimized sequences wastes statistical power that could be gained at no additional cost

Minimum Reporting Checklist

Based on COBIDAS guidelines (Nichols et al., 2017) and Poldrack et al. (2008):

• Design type (block, event-related, mixed, rapid event-related)
• Block duration (for block designs) or ISI distribution parameters (for event-related)
• Number of conditions and number of trials per condition
• Stimulus duration and response window
• Jittering strategy and ISI range (min, max, mean, distribution)
• Null trial proportion and distribution
• Condition ordering method (optimized, m-sequence, pseudo-random)
• Total scan duration per run and number of runs
• Design optimization tool used (if any) and efficiency metric
• Response mapping and counterbalancing scheme
• Practice procedure (in-scanner or out-of-scanner, duration)
• TR and its relationship to stimulus timing

References

• Bandettini, P. A., Jesmanowicz, A., Wong, E. C., & Hyde, J. S. (1993). Processing strategies for time-course data sets in functional MRI of the human brain. Magnetic Resonance in Medicine, 30(2), 161-173.
• Buracas, G. T., & Boynton, G. M. (2002). Efficient design of event-related fMRI experiments using m-sequences. NeuroImage, 16(3), 801-813.
• Dale, A. M. (1999). Optimal experimental design for event-related fMRI. Human Brain Mapping, 8(2-3), 109-114.
• Desmond, J. E., & Glover, G. H. (2002). Estimating sample size in functional MRI (fMRI) neuroimaging studies: Statistical power analyses. Journal of Neuroscience Methods, 118(2), 115-128.
• Durnez, J., Blair, R., & Poldrack, R. A. (2017). NeuroDesign: Optimal experimental designs for task fMRI. bioRxiv, 119594.
• Friston, K. J., Zarahn, E., Josephs, O., Henson, R. N. A., & Dale, A. M. (1999). Stochastic designs in event-related fMRI. NeuroImage, 10(5), 607-619.
• Glover, G. H. (1999). Deconvolution of impulse response in event-related BOLD fMRI. NeuroImage, 9(4), 416-429.
• Haatveit, B. C., Sundet, K., Hugdahl, K., et al. (2010). The validity of d prime as a working memory index. Neuropsychology, 24(5), 629-640.
• Hagberg, G. E., Zito, G., Patria, F., & Sanes, J. N. (2001). Improved detection of event-related functional MRI signals using probability functions. NeuroImage, 14(5), 1193-1205.
• Liu, T. T., Frank, L. R., Wong, E. C., & Buxton, R. B. (2001). Detection power, estimation efficiency, and predictability in event-related fMRI. NeuroImage, 13(4), 759-773.
• Maus, B., van Breukelen, G. J. P., Goebel, R., & Berger, M. P. F. (2010). Optimal design of multi-subject blocked fMRI experiments. NeuroImage, 51(3), 1338-1348.
• Mumford, J. A., & Nichols, T. E. (2008). Power calculation for group fMRI studies accounting for arbitrary design and temporal autocorrelation. NeuroImage, 39(1), 261-268.
• Murphy, K., & Garavan, H. (2005). Deriving the optimal number of events for an event-related fMRI study based on the spatial extent of activation. NeuroImage, 27(4), 771-777.
• Nichols, T. E., Das, S., Eickhoff, S. B., et al. (2017). Best practices in data analysis and sharing in neuroimaging using MRI (COBIDAS). Nature Neuroscience, 20(3), 299-303.
• Petersen, S. E., & Dubis, J. W. (2012). The mixed block/event-related design. NeuroImage, 62(2), 1177-1184.
• Poldrack, R. A., Fletcher, P. C., Henson, R. N., et al. (2008). Guidelines for reporting an fMRI study. NeuroImage, 40(2), 409-414.
• Poldrack, R. A., Mumford, J. A., & Nichols, T. E. (2011). Handbook of Functional MRI Data Analysis. Cambridge University Press.
• Wager, T. D., & Nichols, T. E. (2003). Optimization of experimental design in fMRI: A general framework using a genetic algorithm. NeuroImage, 18(2), 293-309.

See references/ for detailed design optimization examples and parameter lookup tables.