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
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
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*
*[email protected] <mailto:[email protected]>*
*518.577.5881*
--
Andrew Cristinacce
Scientist, Wright Dose Ltd
