SVMCrossVal.git

classify.m          1) function classify(varargin)

classify.m          2) 

classify.m          3) switch nargin
classify.m          4)     case 1
classify.m          5)         paramModel = varargin{1};

classify.m          6)         %PROJECT_BASE_PATH = 'D:\Analyze\Stimolos';
classify.m          7)         PROJECT_BASE_PATH = 'D:\Analyze\Choice\24pilot';

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

classify.m          9)     otherwise
classify.m         10)         error('spmtoolbox:SVMCrossVal:arginError','Please Specify action and parameter model');

classify.m         11) end

private/classify.m 12)        

classify.m         13)         % common params

classify.m         14)         calculateParams  = struct;

private/classify.m 15)         calculateParams.smoothed        = getChkValue(paramModel.chkSmoothed);

classify.m         16) 

classify.m         17)         calculateParams.frameShiftStart = getDouble(paramModel.txtFrameShiftStart);  % -20;
classify.m         18)         calculateParams.frameShiftEnd   = getDouble(paramModel.txtFrameShiftEnd); %15;
classify.m         19)         calculateParams.decodeDuration  = getDouble(paramModel.txtFrameShiftDur);
classify.m         20)         calculateParams.psthStart       = getDouble(paramModel.txtPSTHStart); % -25;
classify.m         21)         calculateParams.psthEnd         = getDouble(paramModel.txtPSTHEnd); % 20;
classify.m         22)         calculateParams.baselineStart   = getDouble(paramModel.txtBaselineStart); % -22;
classify.m         23)         calculateParams.baselineEnd     = getDouble(paramModel.txtBaselineEnd); % -20;
classify.m         24) 
classify.m         25)         calculateParams.svmargs         = get(paramModel.txtSVMopts,'String');

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

classify.m         27) 

classify.m         28)         classStruct = parseClassDef(paramModel);
classify.m         29)         
classify.m         30)         

classify.m         31)         calculateParams.labelMap        = LabelMap(classStruct.labelCells , classStruct.conditionCells, 'auto'); % LabelMap({'<','>','<+<','>+>','<+>','>+<'},{-2,-1,1,2,3,4}); 0 is autolabel
classify.m         32)         calculateParams.classList       = getClasses(calculateParams.labelMap);
classify.m         33)         calculateParams.eventList       = classStruct.eventMatrix; %[9,11,13; 10,12,14];
classify.m         34) %         calculateParams.eventList       = getPSTEventMatrix(calculateParams.labelMap);
classify.m         35)         

classify.m         36)         subjectSelection = getSubjectIDString(paramModel);
classify.m         37)         decode = struct;
classify.m         38)         decode.decodePerformance = [];
classify.m         39)         decode.rawTimeCourse     = [];

classify.m         40)         

classify.m         41)         for subjectCell = subjectSelection
classify.m         42)             SubjectID = cell2mat(subjectCell);

classify.m         43)             namehelper = strcat('s',SubjectID); %Vars can not start with numbers.

classify.m         44) 

classify.m         45)             display('loading SPM.mat ...');

classify.m         46)             spm = load(fullfile(PROJECT_BASE_PATH,SubjectID,PROJECT_RESULT_PATH));

classify.m         47)             display('... done.');

classify.m         48) 
classify.m         49)             %% calculate

classify.m         50)             calculateParams.(namehelper).des             = spm.SPM;
classify.m         51)             calculateParams.(namehelper).voxelList       = parseVoxelList(paramModel,SubjectID);

classify.m         52)             assignin('base','calculateParams',calculateParams);

classify.m         53) 
classify.m         54)             display(sprintf('calculating cross-validation performance time-shift for Subject %s. Please Wait. ...',SubjectID));

classify.m         55)             display('switching off all warnings');

classify.m         56)             warning_state               = warning('off','all');

classify.m         57)             display('calculating ...');

classify.m         58)             decode.(namehelper)         = calculateDecodePerformance(calculateParams,SubjectID);

private/classify.m 59)             

classify.m         60)             display('... done');

classify.m         61)             display('restoring warnings');
classify.m         62)             warning(warning_state);
classify.m         63)             

classify.m         64)             decode.decodePerformance    = [decode.decodePerformance decode.(namehelper).decodePerformance];
classify.m         65)             decode.rawTimeCourse        = [decode.rawTimeCourse decode.(namehelper).rawTimeCourse];

classify.m         66) 
classify.m         67)             assignin('base','decode',decode);
classify.m         68)         end

classify.m         69) 
classify.m         70)         display('Finished calculations.');

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

classify.m         72) 

classify.m         73)         plotParams                   = struct;
classify.m         74)         plotParams.psthStart         = calculateParams.psthStart;
classify.m         75)         plotParams.psthEnd           = calculateParams.psthEnd;
classify.m         76)         plotParams.nClasses          = length(calculateParams.classList);

classify.m         77)         

classify.m         78)         plotParams.frameShiftStart   = calculateParams.frameShiftStart;
classify.m         79)         plotParams.frameShiftEnd     = calculateParams.frameShiftEnd;
classify.m         80)         plotParams.decodePerformance = decode.decodePerformance;
classify.m         81)         plotParams.rawTimeCourse     = decode.rawTimeCourse;
classify.m         82)         

classify.m         83)         if numel(subjectSelection) == 1
classify.m         84)           plotParams.SubjectID         = SubjectID;
classify.m         85)         else
classify.m         86)           plotParams.SubjectID         = 'Multiple';
classify.m         87)         end
classify.m         88) 
classify.m         89)         plotParams.smoothed          = boolToYesNoString(calculateParams.smoothed);
classify.m         90) 

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

classify.m         94)             

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