NONMEM Users Network Archive

Hosted by Cognigen

RE: [NMusers] BAYES method and initial estimate bounds

From: Bauer, Robert <Robert.Bauer_at_iconplc.com>
Date: Tue, 22 Dec 2015 17:40:56 +0000

Bill:
I just tested this on a simple script (see below), and the boundaries are imposed in BAYES analysis for me.

$PROB RUN# Example 1 (from samp5l)
$INPUT C SET ID JID TIME DV=CONC AMT=DOSE RATE EVID MDV CMT CLX V1X QX V2X SDIX SDSX
$DATA example1.csv IGNORE=C
$SUBROUTINES ADVAN3 TRANS4
$PK
MU_1=THETA(1)
MU_2=THETA(2)
MU_3=THETA(3)
MU_4=THETA(4)
IF (THETA(1)<1.63.OR.THETA(1)>1.8) WRITE(50,*) THETA(1)
CL=DEXP(THETA(1)+ETA(1))
V1=DEXP(THETA(2)+ETA(2))
Q=DEXP(MU_3+ETA(3))
V2=DEXP(MU_4+ETA(4))
S1=V1
$ERROR
Y = F + F*EPS(1)
$THETA (1.63, 1.67,1.8) 2.0 2.0 2.0
$OMEGA BLOCK(4) VALUES(0.15,0.01)
$SIGMA (0.6 )
$PRIOR NWPRI
$THETAP (2.0 FIX) (2.0 FIX) (2.0 FIX) (2.0 FIX)
$THETAPV BLOCK(4) FIX VALUES(10000,0.0)
$OMEGAP BLOCK(4) FIX VALUES(0.2,0.0)
$OMEGAPD (4 FIX)
$EST METHOD=BAYES INTERACTION NITER=1000 PRINT=100 CTYPE=3
$COV MATRIX=R PRINT=E UNCONDITIONAL

Robert J. Bauer, Ph.D.
Vice President, Pharmacometrics R&D
ICON Early Phase
Office: (215) 616-6428
Mobile: (925) 286-0769
Robert.Bauer_at_iconplc.com<mailto:Robert.Bauer_at_iconplc.com>
www.iconplc.com<http://www.iconplc.com>

From: owner-nmusers_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com] On Behalf Of Bill Gillespie
Sent: Tuesday, December 22, 2015 8:14 AM
To: nmusers_at_globomaxnm.com
Subject: [NMusers] BAYES method and initial estimate bounds

Hi all,

If bounds are specified in the $THETA statements for initial estimates, what, if any, effect do they have when the BAYES method is used?

I initially thought they might work in conjunction with the prior specifications to implement truncated priors. However that does not seem to be true because statements like "IF(ABS(THETA(1)) > 10) EXIT 1 14" sometimes get tripped even when THETA(1) is bounded between -10 and 10 in the $THETA statement.

Happy Holidays,
Bill
<br /><br />
ICON plc made the following annotations.
------------------------------------------------------------------------------
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,

ICON plc
South County Business Park
Leopardstown
Dublin 18
Ireland
Registered number: 145835

Received on Tue Dec 22 2015 - 12:40:56 EST

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request_at_iconplc.com. Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.