devel/pheno_reduced.mod

   1 $PROBLEM  PHENOBARB SIMPLE MODEL
   2 $INPUT  ID TIME AMT DV
   3 $DATA  pheno_reduced.dta IGNORE=I
   4 $SUBROUTINE  ADVAN1 TRANS2
   5 $PK
   6
   7
   8       TVCL=THETA(1)
   9       TVV=THETA(2)
  10       CL=TVCL*EXP(ETA(1))
  11       V=TVV*EXP(ETA(2))
  12       S1=V
  13 $ERROR
  14
  15
  16       W=F
  17       Y=F+W*EPS(1)
  18
  19       IPRED=F         ;  individual-specific prediction
  20       IRES=DV-IPRED   ;  individual-specific residual
  21       IWRES=IRES/W    ;  individual-specific weighted residual
  22 $THETA  (0,0.0105)    ; CL
  23
  24 $THETA  (0,1.0500)    ; V
  25
  26 $OMEGA            .4  ; IVCL
  27                   .25 ; IVV
  28
  29 $SIGMA            .04
  30
  31 $ESTIMATION  MAXEVAL=9999 SIGDIGITS=4 POSTHOC
  32 ;$COVARIANCE  PRINT=E
  33 $TABLE ID TVCL TVV NOPRINT ONEHEADER FILE=patab1