SVMCrossVal.git

private/plotDecodePerformance.m   1) function plotDecodePerformance(type,decode,subjectData)

plotDecodePerformance.m           2) 

private/plotDecodePerformance.m   3) header = decode.header;

private/plotDecodePerformance.m   4) timeline = header.timeline;

private/plotDecodePerformance.m   5) frameshift = header.frameShift;

private/plotDecodePerformance.m   6) 

private/plotDecodePerformance.m   7) psthStart         = timeline.psthStart;
private/plotDecodePerformance.m   8) psthEnd           = timeline.psthEnd;

private/plotDecodePerformance.m   9) frameStart        = frameshift.frameShiftStart;
private/plotDecodePerformance.m  10) frameEnd          = frameshift.frameShiftEnd;

plotDecodePerformance.m          11) 

private/plotDecodePerformance.m  12) nClasses          = numel(header.classDef.labelCells);
private/plotDecodePerformance.m  13) decodePerformance = decode.decodePerformance;
private/plotDecodePerformance.m  14) psth              = decode.rawTimeCourse;

private/plotDecodePerformance.m  15) ert               = psth; % Rename to EventRelatedTimeline

private/plotDecodePerformance.m  16) SubjectID         = subjectData;

private/plotDecodePerformance.m  17) 

private/plotDecodePerformance.m  18) nSubjects         = size(SubjectID,2);

plotDecodePerformance.m          19) 

private/plotDecodePerformance.m  20) nTrials           = getNTrials(psth);

plotDecodePerformance.m          21) 
plotDecodePerformance.m          22)     f = figure;

private/plotDecodePerformance.m  23)   

plotDecodePerformance.m          24)     hold on;
plotDecodePerformance.m          25)     chanceLevel = 100/nClasses;

private/plotDecodePerformance.m  26) 

plotDecodePerformance.m          27)     axis([psthStart psthEnd 0 100])

plotDecodePerformance.m          28)     xlabel('time [sec]');
plotDecodePerformance.m          29)     ylabel('decode performance [%]');

private/plotDecodePerformance.m  30)     
private/plotDecodePerformance.m  31)     

private/plotDecodePerformance.m  32) plottime= tic;

private/plotDecodePerformance.m  33)     switch type

private/plotDecodePerformance.m  34)         case 'psth'

private/plotDecodePerformance.m  35)             plotPSTH(ert,psthStart,psthEnd);  
private/plotDecodePerformance.m  36)             
private/plotDecodePerformance.m  37)         case 'trajectory'
private/plotDecodePerformance.m  38)             plotTrajectory(ert,psthStart,psthEnd);
private/plotDecodePerformance.m  39) 
private/plotDecodePerformance.m  40)         case 'trajectory3'
private/plotDecodePerformance.m  41)             plotTrajectory3(ert);

private/plotDecodePerformance.m  42) 

private/plotDecodePerformance.m  43)         case 'simple'
private/plotDecodePerformance.m  44)             plotDecodePerformanceWithSE(frameStart,frameEnd,decodePerformance)

private/plotDecodePerformance.m  45) 
private/plotDecodePerformance.m  46)             plot([psthStart psthEnd],[chanceLevel chanceLevel],'k:');
private/plotDecodePerformance.m  47)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.05,chanceLevel/100,'--');
private/plotDecodePerformance.m  48)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.01,chanceLevel/100,'-.');
private/plotDecodePerformance.m  49)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.001,chanceLevel/100,':');
private/plotDecodePerformance.m  50) 

private/plotDecodePerformance.m  51)         case 'class performance'
private/plotDecodePerformance.m  52)             plotClassPerformance(frameStart,frameEnd,decodePerformance,nClasses)

private/plotDecodePerformance.m  53) 
private/plotDecodePerformance.m  54)             plot([psthStart psthEnd],[chanceLevel chanceLevel],'k:');
private/plotDecodePerformance.m  55)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.05,chanceLevel/100,'--');
private/plotDecodePerformance.m  56)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.01,chanceLevel/100,'-.');
private/plotDecodePerformance.m  57)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.001,chanceLevel/100,':');

