devel/pheno.mod

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