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;

calculateDecodePerformance.m          4) 

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

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

private/calculateDecodePerformance.m  9) 

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

private/calculateDecodePerformance.m 12) 

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

private/calculateDecodePerformance.m 21) 

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

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

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

private/calculateDecodePerformance.m 40)     decodePerformance = [decodePerformance; svm_single_crossval(svmlabel,svmdata,svmopts)];
private/calculateDecodePerformance.m 41)     
private/calculateDecodePerformance.m 42)     

private/calculateDecodePerformance.m 43) end
private/calculateDecodePerformance.m 44) 
private/calculateDecodePerformance.m 45) outputStruct.decodePerformance  = decodePerformance;
private/calculateDecodePerformance.m 46) outputStruct.svmdata            = svmdata;
private/calculateDecodePerformance.m 47) outputStruct.svmlabel           = svmlabel;

private/calculateDecodePerformance.m 48) outputStruct.rawTimeCourse      = subjectStruct.pst;
private/calculateDecodePerformance.m 49) % outputStruct.minPerformance     = minPerformance;
private/calculateDecodePerformance.m 50) % outputStruct.maxPerformance     = maxPerformance;

calculateDecodePerformance.m         51) end
calculateDecodePerformance.m         52)