private/plotDecodePerformance.m  58)         case 'x-subject-val'
private/plotDecodePerformance.m  59)             for c = 1:nSubjects
private/plotDecodePerformance.m  60)                 plot(frameStart:frameEnd, decodePerformance(:,c) ,[colorChooser(mod(c,nSubjects)+3) '-']);
private/plotDecodePerformance.m  61)             end

private/plotDecodePerformance.m  62)             plotDecodePerformanceWithSE(frameStart,frameEnd,decodePerformance)

private/plotDecodePerformance.m  63)             
private/plotDecodePerformance.m  64)             plot([psthStart psthEnd],[chanceLevel chanceLevel],'k:');
private/plotDecodePerformance.m  65)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.05,chanceLevel/100,'--');
private/plotDecodePerformance.m  66)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.01,chanceLevel/100,'-.');
private/plotDecodePerformance.m  67)             plotBinConfidenceIntervall(psthStart,psthEnd,nTrials,0.001,chanceLevel/100,':');

private/plotDecodePerformance.m  68) 

private/plotDecodePerformance.m  69)     end
private/plotDecodePerformance.m  70) toc(plottime);    

plotDecodePerformance.m          71)     hold off;

plotDecodePerformance.m          72)     

private/plotDecodePerformance.m  73) end
private/plotDecodePerformance.m  74) 

private/plotDecodePerformance.m  75) function plotBinConfidenceIntervall(pStart,pEnd,nTrials,alpha,limit,color)
private/plotDecodePerformance.m  76)     [pLevel Z] = rev_ki_bin(nTrials,alpha,limit);
private/plotDecodePerformance.m  77)     [lower upper proz] = ki_bin(Z,nTrials,alpha);
private/plotDecodePerformance.m  78) 

private/plotDecodePerformance.m  79)     plot([pStart pEnd],[pLevel*100 pLevel*100],[color 'k']);
private/plotDecodePerformance.m  80) %     plot([pStart pEnd],[lower*100 lower*100],[color 'k']);
private/plotDecodePerformance.m  81) %     plot([pStart pEnd],[upper*100 upper*100],[color 'k']);

private/plotDecodePerformance.m  82) end
private/plotDecodePerformance.m  83) 
private/plotDecodePerformance.m  84) function n = getNTrials(psth)

private/plotDecodePerformance.m  85) nOverallConditions = size(psth,2);

private/plotDecodePerformance.m  86) 
private/plotDecodePerformance.m  87) n = 0;

private/plotDecodePerformance.m  88) for ns = 1:nOverallConditions

private/plotDecodePerformance.m  89)     nClasses = size(psth{ns},2);
private/plotDecodePerformance.m  90)     for nc = 1:nClasses
private/plotDecodePerformance.m  91)         n = n + size(psth{ns}{nc},1);
private/plotDecodePerformance.m  92)     end
private/plotDecodePerformance.m  93) end

private/plotDecodePerformance.m  94) n=n/(nOverallConditions/nClasses);
private/plotDecodePerformance.m  95) 

private/plotDecodePerformance.m  96) end
private/plotDecodePerformance.m  97) 

private/plotDecodePerformance.m  98) function plotClassPerformance(frameStart,frameEnd,decodePerformance,nClasses)
private/plotDecodePerformance.m  99) for c = 1:size(decodePerformance,2)
private/plotDecodePerformance.m 100)     plot(frameStart:frameEnd, decodePerformance(:,c) ,[colorChooser(mod(c,nClasses)+3) '-']);
private/plotDecodePerformance.m 101) end
private/plotDecodePerformance.m 102) plot(frameStart:frameEnd, mean(decodePerformance,2) ,'b','LineWidth',2);
private/plotDecodePerformance.m 103) end

plotDecodePerformance.m         104) 

private/plotDecodePerformance.m 105) function plotDecodePerformanceWithSE(frameStart,frameEnd,decodePerformance)
private/plotDecodePerformance.m 106) plot(frameStart:frameEnd, mean(decodePerformance,2) ,'b','LineWidth',2);
private/plotDecodePerformance.m 107) 
private/plotDecodePerformance.m 108) se = myStdErr(decodePerformance,2);
private/plotDecodePerformance.m 109) plot(frameStart:frameEnd, mean(decodePerformance,2)+se ,'b:');
private/plotDecodePerformance.m 110) plot(frameStart:frameEnd, mean(decodePerformance,2)-se ,'b:');
private/plotDecodePerformance.m 111) end

