From LIAO@cosmos.prius.jnj.com Tue Feb 13 08:51:00 1996
Subject: FORTRAN error
Hi, I am working on an PK/PD analysis using an indirect inhibition model. When I fixed the gamma to 1 in the Hill function, I was able to get somewhat reasonable fit. However, when I tried to estimate gamma (greater than 1) in order to fit the steeper recovery of the PD data, I got an error msg at the begining of the iteration-
run-time error M6202: MATH - log: SING error
I am running the DOUBLE PRECISION NONMEM VERSION IV LEVEL 2.1 on a PC. The control stream file is shown below. Any suggestion will be greatly appreciated.
=========================================================
$PROB A INDIRECT INHIBITION PK/PD MODEL
$INPUT ID AMT TIME DOSE PD=DV CL V KA K12 K21
$DATA XXXX4.DAT
$SUBROUTINES ADVAN8 TOL=5
$MODEL NPARAMETERS=20
COMP=(CENTRAL)
$PK
CALLFL=1
KEL= CL/V
SUMR=K12+KEL+K21
DISC=DSQRT(SUMR*SUMR-4.0*K21*KEL)
ALPH=(SUMR+DISC)/2.0
BETA=(SUMR-DISC)/2.0
M0=DOSE*KA/V
M1=(K21-ALPH)/(KA-ALPH)/(BETA-ALPH)*M0
M2=(K21-BETA)/(KA-BETA)/(ALPH-BETA)*M0
M3=(K21-KA)/(ALPH-KA)/(BETA-KA)*M0
TR01=THETA(1)
IC50=THETA(2)
EMAX=THETA(3)
GAMMA=THETA(4)
R01=TR01*(1+ETA(1))
K10=R01/100
F1=100 ;BOUNDARY CONDITION FOR EQUATION 1
P(1)=KA
$DES
CP=M1*DEXP(-ALPH*T) + M2*DEXP(-BETA*T) + M3*DEXP(-P(1)*T)
CPG=CP**GAMMA
IC=EMAX*CPG/(IC50**GAMMA + CPG)
DADT(1)= -K10*A(1) + R01*(1 - IC)
$ERROR
YCAL=F
Y=F*(1+EPS(1))
$THETA (5, 30, 100) (50, 800, 3000) 1.0 FIXED (1.1, 3, 10)
$OMEGA .1
$SIGMA .1
$EST POSTHOC MAXEVAL=9999 PRINT=10 MSFO=MSF2
$COV
$TABLE ID TIME CP YCAL K10 IC50 R01 EMAX FILE=XXX1.OUT
NOHEADER NOPRINT