devel/pheno_singular.mod

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