Code covered by the BSD License  

Highlights from
Using SimBiology for Mechanism-Based PK/PD Modeling in Preclinical & Discovery

image thumbnail

Using SimBiology for Mechanism-Based PK/PD Modeling in Preclinical & Discovery

by

 

Examples used in the March 2012 webinar: Using SimBiology for Mechanism-Based PK/PD Modeling

plotResponseCurves(hax, results, fitInfo)
function plotResponseCurves(hax, results, fitInfo)
% Copyright 2011 - 2012 MathWorks, Inc.

xl      = [2e-2 2e4] ;
xtck    = logspace(-1,5, 4);
xsc     = 'log' ;
xlfit   = logspace(-1,5,30) ;
ytckLbl = '' ;
ytck    = [0 3 6 9] ;
  
% Plot
hold(hax, 'on')
hl(1) = plot(hax, results.Dose, results.log10CFU_at_24hr, 'bo:' ) ; 

% Set axes properties
set(hax, 'Xlim'     , xl        ,...
         'XTick'    , xtck      ,...
         'XScale'   , xsc       ,...
         'YTick'    , ytck      ,...
         'Ylim'     , [0 10]    ,...
         'YTickLabel', ytckLbl  ) ;
     
xlabel(hax, 'Dose',  'FontName', 'Calibri',  'FontSize', 12) ;
box on ;


if ~isempty(fitInfo)

    hl(2) = plot(hax, xlfit, fitInfo.model(xlfit), 'r--', 'LineWidth', 2.5);
    axes(hax)
    hl(3) = text( 1.2*xl(1), 6.5, num2str(fitInfo.bstat),...
        'FontWeight', 'bold', 'FontName', 'Calibri', 'FontSize', 10 ) ; 
    hl(4) = text( 1.2*xl(1), 2.5, num2str(fitInfo.bcid),...
        'FontWeight', 'bold', 'FontName', 'Calibri', 'FontSize', 10 ) ;

    hl(5) = line([xl(1), xl(end)], [3 3], 'LineStyle', '--', 'LineWidth', 2.5, 'Color', [0.5 0.5 0.5]) ;
    hl(6) = line([xl(1), xl(end)], [6 6], 'LineStyle', '--', 'LineWidth', 2.5, 'Color', [0.5 0.5 0.5]) ; 

end

legend('boxoff')
if strcmp(get(hax, 'Visible'), 'off')
    set(hl, 'Visible', 'off') 
end
    



Contact us