NONMEM Users Network Archive

Hosted by Cognigen

RE: [External] Re: [NMusers] COMPARTMENT ASSOCIATED WITH THE PREDICTION IS OFF - CMT= 3

From: Singla, Sumeet K <sumeet-singla_at_uiowa.edu>
Date: Wed, 23 Oct 2019 12:48:38 +0000

Thanks everyone for such a quick response. Changing CMT from 3 to 2 in my d=
ataset was enough. Thank you so much!

Regards,
Sumeet Singla

From: Andrew Cristinaccce <andy_at_wright-dose.com>
Sent: Wednesday, October 23, 2019 8:31 AM
To: Singla, Sumeet K <sumeet-singla_at_uiowa.edu>; nmusers_at_globomaxnm.com
Subject: [External] Re: [NMusers] COMPARTMENT ASSOCIATED WITH THE PREDICTIO=
N IS OFF - CMT= 3


Hi Sumeet

You are using the wrong ADVAN. ADVAN3 is for a fixed 2 compartment PK model=
 with no heart rate compartment. There is no 3rd compartment in this model,=
 hence your error message.

You need to use an ADVAN with no fixed model, eg. 6 or 13, and define your =
compartments in the $MODEL and $DES blocks.

Best wishes

Andy
On 23/10/2019 13:01, Singla, Sumeet K wrote:

Hello everyone,

I am getting the following error while running a 2-compartment PK with Emax=
 model for heart rate for simultaneous estimation of parameters. Can anyone=
 please help me with that? Looks like CMT=3 which is associated with hear=
t rate is turned off. How to turn it back on? I have attached error, sample=
 of my dataset and control stream.
"
0DATA REC 3: COMPARTMENT ASSOCIATED WITH THE PREDICTION IS OFF: CMT= =
  3
0DATA REC 5: COMPARTMENT ASSOCIATED WITH THE PREDICTION IS OFF: CMT= =
  3
0DATA REC 7: COMPARTMENT ASSOCIATED WITH THE PREDICTION IS OFF: CMT= =
  3
0DATA REC 9: COMPARTMENT ASSOCIATED WITH THE PREDICTION IS OFF: CMT= =
  3"

ID
TIME
TIME1
AMT
DV
CMT
MDV
EVID
TYPE1
TYPE3
1402
0
9:39
40000
0
1
1
1
0
1
1402
0.183333
9:50
.
271.8
1
0
0
0
1
1402
0.183333
9:50
.
11
3
0
0
1
0
1402
0.483333
10:08
.
24.6
1
0
0
0
1
1402
0.483333
10:08
.
6
3
0
0
1
0



$INPUT C EVENT=DROP ID1=DROP ID TIME TIME1=DROP AMT DV CMT MDV EVID

$SUBROUTINE ADVAN3 TRANS4

$PK
TVV1 = THETA(1) ;Central Volume of distributi=
on in L
V1 = TVV1*EXP(ETA(1))
TVCL = THETA(2)
CL = TVCL*EXP(ETA(2)) ;Clearance
TVQ = THETA(3)
Q = TVQ*EXP(ETA(3)) ;Intercompartment Clearance
TVV2 = THETA(4)
V2 = TVV2*EXP(ETA(4)) ;Peripheral volume in L
S1=V1
EMAX = THETA(5)*EXP(ETA(5))
EC50 = THETA(6)*EXP(ETA(6))

$ERROR
CP = A(1)/S1
EFF= EMAX*CP/(EC50+CP)
IF (CMT.EQ.1) THEN
IPRED=CP
Y = IPRED*(1+ERR(1))
ELSE
IPRED=EFF
Y= IPRED +ERR(2)
ENDIF

$THETA
(0, 10); [V1]
(0, 50); [CL]
(0, 50); [Q]
(0, 100); [V2]
(50) ;
(0.01,100) ;

$OMEGA
0.16; [P] omega(1,1)
0.08; [P] omega(2,2)
0.20; [P] omega(3,3)
0.30; [P] omega(4,4)
(0, 0.2) ;BSVEMAX
(0.4) ; EC50

$SIGMA
0.40 ;ERR1
0.40 ;ERR2

$ESTIMATION METHOD=1 MAXEVAL=9999 SIG=3 PRINT=5 POSTHOC NOABORT MSF=
O=simultaneous.MSF

$TABLE ID TIME DV IPRED CP PRED NOPRINT ONEHEADER FILE=mytab001

Regards,

Sumeet K. Singla
Ph.D. Candidate
Division of Pharmaceutics and Translational Therapeutics
College of Pharmacy | University of Iowa
Iowa City, Iowa
sumeet-singla_at_uiowa.edu<mailto:sumeet-singla_at_uiowa.edu>
518.577.5881


--

Andrew Cristinacce

Scientist, Wright Dose Ltd


Received on Wed Oct 23 2019 - 08:48:38 EDT

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request_at_iconplc.com. Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.