Re: [NMusers] Nonparametric Modelling

Dear Julie,
 
Thank you for your comments! Unfortunately the parametric step runs just =
fine... :-)
 
Cheers,
 
Rik

>>> "Julie ANTIC" <j.antic_at_envt.fr> 04-07-07 3:48 PM >>>
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 $NONPARAMETR=
IC option, and correct eventuel errors...
 
I hope it'll help, good luck,
Julie Antic


----- Original Message -----
From: Rik Schoemaker ( mailto:RSchoemaker_at_chdr.nl )
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 - 11:01:43 EDT