NONMEM Users Network Archive

Hosted by Cognigen

Re: forrtl: error (73): floating divide by zero

From: Leonid Gibiansky <LGibiansky>
Date: Tue, 02 Jun 2009 22:27:54 -0400

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




hoxu8507
> 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: hongmei
>
>
>
>
>
>
> ----------------------------------------------------------------
> This message was sent using IMP, the Internet Messaging Program.
>
>
>
Received on Tue Jun 02 2009 - 22:27:54 EDT

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.