I would try
IF(W.LE.0.000001) W = 1
Leonid
--------------------------------------
Leonid Gibiansky, Ph.D.
President, QuantPharm LLC
web: www.quantpharm.com
e-mail: LGibiansky at quantpharm.com
tel: (301) 767 5566
[email protected] wrote:
Dear Nonmem users,
I have been running a POPPK model with sparse sampling. This drug is
given intravenously and the half life is about 15 mins.
NONMEM VI run gives the following error:
forrtl: error (73): floating divide by zero
Image PC Routine
Line Source
nonmem1.exe 0040141C Unkown Unkown
Unkown
.......
Estimation omitted
The control stream is as following:
$SUBR ADVAN1 TRANS2
$PK
TVCL=THETA(1)
TVV=THETA(2)
CL = TVCL*EXP(ETA(1)) ;CL
V = TVV*EXP(ETA(2)) ;V
S1 = V ;Scaling factor
$ERROR
Y=F*EXP(ERR(1))+ERR(2) ;EXPONENTIAL and additive residual
IPRED=F
IRES = DV - F
W = F
IF(W.EQ.0) W = 1
IWRES = IRES/W
$THETA
(0.1, 50) ;POPCL
(0.01, 13.89) ;POPV
$OMEGA
0.36 ;BSVCL
0.25 ;BSVV
$SIGMA
0.49 ; EXPONENTIAL portion
0.123 ; Additive portion
$ESTIMATION MAX=9990 SIG=3 PRINT=1 NOABORT POSTHOC
METHOD=COND INTER
$COVARIANCE MATRIX=S
Any suggestions would be highly appreciated. Thanks.
Dr. Hongmei Xu
Faculty of Pharmacy (A15)
University of Sydney, NSW 2006
Email: [email protected]
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
