NONMEM Users Network Archive

Hosted by Cognigen

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

From: hoxu8507
Date: Wed, 03 Jun 2009 13:12:28 +1000

Dear Leonid,

It works. Many thanks for your quick response.

Hongmei

Quoting Leonid Gibiansky <LGibiansky

> 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.
>>
>>
>>
>



----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.

Received on Tue Jun 02 2009 - 23:12:28 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.