pmxTools (version 1.5)

get_omega: Extract variability parameter estimates from a NONMEM output object.

Description

Extract variability parameter estimates from a NONMEM output object.

Usage

get_omega(x, output = "est", sigdig = 6, sep = "-", est.step = NULL)

Value

A symmetrical matrix, or a list of symmetrical matrices if all is specified.

Arguments

x

A NONMEM output object generated using read_nm.

output

A flag specifying the matrix or matrices to be output. Valid flag values are est (the default), se, rse, cor, cse, 95ci, or all.

sigdig

Specifies the number of significant digits to be provided (default=6).

sep

Specifies the separator character to use for 95% confidence intervals (default="-").

est.step

Specifies which estimation step to return parameters from (default is the last).

Author

Justin Wilkins,

Details

Output options are as follows:

• est returns the estimated OMEGA variance-covariance matrix.

• se returns the standard errors for the estimated OMEGA variance-covariance matrix.

• rse returns the relative standard errors for the estimated OMEGA variance-covariance matrix (se/est*100).

• cor returns the correlation matrix matrix.

• cse returns the standard errors for the correlation matrix.

• 95ci returns the asymptotic 95% confidence intervals for the elements of the OMEGA variance-covariance matrix (est +/- 1.96*se).

• all returns all available OMEGA matrices.

See Also

NONMEM (https://www.iconplc.com/solutions/technologies/nonmem)

Examples

if (FALSE) {
 nmOutput  <- read_nm("run315.xml")
 omegas    <- get_omega(nmOutput)
 omegaRSEs <- get_omega(nmOutput, "rse")
}