n运输单分布隔间PK simbiology吸收模型
5视图(30天)
显示旧的评论
我已经建立了一个与单一distibutional n-transit室药代动力学使用SimBiology吸收模型。这个模型是基于交通模型种以下出版物:萨维奇,智慧化彼特,D.M.Kerbusch, t . et al . Pharmacokinet Pharmacodyn 34 (2007): 711。
https://doi.org/10.1007/s10928 - 007 - 9066 - 0
我选择先加载一个隔间吸收模型与间隙parmaeterization(宏观参数),然后修改了吸收反应方程n运输舱基于解析解来确定运输舱数量作为连续变量(而不是整数)。
分析解决方案的隔间:
![n分析solution.PNG](http://www.tianjin-qmedu.com/matlabcentral/answers/uploaded_files/214672/n%20analytical%20solution.png)
单一分布间歌唱:
![nTransit 1 comp ODE.PNG](http://www.tianjin-qmedu.com/matlabcentral/answers/uploaded_files/214673/nTransit%201Comp%20ODE.png)
变换的颂歌n很大时数值稳定:
![ODE transform.PNG](http://www.tianjin-qmedu.com/matlabcentral/answers/uploaded_files/214674/ODE%20transform.png)
武官总结报告是我在simbiology定义模型。
我不确定如果我用错了语法,但我看到负面cincentrations第一次模拟时间点由这仿真图如图所示:
![sim.PNG](http://www.tianjin-qmedu.com/matlabcentral/answers/uploaded_files/214675/sim.png)
仿真应该是一个逐渐吸收阶段与典型tlag吸收simbiology PK模型中的模型库。
欢呼,
丹
答案(1)
Wojciech Krzyzanski
2019年4月19日
亲爱的丹,
积极的投入到Aa舱和0初始值(Aa(0) = 0),你的氢氧化钠Aa (t)应该总是正的。你的模拟吸收曲线意味着ODE求解器的问题。既然你提到中转仓的数量很大,我认为这是一个刚度问题,你的ODE求解器不适合激烈的方程。你得到较小的n - Aa吗?大型ka刚度会发生。
问候,
Wojciech
2的评论
杰里米Huard
2019年4月25日
编辑:杰里米Huard
2019年4月26日
嗨,丹尼尔,
定义运输舱的方程只适用于单个丸。
但您可以修改重复计量方程代替“时间”的时候从这里推荐的剂量:
https://dx.doi.org/10.1128%2FAAC.00461-07
在SimBiology您需要添加一个事件存储计量的时间。如果你量量是常数,这是所有你需要做的。如果你的剂量取决于模型参数数量,您可以修改剂量在这个事件。
我实现了这个在你的模型(见附件)。
编辑:这是错误的URL。我已经修正它。