NONMEM Users Network Archive

Hosted by Cognigen

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

From: Andrew Cristinaccce <andy_at_wright-dose.com>
Date: Wed, 23 Oct 2019 13:30:50 +0100

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*
>
> *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:30:50 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_at_globomaxnm.com.