NONMEM Users Network Archive

Hosted by Cognigen

RE: [NMusers] Why should we avoid using micro rate constants?

From: STANDING, Joseph (GREAT ORMOND STREET HOSPITAL FOR CHILDREN NHS FOUNDATION TRU <joseph.standing_at_nhs.net>
Date: Tue, 5 Feb 2019 10:53:21 +0000

Dear Sumeet,

If you are assuming a distribution for your parameters (e.g. log-Normal p =
= theta * exp(eta)) then it might matter if you use rate constants versus=
 clearances and volumes. In general, if you want to make the log-Normal as=
sumption you should use clearances and volumes as there is reasonable biolo=
gical prior knowledge to show these generally follow a log-Normal distribut=
ion (do some reading on the occurrence of log-Normal distributions in biolo=
gy).

The rate constant is a ratio of two (usually) log-Normally distributed vari=
ables (e.g. k = CL/V) and hence may not necessarily be a shape that can i=
tself be described as a log-Normal. Here is some R-code that highlights th=
is:

# Simulate some realistic PK for a water soluble renally cleared drug
vd <- 40 * exp(rnorm(10000, sd = 0.5))
cl <- 6 * exp(rnorm(10000, sd = 0.5))
k <- cl / vd
# Visualise the histograms and use fitdistr function to
# fit a log-Normal
require(MASS)
# Volume:
hist(vd, freq = FALSE)
fit <- fitdistr(vd, "log-normal")$estimate
lines(dlnorm(0:max(vd), fit[1], fit[2]), lwd = 3)
# ...yes
#
# Clearance:
hist(cl, freq = FALSE)
fit <- fitdistr(cl, "log-normal")$estimate
lines(dlnorm(0:max(cl), fit[1], fit[2]), lwd = 3)
# ...yes
#
# K
hist(k, freq = FALSE)
fit <- fitdistr(k, "log-normal")$estimate
lines(dlnorm(0:max(k), fit[1], fit[2]), lwd = 3)
# ...no


People who do not like to make assumptions on distributions of parameters u=
se a nonparametric approach, and in this case it does not matter whether yo=
u use rate constants or clearances and volumes. However, unless you collec=
t rich informative data (to get good individual parameter estimates) and lo=
ts of it (to get a true idea of the distribution of parameters in the popul=
ation) it is usually advised to make a distributional assumption, and the l=
og-Normal is often sensible.

BW,

Joe




Joseph F Standing
MRC Fellow, UCL Institute of Child Health
Antimicrobial Pharmacist, Great Ormond Street Hospital
Honorary Senior Lecturer, St George's University of London
Tel: +44(0)207 905 2370
Mobile: +44(0)7970 572435
________________________________________
From: owner-nmusers_at_globomaxnm.com [owner-nmusers_at_globomaxnm.com] on behalf=
 of janet.wade_at_telia.com [janet.wade_at_telia.com]
Sent: 05 February 2019 06:51
To: 'Leonid Gibiansky'; 'Singla, Sumeet K'
Cc: nmusers_at_globomaxnm.com
Subject: RE: [NMusers] Why should we avoid using micro rate constants?

Hi All,

It could also be the statistical model. If you are estimating 4 parameters =
then different parameterisations should be fairly equivalent if a BLOCK(4) =
structure is used for both parameterisations. If only the diagonal option i=
s used, then this could be why different results/minimisations are obtained=
 for different parameterisations.

Kind regards,
Janet


Janet R Wade, PhD
Occams
Senior Consultant



From: owner-nmusers_at_globomaxnm.com <owner-nmusers_at_globomaxnm.com> On Behalf=
 Of Leonid Gibiansky
Sent: 04 February 2019 07:30
To: Singla, Sumeet K <sumeet-singla_at_uiowa.edu>
Cc: nmusers_at_globomaxnm.com
Subject: Re: [NMusers] Why should we avoid using micro rate constants?

It could be just coding error, could you show the control stream?
Thanks
Leonid

On Feb 3, 2019, at 12:44 PM, Singla, Sumeet K <sumeet-singla_at_uiowa.edu<mail=
to:sumeet-singla_at_uiowa.edu>> wrote:
Hello everyone!

I have a question. I was trying to build a 2-compartment PK model for marij=
uana use in occasional and chronic smokers. Initially, I was using providin=
g rate constants K12 and K21 ­in PK block and it resulted in poor fitting=
. Then, I later changed to CL,V1, V2 , Q and it resulted in proper fitting.=
 I was perplexed as to why I couldn’t get a proper fit by providing rate =
constants? I tried to look online but couldn’t find any proper explanatio=
n about when (or not) we should use micro constants in PK block to define o=
ur model in NONMEM? Does anyone has any useful insights into this?

Regards,
Sumeet Singla
Graduate Student
Dpt. of Pharmaceutics & Translational Therapeutics
College of Pharmacy- University of Iowa



***************************************************************************=
*****************************************

This message may contain confidential information. If you are not the inten=
ded recipient please inform the
sender that you have received the message in error before deleting it.
Please do not disclose, copy or distribute information in this e-mail or ta=
ke any action in relation to its contents. To do so is strictly prohibited =
and may be unlawful. Thank you for your co-operation.

NHSmail is the secure email and directory service available for all NHS sta=
ff in England and Scotland. NHSmail is approved for exchanging patient data=
 and other sensitive information with NHSmail and other accredited email se=
rvices.

For more information and to find out how you can switch, https://portal.nhs=
.net/help/joiningnhsmail


Received on Tue Feb 05 2019 - 05:53:21 EST

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.