Code covered by the BSD License  

Highlights from
Response Surface Models of Drug Interactions with Curve Fitting Toolbox

image thumbnail

Response Surface Models of Drug Interactions with Curve Fitting Toolbox

by

 

20 Mar 2009 (Updated )

"Response Surface Models of Drug Interactions with Curve Fitting Toolbox" webinar demo files

myCreateSurfaceFit(Propofol, Remifentanil, Effect, EffectName)
function [fitresult, gof] = myCreateSurfaceFit(Propofol, Remifentanil, Effect, EffectName)
%CREATESURFACEFIT2(PROPOFOL,REMIFENTANIL,Effect)
%  Fit surface to data.
%
%  Data for 'Model' fit:
%      X Input : Propofol
%      Y Input : Remifentanil
%      Z output: Effect
%      Weights : (none)
%
%  Output:
%      fitresult : an sfit object representing the fit.
%      gof : structure with goodness-of fit info.
%
%  See also FIT, SFIT.

%  Auto-generated by MATLAB on 16-Mar-2009 13:34:42

% Initialize arrays to store fits and goodness-of-fit.
% Copyright 2009, The MathWorks, Inc.

fitresult = cell( 1, 1 );
gof = struct( 'sse', cell( 1, 1 ), ...
    'rsquare', [], 'dfe', [], 'adjrsquare', [], 'rmse', [] );

%% Fit: 'Model'.
ft = fittype( 'combinedEffect(x,y, IC50A, IC50B, alpha, n)', 'indep', {'x', 'y'}, 'depend', 'z' );
opts = fitoptions( ft );
opts.Display = 'Off';
opts.Lower = [0 0 -5 -0];
opts.Robust = 'LAR';
opts.StartPoint = [0.00893838724332152 0.706165672266879 1 0.746030748284422];
opts.Upper = [Inf Inf Inf Inf];
opts.Weights = zeros(1,0);
[fitresult{1}, gof(1)] = fit( [Propofol, Remifentanil], Effect, ft, opts ) ; 

% Create a figure for the plots.
f = figure( 'Name', 'Model', 'Color', [1 1 1]);

% Plot fit with data.
subplot1 = subplot( 2, 1, 1 );
h = plot( fitresult{1}, [Propofol, Remifentanil], Effect );
grid on
view(subplot1 , [-60 20]);
% Label axes
xlabel( 'Propofol' );
ylabel( 'Remifentanil' );
zlabel( EffectName );
% legend( h, 'Model Fit', 'Data', 'Location', 'NorthEast' );
set(h(1), 'FaceColor',[0.702 0.7804 1],...
          'FaceAlpha',0.6,...
          'EdgeAlpha',0.5);
      
set(h(2), 'MarkerFaceColor',[0, 0 ,1],...
          'MarkerEdgeColor',[0 0 1]  ,...
          'MarkerSize'     , 5);
       

% Plot residuals.
subplot2 = subplot( 2, 1, 2 );
h = plot( fitresult{1}, [Propofol, Remifentanil], Effect, 'Style', 'Residual' );
grid on
view(subplot2,[-60 20]);
% Label axes
xlabel( 'Propofol' );
ylabel( 'Remifentanil' );
zlabel( ['Residuals ', EffectName] );
% legend( h, 'Model - residuals', 'Location', 'NorthEast' );

set(h, 'MarkerFaceColor',[0, 0 ,1],...
        'MarkerEdgeColor',[0 0 1]  ,...
        'MarkerSize'     , 5);
      

Contact us