runSvdprojLda

%% Making data

0001 %%% Making data
0002 clear trials classes parameterSets;
0003 
0004 load_eegdata;
0005 
0006 subject = 1
0007 tasks = [1 2 3 4 5]
0008 %tasks = [2 3]
0009 trials = 1:5
0010 
0011 trials = squeeze(eegdata(subject,tasks,trials));
0012 % Remove EOG channel
0013 [nr,nc] = size(trials);
0014 for r = 1:nr
0015   for c = 1:nc
0016     trials{r,c} = trials{r,c}(1:6,:);
0017   end
0018 end
0019 
0020 classes = repmat((1:length(tasks))',1,length(trials));
0021 
0022 %%% Setting up matrix of parameter values
0023 lags = 0:4;
0024 nchannels = 6;
0025 winShift = 32;
0026 k = 0;
0027 for nl = lags
0028   for nModes = 1:nchannels*(nl+1)
0029     for firstMode = 1:(nchannels*(nl+1))-nModes+1
0030       k = k + 1;
0031       parameterSets(k).nlags = nl;
0032       parameterSets(k).firstMode = firstMode;
0033       parameterSets(k).nModes = nModes;
0034     end
0035   end
0036 end
0037 
0038 %%% Run all experiments
0039 valreps = 5;
0040 CV = general(trials,classes,parameterSets,0.2,valreps,...
0041     @makeSVDProj,@LDA,@selectSamplesRows);
0042   
0043 %%% Display results
0044 
0045 summarizeTestResults(CV);