RE : Con't: fast and slow absorption
From: Sam Liao sliao@pharmaxresearch.com
Subject: RE : [NMusers] Con't: fast and slow absorption
Date: Wed, 15 Nov 2006 15:14:08 -0500
Jian:
Here is an example NONMEM control stream for the pk model you need. Please make sure to
specify CMT=3 for your PK data. You also need one dose event record with CMT=1 and the second
dose event record with CMT=2. The AMT will be the same in both records.
Best regards,
Sam Liao
Pharmax Research
=========================================================
$PROB PHARMACOKINETIC MODEL 1-CMPT WITH ONE SLOW AND ONE FAST ABSORPTION
$INPUT ID TIME DAY AMT CONC=DV CMT AGE WT SEX RACE
$DATA nm.prn IGNORE=#
$SUBROUTINE ADVAN6 TOL=6
$MODEL
NCOMP=3
COMP=(DEPOT1,DEFDOSE)
COMP=(DEPOT2)
COMP=(CENT,DEFOBS)
$PK
TVCL= THETA(1)
TVVD= THETA(2)
TKA1= THETA(3)
TKA2= THETA(4)
TF1 =THETA(5)
TLAG= THETA(6)
CL = TVCL*EXP(ETA(1))
V = TVVD*EXP(ETA(2))
KA1 = TKA1*EXP(ETA(3))
KA2 = TKA2*EXP(ETA(4))
F1 = TF1 *EXP(ETA(5))
IF (F1.GT.1) F1=0.999
F2 = 1 - F1
ALAG1=TLAG + ETA(6)
IF (ALAG1.LE.0) ALAG1=0
ALAG2 = ALAG1
K30 = CL/V
S3 = V
;ODE OF THE DECONJUGATION MODEL
$DES
DADT(1)= -KA1*A(1) ; KA1
DADT(2)= -KA2*A(2) ; KA2
DADT(3)= KA1*A(1) +KA2*A(2)-K30*A(3) ; CP
$ERROR
Y=F*(1+ERR(1)) + ERR(2)
IPRE= F
IRES= DV - IPRE
$THETA
(1,25,100) ;1 CL
(100,400,900) ;2 V
(0.1,.3,5) ;3 KA1
(0.001,.1,.3) ;4 KA2
(0.1,.5,1) ;5 F1
0 FIX ;6 TLAG
$OMEGA
0.1 ;1 VAR IN CL
0.1 ;2 VAR IN V
0.1 ;3 VAR IN KA1
0.1 ;4 VAR IN KA2
0.1 ;5 VAR IN F1
0. FIX ;6 VAR IN TLAG
$SIGMA 0.1 1.0
$EST SIG=3 MAXEVAL=9999 PRINT=5 METHOD=0 POSTHOC NOABORT
$COV
$TABLE ID TIME MDV K30 KA1 KA2 CL V ALAG1 F1 IPRE IRES
ONEHEADER NOPRINT FILE=fit