call_create_hrfpl_maps.m

% Script to call function 'create_hrfpl_maps.m'. Selects betas, names
% conditions and loops over subjects. Example for original sample

clear all

Root        = '/Users/immanuel_elbau/Desktop/PNAS_Repository/FirstLevels_Original/'; % Adapt path 

Conditions  = {'PreStress', 'Stress', 'PostStress'};

Subs        = dir([Root,'IST*']);

%% User-Area
% Enter file names of canonical and derivative beta-images per condition/
% session{n}. Note that the relevant betas correspond to the
% mental-arithmetic regressor.

% For Experimental group (IST0001 - IST0040)

Canonical{1}     = ['beta_0001.img'];
Canonical{2}     = ['beta_0013.img'];
Canonical{3}     = ['beta_0025.img'];

Derivative{1}    = ['beta_0002.img'];
Derivative{2}    = ['beta_0014.img'];
Derivative{3}    = ['beta_0026.img'];

% Cave: For Control group (IST0041 - IST0063) - no verbal-feedback
% regressors:

% Canonical{1}     = ['beta_0001.img'];
% Canonical{2}     = ['beta_0011.img'];
% Canonical{3}     = ['beta_0021.img'];
%
% Derivative{1}    = ['beta_0002.img'];
% Derivative{2}    = ['beta_0012.img'];
% Derivative{3}    = ['beta_0022.img'];

%% Loop over subjects and sessions and call of function

for master = 1:numel(Subs)
    
    Imagefolder    = [Root,Subs(master).name,'/IST/residualized/firstlevel_multi_16_tempder_new/'];
        
    disp(['### Working on ', Imagefolder, ' ###']);
    
    for z = 1:3
        
        create_hrfpl_maps([Imagefolder,Canonical{z}], [Imagefolder,Derivative{z}], Conditions{z});
        
    end
    
end