Re: Simulate simultaneous zero and first order absorption
Dear Justus,
Quoted reply history
On the second dosing record, code SS = 2 for superposition of a second SS
dose. SS = 1 will ignore prior dosing records.
Best regards,
Marc
On Wed, Mar 5, 2014 at 9:51 PM, Justus Bingham <
[email protected]> wrote:
> 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
>
> CIDCMTAMTEVIDTIMERATEIISSCLIVIKAIDV.1110001002410.050.30.05..12100010-2241
> 0.050.30.05..12.00.01...0.050.30.05..12.00.5...0.050.30.05..12.01...0.05
> 0.30.05..12.01.5...0.050.30.05..12.02...0.050.30.05..12.02.5...0.050.30.05
> ..12.03...0.050.30.05..12.04...0.050.30.05..12.05...0.050.30.05..12.06...
> 0.050.30.05..12.07...0.050.30.05..12.08...0.050.30.05..12.09...0.050.30.05
> ..12.010...0.050.30.05..12.011...0.050.30.05..12.012...0.050.30.05..12.014
> ...0.050.30.05..12.016...0.050.30.05..12.018...0.050.30.05..12.020...0.05
> 0.30.05..12.022...0.050.30.05..12.024...0.050.30.05.
>
> Any help would be appreciated.
>
> Regards,
>
> Justus
>