NMproject (version 0.6.9)

rr: Run record

Description

Displays the transformed parameters of a completed or running model. Normally used inside of a diagnostic template, but can be useful for quickly seeing parameter estimates of several models.

Usage

rr(m, trans = TRUE)

Value

A tibble with NONMEM run results.

Arguments

m

An nm object.

trans

Logical. If TRUE (default) will transform using control file $THETA/OMEGA conventions.

NONMEM coding conventions used by NMproject

The convention for $THETA comments used by NMproject is value ; name ; unit ; transformation

e.g. $THETA 0.1 ; KA ; h-1 ; LOG

The options for THETA transformations are: LOG, LOGIT, RATIO and missing. LOG and LOGIT refer to log and logit transformed THETAs, respectively where the parameters should be back-transformed for reporting. RATIO refers to ratio data types, i.e. parameters that are positive and have a meaningful zero. Most parameters like KA, CL, EMAX fall into this category, but covariates effects which can go negative do not. RSEs are calculated for ratio data. Missing transformations are suitable for all other parameters, here no RSEs will be calculated, only raw SE values will be reported.

The convention for $OMEGA is similar but without a unit item: value ; name ; transformation

e.g. $OMEGA 0.1 ; IIV_KA ; LOG

The options for OMEGA are either LOG or missing. LOG indicating that the individual parameter distribution is log normally distributions and should be reported as a CV% (and associated RSE%) rather than as the raw NONMEM estimate.

The convention for $OMEGA is just : value ; name.

THETA transformations using <code>trans = TRUE</code>

The value of FINAL and RSE% (always accompanied with a % symbol in outputs) in the returned tibble is the reported standard error (where applicable) where \(\theta\) and \(se(\theta)\) are the NONMEM reported values of parameters and standard errors, respectively:

LOG

\(FINAL = exp(\theta), RSE = 100\sqrt(exp(se(\theta)^2) - 1)\)

RATIO

\(FINAL = \theta, RSE = 100se(\theta)/\theta\)

LOGIT

\(FINAL = 100/(1 + exp(-\theta)), SE = se(\theta)\)

missing

\(FINAL = \theta, SE = se(\theta)\)

OMEGA transformations using <code>trans = TRUE</code>

The value of FINAL and RSE% (always accompanied with a % symbol in outputs) in the returned tibble is the reported standard error (where applicable) where \(\omega^2\) and \(se(\omega^2)\) are the NONMEM reported values of parameters and standard errors, respectively

LOG

\(FINAL = 100\sqrt(exp(\omega^2) - 1), RSE = 100(se(\omega^2)/\omega^2)/2\)

missing

\(FINAL = \omega^2, SE = se(\omega^2)\)

SIGMA transformations using <code>trans = TRUE</code>

The value of FINAL and RSE% (always accompanied with a % symbol in outputs) in the returned tibble is the reported standard error (where applicable) where \(\sigma^2\) and \(se(\sigma^2)\) are the NONMEM reported values of parameters and standard errors, respectively. All sigmas are reported as standard deviations.

all sigmas

\(FINAL = \sqrt(\sigma^2), RSE = 100se(\sigma^2) / \sigma^2\)

See Also

nm_render()

Examples

## requires NONMEM to be installed
if (FALSE) {

rr(m1)

## compare m1 and m2

rr(c(m1, m2))
}