function ppsth = runCoordTable(baseDir,subjects,timeline,coordinates,volumeFileList,classDef) % paramModel = varargin{1}; disp('run coord table') disp('does not work yet.') return PROJECT_BASE_PATH = baseDir; PROJECT_RESULT_PATH = 'results\SPM.mat'; timelineParams = timeline; % common params calculateParams = struct; calculateParams.PROJECT_BASE_PATH = PROJECT_BASE_PATH; calculateParams.PROJECT_RESULT_PATH = PROJECT_RESULT_PATH; calculateParams.RANDOMIZE = 0; calculateParams.sessionList = 1:3; calculateParams.eventList = classStruct.eventMatrix; %[9,11,13; 10,12,14]; subjectSelection = subjects; calculateParams.subjectSelection = subjectSelection; decode = claculateMultiSubjectDecodePerformance(timelineParams,calculateParams,paramModel); display('Finished calculations.'); display('Plotting...'); plotParams = struct; % plotParams.SVMCROSSVAL_CROSSVAL_METHOD_DEF = SVMCROSSVAL_CROSSVAL_METHOD_DEF; plotParams.CROSSVAL_METHOD = calculateParams.CROSSVAL_METHOD; plotParams.nClasses = length(calculateParams.classList); plotParams.decodePerformance = decode.decodePerformance; plotParams.rawTimeCourse = decode.rawTimeCourse; plotParams.SubjectID = subjectSelection; plotParams.smoothed = boolToYesNoString(calculateParams.smoothed); assignin('base','plotParams',plotParams); % plotDecodePerformance(params.psthStart,params.psthEnd,params.nClasses,decode.decodeTable,params.frameShiftStart,params.frameShiftEnd,decode.rawTimeCourse); plotDecodePerformance(timelineParams,plotParams); display('all done.'); end %% subject loop function decode = claculateMultiSubjectDecodePerformance(timelineParams,calculateParams,paramModel) decode = struct; decode.decodePerformance = []; decode.rawTimeCourse = []; for subjectCell = calculateParams.subjectSelection SubjectID = cell2mat(subjectCell); namehelper = strcat('s',SubjectID); %Vars can not start with numbers. display('loading SPM.mat ...'); spm = load(fullfile(calculateParams.PROJECT_BASE_PATH,SubjectID,calculateParams.PROJECT_RESULT_PATH)); display('... done.'); % per subject params.. subjectParams = struct; subjectParams.des = spm.SPM; subjectParams.voxelList = mapVoxelList(voxelList,SubjectID); subjectParams.SubjectID = SubjectID; subjectParams.namehelper = namehelper; % assignin('base','calculateParams',calculateParams); display(sprintf('calculating cross-validation performance time-shift for Subject %s. Please Wait. ...',SubjectID)); display('switching off all warnings'); warning_state = warning('off','all'); display('calculating ...'); decode.(namehelper) = calculateDecodePerformance(timelineParams,calculateParams,subjectParams); display('... done'); display('restoring warnings'); warning(warning_state); decode.decodePerformance = [decode.decodePerformance decode.(namehelper).decodePerformance]; decode.rawTimeCourse = [decode.rawTimeCourse decode.(namehelper).rawTimeCourse]; assignin('base','decode',decode); end end