$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