ROUNDING ERRORS (ERROR=134)
From: "Jerry Zhao"
Subject: [NMusers] ROUNDING ERRORS (ERROR=134)
Date:Tue, 13 May 2003 10:56:55 -0800
Dear All,
I am trying to simultaneously model parent and metabolite (control file attached to the
end of the message). I used log transformed data. But I kept getting Rounding Error
message:
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
0MINIMIZATION TERMINATED
DUE TO ROUNDING ERRORS (ERROR=134)
NO. OF FUNCTION EVALUATIONS USED: 958
NO. OF SIG. DIGITS UNREPORTABLE
ETABAR IS THE ARITHMETIC MEAN OF THE ETA-ESTIMATES,
AND THE P-VALUE IS GIVEN FOR THE NULL HYPOTHESIS THAT THE TRUE MEAN IS 0.
ETABAR: 0.13E-01 -0.62E-02 0.66E-02 -0.39E-02 0.49E-01 0.44E-01
P VAL.: 0.86E+00 0.93E+00 0.94E+00 0.96E+00 0.74E+00 0.67E+00
1
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Could someone please shed some light on what is causing this error and what I should do
to get rid of it? The model fit at the point of termination looks quite good, and the
parameter estimates are quite good too.
Thank you.
Jeri Sottos
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
$PROBLEM 2Comp Parent Drug 1Comp Metabolite Model -- Log Transformed DV
$DATA DATA_LOG.PRN
$INPUT ID AMT RATE TIME DV CMT EVID
$SUBROUTINES ADVAN6 TRANS1 TOL=5
$MODEL NPAR=6 NCOMP=3
COMP=(PARCEN DEFDOSE DEFOBS) ;CENTRAL COMPARTMENT FOR PARENT DRUG
COMP=(PARPER NODOSE) ;PERIPHERAL COMPARTMENT FOR PARENT DRUG
COMP=(METCEN NODOSE) ;CENTRAL COMPARTMENT FOR METABOLITE
$PK
CLt=THETA(1)*EXP(ETA(1)) ;PARENT DRUG CENTRAL CLEARANCE
Vc=THETA(2)*EXP(ETA(2)) ;PARENT DRUG CENTRAL COMPARTMENT VOLUME
CLd=THETA(3)*EXP(ETA(3)) ;PARENT DRUG DISTRIBUTION CLEARANCE
Vp=THETA(4)*EXP(ETA(4)) ;PARENT DRUG PERIPHERAL COMPARTMENT VOLUME
CLtm=THETA(5)*EXP(ETA(5)) ;METABOLITE CENTRAL CLEARANCE
Vcm=THETA(6)*EXP(ETA(6)) ;METABOLITE CENTRAL COMPARTMENT VOLUME
S1=Vc/1000
S3=Vcm/1000
$DES
DADT(1)=-(CLt+CLd)/Vc*A(1)+CLd/Vp*A(2)
DADT(2)=CLd/Vc*A(1)-CLd/Vp*A(2)
DADT(3)=-CLtm/Vcm*A(3)+.4045*CLt/Vc*A(1)
;40.45% OF PARENT DRUG IS COVERTED INTO METABOLITE
$ERROR
FLAG=0
IF(AMT.NE.0)FLAG=1 ;dosing records only
IPRED=LOG(F+FLAG) ;transform the prediction to the log of the
;prediction
; IPRED=log(f) for concentration records and
; IPRED=log(f+1) for dose records
R1=0
IF (CMT.EQ.1) R1=1
R2=0
IF (CMT.EQ.3) R2=1
Y1=IPRED+EPS(1)
Y3=IPRED+EPS(2)
Y=R1*Y1+R2*Y3
YORI=EXP(Y) ; y in normal scale
DVORI=EXP(DV) ; DV in normal scale
IRES=EXP(DV)-EXP(IPRED)
$THETA (10,33) (10,27) (30,60) (30,86) (10,17) (10,36)
$OMEGA BLOCK(6)
.24
.01 .28
.27 .01 .36
.20 .16 .25 .32
.24 .05 .24 .10 0.77
.17 .03 .09 .02 0.43 0.47
$SIGMA BLOCK (2)
.06
.09 .18
$ESTIMATION METHOD=1 INTERACTION PRINT=10 POSTHOC REPEAT MAXEVAL=9999
$TABLE ID TIME CMT IPRED IRES CLt Vc CLd Vp CLtm Vcm NOPRINT
FILE=mTBL.TXT
$COV
$SCAT DVORI VS TIME BY CMT
$SCAT YORI VS TIME BY CMT
$SCAT IRES VS TIME BY CMT
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^