NONMEM Users Network Archive

Hosted by Cognigen

Effect-compartment model

From: Singla, Sumeet K <sumeet-singla>
Date: Tue, 3 Dec 2019 03:49:48 +0000

Hi Everyone,

I am trying to relate my compound, THC and metabolite, THC-OH, concentratio=
ns to psychological highness (on scale of 0-10) using effect compartment ap=
proach. This effect compartment model is well validated in literature for t=
hese two compounds. When I run my model, I get the following error. I have =
tried ADVAN 6,8 AND 13 but no success so far. I have attached sample of my =
dataset and the control stream. Both THC and THC-OH has peripheral compartm=
ents. TYPE 3 equates to THC conc, TYPE2 equates to metabolite conc, and TYP=
E 1 equates to highness. CMT 1 and CMT4 denote THC and THCOH concentration =
in central compartment, respectively.

Please let me know of any possible solutions.

"OCCURS DURING SEARCH FOR ETA AT A NONZERO VALUE OF ETA
NUMERICAL DIFFICULTIES WITH INTEGRATION ROUTINE.
 NO. OF REQUIRED SIGNIFICANT DIGITS IN SOLUTION VECTOR
 TO DIFFERENTIAL EQUATIONS, 3, MAY BE TOO LARGE. "

          ID TIME =
    DV CMT TYPE1 TYPE2 TYPE3 =
    MDV
1204
9:38
21138
79.9
1
0
0
1
1
1204
9:38
.
7
             .
1
0
0
0
1204
9:38
.
10.1
4
0
1
0
0


$SUBROUTINE ADVAN6 TRANS1 TOL=3

$MODEL NCOMP = 6
COMP = (CENTRAL)
COMP = (PERIPH)
COMP = (EFFTHC)
COMP = (METABOL)
COMP = (METABOL2)
COMP= (EFFTHCOH)

$PK
Q1 =THETA(1)*EXP(ETA(1)) ; Inter-compartment clea=
rance for THC
V1 =THETA(2)*EXP(ETA(2)) ; V - THC - Central Com=
p
V2 =THETA(3)*EXP(ETA(3)) ; V - THC - Peripheral =
Comp
K14 =THETA(4)*EXP(ETA(4)) ; Metabolite conversion
Q2 =THETA(5)*EXP(ETA(5)) ; Inter-compartment cle=
arance for THC-OH
V4 =THETA(6)*EXP(ETA(6)) ; V - THCOH - Central C=
omp
V5 =THETA(7)*EXP(ETA(7)) ; V - THCOH - Periphera=
l Comp
CL =THETA(8)*EXP(ETA(8)) ; Clearance
KE01 =THETA(9)*EXP(ETA(9)) ; rate constant for THC e=
ffect compartment
KE02 =THETA(10)*EXP(ETA(10)) ; rate constant for THCOH =
effect compartment
S1=V1
S4=V4
A_0(1)=0
A_0(3)=0
A_0(4)=0
A_0(6)=0

$DES
C1=A(1)/V1
C3=A(3)
C4=A(4)/V4
C6=A(6)
DADT(1) = (Q1/V2)*A(2) - (Q1/V1)*A(1) - K14*A(1)
DADT(2) = (Q1/V1)*A(1) - (Q1/V2)*A(2)
DADT(3) = KE01*C1 - KE01*A(3)
DADT(4) = K14*A(1) - (CL/V4)*A(4) - (Q2/V4)*A(4) + (Q2/V5)*A(5)
DADT(5) = (Q2/V4)*A(4) - (Q2/V5)*A(5)
DADT(6) = KE02*C4 - KE02*A(6)

$ERROR
CE1=A(3)
CP=A(1)/V1
CE2=A(6)
CM=A(4)/V4
EMAX = THETA(11)*EXP(ETA(11))
EC50 = THETA(12)*EXP(ETA(12))
HILL = THETA(13)*EXP(ETA(13))
CONC = CP*(1 + ERR(1)) + ERR(2)
CONCMET = CM + ERR(3)
E = EMAX*(((CE1**HILL)+(CE2**HILL))/((EC50**HILL)+((CE1**HILL)+(CE2**HILL=
))))
EFF = E + ERR(4)
IF(TYPE3.EQ.1) IPRED = CP
IF(TYPE2.EQ.1) IPRED = CM
IF(TYPE1.EQ.1) IPRED = E
Y = (CONC*TYPE3) + (EFF*TYPE1) + (CONCMET*TYPE2)

$THETA
(0, 6.964) ;[Q1]
(0, 20) ;[V1]
(0, 50) ;[V2]
(0, 10) ;[K14]
(0, 10) ;[Q2]
(0, 40) ;[V4]
(0, 50) ;[V5]
(0, 102.95) ;[CL]
(0,5,50) ;[KE01]
(0,5,50) ;[KE02]
10 FIX ;[EMAX]
(1, 50, 400) ;[EC50]
(0.1, 1,5) ;[HILL]

$OMEGA
0.0152 ; [P] omega(1,1)
0.005 ; [P] omega(2,2)
0.06 ; [P] omega(3,3)
0.04 ; [P] omega(4,4)
0.03 ; [P] omega(5,5)
0.06 ; [P] omega(6,6)
0.03 ; [P] omega(7,7)
0.07 ; [P] omega(8,8)
0.06 ; [P] omega(9,9)
0.05 ; [P] omega(10,10)
0.05 ; [P] omega(11,11)
0.06 ; [P] omega(12,12)
0.06 ; [P] omega(13,13)

$SIGMA
0.5 ; [ERR1]
0.5 ; [ERR3]
0.5 ; [ERR4]

$SIGMA
0.00001 FIX ;[ERR2]


$EST METHOD=1 MAXEVAL=9999 PRINT=5 SIG=3 NOABORT MSFO=2010-2013.M=
SF
$COV

Regards,

Sumeet K. Singla
Ph.D. Candidate
Division of Pharmaceutics and Translational Therapeutics
College of Pharmacy | University of Iowa
Iowa City, Iowa
sumeet-singla
518.577.5881


Received on Mon Dec 02 2019 - 22:49:48 EST

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.