Problem in Simulation
Dear NMUSERS,
I am trying to simulate ANC using Transit compartment Neutrapenia model. I am
simulating for different schedule. I have used this for another drug and it
worked fine. But rightnow it gives me error as below. If someone can guide me
regarding this problem, it would be helpful.
PROBLEM NO.: 1 SUBPROBLEM NO.: 1
0PRED EXIT CODE = 10
INDIVIDUAL NO. 61 ID= 6.10000000000000E+01 (WITHIN-INDIVIDUAL) DATA
REC NO. 207 THETA= 5.45E+00 1.35E+02 1.74E-01 2.14E-01 0.00E+00
DES SUBROUTINE: ERROR IN
COMPUTATION
ATTEMPT TO COMPUTE BASE**POWER WITH
BASE<0.
MESSAGE ISSUED FROM SIMULATION STEP==== END TIME ==== Fri 06/01/2012 02:09 PM
Please find below control stream. If anyone can take a look and advice on that
it would be helpful.
$PROB ANC Predictions
$INPUT C=DROP ID EVID CMT AMT RATE DV TRET TIME DAYS
$DATA all_schedules_grid.csv IGN=C
$SUBS ADVAN8 TOL=5
$MODEL
COMP=(STEM) ;1
COMP=(TRANSIT1) ;2
COMP=(TRANSIT2) ;3
COMP=(TRANSIT3) ;4
COMP=(TRANSIT4) ;5
COMP=(Central) ;6
COMP=(Peri) ;7
$PK
"FIRST" USE PRCOM_INT,ONLY:IMAX
" MAIN
" IMAX=9900000000
CL = 1.14
V1 = 6
K12 = 0.14
K21 = 0.06
K10 = CL/V1
CIRC0 = THETA(1)*EXP(ETA(1))
MTT = THETA(2)*EXP(ETA(2))
KTR = 4/MTT
GAM = THETA(3)
SLOPE = THETA(4)*EXP(ETA(3))
A_0(1) = CIRC0
A_0(2) = CIRC0
A_0(3) = CIRC0
A_0(4) = CIRC0
A_0(5) = CIRC0
$DES
EDRUG = 0
CP=A(6)/V1
IF(TRET.EQ.1.OR.TRET.EQ.3) EDRUG = SLOPE*CP
DRUG = EDRUG
DADT(1) = -KTR*A(1) + KTR*A(1)*(1-DRUG)*((CIRC0/A(5))**GAM)
DADT(2) = -KTR*A(2) + KTR*A(1)
DADT(3) = -KTR*A(3) + KTR*A(2)
DADT(4) = -KTR*A(4) + KTR*A(3)
DADT(5) = -KTR*A(5) + KTR*A(4) ; CIRCULATING CELLS
DADT(6) = -K10*A(6) -K12*A(6) +K21*A(7)
DADT(7) = -K21*A(7) +K12*A(6)
$ERROR
IPRED = 0.0001
IF(A(5).GT.0) IPRED = A(5)
W = THETA(5)*IPRED
IRES = DV - IPRED
IWRES = IRES/W
Y = IPRED+W*EPS(1) REP=IREP
AA9 = A(6)/V1
$THETA 5.45 FIX ;1 BASE
$THETA 135 FIX ;2 MTT (h)
$THETA 0.174 FIX ;3 POWER
$THETA 0.2141 FIX ;4 SLOPE
$THETA 0 FIX; (.554 ) ;7 Res err
$OMEGA 0.168 FIX ;1 IIV CIRC0
$OMEGA 0.0256 FIX ;2 IIV MTT
$OMEGA 0.36 FIX ;3 IIV SLOPE
$SIGMA 0 FIX ;1 FIX
$SIM (889215690) ONLYSIM SUBPROBLEM = 5
$TABLE REP ID TIME TRET DAYS CP EVID IPRED CIRC0 MTT GAM SLOPE NOPRINT NOHEADER
FILE=ResultANC.tab
Schedule file:
C 24.7143mg_1on-0off
C ID Evid CMT AMT RATE DV TRET TIME DAYS
1 0 5 0 0 . 3 0 0
1 1 6 24.71 24.71 . 3 0 0
1 0 5 0 0 . 3 12 0.5
1 0 5 0 0 . 3 24 1
1 1 6 24.71 24.71 . 3 24 1
1 0 5 0 0 . 3 36 1.5
Mayank PatelPh.D Student,
Long Island University, NY