Analysis-code/average_traces.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
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
%% average traces: DNA       **********************************************

%% raw trace
load('C:\Users\Rafael\Documents\ITQB estágio\DNA\DNA4_8');
plot(mean(acc,1));

%% step trace
load('C:\Users\Rafael\Documents\ITQB estágio\code and pdf\DNA4_8 (stepSizes)');

plot(mean(stepFit,1));

%% num trace

load('C:\Users\Rafael\Documents\ITQB estágio\code and pdf\DNA4_8 (stepSizes)');

plot(mean(numFit,1));


%% average traces: LatA      **********************************************

%% raw trace
files = {...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw1_1uMactin_LatA_nobump_thr150_Salapaka50',...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw5_1uM_LatA_cut_thr150_Salapaka50'};

raw = [];

for file = 1:length(files)
    load(files{file});

    raw(end+1:end+size(Raw,1),1:3968) = Raw(:,1:3968);
end

plot(mean(raw,1));

%% step trace
files = {...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw1_1uMactin_LatA_nobump_thr150_Salapaka50',...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw5_1uM_LatA_cut_thr150_Salapaka50'};

step = [];

for file = 1:length(files)
    load(files{file});

    step(end+1:end+size(Step,1),1:3968) = Step(:,1:3968);
end

plot(mean(step,1));


%% num trace
files = {...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw1_1uMactin_LatA_nobump_thr150_Salapaka50',...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw5_1uM_LatA_cut_thr150_Salapaka50'};

num = [];

for file = 1:length(files)
    load(files{file});
    for i = 1:size(Step,1)
        fit = getNumSteps(Step(i,:));
        num(end+1,1:3968) = fit(:,1:3968);
    end
end

plot(mean(num,1));


%% average traces: Actin      **********************************************

%% raw trace
files = {...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw1_1uMactin_nobump_thr300_Salapaka50',...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw2_BSA_nobump_thr300_Salapaka50'};

raw = [];

for file = 1:length(files)
    load(files{file});

    raw(end+1:end+size(Raw,1),1:3909) = Raw(:,1:3909);
end

plot(mean(raw,1));

%% step trace
files = {...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw1_1uMactin_nobump_thr300_Salapaka50',...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw2_BSA_nobump_thr300_Salapaka50'};

step = [];

for file = 1:length(files)
    load(files{file});

    step(end+1:end+size(Step,1),1:3909) = Step(:,1:3909);
end

plot(mean(step,1));


%% num trace
files = {...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw1_1uMactin_nobump_thr300_Salapaka50',...
    'C:\Users\Rafael\Documents\ITQB estágio\code and pdf\zmw2_BSA_nobump_thr300_Salapaka50'};

num = [];

for file = 1:length(files)
    load(files{file});
    for i = 1:size(Step,1)
        fit = getNumSteps(Step(i,:));
        num(end+1,1:3909) = fit(:,1:3909);
    end
end

plot(mean(num,1));