生物制药数据表面与自定义拟合方程
这个例子展示了如何使用曲线拟合工具箱™适合反应表面麻醉数据来分析药物交互作用。响应面模型提供一个好的方法对于理解药物的药效学相互作用行为的组合。
这个数据是基于结果摘要:克恩,谢G,白色杰,伊根道明。Opioid-hypnotic协同效应:propofol-remifentanil药效学相互作用的响应面分析志愿者。麻醉学2004;100:1373 - 81。
麻醉通常至少两种药物的过程,包括阿片类和镇静催眠。这个示例使用异丙酚和Reminfentanil药物类原型。他们的交互是由四种不同的镇痛和镇静措施应对药物组合。痛觉计,手足抽搐、镇静和Laryingoscopy组成的四个措施替代药物效应在不同浓度的异丙酚和Reminfentanil组合。
下面的代码,使用曲线拟合工具箱方法,再现了互动与曲线拟合表面构建工具中描述表面生物制药数据拟合。
加载数据
从文件加载数据。
数据= importdata (“OpioidHypnoticSynergy.txt”);异丙酚= data.data (: 1);Remifentanil = data.data (:, 2);痛觉计= data.data (: 3);强直= data.data (: 4);镇静= data.data (: 5);Laryingoscopy = data.data (: 6);
创建模型适合类型
您可以使用fittype函数定义的模型,在哪里CA和CB药物浓度,IC50A,IC50B,α,n估计的系数。创建模型适合类型。
英国《金融时报》= fittype (“Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1) ',…“独立”,{“CA”,“CB”},“依赖”,“z”,“问题”,“Emax”)
英尺=一般模型:英尺(IC50A IC50B,α,n, Emax, CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)
假设Emax = 1因为输出归一化的影响。
Emax = 1;
设置合适的选项
设置合适的选项为健壮的拟合,界限,开始点。
选择= fitoptions(英尺);选择。低= [0 0 5 0];选择。健壮的=“守护神”;选择。曾经繁荣= (0.0089,0.706,1.0,0.746);
适合和情节痛觉计的表面
(f, gof) =适合([异丙酚,Remifentanil],痛觉计,英国《金融时报》,…选择,“问题”Emax)
成功,但拟合停止因为残差变化量小于公差(TolFun)。
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 4.149 (4.123, 4.174) IC50B = 9.044(8.971, 9.118)α= 8.502 (8.316,8.687)n = 8.289(8.131, 8.446)问题参数:Emax = 1
gof =结构体字段:上交所:0.0842 rsquare: 0.9991教育部:393 adjrsquare: 0.9991 rmse: 0.0146
情节(f,异丙酚,Remifentanil],痛觉计);
适合表面手足抽搐
重用相同的fittype创建一个响应面手足抽搐。
(f, gof) =适合([异丙酚,Remifentanil],手足抽搐,英国《金融时报》,选择,“问题”Emax)
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 4.544 (4.522, 4.567) IC50B = 21.22(21.04, 21.4)α= 14.94 (14.67,15.21)n = 6.132(6.055, 6.209)问题参数:Emax = 1
gof =结构体字段:上交所:0.0537 rsquare: 0.9993教育部:393 adjrsquare: 0.9993 rmse: 0.0117
情节(f,异丙酚,Remifentanil],强直);
适合表面镇静
(f, gof) =适合(异丙酚,Remifentanil,镇静,英国《金融时报》,选择,“问题”Emax)
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 1.843 (1.838, 1.847) IC50B = 13.7(13.67, 13.74)α= 1.986 (1.957,2.015)n = 44.27(42.56, 45.98)问题参数:Emax = 1
gof =结构体字段:上交所:0.0574 rsquare: 0.9994教育部:393 adjrsquare: 0.9994 rmse: 0.0121
情节(f,异丙酚,Remifentanil],镇静);
适合Laryingoscopy表面
(f, gof) =适合([异丙酚,Remifentanil], Laryingoscopy,英国《金融时报》,选择,“问题”Emax)
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 5.192 (5.177, 5.207) IC50B = 37.77(37.58, 37.97)α= 19.67 (19.48,19.86)n = 37(35.12, 38.87)问题参数:Emax = 1
gof =结构体字段:上交所:0.1555 rsquare: 0.9982教育部:393 adjrsquare: 0.9982 rmse: 0.0199
情节(f,异丙酚,Remifentanil], Laryingoscopy);
