How to translate Monolix Combined 1 ( additive + proportional) error model in Phoenix

General
1

Hi, I am trying to replicate a published PopPK model built in Monolix with Combined1 error model:

image
image819×108 17.6 KB

There we get additive and proportional error value. However, in Phoenix we have addi+multi error model as:
observe(CObs = C + CEps * sqrt(1 + C^2 * (CMultStdev/sigma())^2)). So, how to define a and b in this?
Currently, i understand is
CEps = 1 and
observe(CObs = C + CEps * sqrt(a^2 + (b * (C)^2)). Is this the right way to represent in Phoenix NLME, please guide.

2

Hi Bhim, what I would suggest is use the built-in or graphical models to learn some of the syntax conventions; as you change the options you will see the generated code and statements change for you which is a great way to learn e.g.
(Model text and also Warnings tab can be informative too !)

image
image816×356 16.6 KB

so typically we represent as you said as
observe(CObs = C + CEps * sqrt(1 + C^2 * (CMultStdev/sigma())^2))
I would recommend to look at the PML guide, particularly page 21 found under https://help.certara.net/en/articles/10769710-phoenix-8-6-x-documentation which explains that

Only a single error variable can be used in an observe statement such as the one on line 23.
Compound residual error models for any given observed variable, such as mixed additive
and proportional, must be built using a combination of fixed effects and a single error vari-
able rather than multiple error variables. For example, the following statements are correct:
stparm(CMixRatio = tvCMixRatio)
fixef(tvCMixRatio = c(, 0.1, ))
error(eps1 = 0.01)
observe(cObs = c + eps1(1 + c*MixRatio))*
The following statements are incorrect as the observe statement contains two error vari-
ables, eps1 and eps2:
error(eps1 = 0.01)
error(eps2 = 0.001)
observe(cObs = c*(1 + eps2) + eps1)

For more information on the observe statement, see the “Observations” section.

3

Thank you SimonDavis.