SVMCrossVal.git

private/runCoordTable.m  1) function subjectData = runCoordTable(args)

private/runCoordTable.m  2) 

private/runCoordTable.m  3) % global SVMCROSSVAL_SUBJECTSTRUCT_NAME;

private/runCoordTable.m  4) 

private/runCoordTable.m  5)     disp('run coord table')
private/runCoordTable.m  6)     

private/runCoordTable.m  7)     subjects = args.subjects;

private/runCoordTable.m  8)     nSubjects = numel(subjects);

private/runCoordTable.m  9)     sessionlist = args.sessionList;

private/runCoordTable.m 10)     pstopts = args.psthOpts;

private/runCoordTable.m 11)     

private/runCoordTable.m 12)     disp(sprintf('batch processing %g subjects.',nSubjects));

private/classify.m      13) 

private/runCoordTable.m 14)     for s = 1:nSubjects
private/runCoordTable.m 15)         subjectStruct{s}.dir = fullfile(args.basedir,cell2mat(subjects(s)));
private/runCoordTable.m 16)         d = load(fullfile(subjectStruct{s}.dir,'results','SPM.mat'));
private/runCoordTable.m 17)         subjectStruct{s}.des = d.SPM;
private/runCoordTable.m 18)         subjectStruct{s}.name = cell2mat(subjects(s));

classify.m              19)         

private/runCoordTable.m 20)         map = load(fullfile(subjectStruct{s}.dir,'results','roi','coord_map.mat'));
private/runCoordTable.m 21)         subjectStruct{s}.coords = getSubjectCoordinates(args.coords,map);

classify.m              22)         

private/runCoordTable.m 23) %         nVoxel = size(subjectStruct{s}.coords,1);

private/runCoordTable.m 24)         

private/runCoordTable.m 25)         disp('fetching volume definitions, please wait');
private/runCoordTable.m 26)         subjectStruct{s}.volumes = spm_vol(getImageFileList(subjectStruct{s}.dir,sessionlist,args.mask));

private/runCoordTable.m 27)       

private/runCoordTable.m 28)         [extr nExtractedVoxel] = calculateImageData(subjectStruct{s}.volumes,subjectStruct{s}.coords);

private/runCoordTable.m 29)          
private/runCoordTable.m 30)         pstopts.des = subjectStruct{s}.des;
private/runCoordTable.m 31)         pstopts.eventList = args.eventList;
private/runCoordTable.m 32)         pstopts.sessionList = sessionlist;

private/classify.m      33)         

private/runCoordTable.m 34)         disp(sprintf('computing psth for %g voxel',nExtractedVoxel));
private/runCoordTable.m 35)         for iVoxel = 1:nExtractedVoxel

private/runCoordTable.m 36)             rawdata = [];
private/runCoordTable.m 37)             for iImage = 1:length(extr);
private/runCoordTable.m 38)                 tmp = extr(iImage);

private/runCoordTable.m 39)                 rawdata = [rawdata tmp.dat(iVoxel)];

private/runCoordTable.m 40)             end

private/runCoordTable.m 41)             subjectStruct{s}.rawData{iVoxel} = rawdata;

private/runCoordTable.m 42)             subjectStruct{s}.pst{iVoxel} = calculatePST(args.timeline,pstopts,rawdata);
private/runCoordTable.m 43)         end

private/runCoordTable.m 44)         

private/runCoordTable.m 45)        disp(sprintf('done %g / %g',s,nSubjects));

private/runCoordTable.m 46)     end

private/runCoordTable.m 47)     

private/runCoordTable.m 48)     subjectData = subjectStruct;
private/runCoordTable.m 49) %     assignin('base',SVMCROSSVAL_SUBJECTSTRUCT_NAME,subjectStruct);