Re: Parallel first order and Michaelis-Menten elimination
ADVAN6 ADVAN8 or (nm7) ADVAN13
The code is below
Leonid
-------------------
$SUBROUTINE ADVAN6 TOL=9
$MODEL
NCOMP = 2
COMP = (CENTRAL) ;1
COMP = (PERIPH) ;2
$PK
CL= THETA(1)*EXP(ETA(1))
V1= THETA(2)*EXP(ETA(2))
Q = THETA(3)*EXP(ETA(3))
V2= THETA(4)*EXP(ETA(4))
VM= THETA(5)*EXP(ETA(5))
KM= THETA(6)
K10 = CL/V1
K12 = Q/V1
K21 = Q/V2
S1 = V1
S2 = V2
$DES
C1 = A(1)/S1
DADT(1) = -K10*A(1)-A(1)*VM/(KM+C1)-K12*A(1)+K21*A(2)
DADT(2) = K12*A(1)-K21*A(2)
$ERROR
TY = A(1)/V1
IPRED=TY
W = SQRT(THETA(7)**2/TY**2+THETA(8)**2)
Y = IPRED*EXP(W*ERR(1))
$THETA
.....
$OMEGA
.....
$SIGMA
1 FIX ; ~ERR
--------------------------------------
Leonid Gibiansky, Ph.D.
President, QuantPharm LLC
web: www.quantpharm.com
e-mail: LGibiansky at quantpharm.com
tel: (301) 767 5566
[email protected] wrote:
> Dear All,
>
> I am working with a Biologic and would like to have a PK model with parallel first order and Michaelis-Menten elimination. Any suggestion about which subroutine I am supposed to use? If you can share an example for the control stream, that will be a great help.
>
> Thanks,
>
> Yuhong