SVMCrossVal.git

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

calculateDecodePerformance.m          2) outputStruct = struct;

private/calculateDecodePerformance.m  3) RANDOMIZE_DATAPOINTS = 1;

calculateDecodePerformance.m          4) 

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

private/calculateDecodePerformance.m  6) 

private/calculateDecodePerformance.m  7) timeLineStart   = timeline.frameShiftStart;
private/calculateDecodePerformance.m  8) timeLineEnd     = timeline.frameShiftEnd;

private/calculateDecodePerformance.m  9) 

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

private/calculateDecodePerformance.m 18) 

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

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

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

private/calculateDecodePerformance.m 37)     decodePerformance = [decodePerformance; svm_single_crossval(svmlabel,svmdata,svmopts)];
private/calculateDecodePerformance.m 38)     
private/calculateDecodePerformance.m 39)     

private/calculateDecodePerformance.m 40) end
private/calculateDecodePerformance.m 41) 
private/calculateDecodePerformance.m 42) outputStruct.decodePerformance  = decodePerformance;
private/calculateDecodePerformance.m 43) outputStruct.svmdata            = svmdata;
private/calculateDecodePerformance.m 44) outputStruct.svmlabel           = svmlabel;

private/calculateDecodePerformance.m 45) outputStruct.rawTimeCourse      = subjectStruct.pst;
private/calculateDecodePerformance.m 46) % outputStruct.minPerformance     = minPerformance;
private/calculateDecodePerformance.m 47) % outputStruct.maxPerformance     = maxPerformance;

calculateDecodePerformance.m         48) end
calculateDecodePerformance.m         49)