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

From: Andrew Cristinaccce <[email protected]> Sent: Wednesday, October 23, 2019 8:31 AM To: Singla, Sumeet K <[email protected]>; [email protected] 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 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