Concentration dependent volume of distribution
From: Johan.Rosenborg@astrazeneca.com
Subject: [NMusers] Concentration dependent volume of distribution
Date: Monday, 18 March 2002 17:47
All,
I would like to model concentration or possibly amount dependent volume of
distribution, using ADVAN8. Does anyone have a suggestion on how to deal
with this problem? I have tried to introduce VD, VMAX and A50 as shown
below, but it did not work.
/ Johan
$PROBLEM
$INPUT ID AMT SS II RATE TIME EVID CMT DV FLAG STUD
$DATA dat.prn IGNORE=#
$SUBROUTINES ADVAN8 TOL=4
$MODEL
COMP=(CENTRAL)
COMP=(PERIPHL)
COMP=(PERIPHC)
$PK
K10 =THETA(1)*(1+ETA(1))
K12 =THETA(2)
K21 =THETA(3)
K13 =THETA(4)*(1+ETA(2))
K31 =THETA(5)
;------------ Volume-------------------------
VD =THETA(6)*(1+ETA(3))
VMAX=THETA(9)*(1+ETA(5))
A50 =THETA(10)*(1+ETA(6))
;--------------------------------------------
F0 =THETA(7)*(1+ETA(4))
$DES
S1 =VD+VMAX*A(1)/(A50+A(1))
C1 =A(1)/S1
DADT(1)=K41*A(4)-(K10+K12+K13)*C1*S1+K21*A(2)+K31*A(3)
DADT(2)=K12*C1*S1-K21*A(2)
DADT(3)=K13*C1*S1-K31*A(3)
$ERROR
IF (EVID.EQ.0.AND.F.GT.0) THEN
IPRED=LOG(F)
ELSE
IPRED=LOG(F+0.001)
ENDIF
IRES=DV-IPRED
IWRES=IRES
X6=0
X7=0
IF(CMT.EQ.1)X6=1
IF(CMT.EQ.(-5))X7=1
Y=IPRED+X6*ERR(1)+X7*ERR(2)