NONMEM Users Network Archive

Hosted by Cognigen

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

From: Saeheum Song <ss.pkpdmodel_at_gmail.com>
Date: Wed, 23 Oct 2019 08:50:16 -0400

His model is 2 CMT model. Prediction is not compartment but parameter. So
changing data column name to measurement should be enough.

Kind regards

On Wed, Oct 23, 2019, 8:45 AM Andrew Cristinaccce <andy_at_wright-dose.com>
wrote:

> 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 you=
r
> 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 w=
ith
> 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 <sumeet-singla_at_uiowa.edu>*
>
> *518.577.5881*
>
>
>
> --
> Andrew Cristinacce
> Scientist, Wright Dose Ltd
>
>


Received on Wed Oct 23 2019 - 08:50:16 EDT

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com. Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.