此示例显示如何使用曲线拟合工具箱™将响应曲面适合某些麻醉数据以分析药物交互效果。响应面模型提供了理解药物组合的药效互动行为的良好方法。
此数据基于本文的结果:Kern SE,Xie G,White JL,Egan Td。阿片类药物催眠协同:志愿者在志愿者中Photofol-Remifentanil药效互动的响应表面分析。麻醉学到2004;100:1373-81。
麻醉通常是至少一种双药方法,由阿片类药物和镇静催眠组成。此示例使用异丙酚和Reminfentanil作为药物类原型。它们的相互作用是通过对药物组合的四种不同衡量标和镇静反应的四种不同措施来测量。Algometry,Tetany,镇静和喉镜检查包括异丙酚和Reminfentanil各种浓度组合的四种替代药物效果的衡量标准。
以下代码使用曲线拟合工具箱方法,再现交互式表面构建与曲线配件工具表面适合生物制药数据。
从文件加载数据。
data = importData('opioidhypnoticsynergy.txt');propofol = data.data(:,1);remifentanil = data.data(:,2);Algometry = data.data(:,3);tetany = data.data(:,4);镇静= data.data(:,5);learyingoscopy = data.data(:,6);
你可以使用fittype.
从纸张中定义模型的功能,在哪里加利福尼亚州
和CB.
是药物浓度,还有药物浓度IC50A
那IC50B.
那α
, 和N.
是要估计的系数。创建模型适合类型。
ft = fittype('Emax *(CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N /((CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB /IC50B)))^ n + 1)'那......'独立的',{'CA'那'cb'},'依赖'那'z'那'问题'那'emax')
ft =一般型号:ft(IC50a,IC50b,alpha,n,Emax,CA,Cb)= Emax *(CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N /((CA / IC50A + CB / IC50B + ALPHA *(CA / IC50A)*(CB / IC50B))^ N + 1)
认为Emax = 1
因为效果输出归一化。
Emax = 1;
为强大的拟合,界限和起点设置适合选项。
选择= fitoptions(ft);opts.lower = [0,0,-5,-0];opts.robust ='lar';opts.startpoint = [0.0089,0.706,1.0,0.746];
[F,GOF] = Fit([Pronofol,雷芬丹],Algometry,Ft,......选择,'问题',Emax)
成功,但拟合停止,因为剩余的变化小于容忍度(Tolfun)。
一般型号:F(CA,CB)= EMAX *(CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N /((CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N + 1)系数(具有95%的置信度):IC50A = 4.148(4.123,4.173)IC50B = 9.043(8.97,9.117)Alpha = 8.499(8.314,8.685)n = 8.294(8.136,8.452)问题参数:Emax = 1
GOF =结构与字段:SSE:0.0842 RSQUARE:0.9991 DFE:393 Addrsquare:0.9991 RMSE:0.0146
图(F,[PROPOFOL,REMIFENTANIL],ALPORERY);
重用相同的fittype.
为Tetany创建响应面。
[F,GOF] = FIT(【Pronofol,雷芬丹],Tetany,Ft,Opts,'问题',Emax)
一般型号:F(CA,CB)= EMAX *(CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N /((CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N + 1)系数(具有95%的置信度):IC50A = 4.544(4.522,4.567)IC50B = 21.22(21.04,21.4)Alpha = 14.94(14.67,15.21)n = 6.132(6.055,6.209)问题参数:Emax = 1
GOF =结构与字段:SSE:0.0537 RSQUARE:0.9993 DFE:393 Adjrsquare:0.9993 RMSE:0.0117
图(F,[PROPOFOL,REMIFENTANIL],TETANY);
[F,GOF] =配合([异丙酚,雷芬丹],镇静,FT,OPTS,'问题',Emax)
一般型号:F(CA,CB)= EMAX *(CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N /((CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N + 1)系数(具有95%的置信度):IC50A = 1.843(1.838,1.847)IC50B = 13.7(13.67,13.74)Alpha = 1.986(1.957,2.015)n = 44.27(42.56,45.98)问题参数:Emax = 1
GOF =结构与字段:SSE:0.0574 RSQUARE:0.9994 DFE:393 ADJRSQUARE:0.0994 RMSE:0.0121
图(F,[异丙酚,雷芬丹],镇静);
[F,GOF] = Fit(【PrOPOFOL,雷芬丹],喉镜,FT,OPTS,'问题',Emax)
一般型号:F(CA,CB)= EMAX *(CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N /((CA / IC50A + CB / IC50B + Alpha *(CA / IC50A)*(CB / IC50B))^ N + 1)系数(具有95%的置信度):IC50A = 5.192(5.177,5.207)IC50B = 37.77(37.58,37.97)Alpha = 19.67(19.48,19.86)n = 37(35.12,38.87)问题参数:Emax = 1
GOF =结构与字段:SSE:0.1555 RSQUARE:0.9982 DFE:393 Adjrsquare:0.0199
图(F,[异丙酚,雷芬丹],喉镜检查);