Re: [NMusers] batches with different half-life - linear PK

Please see attached remarks and files.

----------
On Thu, 19 Apr 2007 14:54:11 +0000 (GMT), nonmem_at_optonline.net said:
> Hello NONMEM Users,
>
> Sometimes the same subjects receive the same drug from different batches.
> Potentially, half-life can be different accross batches. Theoretically,
> it is easy to model if PKis linear because NONMEM superimposes
> single-dose processes to obtain predicted concentrations for multiple
> doses. Does enyone know how to implement it in NONMEM? Are there any
> classes or publications which can help?
>
> Thanks!
> Pavel
--
  Alison Boeckmann
  alisonboeckmann_at_fastmail.fm

$PROB batch1
$DATA dat1.txt IGNORE=#
$INPUT ID TIME AMT RATE BAT CMT DV EVID PCMT
$BIND - - - - DOSE
$SUBR ADVAN1
$PK (NONEVENT)
IF (BAT.EQ.1) THEN
  F1=1
  D1=1
  ALAG1=1
ENDIF
IF (BAT.EQ.2) THEN
  F1=.5
  D1=.5
  ALAG1=.5
ENDIF
K=1
S1=1+ETA(1)
$ERROR Y=F+ERR(1)
$OMEGA 1
$SIGMA 1
$TABLE TIME AMT BAT CMT FILE=tab1.txt

$PROB batch2
$DATA dat2.txt IGNORE=#
$INPUT ID TIME AMT RATE BAT CMT DV EVID PCMT
$BIND - - - - DOSE
$SUBR ADVAN1
$PK (NONEVENT)
IF (BAT.EQ.1) THEN
  F1=1
  D1=1
  ALAG1=1
ENDIF
IF (BAT.EQ.2) THEN
  F1=.5
  D1=.5
  ALAG1=.5
ENDIF
K=1
S1=1+ETA(1)
$ERROR Y=F+ERR(1)
$OMEGA 1
$SIGMA 1
$TABLE TIME AMT BAT CMT FILE=tab2.txt

$PROB batch1+batch2
$DATA dat3.txt IGNORE=#
$INPUT ID TIME AMT RATE BAT CMT DV EVID PCMT
$BIND - - - - DOSE
$SUBR ADVAN1
$PK (NONEVENT)
IF (BAT.EQ.1) THEN
  F1=1
  D1=1
  ALAG1=1
ENDIF
IF (BAT.EQ.2) THEN
  F1=.5
  D1=.5
  ALAG1=.5
ENDIF
K=1
S1=1+ETA(1)
$ERROR Y=F+ERR(1)
$OMEGA 1
$SIGMA 1
$TABLE TIME AMT BAT CMT FILE=tab3.txt

$PROB batch1
$DATA datp.txt IGNORE=#
$INPUT ID TIME AMT RATE CMT DV EVID PCMT
$SUBR ADVAN6 TOL=6
$MODEL COMP=(B1) COMP=(B2)
$PK
 F1=1
 D1=1
 ALAG1=1
 F2=.5
 D2=.5
 ALAG2=.5
 K13=1
 K23=1
 S1=1+ETA(1)
 S2=1+ETA(1)
$ERROR
Y=A(1)/S1+A(2)/S2+ERR(1)
IF (PCMT.EQ.3) Y=A(3)
$DES
DADT(1)=-K13*A(1)
DADT(2)=-K23*A(2)
$OMEGA 1
$SIGMA 1
$TABLE TIME AMT CMT FILE=tabd.txt

$PROB batch1
$DATA datp.txt IGNORE=#
$INPUT ID TIME AMT RATE CMT DV EVID PCMT
$SUBR ADVAN7
$MODEL COMP=(B1) COMP=(B2)
$PK
 F1=1
 D1=1
 ALAG1=1
 F2=.5
 D2=.5
 ALAG2=.5
 K13=1
 K23=1
 S1=1+ETA(1)
 S2=1+ETA(1)
$ERROR
Y=A(1)/S1+A(2)/S2+ERR(1)
IF (PCMT.EQ.3) Y=A(3)
$OMEGA 1
$SIGMA 1
$TABLE TIME AMT CMT FILE=tabp.txt

# ID TIME AMT RATE BAT CMT DV EVID PCMT
1 0 1000 -2 1 1 . 1 .
1 0 . . 0 2 . 2 .
1 1 0 0 0 0 10 0 .
1 2 0 0 0 0 10 0 .
1 3 0 0 0 0 10 0 .
1 10 0 0 0 0 10 0 .
1 999 . . . 2 0 2 2

# ID TIME AMT RATE CMT DV EVID PCMT
1 0 1000 -2 1 . 1 .
1 0 1000 -2 2 . 1 .
1 0 . . 3 . 2 .
1 1 0 0 0 10 0 .
1 2 0 0 0 10 0 .
1 3 0 0 0 10 0 .
1 10 0 0 0 10 0 .
1 999 . . 3 0 2 3

# ID TIME AMT RATE BAT CMT DV EVID PCMT
1 0 1000 -2 2 1 . 1 .
1 0 . . 0 2 . 2 .
1 1 0 0 0 0 10 0 .
1 2 0 0 0 0 10 0 .
1 3 0 0 0 0 10 0 .
1 10 0 0 0 0 10 0 .
1 999 . . . 2 0 2 2

# ID TIME AMT RATE BAT CMT DV EVID PCMT
1 0 1000 -2 1 1 . 1 .
1 0 1000 -2 2 1 . 1 .
1 0 . . 0 2 . 2 .
1 1 0 0 0 0 10 0 .
1 2 0 0 0 0 10 0 .
1 3 0 0 0 0 10 0 .
1 10 0 0 0 0 10 0 .
1 999 . . . 2 0 2 2

# ID TIME AMT RATE CMT DV EVID PCMT
1 0 1000 -2 1 . 1 .
1 0 1000 -2 2 . 1 .
1 0 . . 3 . 2 .
1 1 0 0 0 10 0 .
1 2 0 0 0 10 0 .
1 3 0 0 0 10 0 .
1 10 0 0 0 10 0 .
1 999 . . 3 0 2 3

Received on Sat Apr 21 2007 - 13:51:40 EDT