Analysis-code/BleachUpTrace.m at master · BioMolecularSelfOrganizationLab/Analysis-code · 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
function BleachUpTrace

Brate = 0.009;
Nc = 200;
cycles = 1;
C = 1;
N = 200;
DiffDown = zeros(1,N);
MonDiffBleachTrace = zeros(4000,N);
h = waitbar(0,'Simulating, please wait...');

for j = 1:N
    waitbar(j / N)
    [simZMW] = SIM_Gelsolin(C,cycles,Brate);
    MonTrace = cell2mat(simZMW.bleach);

    trace_d = diff(MonTrace);
    up = find(trace_d>0);
    UpTrac = zeros(1,length(MonTrace));
    UpTrac(up) = 1;
    UpTrace = cumsum(UpTrac);
    %up_c = [0; up];
    %du = diff(up_c);
    %     if i==1
    %         dup = du;
    %     else
    %         dup = [dup;du];
    %     end
    %

%     Bleach = zeros(4000,1);
    BleachUP = zeros(4000,Nc);
%     Bleach(:,1) = ZMW_sim_Bleaching(MonTrace,Brate);


    for i=1:Nc

        %     Bleach(:,i) = ZMW_sim_Bleaching(MonTrace,Brate);
        BleachUP(:,i) = ZMW_sim_Bleaching(UpTrace,Brate);

        Ndo(i) = length(find(diff(BleachUP(:,i))<0));

    end

    MonDiffBleachTrace(:,j) = MonTrace - mean(BleachUP,2);

    NdoMT = length(find(diff(MonTrace)<0));
    DiffDown(j) = NdoMT-mean(Ndo);

end

% figure;
% [n,x] = hist(DiffDown,calcnbins(DiffDown,'fd'));
% plot(x,n,'')

MonDBT = mean(MonDiffBleachTrace,2);

figure;
plot(MonDBT,'')

% figure;
% [n,x] = hist(MonDBT,calcnbins(MonDBT,'fd'));
% plot(x,n,'')

close(h)
% figure;
% plot(UpTrace, 'r');
% hold on;
% plot(MonTrace, 'k');
% hold off;
%
% figure;
% plot(BleachUP,'r')
% hold on
% plot(Bleach,'b')
% hold off