Team:DTU-Denmark/plot mutant2

From 2013.igem.org

(Difference between revisions)
(Created page with "{{:Team:DTU-Denmark/Templates/StartPage|}} ==Timer series plots of Mutant2== <pre> code </pre> {{:Team:DTU-Denmark/Templates/EndPage}}")
(Timer series plots of Mutant2)
 
Line 4: Line 4:
<pre>
<pre>
-
code
+
clear all
 +
% create model
 +
model=SBmodel('mutant_2.txt')
 +
SBPDmakeMEXmodel(model,'modelmex');
 +
 
 +
%calculate steady state
 +
SBsteadystate(model)
 +
 
 +
conc = [5.872e-2 5.872e-1 5.872 29.358];
 +
enz = [2.115e-4 1.0575e-03 2.1115e-3 2.115e-1];
 +
count = 0;
 +
for i=1:4
 +
    count = 0;
 +
    for j=1:4
 +
        count = count+1;
 +
        % make parametervector
 +
        parametervector = makeparamvecSBPD('modelmex',{'NO20','E'},[conc(i); enz(j)]);
 +
        output=modelmex([0:0.001:10],[],parametervector);
 +
        % save results in separate vectors
 +
        % check in which order they are (double click on output in workspace)
 +
        NO(:,count) = output.statevalues(:,1);
 +
        N2O(:,count) = output.statevalues(:,2);
 +
    end
 +
    t = output.time';
 +
    % plot concentrations over time:
 +
    figure
 +
    plot(t,N2O(:,1),t,N2O(:,2),'--',t,N2O(:,3),':',t,N2O(:,4),'-.','LineWidth',2)
 +
    if i == 1
 +
        title('1 mg/l nitrite (58.72 uM)')
 +
    elseif i == 2
 +
        title('10 mg/l nitrite (587.2 uM)')
 +
    elseif i == 3
 +
        title('100 mg/l nitrite (5.872 mM)')
 +
    elseif i == 4
 +
        title('500 mg/l nitrite (29.358 mM)')
 +
    end
 +
   
 +
    xlabel('Time in [min]','FontSize',16)
 +
    ylabel('N2O concentration in [mM]','FontSize',16)
 +
    hleg1 = legend('100 enzymes per cell','500 enzymes per cell','1000 enzymes per cell','100 000 enzymes per cell');
 +
    set(hleg1,'Location','SouthEast','FontSize',12)
 +
end
</pre>
</pre>
{{:Team:DTU-Denmark/Templates/EndPage}}
{{:Team:DTU-Denmark/Templates/EndPage}}

Latest revision as of 18:44, 4 October 2013

Timer series plots of Mutant2

clear all
% create model
model=SBmodel('mutant_2.txt')
SBPDmakeMEXmodel(model,'modelmex');

%calculate steady state
SBsteadystate(model)

conc = [5.872e-2 5.872e-1 5.872 29.358];
enz = [2.115e-4 1.0575e-03 2.1115e-3 2.115e-1];
count = 0;
for i=1:4
    count = 0;
    for j=1:4
        count = count+1;
        % make parametervector
        parametervector = makeparamvecSBPD('modelmex',{'NO20','E'},[conc(i); enz(j)]);
        output=modelmex([0:0.001:10],[],parametervector);
        % save results in separate vectors
        % check in which order they are (double click on output in workspace)
        NO(:,count) = output.statevalues(:,1);
        N2O(:,count) = output.statevalues(:,2);
    end
    t = output.time';
    % plot concentrations over time:
    figure
    plot(t,N2O(:,1),t,N2O(:,2),'--',t,N2O(:,3),':',t,N2O(:,4),'-.','LineWidth',2)
    if i == 1
        title('1 mg/l nitrite (58.72 uM)')
    elseif i == 2
        title('10 mg/l nitrite (587.2 uM)')
    elseif i == 3
        title('100 mg/l nitrite (5.872 mM)')
    elseif i == 4
        title('500 mg/l nitrite (29.358 mM)')
    end
    
    xlabel('Time in [min]','FontSize',16)
    ylabel('N2O concentration in [mM]','FontSize',16)
    hleg1 = legend('100 enzymes per cell','500 enzymes per cell','1000 enzymes per cell','100 000 enzymes per cell');
    set(hleg1,'Location','SouthEast','FontSize',12)
end