Behavioral-and-Neural-Analysis-Toolkit/getConsecutiveTrialIndex.m at master · CPapadim/Behavioral-and-Neural-Analysis-Toolkit · 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
function [trialscurr, trialsprev, trialsIDcurr, trialsIDprev, keepidxCurr, keepidxPrev] = getConsecutiveTrialIndex(filename, previoustrialsID, currenttrialsID)
% Generate indices for current and previous trials to keep that are consecutive
%
% Inputs:
%		filename - name of file to operate on
%		previoustrialsID - previous trial IDs
%		currenttrialsID - current trial IDs
%
% Output:
%		keepidxCurr - index of current trials to keep
%		keepidxPrev - index of previous trials to keep
%
% NOTE:  Indices can be used like this:  previoustrials=previoustrials(keepidxPrev), etc

	[trialsID trialNums stacks classes success] = getTrialData(filename, '', '');
	trialNumsPadded=[trialNums; -99999]; % Pad to make diff vector the same size as the rest

	trialNumsConsecutive=trialNums(diff(trialNumsPadded)==1 & ismember(trialsID, previoustrialsID) & ismember(trialsID+1,currenttrialsID));
	trialNumsConsecutiveID=trialsID(diff(trialNumsPadded)==1 & ismember(trialsID, previoustrialsID) & ismember(trialsID+1,currenttrialsID));

	keepidxCurr=ismember(trialsID,trialNumsConsecutiveID+1);
	keepidxPrev=ismember(trialsID+1,trialsID(keepidxCurr));

	trialsprev=trialNums(keepidxPrev);
	trialscurr=trialNums(keepidxCurr);
	trialsIDprev=trialsID(keepidxPrev);
	trialsIDcurr=trialsID(keepidxCurr);


	%previoustrials=previoustrials(1:end-1);

	%trialsconsecutive=previoustrials(diff(previoustrials)==1);
	%trialsconsecutiveID=previoustrialsID(diff(previoustrials)==1);

	%keepidxCurr=ismember(currenttrialsID,trialsconsecutiveID+1);
	%keepidxPrev=ismember(previoustrialsID+1,currenttrialsID(keepidxCurr));


end