ADVAN5 code
From: "Nele Plock" nplock@zedat.fu-berlin.de
Subject: [NMusers] ADVAN5 code
Date: Thu, April 28, 2005 8:22 am
Dear all,
I have been trying to write a model file using ADVAN5. I have attached the
code below. The problem is I cant get it to run. I get the following error message:
PRED EXIT CODE = 1
0INDIVIDUAL NO. 1 ID=0.10000000E+01 (WITHIN-INDIVIDUAL) DATA REC NO. 22
THETA=
7.78E+00 6.00E+00 1.49E+01 4.90E+01 2.24E+02 3.52E+01 7.59E-01 1.26E+00 7.00E-01
OCCURS DURING SEARCH FOR ETA AT INITIAL VALUE, ETA=0
NUMERICAL DIFFICULTIES OBTAINING THE SOLUTION.
THE COEFFICIENT MATRIX IS ALGORITHMICALLY SINGULAR.
0PROGRAM TERMINATED BY OBJ
MESSAGE ISSUED FROM ESTIMATION STEP
AT INITIAL OBJ. FUNCTION EVALUATION
I have a dataset with the first dose given intravenously (OCC=1), then some
observations, and afterwards only po doses (OCC>1) that I want to model with
a zero order absorption. Data record 22 is the second of those po doses. From
earlier seperate modelling of the data I know that the iv-data can be best
described with a 2-CMT-model while the po-data is only supported by a
1-CMT-model. Can anyone tell me if there is still an error in the code?
Any help would be highly appreciated!
Thanks,
Nele
$SUBROUTINES ADVAN5
$MODEL
NCOMPS=2
COMP=(CENTRAL, DEFDOSE, DEFOBS)
COMP=(PERIPH)
$PK
CL=THETA(1)*EXP(ETA(1))
IF (OCC.GT.1) CL=THETA(2)*EXP(ETA(2))
V1=THETA(3)
IF (OCC.GT.1) V1=THETA(4)*EXP(ETA(3))
Q=THETA(5)*EXP(ETA(4))
V2=THETA(6)*EXP(ETA(5))
D1=THETA(7)*EXP(ETA(6))
ALAG1=0
IF (ID.EQ.3.AND.OCC.GT.1) ALAG1=THETA(8)
IF (ID.EQ.6.AND.OCC.GT.1) ALAG1=THETA(8)
IF (ID.EQ.9.AND.OCC.GT.1) ALAG1=THETA(8)
PHI=LOG(THETA(9)/(1-THETA(9)))
F1=EXP(PHI+ETA(7))/(1+EXP(PHI+ETA(7)))
S1=V1
K10=CL/V1
K12=Q/V1
IF(OCC.GT.1) K12=0
K21=Q/V2
IF(OCC.GT.1) K21=0
$ERROR
IPRED=F
DEL=0
IF (IPRED.EQ.0) DEL=0.0001
W=F
IRES=DV-IPRED
IWRES=IRES/(W+DEL)
Y=IPRED+W*EPS(1)+EPS(2)
;---------------------------------------INITIAL ESTIMATES-----------------------------------
$THETA 7.78 FIX ; 1 CLiv
$THETA (0,6) ; 2 CLpo
$THETA 14.9 FIX ; 3 V1iv
$THETA (0,49) ; 4 V1po
$THETA 224 FIX ; 5 Q
$THETA 35.2 FIX ; 6 V2
$THETA 0.759 FIX ; 7 D1
$THETA 1.26 FIX ; 8 ALAG1
$THETA (0,0.7,1) ; 9 F1
$OMEGA 0.132 FIX ; 1 IIV_CLiv
$OMEGA 0.1 ; 2 IIV_CLpo
$OMEGA 0.1 ; 3 IIV_V1po
$OMEGA 1.66 FIX ; 4 IIV_Q
$OMEGA 0.0451 FIX ; 5 IIV_V2
$OMEGA 0.71 FIX ; 6 IIV_D1
$OMEGA 0.1 ; 7 IIV_F1
$SIGMA 0.1 ; prop.error
$SIGMA 1 ; add. error
$ESTIMATION PRINT=5 METHOD=1 INTERACTION MSFO=msf MAXEVAL=9999 NOABORT
$COV
$TABLE ID TIME IPRED IWRES ONEHEADER NOPRINT FILE=sdtab
$TABLE ID ETA1 ETA2 ETA3 ETA4 ETA5 ETA6 ETA7 CL V1 Q V2 D1 ALAG1 F1 ONEHEADER
NOPRINT FILE=patab