Error: Initial Values for ATOL / LSODI1 CODE 205 ?
From:"Kyun-Seop Bae"
Subject:[NMusers] Error: Initial Values for ATOL / LSODI1 CODE 205 ?
Date:Sun, February 10, 2002 8:30 pm
Dear Juan Jose and Sam,
Thank you very much for your immediate help.
CMT=4 was typographical error, and with the correction of this, the problem still
remains.
TOL=5 or 6 was also failed, I increased it.
Fixing eta as 0 also leads the same error.
For Sam's advice, I don't have your error message without estimation step.
But, similar to the previous one : 'ERROR IN LSODI1: CODE 205'
Maybe this is more fundamental cause than previous mail.
I found the fortran code in the file of 'LSODI1.for', but it is not easy to figure out.
(Just I found that state of KFLAG variable determines the error code.)
Sincerely,
Kyun-Seop Bae MD
Clinical Trial Center
Seoul National University Hospital
28 Yongon-dong Chongno-gu
Seoul, 110-744, Korea
Tel: +82-2-760-1666
Email: ksbae@snu.ac.kr
===========================
// OUTPUT FILE WITH ERROR MESSAGES //
0PRED EXIT CODE = 1
0INDIVIDUAL NO. 1 ID=0.10000000E+01 (WITHIN-INDIVIDUAL) DATA REC NO. 5
THETA= 5.00E-01 2.00E+00 1.00E-02
OCCURS DURING SEARCH FOR ETA AT INITIAL VALUE, ETA=0
ERROR IN LSODI1: CODE 205
ERROR OCURRED WHILE ATTEMPTING TO OBTAIN INITIAL VALUES FOR ATOL
0PROGRAM TERMINATED BY OBJ
MESSAGE ISSUED FROM ESTIMATION STEP
AT INITIAL OBJ. FUNCTION EVALUATION
1THERE ARE ERROR MESSAGES IN FILE PRDERR
Stop - Program terminated.
===========================
// CONTROL FILE //
$PROBLEM MODEL XXX
$INPUT ID TIME AMT DV CMT
$DATA A9_2.CSV IGNORE=#
$SUBROUTINE ADVAN9 TOL=10 ; TOL=5 or 6 shows the same error
$MODEL
NCOMP = 5 ;
COMP = (WATER) ; A(1), W
COMP = (ACH) ; A(2), X
COMP = (AC) ; A(3), YY
COMP = (HCL) ; A(4), Z
COMP = (H, DEFOBS, EQU) ; A(5), H
$PK
PKA = 5.5;
KA = 10 ** (-1 * PKA); ; Cause of error is too frequent use of
exponentials ?
VMAX = THETA(1)*EXP(ETA(1)) ; Error does not change with additive eta
KM = THETA(2);
KGE = THETA(3);
$AESINIT
INIT = 0; 0=APPROXIMATE, 1=EXACT
A(5) = 0.2;
$AES
E(5) = A(5)*(A(3)-A(4)+A(5))/(A(2)+A(4)-A(5)) - KA * A(1) ; Is this too complex?
$DES
PH = -1 * LOG (A(5)/A(1))
DADT(1) = -1 * KGE * A(1)
DADT(2) = -1 * KGE * A(2)
DADT(3) = -1 * KGE * A(3)
DADT(4) = VMAX * PH / (KM + PH) - KGE * A(4)
$ERROR
Y = F * EXP(EPS(1))
$THETA
(0, 0.5) ;THETA(2), VMAX - Is my initial guesses not adequate?
(0, 2) ;THETA(3), KM
(0, 0.01) ;THETA(4), KGE
$OMEGA
0.2
$SIGMA
0.02
$ESTIMATION METHOD=1 ; No matter with estimation method. the same error
message
===============================
// DATA FILE, CSV FORMAT IN ORIGINAL// -> TYPO CORRECTED(CMT 4 -> 5)
#ID TIME AMT DV CMT
1 0 0.9 . 1
1 0 0.9 . 2
1 0 0.11 . 3
1 0 0.4 . 4
1 1 . 5.37E-06 5
1 2 . 1.13E-05 5
1 3 . 1.43E-05 5
1 4 . 3.87E-06 5
1 5 . 4.92E-06 5
1 6 . 6.50E-06 5
1 7 . 9.16E-06 5
1 8 . 1.45E-05 5
1 9 . 3.08E-05 5
1 10 . 0.003155421 5
1 11 . 0.286965292 5
1 12 . 0.123910793 5
1 13 . 0.2 5
2 0 1 . 1
2 0 1 . 2
2 0 0.1 . 3
2 0 0.5 . 4
2 1 . 0.00001 5
2 2 . 0.00001 5
2 3 . 1.58E-05 5
2 4 . 3.16E-06 5
2 5 . 5.01E-06 5
2 6 . 6.31E-06 5
2 7 . 0.00001 5
2 8 . 1.58E-05 5
2 9 . 3.16E-05 5
2 10 . 0.003162278 5
2 11 . 0.316227766 5
2 12 . 0.1 5
2 13 . 0.199526231 4