SVMCrossVal.git

1) function subjectStruct  = runROIImageMaskMode(args)

2) 

3) % global SVMCROSSVAL_SUBJECTSTRUCT_NAME;

4) 

5) subjects = args.subjects;

6)    

7) nSubjects = numel(subjects);

8) sessionlist = args.sessionList;

9) 

10) disp(sprintf('we have %g subjects. Press ANY-Key to continue.\n Use Retrun if your Keyboard lacks the ANY-Key.',nSubjects));
11) pause

12) 
13) for s = 1:nSubjects
14)     subjectStruct{s}.dir = fullfile(args.basedir,cell2mat(subjects(s)));

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

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

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

30) 
31)     disp('computing volume values, please wait');

32)     

33)     extr = calculateRoiImageData(subjectStruct{s}.volumes,subjectStruct{s}.roiFile);
34)     
35)     nVoxel = size(extr(1).dat,1);
36) 

37)     % calculate psth

38)     pstopts.des = subjectStruct{s}.des;
39)     pstopts.eventList = args.eventList;
40)     pstopts.sessionList = sessionlist;
41) 
42)     disp('computing psth');

43)     for iVoxel = 1:nVoxel
44)         rawdata = [];
45)         for iImage = 1:length(extr);
46)             tmp = extr(iImage);
47)             rawdata = [rawdata tmp.dat(iVoxel)];
48)         end
49)         subjectStruct{s}.pst{iVoxel} = calculatePST(args.timeline,pstopts,rawdata);
50)     end

51)     disp('done');

52) end
53) 

54) % assignin('base',SVMCROSSVAL_SUBJECTSTRUCT_NAME,subjectStruct);

55)