clear all
% Data extaction:
directory = 'C:\Users\Usuario\Documents\Carlos AV\UPV\iGEM\Modelling\Experimentos\Comparison_GFP_YFP\';
sheet = 'Plate 1 - Sheet1';
fname = 'Experiment_GFP_YFP_data';

FOD = xlsread([directory fname '.xlsx'],sheet,'L438:AA507');
% Time vector:
    % Total time: 06:00:00 (HH:MM:SS)
    % Interval: 05:00 (MM:SS)
time = 5:5:(6*60-10); % We have detected a NaN value at the end of FOD vector
FOD_GFP = FOD(:,1:8);
FOD_YFP = FOD(:,9:16);

save([directory fname '.mat'],'FOD_GFP','FOD_YFP','FOD','time');
%% GRAPHICS
clear all
directory = 'C:\Users\Usuario\Documents\Carlos AV\UPV\iGEM\Modelling\Experimentos\Comparison_GFP_YFP\';
fname = 'Experiment_GFP_YFP_data';
load([directory fname '.mat'])

% Corrected FOD with excitation/emission wavelengths efficiency 
% (490/530 nm) and gain
Gain = 60;
GFP_Exc_eff = 0.6054;
GFP_Em_eff = 0.585;
YFP_Exc_eff = 0.3246;
YFP_Em_eff = 0.9952;

MEAN_FOD_GFP = mean(FOD_GFP,2)/(GFP_Exc_eff*GFP_Em_eff*Gain);
MEAN_FOD_YFP = mean(FOD_YFP,2)/(YFP_Exc_eff*YFP_Em_eff*Gain);

% We select a stationary period of time
fig_comp= figure('Color','w');
% Create axis
axes = axes('Parent',fig_comp);
hold(axes,'on');
plot(time,MEAN_FOD_GFP,'r',time,MEAN_FOD_YFP,'LineWidth',1)
title('Comparison between relative intensity of YFP and GFP reporter proteins','FontName','Lato','FontSize',12);
ylabel('Fluorescence (RFU)','FontName','Lato')
xlabel('Time (min)','FontName','Lato')
legend({'GFP','YFP'},'FontName','Lato','Location','east');
xlim([time(1) time(end)])

box(axes,'on');
set(axes,'FontName','Lato','XGrid','on','YGrid','on');


% Number of equivalent fluorescein molecules
% 273.7 MEFL/RFU(GFP)
MOL = 273.7*MEAN_FOD_GFP(40:end);

% Number of MEFL = Number of GFP molecules = Number of YFP molecules

% MEFL/RFU factor
K = mean(MOL./MEAN_FOD_YFP(40:end))