From: Singla, Sumeet K <*sumeet-singla*>

Date: Wed, 23 Oct 2019 12:48:38 +0000

Thanks everyone for such a quick response. Changing CMT from 3 to 2 in my d=

ataset was enough. Thank you so much!

Regards,

Sumeet Singla

From: Andrew Cristinaccce <andy

Sent: Wednesday, October 23, 2019 8:31 AM

To: Singla, Sumeet K <sumeet-singla

Subject: [External] Re: [NMusers] COMPARTMENT ASSOCIATED WITH THE PREDICTIO=

N 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 hear=

t 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 distributi=

on 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 MSF=

O=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

518.577.5881

--

Andrew Cristinacce

Scientist, Wright Dose Ltd

Received on Wed Oct 23 2019 - 08:48:38 EDT

Date: Wed, 23 Oct 2019 12:48:38 +0000

Thanks everyone for such a quick response. Changing CMT from 3 to 2 in my d=

ataset was enough. Thank you so much!

Regards,

Sumeet Singla

From: Andrew Cristinaccce <andy

Sent: Wednesday, October 23, 2019 8:31 AM

To: Singla, Sumeet K <sumeet-singla

Subject: [External] Re: [NMusers] COMPARTMENT ASSOCIATED WITH THE PREDICTIO=

N 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 hear=

t 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 distributi=

on 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 MSF=

O=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

518.577.5881

--

Andrew Cristinacce

Scientist, Wright Dose Ltd

Received on Wed Oct 23 2019 - 08:48:38 EDT