SVMCrossVal.git

1) %% subject loop
2) function decode = xsvm_subject_loop(header,subjectdata,svmopts)
3) 

4) addpath 'libsvm-mat-2.88-1';
5) 

6) nSubjects = numel(subjectdata);
7) 

8) RANDOMIZE_DATAPOINTS = 0;
9) 

10) decode = struct;
11) decode.decodePerformance = [];
12) decode.rawTimeCourse     = [];
13) 

14) disp(sprintf('computinig additional datastructs for %u subjects',nSubjects));

15) 
16) timeline = header.timeline;

17) timeline.frameShiftStart = header.frameShift.frameShiftStart;
18) timeline.frameShiftEnd   = header.frameShift.frameShiftEnd;
19) timeline.decodeDuration  = header.frameShift.decodeDuration;

20) 

21) % TimePointMatrix

22) for subjectDataID = 1:nSubjects
23)     currentSubject = subjectdata{subjectDataID};
24)     timePointArgs.pst           = currentSubject.pst;
25)     timePointArgs.labelMap      = LabelMap(header.classDef.labelCells,header.classDef.conditionCells);
26)     timePointArgs.eventList     = header.classDef.eventMatrix;
27) 
28)     timePointMatrix{subjectDataID} = buildTimePointMatrix(timeline,timePointArgs);

29)     
30)     decode.rawTimeCourse = [decode.rawTimeCourse currentSubject.pst];

31) end
32) 

33) % timeframe x-subject validation

34) timeLineStart   = timeline.frameShiftStart;
35) timeLineEnd     = timeline.frameShiftEnd;
36) 

37) display(sprintf('%u -fold cross validation for %u timeslices.\n',nSubjects,size(1:timeLineEnd-timeLineStart+1,2)));

38) % disp(sprintf('Press ANY-Key to continue.\n Use Retrun if your Keyboard lacks the ANY-Key.'));
39) % pause
40) 

41) for timeIndex = 1:timeLineEnd-timeLineStart+1

42)     cross_value = [];

43)     for validationSubjectID = 1:nSubjects

44)         svm_train_label = [];
45)         svm_train_data  = [];
46)         svm_validation_label = [];
47)         svm_validation_data  = [];

48)         for subjectDataID = 1:nSubjects
49)             svmstruct = calculateSVMTables(timePointMatrix{subjectDataID},timeIndex);
50)             if subjectDataID == validationSubjectID
51)                 svm_validation_label = svmstruct.svmlabel;
52)                 svm_validation_data  = svmstruct.svmdata;
53)             else
54)                 svm_train_label = [svm_train_label; svmstruct.svmlabel];
55)                 svm_train_data  = [svm_train_data;  svmstruct.svmdata];
56)             end
57)         end
58)         

59)         if RANDOMIZE_DATAPOINTS
60)             rndindex  = randperm(length(svm_train_label));
61)             svm_train_data   = svm_train_data(rndindex,:);
62)             svm_train_label  = svm_train_label(rndindex);
63)         end
64) 

65)         svmmodel = svmtrain(svm_train_label,svm_train_data,svmopts);
66)         
67)         [plabel accuracy dvalue] = svmpredict(svm_validation_label,svm_validation_data,svmmodel,'');

68)         cross_value = [cross_value accuracy(1)];

69)         

70)     end

71)     decode.decodePerformance = [decode.decodePerformance; cross_value];

72)     

73) end

74) 
75) disp('decode done');