fitKernel/getChoice.m at main · SysNeuroHub/fitKernel · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40

 function [outcome, time] = getChoice(d)
 %[outcome, time] = getChoice(d)
 % returns true for trials ending in SUCCESS, false for trials ending in FAIL
 % time: time when response is registered as SUCCESS or FAIL from the start of a trial

 %cf. cuesaccade.get.cOnset (only return SUCCESS trials)

      [t,trial,~,state] = d.meta.choice.state;

      stateNEW = state;
      stateNEW = strjoin(stateNEW);
      fixedStates = strrep(stateNEW, 'FIXATING FREEVIEWING FIXATING SUCCESS', 'FIXATING FREEVIEWING FIXATING FAIL');
      fixedStates = strrep(fixedStates, '  ', ' 0 ');
      fixedStates = strsplit(fixedStates);

      % we want the start of the FIXATING state, but only for SUCCESSful trials...
      success = strcmp(fixedStates,'SUCCESS');
      fixating = strcmp(fixedStates,'FIXATING');

      outcome = nan([d.numTrials,1]);
      time = nan([d.numTrials,1]);
      %ix_s = strcmpi(state,'SUCCESS');
      ix_s = ismember(trial,trial(success)) & fixating;

      outcome(trial(ix_s)) = true;
      time(trial(ix_s)) = t(ix_s);

      %       ix_f = strcmpi(fixedStates,'FAIL') & fixating;
      ix_f = strcmpi(fixedStates,'FAIL');
      outcome(trial(ix_f)) = false;
      time(trial(ix_f)) = t(ix_f);

      firstFrames = d.meta.cic.firstFrame.time(1)';
      time = time - firstFrames;

      tOnset = d.tOnset;
      time(isnan(tOnset)) = nan;
      outcome(isnan(tOnset)) = nan;
 end