RE: Enterohepatic Circulation!!!
From:VPIOTROV@PRDBE.jnj.com
Subject:RE: [NMusers] Enterohepatic Circulation!!!
Date:Mon, 10 Jun 2002 21:51:10 +0200
Johan,
The control used by Sreenivasa was proposed by Stu Beal in his Intermediate
level NONMEM workshop. I recommended to Sreenivasa to include one line of
code because it was in the original control. I did not test it myself. Below you
will find my own control stream which may help you with your problem. It needs
some tuning of course. The idea is to use 2 dummy compartments for switching
EHC on and off. The data set for simulations is attached to this mail.
Best regards,
Vladimir
$PROBLEM EHC
$INPUT ID TIME AMT DV CMT
$DATA nmd.txt
$SUB ADVAN9 TOL=5
$MODEL
COMP = (DEPOT,DEFDOSE) ; 1
COMP = (CENTRAL,DEFOBS) ; 2
COMP = (GALLBLD) ; 3
COMP = (DUMMY1) ; 4
COMP = (DUMMY2) ; 5
$PK
CL = THETA(1)
V = THETA(2)
CLBI = THETA(3) ; CL for bile excretion
KA = THETA(4)
KREL = THETA(5) ; K for release from gall bladder
ALAG4= THETA(6) ; Delay in release
ALAG5= ALAG4 + THETA(7) ; Release stop
K12 = KA
K20 = CL/V
K23 = CLBI/V
K31 = KREL
S2 = V
$DES
DADT(1) = -K12 * A(1) + K31 * A(4) * (1-A(5)) * A(3)
DADT(2) = K12 * A(1) - K20 * A(2) - K23 * A(2)
DADT(3) = K23 * A(2) - K31 * A(4) * (1-A(5)) * A(3)
DADT(4) = 0
DADT(5) = 0
$ERROR
Y = F*(1+ERR(1))
$THETA
(0 2) ; CL
(0 50) ; V
(0 2) ; CLBI
(0 5) ; KA
(100 FIX) ; KREL
(0 2) ; ALAG4
(0 .1) ; ALAG5
$OMEGA .001
$SIM (97834)
$EST
Data Set:
1 0 10 0 1
1 0 1 0 4
1 0 1 0 5
1 0.4 0 1 2
1 0.8 0 1 2
1 1.2 0 1 2
1 1.6 0 1 2
1 2 0 1 2
1 2.4 0 1 2
1 2.8 0 1 2
1 3.2 0 1 2
1 3.6 0 1 2
1 4 0 1 2
1 4.4 0 1 2
1 4.8 0 1 2
1 5.2 0 1 2
1 5.6 0 1 2
1 6 0 1 2
1 6.4 0 1 2
1 6.8 0 1 2
1 7.2 0 1 2
1 7.6 0 1 2
1 8 0 1 2