SVMCrossVal.git

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
function subjectStruct  = runROIImageMaskMode(args)
subjects = args.subjects;
   
nSubjects = numel(subjects);
sessionlist = args.sessionList;
pstopts = args.psthOpts;
pstopts.eventList = args.eventList;
pstopts.sessionList = sessionlist;
disp(sprintf('batch processing  %g subjects.',nSubjects));
for s = 1:nSubjects
    subjectStruct{s}.dir = fullfile(args.basedir,cell2mat(subjects(s)));
    d = load(fullfile(subjectStruct{s}.dir,'results','SPM.mat'));
    subjectStruct{s}.des = d.SPM;
    subjectStruct{s}.name = cell2mat(subjects(s));
    subjectStruct{s}.roiFile = ui_selectRoiImage(...
        sprintf('Select ROI Files for %s',subjectStruct{s}.name),...
        fullfile(subjectStruct{s}.dir,'results','roi'));
end
for s = 1:nSubjects
    % load image data
   
    disp('fetching volume definitions, please wait');
    subjectStruct{s}.volumes = spm_vol(getImageFileList(subjectStruct{s}.dir,sessionlist,args.mask));
    disp('computing volume values, please wait');
    
    [extr x y z] = calculateRoiImageData(subjectStruct{s}.volumes,subjectStruct{s}.roiFile);
    
    nVoxel = size(extr(1).dat,1);
    % calculate psth
    pstopts.des = subjectStruct{s}.des;
    disp('computing psth');
    for iVoxel = 1:nVoxel
        [ x(iVoxel) y(iVoxel) z(iVoxel)]
        rawdata = [];
        for iImage = 1:length(extr);
            tmp = extr(iImage);
            rawdata = [rawdata tmp.dat(iVoxel)];
        end
        subjectStruct{s}.pst{iVoxel} = calculatePST(args.timeline,pstopts,rawdata);
    end
    disp('done');
end
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX