PK e

1 cpt

$PROB PK 1 cpt voie orale
$INPUT ID TIME DV AMT AGE WT CREA
$DATA ..\data.csv IGNORE=#
$SUBROUTINE ADVAN2 TRANS2 ; 1 cpt voie orale
$PK
CL =THETA(1) * EXP(ETA(1))
V =THETA(2) * EXP(ETA(2))
KA =THETA(3) * EXP(ETA(3))

K=CL/V ; REPARAMETRISATION
S2=V ; SCALE

$ERROR
DEL=0
IF (F.EQ.0) DEL=1
W=F+DEL
Y=F+W*EPS(1)
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W

$THETA (0, 5) ; CL
$THETA (0, 50) ; VC
$THETA (0, 1) ; KA

$OMEGA 0.5; OM_CL
$OMEGA 0.1; OM_VC
$OMEGA 0.1; OM_KA

$SIGMA 0.1; SIM_P

$EST METHOD=1 INT MAXEVAL=8000 PRINT=5 NOABORT POSTHOC NSIG=3
$COVARIANCE
$TABLE ID TIME DV PRED WRES RES IPRED IRES IWRES FILE=base.fit
$TABLE ID AMT KA CL V AGE WT CREA ETA1 ETA2 ETA3
FIRSTONLY NOPRINT NOAPPEND
FILE=parameter.fit

---

2 cpt:

[quote]$PROB PK 2 cpts voie orale
$INPUT ID TIME DV AMT AGE WT CREA
$DATA ..\data.csv IGNORE=#
$SUBROUTINE ADVAN2 TRANS2 ; 1 cpt voie orale
$PK
CL =THETA(1) * EXP(ETA(1))
V2 =THETA(2) * EXP(ETA(2))
Q =THETA(3) * EXP(ETA(3))
V3 =THETA(4) * EXP(ETA(4))
KA =THETA(5) * EXP(ETA(5))

K=CL/V ; REPARAMETRISATION
S2=V2 ; SCALE

$ERROR
DEL=0
IF (F.EQ.0) DEL=1
W=F+DEL
Y=F+W*EPS(1)
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W

$THETA (0, 5) ; CL
$THETA (0, 50) ; V2
$THETA (0, 1) ; Q
$THETA (0, 30) ; V3
$THETA (0, 1) ; KA

$OMEGA 0.5; OM_CL
$OMEGA 0.1; OM_V2
$OMEGA 0.1; OM_Q
$OMEGA 0.1; OM_V3
$OMEGA 0.1; OM_KA

$SIGMA 0.1; SIM_P

$EST METHOD=1 INT MAXEVAL=8000 PRINT=5 NOABORT POSTHOC NSIG=3
$COVARIANCE
$TABLE ID TIME DV PRED WRES RES IPRED IRES IWRES FILE=base.fit
$TABLE ID AMT KA CL V AGE WT CREA ETA1 ETA2 ETA3
FIRSTONLY NOPRINT NOAPPEND
FILE=parameter.fit

---

3 cpt

$PROB PK 3 cpts voie orale
$INPUT ID TIME DV AMT AGE WT CREA
$DATA ..\data.csv IGNORE=#
$SUBROUTINE ADVAN12 TRANS4
$PK

CL = THETA(1) * EXP(ETA(1))
V2 = THETA(2) * EXP(ETA(2))
KA = THETA(3) * EXP(ETA(3))
Q3 = THETA(4)
V3 = THETA(5)
Q4 = THETA(6)
V4 = THETA(7)

S2 = V2
K = CL/V2

$ERROR
DEL=0
IF (F.EQ.0) DEL=1
W=F+DEL
Y=F+W*EPS(1)
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W

$THETA
(0, 10) ; CL
(0, 30) ; V2
(0, 0.5) ; KA
(0, 30) ; Q3
(0, 30) ; V3
(0, 30) ; Q4
(0, 500) ; V4

$OMEGA 0.5; OM_CL
$OMEGA 0.1; OM_V2
$OMEGA 0.1; OM_KA

$SIGMA 0.5; SIM_P
$EST METHOD=1 INT MAXEVAL=8000 PRINT=5 NOABORT POSTHOC NSIG=3
$COVARIANCE
$TABLE ID TIME DV PRED WRES RES IPRED IRES IWRES FILE=base.fit
$TABLE ID AMT KA CL V2 Q3 V3 Q4 V4 AGE WT CREA ETA1 ETA2 ETA3
FIRSTONLY NOPRINT NOAPPEND
FILE=parameter.fit