rr: Run record

In NMproject: Script Based 'NONMEM' Model Development

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)

Arguments

m	An nm object.
trans	Logical. If TRUE (default) will transform using control file $THETA/OMEGA conventions.

Value

A tibble with NONMEM run results.

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 trans = TRUE

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 θ and se(θ) are the NONMEM reported values of parameters and standard errors, respectively:

LOG

FINAL = exp(θ), RSE = 100√(exp(se(θ)^2) - 1)

RATIO

FINAL = θ, RSE = 100se(θ)/θ

LOGIT

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

missing

FINAL = θ, SE = se(θ)

OMEGA transformations using trans = TRUE

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 ω^2 and se(ω^2) are the NONMEM reported values of parameters and standard errors, respectively

LOG

FINAL = 100√(exp(ω^2) - 1), RSE = 100(se(ω^2)/ω^2)/2

missing

FINAL = ω^2, SE = se(ω^2)

SIGMA transformations using trans = TRUE

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 σ^2 and se(σ^2) are the NONMEM reported values of parameters and standard errors, respectively. All sigmas are reported as standard deviations.

all sigmas

FINAL = √(σ^2), RSE = 100se(σ^2) / σ^2

See Also

nm_render()

Examples

## requires NONMEM to be installed
## Not run: 

rr(m1)

## compare m1 and m2

rr(c(m1, m2))

## End(Not run)

NMproject documentation built on Sept. 30, 2022, 1:06 a.m.