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
Quoted reply history
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*
>
> *sumeet-singla_at_uiowa.edu <mailto:sumeet-singla_at_uiowa.edu>*
>
> *518.577.5881*
>
--
Andrew Cristinacce
Scientist, Wright Dose Ltd