NONMEM Users Network Archive

Hosted by Cognigen

[NMusers] RE: Link different thetas to same omega using mu referencing

From: Michael Fossler <mfossler_at_trevenainc.com>
Date: Wed, 13 Apr 2016 11:56:20 +0000

This looks like it would work - it seems to obey the rules around mu-refere=
ncing. You can simplify the $OMEGA by using SAME(3)


Mike
Michael J. Fossler, Pharm. D., Ph. D., F.C.P.
VP, Quantitative Sciences
Trevena, Inc
mfossler_at_trevenainc.com<mailto:mfossler_at_trevenainc.com>
Office: 610-354-8840, ext. 249
Cell: 610-329-6636

From: owner-nmusers_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com] On=
 Behalf Of Leander, Jacob
Sent: Wednesday, April 13, 2016 7:01 AM
To: nmusers_at_globomaxnm.com
Subject: [NMusers] Link different thetas to same omega using mu referencing

Hi

Consider the case where we have a categorical (potentially time-varying) co=
variate named CAT that can take several values (e.g. 0,1,2,3).
Let's say I want to model CL depending on this categorical covariate (where=
 parameters are estimated on log-scale).
A simple solution to this would be

IF(CAT.EQ.1) TVCL = THETA(1)
IF(CAT.EQ.2) TVCL = THETA(2)
IF(CAT.EQ.3) TVCL = THETA(3)

I let CL be described by a log-normal distribution with the same omega elem=
ent for all CAT values.

CL = EXP(TVCL + ETA(1))

With an omega statement:

$OMEGA 0.1

I now want to switch to the mu referencing framework in NONMEM. Each theta =
needs to be reference to a MU value, and to force it to be described by the=
 same eta I apply the following approach.

MU_1 = THETA(1)
MU_2 = THETA(2)
MU_3 = THETA(3)

IF(CAT.EQ.1) CL = EXP(MU_1 + ETA(1))
IF(CAT.EQ.2) CL = EXP(MU_2 + ETA(2))
IF(CAT.EQ.3) CL = EXP(MU_3 + ETA(3))

With an omega statement:

$OMEGA BLOCK(1) 0.1
$OMEGA BLOCK(1) SAME
$OMEGA BLOCK(1) SAME

where I force the omega element corresponding to ETA1, ETA2 and ETA3 to be=
 the same.

Is this a valid approach to obtain what I need?
Or is there a simpler way in the mu-referencing framework to link different=
 thetas to the same omega element?

All the best,
Jacob

________________________________
Confidentiality Notice: This message is private and may contain confidentia=
l and proprietary information. If you have received this message in error, =
please notify us and remove it from your system and note that you must not =
copy, distribute or take any action in reliance on it. Any unauthorized use=
 or disclosure of the contents of this message is not permitted and may be =
unlawful.

________________________________

Notice: This e-mail message, together with any attachments, contains inform=
ation of Trevena, Inc., 1018 West 8th Avenue, King of Prussia, PA 19406, US=
A. This information may be confidential, proprietary, copyrighted and/or le=
gally privileged.
It is intended solely for use by the individual or entity named on this mes=
sage. If you are not the intended recipient, and have received this message=
 in error, please notify us immediately and delete it and any attachments f=
rom your system.

Received on Wed Apr 13 2016 - 07:56:20 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@globomaxnm.com.