SVMCrossVal.git

1) function runROIImageMaskMode(args)

2) 

3) subjects = args.subjects;

4)    

5) nSubjects = size(subjects);

6) sessionlist = args.sessionList;

7) 
8) 
9) for s = 1:nSubjects
10)     subjectStruct{s}.dir = fullfile(args.basedir,cell2mat(subjects(s)));

11)     d = load(fullfile(subjectStruct{s}.dir,'results','SPM.mat'));
12)     subjectStruct{s}.des = d.SPM;

13)     subjectStruct{s}.name = cell2mat(subjects(s));
14)     subjectStruct{s}.roiFile = ui_selectRoiImage(...
15)         sprintf('Select ROI Files for %s',subjectStruct{s}.name),...
16)         fullfile(subjectStruct{s}.dir,'results','roi'));
17) end
18) 
19) 
20) 
21) for s = 1:nSubjects
22)     % load image data
23)    
24)     disp('fetching volume definitions, please wait');

25)     subjectStruct{s}.volumes = spm_vol(getImageFileList(subjectStruct{s}.dir,sessionlist,args.mask));

26) 
27)     disp('computing volume values, please wait');

28)     
29)     rawData = calculateRoiImageData(subjectStruct{s}.volumes,subjectStruct{s}.roiFile);

30)     % calculate psth

31)     pstopts.des = subjectStruct{s}.des;
32)     pstopts.eventList = args.eventList;
33)     pstopts.sessionList = sessionlist;
34)     
35)     assignin('base','pstopts',pstopts);
36)     assignin('base','rawData',rawData);
37) 
38)     disp('computing psth');
39)     subjectStruct{s}.psth    = calculatePST(args.timeline,pstopts,rawData);
40)     disp('done');

41) end
42) 

43) assignin('base','subjectStruct',subjectStruct);
44)