RE: [NMusers] Nonparametric method

This appears to be a different problem. In one example, an error
message is issued which indictates that none of the subjects have an
active eta, which is clearly not the case. Also, behavior/results of the
NONPARAMETRIC option may change markedly based on 1) whether the last
individual in the data set has nonactive etas or not or 2) whether the
eta's are ordered such that all eta's in common make up the initial part
of the eta vector. A modified NP subroutine that eliminates these
problems is being reviewed and subjected to further testing.

Tom Ludden

-----Original Message-----
From: owner-nmusers_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com]
On Behalf Of Leonid Gibiansky
Sent: Tuesday, July 24, 2007 2:11 PM
To: Brian M. Sadler
Cc: nmusers_at_globomaxnm.com
Subject: Re: [NMusers] Nonparametric method

If it is indeed similar, one can try to solve is with the similar trick,
forcing all etas to work
for all subjects:

http://www.metrumrg.com/publications/aaps_2004_etabug.pdf

Leonid


Brian M. Sadler wrote:
>
>
> Dear Tom,
>
>
>
> How does it differ from the parametric problem of non-influential
> correlated etas addressed by Stuart in article 99nov142002.html
> (http://www.cognigencorp.com/nonmem/nm/99nov142002.html)?
>
>
>
> Cheers... Brian
>
>
>
> * From: * owner-nmusers_at_globomaxnm.com
> [mailto:owner-nmusers_at_globomaxnm.com] *On Behalf Of *Ludden, Thomas
(MYD)
> *Sent:* Tuesday, July 24, 2007 12:02 PM
> *To:* nmusers_at_globomaxnm.com
> *Subject:* [NMusers] Nonparametric method
>
>
>
> Dear nmusers,
>
>
>
> Since Rik Schoemaker's original e-mail a few weeks ago, the code
> associated with the NONPARAMETRIC option has been reviewed. There are

> problems (bugs) that may reveal themselves when the data set contains

> individuals with different active etas. For example, a data set with
> some individuals receiving only IV doses while others receive only
oral
> doses may exhibit this bug. For the individuals receiving only IV
> doses, the eta's associated with absorption-related parameter such
> as F1, ALAG1, KA will be nonactive. At this time it appears, but it
is
> not certain, that the bugs will only affect the resulting estimates
when
> the data set contains individuals with different active etas. Code
> changes to address these bugs are being reviewed and tested. Anyone
who
> has a control stream and data set for which they have applied the
> NONPARAMETRIC option is welcome to send these to me
(luddent_at_iconus.com
> <mailto:luddent_at_iconus.com>) for testing with the modified routine.
We
> hope to provide fixes in the near future.
>
>
>
> Tom Ludden
>
>
>
> * From: * owner-nmusers_at_globomaxnm.com
> [mailto:owner-nmusers_at_globomaxnm.com] *On Behalf Of *Rik Schoemaker
> *Sent:* Wednesday, July 04, 2007 11:02 AM
> *To:* nmusers_at_globomaxnm.com
> *Subject:* 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
> $NONPARAMETRIC 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 <mailto: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
>
>
________________________________________________________________________
_______________________________________________________
>
> This e-mail transmission may contain confidential or legally
privileged
> information that is intended only for the individual or entity named
in
> the e-mail address. If you
> are not the intended recipient, you are hereby notified that any
> disclosure, copying, distribution, or reliance upon the contents of
this
> e-mail is strictly prohibited.
> If you have received this e-mail transmission in error, please reply
to
> the sender, so that ICON plc can arrange for proper delivery and then
> please delete the message. Thank You.
>
________________________________________________________________________
_______________________________________________________
>
>
>
>
>
>



___________________________________________________________________________=
____________________________________________________
This e-mail transmission may contain confidential or legally privileged =
information that is intended only for the individual or entity named in the=
 =
e-mail address. If you
are not the intended recipient, you are hereby notified that any disclosure=
,=
 copying, distribution, or reliance upon the contents of this e-mail is =
strictly prohibited.
If you have received this e-mail transmission in error, please reply to the=
 =
sender, so that ICON plc can arrange for proper delivery and then
please delete the message. Thank You.
___________________________________________________________________________=
____________________________________________________



Received on Tue Jul 24 2007 - 15:19:37 EDT