Skip to content

Commit

Permalink
Update tmfc_LSS_after_FIR.m
Browse files Browse the repository at this point in the history
  • Loading branch information
@Masharipov
Masharipov committed Dec 12, 2024
1 parent 7fd46f5 commit bfdf67b
Showing 1 changed file with 37 additions and 6 deletions.
43 changes: 37 additions & 6 deletions functions/tmfc_LSS_after_FIR.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -147,6 +147,15 @@
for iSub = start_sub:nSub
tic
SPM = load(tmfc.subjects(iSub).path);

% Check if SPM.mat has concatenated sessions
% (if spm_fmri_concatenate.m sript was used)
if size(SPM.SPM.nscan,2) == size(SPM.SPM.Sess,2)
SPM_concat(iSub) = 0;
else
SPM_concat(iSub) = 1;
end
concat(iSub).scans = SPM.SPM.nscan;

if isdir(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')]))
rmdir(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')]),'s');
Expand Down Expand Up @@ -204,8 +213,14 @@
matlabbatch{1}.spm.stats.fmri_spec.timing.fmri_t0 = SPM.SPM.xBF.T0;

% Functional images
for image = 1:SPM.SPM.nscan(sess_num(jSess))
matlabbatch{1}.spm.stats.fmri_spec.sess.scans{image,1} = fullfile(tmfc.project_path,'FIR_regression',['Subject_' num2str(iSub,'%04.f')],['Res_' num2str(SPM.SPM.Sess(sess_num(jSess)).row(image),'%.4d') '.nii,1']);
if SPM_concat(iSub) == 0
for image = 1:SPM.SPM.nscan(sess_num(jSess))
matlabbatch{1}.spm.stats.fmri_spec.sess.scans{image,1} = fullfile(tmfc.project_path,'FIR_regression',['Subject_' num2str(iSub,'%04.f')],['Res_' num2str(SPM.SPM.Sess(sess_num(jSess)).row(image),'%.4d') '.nii,1']);
end
else
for image = 1:size(SPM.SPM.xY.VY,1)
matlabbatch{1}.spm.stats.fmri_spec.sess.scans{image,1} = fullfile(tmfc.project_path,'FIR_regression',['Subject_' num2str(iSub,'%04.f')],['Res_' num2str(SPM.SPM.Sess(jSess).row(image),'%.4d') '.nii,1']);
end
end

% Current trial vs all other trials (of interest and no interrest)
Expand Down Expand Up @@ -268,9 +283,17 @@
spm_jobman('run',batch{kTrial});

% Use explicit mask
SPM_LSS = load(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')],['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'));
SPM_LSS = load(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')], ...
['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'));
SPM_LSS.SPM.xM.TH = -Inf(size(SPM_LSS.SPM.xM.TH));
tmfc_parsave_SPM(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')],['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'),SPM_LSS.SPM);
tmfc_parsave_SPM(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')], ...
['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'),SPM_LSS.SPM);

% If concatenated
if SPM_concat(iSub) == 1
spm_fmri_concatenate(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')], ...
['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'),concat(iSub).scans);
end

% Estimate LSS GLM
spm_jobman('run',batch_2{kTrial});
Expand Down Expand Up @@ -303,9 +326,17 @@
spm_jobman('run',batch{kTrial});

% Use explicit mask
SPM_LSS = load(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')],['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'));
SPM_LSS = load(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')], ...
['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'));
SPM_LSS.SPM.xM.TH = -Inf(size(SPM_LSS.SPM.xM.TH));
tmfc_parsave_SPM(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')],['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'),SPM_LSS.SPM);
tmfc_parsave_SPM(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')], ...
['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'),SPM_LSS.SPM);

% If concatenated
if SPM_concat(iSub) == 1
spm_fmri_concatenate(fullfile(tmfc.project_path,'LSS_regression_after_FIR',['Subject_' num2str(iSub,'%04.f')], ...
['LSS_Sess_' num2str(sess_num(jSess)) '_Trial_' num2str(kTrial)],'SPM.mat'),concat(iSub).scans);
end

% Estimate LSS GLM
spm_jobman('run',batch_2{kTrial});
Expand Down

0 comments on commit bfdf67b

Please sign in to comment.