Re: [NMusers] Nonparametric Modelling

From: Julie ANTIC <j.antic_at_envt.fr>
Date: Wed, 4 Jul 2007 15:48:48 +0200

The nonparametric option will only estimate the distribution of ETA. The =
fixed parameters of your model (THETA(1), THETA(2),...,THETA(8)) should =
be estimated before the non parametric step with a parametric procedure =
($EST...).
By default, if the (parametric) estimation ($EST) step failed, the =
nonparametric step is not run. I think, it is your problem because the =
UNCONDITIONAL option just permits to run the nonparametric step when the =
$EST step has failed.

If I were you, I would try to run the command file without the =
$NONPARAMETRIC option, and correct eventuel errors...

I hope it'll help, good luck,
Julie Antic
  ----- Original Message -----
  From: Rik Schoemaker
  To: nmusers_at_globomaxnm.com
  Sent: Wednesday, July 04, 2007 2:30 PM
  Subject: [NMusers] Nonparametric Modelling


  Dear all,

  I'm trying to implement nonparametric population parameter estimation =
in NONMEM VI but I run into the following error message:

  NONPARAMETRIC ESTIMATE REQUESTED. BUT NO INDIVIDUAL RECORDS EXIST =
WITH DATA MODELED WITH ACTIVE ETA'S

  As a result no tables are output at all. However, if I leave out the
  $NONP UNCONDITIONAL
  the model runs fine and I do get my tables...

  Any clues?

  Thanks in advance!

  Rik Schoemaker


  $PROB PK C1
  $INPUT ID TRT TIME NDV DV CMT AMT SS RATE II MDV WT
  $DATA DRUG4.NM
  $SUBROUTINES ADVAN4 TRANS4
  $PK
  TABS = THETA(1)*EXP(ETA(1))
  KA = 0.692/TABS
  S2 = THETA(2)*EXP(ETA(2))
  V2 = S2
  V3 = THETA(3)*EXP(ETA(3))
  CL = THETA(4)*EXP(ETA(4))
  Q = THETA(5)*EXP(ETA(5))
  ALAG1 = THETA(6)*EXP(ETA(6))
  F1 = THETA(7)*EXP(ETA(7))
  $ERROR
  IPRED = LOG(.025)
  W = THETA(8)
  IF(F.GT.0) IPRED = LOG(F)
  IRES = IPRED-DV
  IWRES = IRES/W
  Y = IPRED+ERR(1)*W
  $THETA 1 15 60 12 15 0.6 0.45 0.3
  $OMEGA .1 .1 .1 .1 .1 .1 .1
  $SIGMA 1 FIXED
  $EST PRINT=1 MAXEVAL 9999 POSTHOC NOABORT METHOD=1
  MSFO=msfb1
  $NONP UNCONDITIONAL
  $COV PRINT=E
  $TABLE ID TIME NDV DV CMT AMT SS RATE II MDV WT
  IPRED IRES IWRES FILE = SDTAB1
  ONEHEADER NOPRINT
  $TABLE ID ALAG1 TABS S2 V3 CL Q F1 FILE=PATAB1
  ONEHEADER NOPRINT


  First few data lines:

  1 0 0 . . 1 175000 1 . 12 1 58
  1 0 0.05 . . 2 . . . . 1 58
  1 0 0.2166666667 . . 2 . . . . 1 58
  1 0 0.3833333333 . . 2 . . . . 1 58
  1 0 0.55 . . 2 . . . . 1 58
  1 0 0.7166666667 . . 2 . . . . 1 58
  1 0 0.8833333333 . . 2 . . . . 1 58
  1 0 1 1206 7.0950643773 2 . . . . 0 58
  1 0 1.05 . . 2 . . . . 1 58
  1 0 1.2166666667 . . 2 . . . . 1 58
  1 0 1.3833333333 . . 2 . . . . 1 58
  1 0 1.55 . . 2 . . . . 1 58
  1 0 1.7166666667 . . 2 . . . . 1 58
  1 0 1.8833333333 . . 2 . . . . 1 58
  1 0 2 1206 7.0950643773 2 . . . . 0 58
  1 0 2.05 . . 2 . . . . 1 58
  1 0 2.2166666667 . . 2 . . . . 1 58
  1 0 2.3833333333 . . 2 . . . . 1 58
  1 0 2.55 . . 2 . . . . 1 58
  1 0 2.7166666667 . . 2 . . . . 1 58
  1 0 2.8833333333 . . 2 . . . . 1 58
  1 0 3 975 6.882437471 2 . . . . 0 58


Received on Wed Jul 04 2007 - 09:48:48 EDT

This archive was generated by hypermail 2.2.0 : Tue Nov 06 2007 - 15:06:24 EST