Thanks everyone for such a quick response. Changing CMT from 3 to 2 in my dataset was enough. Thank you so much!
Regards,
Sumeet Singla
Quoted reply history
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 PREDICTION 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 heart 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 distribution 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
IPREDF
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 MAXEVAL99 SIG=3 PRINT=5 POSTHOC NOABORT MSFO=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