Hello,
To piggyback, I wanted to confirm that you can do resting state DCM (cross-spectral) using the spm_dcm_specify for batch scripting.
It is unclear to me if I need to run spm_dcm_fmri_csd at some point? And if so I would presume it would be during 1st level estimation.
Currently my set-up is an adapted from Peter's PEB tutorial.
Do I need to change spm_dcm_specify or change some other option?
for subject = 1:nsubjects
% Load SPM
sub = project.subjects{subject,1};
subDCMdir = fullfile(project.rootdir,'DCM','RS','GLM',sub,'RSN_Str_SPMmat');
glm_dir = subDCMdir;
SPM = load(fullfile(glm_dir,'SPM.mat'));
SPM = SPM.SPM;
% Load ROIs
f = {fullfile(glm_dir,'VOI_DMN_MPFC_1.mat');
fullfile(glm_dir,'VOI_DMN_PCC_1.mat');
fullfile(glm_dir,'VOI_DMN_rLP_1.mat');
fullfile(glm_dir,'VOI_DMN_lLP_1.mat');
fullfile(glm_dir,'VOI_Limbic-Str_1.mat');
fullfile(glm_dir,'VOI_Exec-Str_1.mat')};
for r = 1:length(f)
XY = load(f{r});
xY(r) = XY.xY;
end
% Move to output directory
cd(glm_dir);
% Select whether to include each condition from the design matrix
% (Task, Pictures, Words)
include = [0 0 0]';
% Specify. Corresponds to the series of questions in the GUI.
s = struct();
s.name = 'DMN_Str_full';
s.u = include; % Conditions
s.delays = repmat(TR,1,nregions); % Slice timing for each region
s.TE = TE;
s.nonlinear = false;
s.two_state = false;
s.stochastic = false;
s.centre = true;
s.induced = 0;
s.a = a;
s.b = b;
s.c = c;
s.d = d;
DCM = spm_dcm_specify(SPM,xY,s);
% Return to script directory
cd(start_dir);
end
Best,
Victoria Okuneye
|