用NONMEM实现IV和IG数据的同时分析

学过药动的人都知道生物利用度（F）是一个非常重要的参数，其反映了机体对药物的利用程度。Wikipedia上对生物利用度（Absolute bioavailability）的解释为：

In pharmacology, in order to determine absolute bioavailability of a drug, a pharmacokinetic study must be done to obtain a plasma drug concentration vs time plot for the drug after both intravenous (iv) and extravascular (non-intravenous, i.e., oral) administration. The absolute bioavailability is the dose-corrected area under curve (AUC) non-intravenous divided by AUC intravenous. For example, the formula for calculating F for a drug administered by the oral route (po) is given below.

那个这个生物利用度F是否可以认为是当药物口服的时候，进入体内的药物的百分比呢？即当生物利用度为10%，给药剂量为100mg时，那么进入体内的总药物量就是10mg。

后来又在百度百科上看到对生物利用度的解释：

IG: dose = 10000 , F = 0.5, CL = 50, V = 200, KA = 1
IV: dose = 1000, CL = 50, V = 200
RTE 1 = IV, RTE 2 = IG
注：模拟时phoenix不支持F参数的设定，把其给药剂量改成dose*F = 5000

第二步，IV和IG的模拟数据分别加上N(0,0.1)的噪音，总共模拟10个subject，IV和IG各5个。数据文件见下面附件(iv-oral-data)。

第三步，用Phoenix的房室模型计算IV和IG模拟数据的AUC。因为我们的数据就是来自房室模型，所以肯定比用NCA计算AUC更准确一些，计算结果如下：

第四步，用NONMEM建立同时拟合IV和IG数据的模型，然后对同一批数据进行分析。分别用ADVAN2口服一室模型（run100.mod）和ADVAN6自定义模型(run101.mod)实现：

总结：

实现了一个模型就可以同时分析IV和IG数据，与传统的AUC比值方法结果一致。但NONMEM这种方法的优势在于其可扩展性，即便药物在吸收过程中发生代谢，同时肝脏中也会代谢等复杂情况，仍可用这种方法计算生物利用度，而这时使用AUC比值的方法就不可靠了。