SVMCrossVal.git

private/runCoordTable.m  1) function ppsth = runCoordTable(baseDir,subjects,timeline,coordinates,volumeFileList,classDef)
private/runCoordTable.m  2) %         paramModel = varargin{1};
private/runCoordTable.m  3)     disp('run coord table')
private/runCoordTable.m  4)     
private/runCoordTable.m  5)     disp('does not work yet.')
private/runCoordTable.m  6)     return
private/runCoordTable.m  7)         
private/runCoordTable.m  8)         PROJECT_BASE_PATH = baseDir;

classify.m               9)         PROJECT_RESULT_PATH = 'results\SPM.mat';

private/runCoordTable.m 10) 

private/classify.m      11) 

private/runCoordTable.m 12)         timelineParams = timeline;

private/classify.m      13)         

classify.m              14)         % common params

classify.m              15)         calculateParams  = struct;

private/classify.m      16)         

private/classify.m      17)         calculateParams.PROJECT_BASE_PATH    = PROJECT_BASE_PATH;
private/classify.m      18)         calculateParams.PROJECT_RESULT_PATH  = PROJECT_RESULT_PATH;

private/classify.m      19)         
private/classify.m      20)         calculateParams.RANDOMIZE       = 0;

classify.m              21)         calculateParams.sessionList     = 1:3;

classify.m              22)         

classify.m              23)         calculateParams.eventList       = classStruct.eventMatrix; %[9,11,13; 10,12,14];
classify.m              24)         

private/runCoordTable.m 25)         subjectSelection = subjects;

private/classify.m      26)         calculateParams.subjectSelection = subjectSelection;

classify.m              27)         

private/runCoordTable.m 28)         decode = claculateMultiSubjectDecodePerformance(timelineParams,calculateParams,paramModel);

classify.m              29) 
classify.m              30)         display('Finished calculations.');

classify.m              31)         display('Plotting...');

classify.m              32) 

classify.m              33)         plotParams                   = struct;

private/classify.m      34)         

private/runCoordTable.m 35) %         plotParams.SVMCROSSVAL_CROSSVAL_METHOD_DEF = SVMCROSSVAL_CROSSVAL_METHOD_DEF;

private/classify.m      36)         plotParams.CROSSVAL_METHOD     = calculateParams.CROSSVAL_METHOD;
private/classify.m      37)         

classify.m              38)         plotParams.nClasses          = length(calculateParams.classList);

private/classify.m      39) 

classify.m              40)         plotParams.decodePerformance = decode.decodePerformance;
classify.m              41)         plotParams.rawTimeCourse     = decode.rawTimeCourse;

private/classify.m      42)         plotParams.SubjectID         = subjectSelection;

classify.m              43)         plotParams.smoothed          = boolToYesNoString(calculateParams.smoothed);
classify.m              44) 

classify.m              45)         assignin('base','plotParams',plotParams);
classify.m              46) %         plotDecodePerformance(params.psthStart,params.psthEnd,params.nClasses,decode.decodeTable,params.frameShiftStart,params.frameShiftEnd,decode.rawTimeCourse);

private/classify.m      47)         plotDecodePerformance(timelineParams,plotParams);

classify.m              48)             

classify.m              49)         display('all done.');

private/classify.m      50) end
private/classify.m      51)     
private/classify.m      52) 
private/classify.m      53) %% subject loop

private/runCoordTable.m 54) function decode = claculateMultiSubjectDecodePerformance(timelineParams,calculateParams,paramModel)

private/classify.m      55) 

private/classify.m      56) decode = struct;
private/classify.m      57) decode.decodePerformance = [];
private/classify.m      58) decode.rawTimeCourse     = [];
private/classify.m      59) 

private/classify.m      60) for subjectCell = calculateParams.subjectSelection

private/classify.m      61)     SubjectID = cell2mat(subjectCell);
private/classify.m      62)     namehelper = strcat('s',SubjectID); %Vars can not start with numbers.
private/classify.m      63) 
private/classify.m      64)     display('loading SPM.mat ...');
private/classify.m      65)     spm = load(fullfile(calculateParams.PROJECT_BASE_PATH,SubjectID,calculateParams.PROJECT_RESULT_PATH));
private/classify.m      66)     display('... done.');
private/classify.m      67) 

private/classify.m      68)     % per subject params..
private/classify.m      69)     subjectParams = struct;
private/classify.m      70)     
private/classify.m      71)     subjectParams.des             = spm.SPM;

private/runCoordTable.m 72) 
private/runCoordTable.m 73) 

private/runCoordTable.m 74)     subjectParams.voxelList       = mapVoxelList(voxelList,SubjectID);

private/runCoordTable.m 75) 

private/classify.m      76)     
private/classify.m      77)     subjectParams.SubjectID       = SubjectID;
private/classify.m      78)     subjectParams.namehelper      = namehelper;
private/classify.m      79) %     assignin('base','calculateParams',calculateParams);

private/classify.m      80) 
private/classify.m      81)     display(sprintf('calculating cross-validation performance time-shift for Subject %s. Please Wait. ...',SubjectID));
private/classify.m      82)     display('switching off all warnings');
private/classify.m      83)     warning_state               = warning('off','all');
private/classify.m      84)     display('calculating ...');

private/classify.m      85)     decode.(namehelper)         = calculateDecodePerformance(timelineParams,calculateParams,subjectParams);

private/classify.m      86) 
private/classify.m      87)     display('... done');
private/classify.m      88)     display('restoring warnings');
private/classify.m      89)     warning(warning_state);
private/classify.m      90) 
private/classify.m      91)     decode.decodePerformance    = [decode.decodePerformance decode.(namehelper).decodePerformance];
private/classify.m      92)     decode.rawTimeCourse        = [decode.rawTimeCourse decode.(namehelper).rawTimeCourse];
private/classify.m      93) 
private/classify.m      94)     assignin('base','decode',decode);
private/classify.m      95) end
private/classify.m      96)