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