plotDecodePerformance.m         112) 

private/plotDecodePerformance.m 113) 
private/plotDecodePerformance.m 114) function plotPSTH(psth,psthStart,psthEnd)
private/plotDecodePerformance.m 115) 
private/plotDecodePerformance.m 116) PSTH_AXIS_MIN = -2;
private/plotDecodePerformance.m 117) PSTH_AXIS_MAX = 2;
private/plotDecodePerformance.m 118) 
private/plotDecodePerformance.m 119)     hold on;
private/plotDecodePerformance.m 120)     if (size(psth) > 0)
private/plotDecodePerformance.m 121)       for voxel = 1:size(psth,2)
private/plotDecodePerformance.m 122)           for label = 1:size(psth{voxel},2)
private/plotDecodePerformance.m 123)               psthData = [];
private/plotDecodePerformance.m 124)               for timepoint = 1:size(psth{voxel}{label},2)
private/plotDecodePerformance.m 125)                   psthData = nanmean(psth{voxel}{label});
private/plotDecodePerformance.m 126)               end

private/plotDecodePerformance.m 127)               PSTH_AXIS_MAX = max(PSTH_AXIS_MAX,nanmax(psthData));

private/plotDecodePerformance.m 128)               PSTH_AXIS_MIN = min(PSTH_AXIS_MIN,nanmin(psthData));

private/plotDecodePerformance.m 129)               plot(psthStart:psthEnd,psthData,[colorChooser(voxel), lineStyleChooser(label)]);

private/plotDecodePerformance.m 130) %               plot(psthStart:psthEnd,psthData,[lineStyleChooser(voxel), colorChooser(label)]);

private/plotDecodePerformance.m 131)           end
private/plotDecodePerformance.m 132)       end
private/plotDecodePerformance.m 133)     end
private/plotDecodePerformance.m 134)     axis([psthStart psthEnd PSTH_AXIS_MIN PSTH_AXIS_MAX])
private/plotDecodePerformance.m 135)     xlabel('time [sec]');
private/plotDecodePerformance.m 136)     ylabel('fMRI-signal change [%]');
private/plotDecodePerformance.m 137)     hold off

plotDecodePerformance.m         138) end
plotDecodePerformance.m         139) 

plotDecodePerformance.m         140) 

