forrtl: error (73): floating divide by zero

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.

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

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.

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.