-
Notifications
You must be signed in to change notification settings - Fork 0
/
purifyDwell.m
58 lines (47 loc) · 2.03 KB
/
purifyDwell.m
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
function dwellData_s = purifyDwell(dwellData, signif_trans)
%transition magnitude vector
dif = diff([dwellData(:, 2)'; dwellData(:,1)']);
if isempty(find(abs(dif) <= signif_trans, 1))
dwellData_s=dwellData;
disp('No unsignificant transition')
end
%repeat until there are no more unsignificant transitions
while length(find(abs(dif) <= signif_trans))
ms = size(dwellData);
%transition magnitude vector
dif = diff([dwellData(:, 2)'; dwellData(:,1)']);
%locate significant & unsignificant transitions
signif = find(abs(dif) > signif_trans | isnan(dif));
unsignif = find(abs(dif) <= signif_trans);
%for all unsignificant transitions,
%add dwelltime to next transition
%and change the start FRET of next
for n = unsignif
if n < ms(1)
%same trace
if dwellData(n+1, 1) == dwellData(n, 2)
dwellData(n+1, 3) = dwellData(n+1, 3) + dwellData(n, 3);
dwellData(n+1, 1) = dwellData(n, 1);
%new trace - add transition to NaN
else
dwellData = [dwellData; dwellData(n, 1) NaN dwellData(n, 3)];
ms = size(dwellData);
signif = [signif ms(1)];
%last trace - add transition to NaN
end
else
dwellData = [dwellData; dwellData(n, 1) NaN dwellData(n, 3)];
ms = size(dwellData);
signif = [signif ms(1)];
end
end
%delete unsignificant transitions
i = 1;
for n = signif
dwellData_s(i, :) = dwellData(n, :);
i = i + 1;
end
disp(['deleted ' int2str(length(unsignif)) ' unsignificant transitions'])
dwellData = dwellData_s;
dif = diff([dwellData(:, 2)'; dwellData(:,1)']);
end