SVMCrossVal.git

private/calculateDecodePerformance.m  1) function outputStruct = calculateDecodePerformance(header,subjectStruct,svmopts)

calculateDecodePerformance.m          2) outputStruct = struct;

private/calculateDecodePerformance.m  3) RANDOMIZE_DATAPOINTS = header.svmrnd;

private/calculateDecodePerformance.m  4) NAN_AS_ZERO = header.nantozero;

calculateDecodePerformance.m          5) 

private/calculateDecodePerformance.m  6) timeline = header.timeline;

private/calculateDecodePerformance.m  7) timeline.frameShiftStart = header.frameShift.frameShiftStart;
private/calculateDecodePerformance.m  8) timeline.frameShiftEnd   = header.frameShift.frameShiftEnd;
private/calculateDecodePerformance.m  9) timeline.decodeDuration  = header.frameShift.decodeDuration;

private/calculateDecodePerformance.m 10) 

private/calculateDecodePerformance.m 11) timeLineStart  = timeline.frameShiftStart;
private/calculateDecodePerformance.m 12) timeLineEnd    = timeline.frameShiftEnd;

private/calculateDecodePerformance.m 13) 

private/calculateDecodePerformance.m 14) % for iVoxel = 1:nVoxel
private/calculateDecodePerformance.m 15) %     rawdata = [];
private/calculateDecodePerformance.m 16) %     for iImage = 1:length(extr);
private/calculateDecodePerformance.m 17) %         tmp = extr(iImage);
private/calculateDecodePerformance.m 18) %         rawdata = [rawdata tmp.dat(iVoxel)];
private/calculateDecodePerformance.m 19) %     end
private/calculateDecodePerformance.m 20) %     pst{iVoxel} = calculatePST(timeline,calculatePstOpts,rawdata);
private/calculateDecodePerformance.m 21) % end

private/calculateDecodePerformance.m 22) 

private/calculateDecodePerformance.m 23) timePointArgs.pst           = subjectStruct.pst;

private/calculateDecodePerformance.m 24) timePointArgs.labelMap      = LabelMap(header.classDef.labelCells,header.classDef.conditionCells);
private/calculateDecodePerformance.m 25) timePointArgs.eventList     = header.classDef.eventMatrix;

private/calculateDecodePerformance.m 26) 
private/calculateDecodePerformance.m 27) timePointMatrix = buildTimePointMatrix(timeline,timePointArgs);
private/calculateDecodePerformance.m 28) 
private/calculateDecodePerformance.m 29) decodePerformance = [];
private/calculateDecodePerformance.m 30) for index = 1:timeLineEnd-timeLineStart+1
private/calculateDecodePerformance.m 31) 
private/calculateDecodePerformance.m 32)     svmdata      = timePointMatrix{index}(:,2:size(timePointMatrix{index},2));
private/calculateDecodePerformance.m 33)     svmlabel     = timePointMatrix{index}(:,1);
private/calculateDecodePerformance.m 34) 
private/calculateDecodePerformance.m 35)     if RANDOMIZE_DATAPOINTS
private/calculateDecodePerformance.m 36)         rndindex  = randperm(length(svmlabel));
private/calculateDecodePerformance.m 37)         svmdata   = svmdata(rndindex,:);
private/calculateDecodePerformance.m 38)         svmlabel  = svmlabel(rndindex);
private/calculateDecodePerformance.m 39)     end
private/calculateDecodePerformance.m 40) 

private/calculateDecodePerformance.m 41)     if NAN_AS_ZERO
private/calculateDecodePerformance.m 42)         svmdata(isnan(svmdata))=0;
private/calculateDecodePerformance.m 43)     end
private/calculateDecodePerformance.m 44)         
private/calculateDecodePerformance.m 45)     

private/calculateDecodePerformance.m 46)     decodePerformance = [decodePerformance; svm_single_crossval(svmlabel,svmdata,svmopts)];
private/calculateDecodePerformance.m 47)     
private/calculateDecodePerformance.m 48)     

private/calculateDecodePerformance.m 49) end
private/calculateDecodePerformance.m 50) 
private/calculateDecodePerformance.m 51) outputStruct.decodePerformance  = decodePerformance;
private/calculateDecodePerformance.m 52) outputStruct.svmdata            = svmdata;
private/calculateDecodePerformance.m 53) outputStruct.svmlabel           = svmlabel;

private/calculateDecodePerformance.m 54) outputStruct.rawTimeCourse      = subjectStruct.pst;
private/calculateDecodePerformance.m 55) % outputStruct.minPerformance     = minPerformance;
private/calculateDecodePerformance.m 56) % outputStruct.maxPerformance     = maxPerformance;

calculateDecodePerformance.m         57) end
calculateDecodePerformance.m         58)