在使用 NONMEM 时注意确保scaling factor对应体积的正确 (尤其当ADVAN 6等自定义模型时)
记录一次NONMEM bug shooting:
场景: 复写文献里的一个 ** 2房室,Weibull吸收,一级代谢的 imatinib popPK。**(模型来自Park, G. J., Park, W. S., Bae, S., Park, S. M., Han, S., & Yim, D. S. (2016). Population pharmacokinetics of imatinib mesylate in healthy Korean subjects. Translational and Clinical Pharmacology, 24(2), 96–104.)
模型结构:
模型参数 (去掉Age covariate):
parame | 值 |
---|---|
CL (L/h) | 13.6 |
Vc (L) | 153 |
Vp (L) | 64 |
Q (L/h) | 8.64 |
K1 | 0.998 |
GAMMA | 2.24 |
模拟数据集中DOSE是mg,在无scaling factor情况下模型得到的浓度为 mg/L, 1mg/L=1000ng/mL
代码1 (错误示范)
细心的小伙伴一看代码就知道哪里错了
$PROB ;;2-comp Weibull abosorption
$DATA imatinib_sim_children_SS2_form_vpcPLOT.csv IGNORE=@
$INPUT ID,TIME,AMT,RATE=DROP,DV,ROUT