fitKernel/concatenate_eye.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
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
function [eyeData_cat, onsets_cat, meta_cat] = ...
    concatenate_eye(eyeData, dd)
%[eyeData_cat] = concatenate_eye(dd)
% returns mamordata.eye concatenated across trials

%[eyeData_cat, cOnset_cat] = concatenate_eye(eyeData, dd)
% additionally returns Onset data concatenated across trials

%[eyeData_cat, onsets_cat, meta_cat] = concatenate_eye(dd)
% additionally reutns saccade and blink times computed in eyelink

%TODO:
%make the function faster with functions like cellfun
%return cuesaccade structure concatenated across trials?


%eyeData = dd.eye;
nTrials = length(eyeData);

if nargout>1
    fdnames = fieldnames(dd)';
    fdnames = fdnames(contains(fdnames, 'Onset') & ~contains(fdnames, 'cOnset')...
        & ~contains(fdnames, 'fOnset'));
    c = cell(length(fdnames),1);
    onsets_cat = cell2struct(c,fdnames);
end
if nargout>2
    edfnames = {'STARTSACC','ENDSACC','STARTBLINK','ENDBLINK'};
    c = cell(length(edfnames),1);
    meta_cat = cell2struct(c,edfnames);
end

%% concatenate across trials as in my Cell Rep paper 2018
%does parea=0 mean blink?
%is there temporal gap between trials?
parea_cat = [];
pwdth_cat = [];
phght_cat = [];
x_cat = [];
y_cat = [];
t_cat = [];
%dd_cat = marmodata.cuesaccade;

if nargout > 1  || ~isempty(varargout(2))
    [~, cOnset] = getChoice(dd); %cOnset time including FAIL
    [outcome, cueOnset] = getCueOnset(dd);
    fOnset = getFOnset(dd);
end

for itr = 1:nTrials
    x_cat = cat(1, x_cat, eyeData(itr).x);
    y_cat = cat(1, y_cat, eyeData(itr).y);
    parea_cat = cat(1, parea_cat, eyeData(itr).parea);
    phght_cat = cat(1, phght_cat, eyeData(itr).phght);
    pwdth_cat = cat(1, pwdth_cat, eyeData(itr).pwdth);

    %     if isempty(t_cat)
    %         t_cat = eyeData(itr).t;
    %     else
    %         t_cat = cat(1, t_cat, max(t_cat)-eyeData(itr).t(1)+eyeData(itr).dt+eyeData(itr).t);
    %     end
    if isempty(t_cat)
        t0 = eyeData(itr).t(1);
    elseif ~isempty(eyeData(itr).t)
        t0 =  max(t_cat)-eyeData(itr).t(1)+eyeData(itr).dt;
    end
    if ~isempty(eyeData(itr).t)
        t_cat = cat(1, t_cat, eyeData(itr).t+t0);
    end

    if nargout > 1  || ~isempty(varargout(2))
        %         for idd = 1:length(ddnames)
        %             dd_cat.(ddnames{idd}) = cat(1, dd_cat.(ddnames{idd}), dd.(ddnames)(itr)+t0);
        %         end
        for ifd = 1:length(fdnames)
            onsets_cat.(fdnames{ifd}) = cat(1, onsets_cat.(fdnames{ifd}), ...
                dd.(fdnames{ifd})(itr)+t0);
        end
        %onsets_cat = structfun(@(x,y)(cat(1,x,y(itr)+t0)), onsets_cat,dd); NG
        onsets_cat.cOnset(itr,1) = cOnset(itr) + t0;
        onsets_cat.cueOnset(itr,1) = cueOnset(itr) + t0;
        onsets_cat.fOnset(itr,1) = fOnset(itr) + t0;

    end
    if nargout > 2 || ~isempty(varargout(3))
        for iedf = 1:length(edfnames)
            meta_cat.(edfnames{iedf}) = cat(1, meta_cat.(edfnames{iedf}), ...
                (dd.meta.edf.(edfnames{iedf})('trial',itr).time)' ...
                - dd.meta.cic.firstFrame('trial',itr).time+t0);
        end
    end
end
t_cat = t_cat(~isnan(t_cat));
[t_cat, ix] = unique(t_cat);
eyeData_cat = marmodata.eye(t_cat, x_cat(ix), y_cat(ix), pwdth_cat(ix), phght_cat(ix));


%% check STARTBLINK/ENDBLINK
if length(meta_cat.STARTBLINK) >  length(meta_cat.ENDBLINK)
    meta_cat.ENDBLINK = [meta_cat.ENDBLINK; eyeData_cat.t(end)];
elseif length(meta_cat.STARTBLINK) <  length(meta_cat.ENDBLINK)
    meta_cat.STARTBLINK = [meta_cat.STARTBLINK; eyeData_cat.t(1)];
end


if ~isempty(find(meta_cat.ENDBLINK - meta_cat.STARTBLINK < 0))
    error('ENDBLINK comes earlier than STARTBLINK');
end

%% check STARTSACC/ENDSACC 28/1/22
if length(meta_cat.STARTSACC) ~= length(meta_cat.ENDSACC)%41
    nSaccs = min(length(meta_cat.STARTSACC), length(meta_cat.ENDSACC));
    ngIdx=find(meta_cat.ENDSACC(1:nSaccs)-meta_cat.STARTSACC(1:nSaccs)<0);

    if isempty(ngIdx)
        meta_cat.STARTSACC = meta_cat.STARTSACC(1:nSaccs);
        meta_cat.ENDSACC = meta_cat.ENDSACC(1:nSaccs);
        %else
        % FILL ME?
    end
end
idx = find(meta_cat.STARTBLINK-meta_cat.ENDBLINK~=0);
meta_cat.STARTBLINK = meta_cat.STARTBLINK(idx);
meta_cat.ENDBLINK = meta_cat.ENDBLINK(idx);

assert(isempty(find(meta_cat.ENDSACC-meta_cat.STARTSACC<0)));
okSacc = find(meta_cat.ENDSACC-meta_cat.STARTSACC>0);
meta_cat.STARTSACC = meta_cat.STARTSACC(okSacc);
meta_cat.ENDSACC = meta_cat.ENDSACC(okSacc);