split_joint_distribution#

pharmpy.modeling.split_joint_distribution(model, rvs=None)[source]#

Splits etas following a joint distribution into separate distributions.

Parameters:

model (Model) – Pharmpy model
rvs (str, list) – Name/names of etas to separate. If None, all etas that are IIVs and non-fixed will become single. None is default.

Returns:

Model – Pharmpy model object

Examples

>>> from pharmpy.modeling import *
>>> model = load_example_model("pheno")
>>> model = create_joint_distribution(model, ['ETA_1', 'ETA_2'])
>>> model.random_variables.etas
⎡ETA₁⎤    ⎧⎡0⎤  ⎡    IVCL      IIV_CL_IIV_V⎤⎫
⎢    ⎥ ~ N⎪⎢ ⎥, ⎢                          ⎥⎪
⎣ETA₂⎦    ⎩⎣0⎦  ⎣IIV_CL_IIV_V      IVV     ⎦⎭
>>> model = split_joint_distribution(model, ['ETA_1', 'ETA_2'])
>>> model.random_variables.etas
ETA₁ ~ N(0, IVCL)
ETA₂ ~ N(0, IVV)