Data Analysis¶
fMRI Preprocessing¶
fmriprep¶
fmriprep
is a pipeline developed by the Poldrack lab at Stanford University for use at the Center for Reproducible
Neuroscience (CRN), as well as for
open-source software distribution.
fmriprep
is designed to provide an easily accessible,
state-of-the-art interface that is robust to variations in scan acquisition
protocols and that requires minimal user input, while providing easily
interpretable and comprehensive error and output reporting.
It performs basic processing steps (coregistration, normalization, unwarping, noise component extraction, segmentation, skullstripping etc.) providing outputs that can be easily submitted to a variety of group level analyses, including task-based or resting-state fMRI, graph theory measures, surface or volume-based statistics, etc.
The fmriprep
workflow takes as principal input the path of the dataset
that is to be processed.
The input dataset is required to be in valid BIDS (Brain Imaging Data
Structure) format, and it must include at least one T1w structural image and
(unless disabled with a flag) a BOLD series.
We highly recommend that you validate your dataset with the free, online
BIDS Validator.
The exact command to run fmriprep
depends on the Installation method.
The common parts of the command follow the BIDS-Apps definition.
Example:
fmriprep data/bids_root/ out/ participant -w work/
GLMs¶
General linear model scripts were run using MATLAB
and SPM 8
.
The model regressors are specificed by the files ending in analyze2
.
The contrasts are calculated in the files that start with contrast2
.
The second level/group analyses are performed by the rfx_par
script.
While the analyze
and contrast
scripts can be run just with the
function, you need to use the following syntax to run the rfx_par
script. Note that you need to provide access to a contrast file.
Example:
f = fullfile('8_pre_liking', preproc_version, 'm8_pre_liking_cons.mat');
load(f);
for con = 1:length(cname)
rfx_par('8_pre_liking',cname(con),good_subjects,preproc_version)
end
Behavioral¶
A Jupyter
notebook (using an R
kernel) for the behavioral results.
DDM¶
We fit both a base model and constant model (with an additional constant parameter that is added to the drift).
The DDM model scripts.
A Jupyter
notebook (using an R
kernel) for the ddm results.
Neural¶
A Jupyter
notebook (using an R
kernel) for the neural results.
Correlations¶
TO BE ADDED…