fitKernel/getTgtBhvLatency.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
function latency_bhv = getTgtBhvLatency(onsets_cat, dd, validEvents, option)
%latency_bhv = getTgtBhvLatency(onsets_cat, dd, validEvents. option)
%
% option = 1: use "srt" stored in dd
% option = 0: cOnset - tOnset

if nargin < 4
    option = 0;
end
if nargin < 3 || isempty(validEvents)
    validEvents = 1:numel(onsets_cat.tOnset);
end

switch option
    case 0
        latency_bhv = onsets_cat.cOnset(validEvents) - onsets_cat.tOnset(validEvents);
    case 1
        latency_bhv = dd.srt(validEvents);
end

% latency_bhv=getTgtBhvLatency(onsets_cat, validEvents);
% srt = dd.srt(validEvents);
% plot(latency_bhv, srt,'.'); hold on;
%
% discrepancy = find(~isnan(latency_bhv)&isnan(srt));
% plot(latency_bhv(discrepancy), .5,'r*');
%
% discrepancy2 = find(isnan(latency_bhv)&~isnan(srt));
% plot(.5,srt(discrepancy2),'r*');
%
% xlabel('cOnset - tOnset');
% ylabel('srt = sacccade.onset - tOnset');
%
% squareplot;
% axis equal padded