clear
addpath(genpath('C:\Users\nh852914\Documents\SPM\spm8\'))
directory = 'E:\Reshma\Task_horizontal\HC\';
onsets_dir = 'E:\Reshma\Beh data\HC\';
output_dir = 'E:\Reshma\Analysis\Contrasts\HC\';
error_dir = 'E:\Reshma\Analysis\Contrasts\errors\';
hard_drive_basepath = 'E:/';
scan_group = 'HC';
cd(directory);
Folders = dir;
spm_get_defaults
global defaults
for ppt = [4]
try
clearvars -except ppt directory error_dir Folders onsets_dir output_dir hard_drive_basepath scan_group
cd(directory);
Folder = Folders(ppt).name;
input_dir = [directory Folder '/'];
pptID = Folder;
output_dir2 = [output_dir pptID '/'];
cd(output_dir2)
load('SPM.mat')
if length(SPM.xCon) == 0
SPM.xCon = spm_FcUtil('Set',cname,'T','c',c(:),SPM.xX.xKXs);
else
SPM.xCon(end + 1) = spm_FcUtil('Set',cname,'T','c',c(:),SPM.xX.xKXs);
end
c = [1 0 zeros(1,29)];
cname = 'Hit_b1';
SPM.xCon(end + 1) = spm_FcUtil('Set',cname,'T','c',c(:),SPM.xX.xKXs);
c = [zeros(1,2) 1 0 zeros(1,27)];
cname = 'Hit_b2';
SPM.xCon(end + 1) = spm_FcUtil('Set',cname,'T','c',c(:),SPM.xX.xKXs);
c = [zeros(1,4) 1 0 zeros(1,25)];
cname = 'Hit_b3';
SPM.xCon(end + 1) = spm_FcUtil('Set',cname,'T','c',c(:),SPM.xX.xKXs);
cd(output_dir2)
participant_number = pptID;
copyfile('con_001.img',strcat(hard_drive_basepath, '/Reshma/Analysis/Contrasts/', scan_group, '/', participant_number, '_Hit_b1.img'));
copyfile('con_001.hdr',strcat(hard_drive_basepath, '/Reshma/Analysis/Contrasts/', scan_group, '/', participant_number, '_Hit_b1.hdr'));
copyfile('con_002.img',strcat(hard_drive_basepath, '/Reshma/Analysis/Contrasts/', scan_group, '/', participant_number, '_Hit_b2.img'));
copyfile('con_002.hdr',strcat(hard_drive_basepath, '/Reshma/Analysis/Contrasts/', scan_group, '/', participant_number, '_Hit_b2.hdr'));
copyfile('con_003.img',strcat(hard_drive_basepath, '/Reshma/Analysis/Contrasts/', scan_group, '/', participant_number, '_Hit_b3.img'));
copyfile('con_003.hdr',strcat(hard_drive_basepath, '/Reshma/Analysis/Contrasts/', scan_group, '/', participant_number, '_Hit_b3.hdr'));
catch ME
save([error_dir pptID '_error_first_level_analysis'], 'ME')
disp(['error for subject ' pptID ]);
continue
end
disp(['successfully ran analysis for subject ' pptID]);
end