Simulating fMRI data: the R package neuRosim

26 downloads 718 Views 2MB Size Report
Summary. Simulating fMRI data: the R package. neuRosim ... Department of Data Analysis .... in-house developed software routines, often not available for the.
neuRosim

Outline Motivation The need for simulation The need for validated simulation The need for software

Simulating fMRI data: the R package neuRosim

Features Goals What can you do with neuRosim?

Marijke Welvaert

How is neuRosim organized?

Example Setting up the design

Department of Data Analysis Ghent University

Simulating the data Exporting and analyzing the data

Summary

Berlin Workshop on Statistics and Neuroimaging 2011

Outline neuRosim

1 Motivation Outline Motivation The need for simulation

The need for simulation The need for validated simulation The need for software

The need for validated simulation The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

2 Features

Goals What can you do with neuRosim? How is neuRosim organized? 3 Example

Setting up the design Simulating the data Exporting and analyzing the data 4 Summary

Knowing the ground truth in MRI neuRosim

Outline Motivation The need for simulation The need for validated simulation The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

Knowing the ground truth in fMRI neuRosim

Outline Motivation The need for simulation The need for validated simulation The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

Reflection in the literature neuRosim

Web of Science publications/year Outline

3000

Motivation The need for simulation

How is neuRosim organized?

Example Setting up the design

1500

What can you do with neuRosim?

Number of publications

Goals

1000

The need for software

Features

2000

2500

The need for validated simulation

Exporting and analyzing the data

500

Simulating the data

0

Summary

1995

2000 Year of publication

2005

2010

fMRI data components neuRosim

Activation Outline

experimentally induced

Motivation

spontaneous

The need for simulation The need for validated simulation The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design

Known artefacts B0 inhomogeneities low-frequency drift Noise system movement

Simulating the data Exporting and analyzing the data

Summary

physiological task-related ... Spatial and temporal correlations

Typical fMRI simulation studies neuRosim

Outline

1

Motivation

known activation combined with real noise e.g. Bianciardi et al. (2004), Lange (1999), Weibull et al. (2008)

The need for simulation The need for validated simulation The need for software

2

white time series known activation combined with white noise i.i.d or AR(1) Gaussian distribution e.g. Lei et al. (2010), Lin et al. (2010), Purdon & Weisskoff (1998), Smith et al. (2011)

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example

hybrid simulation

3

other

Setting up the design Simulating the data Exporting and analyzing the data

Summary

model-based simulation, Bloch equations, noise based on residuals of real data e.g. Drobnjak et al. (2006), Havlicek et al. (2010), Logan & Rowe (2004), Ramsey et al. (2010)

Problems - Discrepancies - Shortcomings neuRosim

Outline Motivation The need for simulation

real noise may contain undesired activity

The need for validated simulation

simulated noise = system noise

The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

beware of the phrase: “. . . simulations under realistic noise conditions. . . ” total ignorance of spatial context no stand-alone simulations often missing (crucial) information while reporting simulation studies

The choice of simulation model matters! neuRosim

Outline Low CNR

0.95

0.30 0.25

0.90

The need for software

Features

Power

0.20

0.75 0.70

Setting up the design

0.15

How is neuRosim organized?

Example

0.10

What can you do with neuRosim?

Power

Goals

0.85

The need for validated simulation

0.80

The need for simulation

1.00

High CNR

0.35

Motivation

White noise AR(1) noise Phys. noise Real noise

0.65

Exporting and analyzing the data

0.05

Simulating the data

0.02

Summary

0.04

0.06 Contrast−to−noise ratio

0.08

0.10

0.20

0.25

0.30

Contrast−to−noise ratio

0.35

Towards a convergence of simulation methods neuRosim

Outline Motivation The need for simulation

in-house developed software routines, often not available for the community

The need for validated simulation

language barrier

The need for software

no widespread software packages

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

But. . . POSSUM (FSL) DCM simulator (SPM) simtb (Matlab Toolbox)

neuRosim wants neuRosim

Outline Motivation The need for simulation The need for validated simulation The need for software

Features

to provide a tool for simulating fMRI data to be a base for more validated simulation studies

Goals What can you do with neuRosim?

to make simulation available for less technical researchers

How is neuRosim organized?

to allow maximum flexibility for the useRs

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

What can you do with neuRosim? neuRosim

Outline Motivation The need for simulation The need for validated simulation

specify your experimental design based on stimulus onsets and durations

The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

specify activated regions using an xyz-coordinate system simulate BOLD activation with the choice of different models simulate resting state activation (still under development) simulate fMRI noise originating from different noise sources generate fMRI data from 1D time series to 4D volume data

Low-level functions neuRosim

Outline Motivation The need for simulation

Building blocks for advanced useRs who want in-depth control over their simulation data

The need for validated simulation The need for software

Activation functions

Noise functions

stimfunction() specifydesign() specifyregion()

systemnoise() temporalnoise() spatialnoise() lowfreqdrift() physnoise() tasknoise()

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

canonicalHRF() gammaHRF() balloon()

High-level functions neuRosim

Outline Motivation The need for simulation The need for validated simulation

Direct simulation of fMRI data

The need for software

Features

Preparation functions

Simulation functions

simprepTemporal() simprepSpatial()

simTSfmri() simVOLfmri() simRestingStatefmri()

Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

Real dataset neuRosim

Consider the data from a repetition priming experiment performed using event-related fMRI (Henson et al., 2002). Outline Motivation The need for simulation The need for validated simulation The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

2 × 2 factorial design famous vs non-famous faces effect of repetition

Setting up the design (1) neuRosim

Temporal Parameters Outline Motivation The need for simulation The need for validated simulation The need for software

Features Goals What can you do with neuRosim? How is neuRosim organized?

Example Setting up the design Simulating the data Exporting and analyzing the data

Summary

R> R> R> R> + + + R> + + + R> + + + R> + + + R> R>

nscan R> R> + R> + R> + R>

region.1A.center

Suggest Documents