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)     radius = args.radius;

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

private/runCoordTable.m 12)     

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

private/classify.m      14) 

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

classify.m              20)         

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

classify.m              23)         

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

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

private/runCoordTable.m 28)       

private/runCoordTable.m 29)         extr = calculateImageData(subjectStruct{s}.volumes,subjectStruct{s}.coords,radius);

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

private/classify.m      34)         

private/runCoordTable.m 35)         disp('computing psth');
private/runCoordTable.m 36)         for iVoxel = 1:nVoxel
private/runCoordTable.m 37)             rawdata = [];
private/runCoordTable.m 38)             for iImage = 1:length(extr);
private/runCoordTable.m 39)                 tmp = extr(iImage);

private/runCoordTable.m 40)                 rawdata = [rawdata tmp.dat(iVoxel,:)];

private/runCoordTable.m 41)             end

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

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

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);