SingleMoleculeAnalysisScripts/MakeSeparateTraces_vbFRET.m at main · NilsSomm/SingleMoleculeAnalysisScripts · 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
clc;
clear;
close all;

% bothFluct = load('000_070 35 mW Cy3 10 mW A488 5Hz 10 min(00002)_bgsub_crop.mat');
cy5Fluct = load('070_070 60mW Cy5 25nW Cy3 10 mW A488 5Hz 10 min(00011)_bgsub_finesave_crop.mat');
% cy5Fluct = load('Cy5fluctuatingCy3high.mat');

ogname = ['070_070.mat'];

%% Normalize All Traces
ind=0;
for datasets = {cy5Fluct}
    datasets{1} = normalizeTraces(datasets{1});
    cy3Traces = [];
    cy3data = [];
    cy5Traces = [];
    cy5data = [];

    %% Generate Complementary traces
    for i = 1:length(datasets{1}.traces)
        ogtraces = datasets{1}.traces;
        cy3Traces(i).acptr = ogtraces(i).donr;
        cy3Traces(i).donr = 1 - ogtraces(i).donr;
        cy3data{i}(:, 2) = ogtraces(i).donr;
        cy3data{i}(:, 1) = 1 - ogtraces(i).donr;

        cy5Traces(i).donr = 1 - ogtraces(i).acptr;
        cy5Traces(i).acptr = ogtraces(i).acptr;
        cy5data{i}(:, 1) = 1 - ogtraces(i).acptr;
        cy5data{i}(:, 2) = ogtraces(i).acptr;
    end

    % if ind == 0
    %     ogname = "both_fluctuating.mat";
    % elseif ind == 1
    %     ogname = "Cy3fluctuatingCy5high.mat";
    % else
    %     ogname = "Cy5fluctuatingCy3high.mat";
    % end


    traces = cy3Traces;
    data = cy3data;
    save(strcat("vb_cy3_", ogname), "traces", "data")

    traces = cy5Traces;
    data = cy5data;
    save(strcat("vb_cy5_", ogname), "traces", "data")

    ind=ind+1;
end


%% Export Both Channels Fluctuating


function output = normalizeTraces(input)
    [donrmin, donrmax, acptrmin, acptrmax] = getMinMax(input);
    for i = 1:size(input.traces,2)
        % trace = input.traces(i);
        input.traces(i).donr = (input.traces(i).donr - donrmin)/(donrmax-donrmin);
        input.traces(i).acptr = (input.traces(i).acptr - acptrmin)/(acptrmax-acptrmin);
    end
    output = input;
end

function [donrmin, donrmax, acptrmin, acptrmax] = getMinMax(sav)
    donrmax = max(arrayfun(@(x) max(x.donr), sav.traces));
    acptrmax = max(arrayfun(@(x) max(x.acptr), sav.traces));
    donrmin = min(arrayfun(@(x) min(x.donr), sav.traces));
    acptrmin = min(arrayfun(@(x) min(x.acptr), sav.traces));
end