private/calculateDecodePerformance.m (d0eb71f) - SVMCrossVal.git

calculateDecodePerformance.m

function outputStruct = calculateDecodePerformance(header,subjectStruct,svmopts)
outputStruct = struct;
RANDOMIZE_DATAPOINTS = 1;

timeline = header.timeline;

timeLineStart   = timeline.frameShiftStart;
timeLineEnd     = timeline.frameShiftEnd;

% for iVoxel = 1:nVoxel
%     rawdata = [];
%     for iImage = 1:length(extr);
%         tmp = extr(iImage);
%         rawdata = [rawdata tmp.dat(iVoxel)];
%     end
%     pst{iVoxel} = calculatePST(timeline,calculatePstOpts,rawdata);
% end

timePointArgs.pst           = subjectStruct.pst;
timePointArgs.labelMap      = LabelMap(header.classDef.labelCells,header.classDef.conditionCells);
timePointArgs.eventList     = header.classDef.eventMatrix;

timePointMatrix = buildTimePointMatrix(timeline,timePointArgs);

decodePerformance = [];
for index = 1:timeLineEnd-timeLineStart+1

svmdata      = timePointMatrix{index}(:,2:size(timePointMatrix{index},2));
    svmlabel     = timePointMatrix{index}(:,1);

if RANDOMIZE_DATAPOINTS
        rndindex  = randperm(length(svmlabel));
        svmdata   = svmdata(rndindex,:);
        svmlabel  = svmlabel(rndindex);
    end

decodePerformance = [decodePerformance; svm_single_crossval(svmlabel,svmdata,svmopts)];
    
    
end

outputStruct.decodePerformance  = decodePerformance;
outputStruct.svmdata            = svmdata;
outputStruct.svmlabel           = svmlabel;
outputStruct.rawTimeCourse      = subjectStruct.pst;
% outputStruct.minPerformance     = minPerformance;
% outputStruct.maxPerformance     = maxPerformance;
end

SVMCrossVal.git

single subject -> svm classification -> plotDecodePerformance working!