PK bolus intra-veineux, bibliothèque ADVAN (1, 2 et 3 cpts)

1 cpt

$PROB PK 1 cpt voie intraveineuse
$INPUT ID TIME DV AMT AGE WT CREA
$DATA ..\data.csv IGNORE=#
$SUBR ADVAN1 TRAN2
$PK

CL=THETA(1) *EXP(ETA(1))
V =THETA(2) *EXP(ETA(2))
K=CL/V
S1=V

$THETA (0,0.3) ; POPCL
$THETA (0,5) ; POPV

$OMEGA 0.2 ; BSVCL
$OMEGA 0.1 ; BSVV

$SIGMA 0.05 ; EPS1

$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

$ESTIMATION METHOD=1 FOCE PRINT=5 INTERACTION SIGDIGIT=3 MAXEVAL=2000 POSTHOC
$COVARIANCE
$TABLE ID AMT DV IPRED IRES IWRES ONEHEADER
NOPRINT FILE=base.fit
$TABLE ID CL V AGE WT CREA ONEHEADER
NOPRINT FILE=covariables.fit


2 cpt

$PROB PK 2 cpts voie intraveineuse
$INPUT ID TIME DV AMT AGE WT CREA
$DATA ..\data.csv IGNORE=#
$SUBROUTINE ADVAN3 TRANS4
$PK

CL =THETA(1) *EXP(ETA(1))
V1 =THETA(2) *EXP(ETA(2))
V2 =THETA(3) *EXP(ETA(3))
Q =THETA(4) *EXP(ETA(4))
S1=V1

$THETA
(0,0.5) ; CL
(0,10) ; V1
(0,8) ; V2
(0,8.4) ; Q

$OMEGA .1 ; ETACL
$OMEGA .1 ; ETAV1
$OMEGA .1 ; ETAV2
$OMEGA .1 ; ETAQ

$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
$SIGMA .05 ; PROP
$ESTIMATION METHOD=1 FOCE PRINT=5 INTERACTION SIGDIGIT=3 MAXEVAL=2000 POSTHOC
$COVARIANCE
$TABLE ID AMT DV IPRED IRES IWRES ONEHEADER
NOPRINT FILE=base.fit
$TABLE ID CL V1 V2 Q AGE WT CREA ONEHEADER
NOPRINT FILE=covariables.fit


3 cpt

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

CL = THETA(1) * EXP(ETA(1))
V1 = THETA(2) * EXP(ETA(2))
Q2 = THETA(3) * EXP(ETA(3))
V2 = THETA(4)
Q3 = THETA(5)
V3 = THETA(6)

S1 = V1
K = CL/V1

$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) ; V1
(0, 0.5) ; Q2
(0, 30) ; V2
(0, 30) ; Q3
(0, 30) ; V3

$OMEGA 0.5; OM_CL
$OMEGA 0.1; OM_V1
$OMEGA 0.1; OM_Q2

$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 CL V1 Q3 V3 Q4 V4 AGE WT CREA ETA1 ETA2 ETA3
FIRSTONLY NOPRINT NOAPPEND
FILE=parameter.fit