Re: forrtl: error (73): floating divide by zero
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.