private/plotDecodePerformance.m 141) function plotTrajectory(psth,psthStart,psthEnd)
private/plotDecodePerformance.m 142) % ERTA_AXIS_MAX = 2;
private/plotDecodePerformance.m 143) % ERTA_AXIS_MIN = 2;
private/plotDecodePerformance.m 144) % ERTB_AXIS_MAX = 2;
private/plotDecodePerformance.m 145) % ERTB_AXIS_MIN = 2;
private/plotDecodePerformance.m 146) 
private/plotDecodePerformance.m 147)     if (mod(size(psth,2),2)==0) %% TODO insufficient test!!!
private/plotDecodePerformance.m 148)     hold on;
private/plotDecodePerformance.m 149) %     grid on;
private/plotDecodePerformance.m 150)     for voxel = 1:2:size(psth,2)
private/plotDecodePerformance.m 151)         for label = 1:size(psth{voxel},2)
private/plotDecodePerformance.m 152)             ertDataA = [];
private/plotDecodePerformance.m 153)             for timepoint = 1:size(psth{voxel}{label},2)
private/plotDecodePerformance.m 154)                 ertDataA = nanmean(psth{voxel}{label}); % mean ueber die trials
private/plotDecodePerformance.m 155)             end
private/plotDecodePerformance.m 156) %               ERTA_AXIS_MAX = max(ERTA_AXIS_MAX,nanmax(ertDataA));
private/plotDecodePerformance.m 157) %               ERTA_AXIS_MIN = min(ERTA_AXIS_MIN,nanmin(ertDataA));
private/plotDecodePerformance.m 158) 
private/plotDecodePerformance.m 159)             voxel2 = voxel+1;
private/plotDecodePerformance.m 160)             ertDataB = [];
private/plotDecodePerformance.m 161)             for timepoint = 1:size(psth{voxel2}{label},2)
private/plotDecodePerformance.m 162)                 ertDataB = nanmean(psth{voxel2}{label}); % mean ueber die trials
private/plotDecodePerformance.m 163)             end
private/plotDecodePerformance.m 164) %               ERTB_AXIS_MAX = max(ERTB_AXIS_MAX,nanmax(ertDataB));
private/plotDecodePerformance.m 165) %               ERTB_AXIS_MIN = min(ERTB_AXIS_MIN,nanmin(ertDataB));
private/plotDecodePerformance.m 166)               
private/plotDecodePerformance.m 167)             plot3(psthStart:psthEnd,ertDataA,ertDataB,['.', colorChooser(label),'-']);
private/plotDecodePerformance.m 168)         end
private/plotDecodePerformance.m 169)     end
private/plotDecodePerformance.m 170) %     axis([psthStart psthEnd ERTA_AXIS_MIN ERTA_AXIS_MAX ERTB_AXIS_MIN ERTB_AXIS_MAX])
private/plotDecodePerformance.m 171)     axis tight;
private/plotDecodePerformance.m 172)     view(90,0);
private/plotDecodePerformance.m 173)     xlabel('time [sec]');
private/plotDecodePerformance.m 174)     ylabel('fMRI-signal change [%] A');
private/plotDecodePerformance.m 175)     zlabel('fMRI-signal change [%] B');
private/plotDecodePerformance.m 176)     
private/plotDecodePerformance.m 177)     hold off
private/plotDecodePerformance.m 178)     
private/plotDecodePerformance.m 179)     
private/plotDecodePerformance.m 180)     else 
private/plotDecodePerformance.m 181)         display('unsuitable # of voxel/roi!');
private/plotDecodePerformance.m 182)     end
private/plotDecodePerformance.m 183) 
private/plotDecodePerformance.m 184) end
private/plotDecodePerformance.m 185) 
private/plotDecodePerformance.m 186) function plotTrajectory3(ert)
private/plotDecodePerformance.m 187)     
private/plotDecodePerformance.m 188)     if (mod(size(ert,2),3)~=0)
private/plotDecodePerformance.m 189)         display('unsuitable # of voxel/roi!');
private/plotDecodePerformance.m 190)         return
private/plotDecodePerformance.m 191)     end
private/plotDecodePerformance.m 192)         
private/plotDecodePerformance.m 193)     nSubjects = size(ert,2)/3;
private/plotDecodePerformance.m 194)     
private/plotDecodePerformance.m 195) %     nLabel = size(ert{1},2);  % this is sane, because the program produces only such things.
private/plotDecodePerformance.m 196)     subplot(2,1,1);
private/plotDecodePerformance.m 197)     hold on;
private/plotDecodePerformance.m 198) %     for label = 1:nLabel 
private/plotDecodePerformance.m 199)         ertDataA = ertDataMatrix(ert,1,1,nSubjects);
private/plotDecodePerformance.m 200)         ertDataB = ertDataMatrix(ert,1,2,nSubjects);
private/plotDecodePerformance.m 201)         ertDataC = ertDataMatrix(ert,1,3,nSubjects);
private/plotDecodePerformance.m 202) 
private/plotDecodePerformance.m 203)         plot3(ertDataA,ertDataB,ertDataC,['.', colorChooser(1),'-']);
private/plotDecodePerformance.m 204) 
private/plotDecodePerformance.m 205)         ertDataA2 = ertDataMatrix(ert,2,1,nSubjects);
private/plotDecodePerformance.m 206)         ertDataB2 = ertDataMatrix(ert,2,2,nSubjects);
private/plotDecodePerformance.m 207)         ertDataC2 = ertDataMatrix(ert,2,3,nSubjects);
private/plotDecodePerformance.m 208) 
private/plotDecodePerformance.m 209)         plot3(ertDataA2,ertDataB2,ertDataC2,['.', colorChooser(2),'-']);
private/plotDecodePerformance.m 210) 
private/plotDecodePerformance.m 211)         %     end
private/plotDecodePerformance.m 212)         
private/plotDecodePerformance.m 213)         for a = 1:length(ertDataA)
private/plotDecodePerformance.m 214)             line([ertDataA(a) ertDataA2(a)],...
private/plotDecodePerformance.m 215)                 [ertDataB(a) ertDataB2(a)],...
private/plotDecodePerformance.m 216)                 [ertDataC(a) ertDataC2(a)],...
private/plotDecodePerformance.m 217)                 'color',[.8/a 0 1-.8/a]);
private/plotDecodePerformance.m 218)             
private/plotDecodePerformance.m 219)             d(a) = sqrt( (ertDataA(a)-ertDataA2(a))^2 +...
private/plotDecodePerformance.m 220)             (ertDataB(a)-ertDataB2(a))^2 + ...
private/plotDecodePerformance.m 221)             (ertDataC(a)-ertDataC2(a))^2);
private/plotDecodePerformance.m 222)         end
private/plotDecodePerformance.m 223) 
private/plotDecodePerformance.m 224)     axis tight;
private/plotDecodePerformance.m 225) %     view(90,0);
private/plotDecodePerformance.m 226)     xlabel('fMRI-signal change [%] A');
private/plotDecodePerformance.m 227)     ylabel('fMRI-signal change [%] B');
private/plotDecodePerformance.m 228)     zlabel('fMRI-signal change [%] C');
private/plotDecodePerformance.m 229)     
private/plotDecodePerformance.m 230)     hold off
private/plotDecodePerformance.m 231) 
private/plotDecodePerformance.m 232)     subplot(2,1,2);
private/plotDecodePerformance.m 233)         
private/plotDecodePerformance.m 234)     plot(d);
private/plotDecodePerformance.m 235)     
private/plotDecodePerformance.m 236) end
private/plotDecodePerformance.m 237) 
private/plotDecodePerformance.m 238) function ertDataMatrix = ertDataMatrix(ert,label,roiID,nSubjects)
private/plotDecodePerformance.m 239)     ertDataMatrix = [];
private/plotDecodePerformance.m 240)     for subject = 1:nSubjects
private/plotDecodePerformance.m 241)         voxel = roiID*subject;
private/plotDecodePerformance.m 242)         for timepoint = 1:size(ert{voxel}{label},2)
private/plotDecodePerformance.m 243)             ertDataMatrix = [ertDataMatrix ; nanmean(ert{voxel}{label})]; % mean ueber die trials
private/plotDecodePerformance.m 244)         end
private/plotDecodePerformance.m 245)     end
private/plotDecodePerformance.m 246)     ertDataMatrix = nanmean(ertDataMatrix,1);
private/plotDecodePerformance.m 247) end
private/plotDecodePerformance.m 248) 
private/plotDecodePerformance.m 249) 

