Dear NMUSERS,
I am having problems simulating simultaneous zero and first order absorption.
Here is my model and data file. Every time it only reads in the second dose
record and simulates based on that record. I have tried this on three different
computers running NONMEM. CLI, VI and KAI are individual nonmem PK parameter
estimates for subjects from sparse data using an established population model
setting MAXEVAL=0. I am trying to simulate a full profile for each subject
based on their PK parameters.
$PROB SIMULATE
$INPUT C ID CMT AMT EVID TIME RATE II SS CLI VI KAI DV
$DATA sim.csv IGNORE=C
$SUBROUTINES ADVAN4 TRANS4
$PK
CL = CLI
V2 = VI
Q = 0.02
V3 = 0.9
KA = KAI
D2 = 6
F1=0.7
F2=0.3
ALAG1=0.125
ALAG2=0.125
$ERROR CALLFL=0
IPRED=0
IF(CMT.EQ.2) IPRED=LOG((A(2)+0.000001)/V2)+ERR(1)
Y = IPRED
IRES = DV - IPRED
IWRES=IRES ;W=1 for additive error model
$OMEGA
0 FIXED
$SIGMA
0 FIXED
$SIMULATION ONLY (12345)
$TABLE ID TIME EVID IPRED
NOPRINT FILE=NCA
C ID CMT AMT EVID TIME RATE II SS CLI
VI KAI DV
. 1 1 1000 1 0 0 24 1 0.05
0.3 0.05 .
. 1 2 1000 1 0 -2 24 1 0.05
0.3 0.05 .
. 1 2 . 0 0.01 . . . 0.05
0.3 0.05 .
. 1 2 . 0 0.5 . . . 0.05
0.3 0.05 .
. 1 2 . 0 1 . . . 0.05
0.3 0.05 .
. 1 2 . 0 1.5 . . . 0.05
0.3 0.05 .
. 1 2 . 0 2 . . . 0.05
0.3 0.05 .
. 1 2 . 0 2.5 . . . 0.05
0.3 0.05 .
. 1 2 . 0 3 . . . 0.05
0.3 0.05 .
. 1 2 . 0 4 . . . 0.05
0.3 0.05 .
. 1 2 . 0 5 . . . 0.05
0.3 0.05 .
. 1 2 . 0 6 . . . 0.05
0.3 0.05 .
. 1 2 . 0 7 . . . 0.05
0.3 0.05 .
. 1 2 . 0 8 . . . 0.05
0.3 0.05 .
. 1 2 . 0 9 . . . 0.05
0.3 0.05 .
. 1 2 . 0 10 . . . 0.05
0.3 0.05 .
. 1 2 . 0 11 . . . 0.05
0.3 0.05 .
. 1 2 . 0 12 . . . 0.05
0.3 0.05 .
. 1 2 . 0 14 . . . 0.05
0.3 0.05 .
. 1 2 . 0 16 . . . 0.05
0.3 0.05 .
. 1 2 . 0 18 . . . 0.05
0.3 0.05 .
. 1 2 . 0 20 . . . 0.05
0.3 0.05 .
. 1 2 . 0 22 . . . 0.05
0.3 0.05 .
. 1 2 . 0 24 . . . 0.05
0.3 0.05 .
Any help would be appreciated.
Regards,
Justus
