$PROB PK 1 cpt voie orale, eq diff, erreur residuelle mixte $INPUT ID TIME AMT DV EVID CMT $DATA ..\data.csv IGNORE=# $SUBROUTINE ADVAN6 TRANS1 TOL=3 $MODEL NCOMP=2 COMP=(DIG, DEFDOSE) COMP=(CENTRAL, DEFOBS) $PK KA = THETA(1) * EXP(ETA(1)) CL = THETA(2) * EXP(ETA(2)) V2 = THETA(3) * EXP(ETA(3)) KE = CL/V2 S2 = V2 $DES DADT(1) = -KA*A(1) DADT(2) = KA*A(1) - KE*A(2) $ERROR IPRED=F W=(F*F*THETA(4)**2+THETA(5)**2)**0.5 Y=F+ERR(1)*W IRES=DV-IPRED IWRES=IRES/W $THETA (0, 2) ; KA (1, 30, 500) ; CL (5, 60, 500) ; V2 (0, .1) ; PROP (0, 500) ; ADD $OMEGA .1 ; ETA_KA .1 ; ETA_CL .1 ; ETA_V2 $SIGMA 1 FIXED $EST METHOD=1 INT MAXEVAL=8000 PRINT=5 NOABORT POSTHOC NSIG=3 $COVARIANCE $TABLE ID TIME DV PRED CMT WRES RES IPRED IWRES IRES FILE=1.fit $TABLE KA CL V2 FIRSTONLY NOPRINT FILE=parameter.fit