plotDecodePerformance.m         250) function color = colorChooser(n)
plotDecodePerformance.m         251)     switch (mod(n,8))
plotDecodePerformance.m         252)     case 0
plotDecodePerformance.m         253)         color = 'r';

plotDecodePerformance.m         254)     case 1
plotDecodePerformance.m         255)         color = 'g';

plotDecodePerformance.m         256)     case 2
plotDecodePerformance.m         257)         color = 'b';
plotDecodePerformance.m         258)     case 3

plotDecodePerformance.m         259)         color = 'c';
plotDecodePerformance.m         260)     case 4
plotDecodePerformance.m         261)         color = 'm';
plotDecodePerformance.m         262)     case 5
plotDecodePerformance.m         263)         color = 'y';

plotDecodePerformance.m         264)     otherwise
plotDecodePerformance.m         265)         color = 'k';
plotDecodePerformance.m         266)     end
plotDecodePerformance.m         267) end
plotDecodePerformance.m         268) 
plotDecodePerformance.m         269) function style = lineStyleChooser(n)
plotDecodePerformance.m         270) switch(mod(n,4))
plotDecodePerformance.m         271)     case 0
plotDecodePerformance.m         272)       style = '--';
plotDecodePerformance.m         273)     case 1
plotDecodePerformance.m         274)         style = '-';
plotDecodePerformance.m         275)     case 2 
plotDecodePerformance.m         276)         style = ':';
plotDecodePerformance.m         277)     case 3

plotDecodePerformance.m         278)         style = '-.';