read_nmext: Extract estimates from NONMEM ext file

In metrumresearchgroup/mrgsolve: Simulate from ODE-Based Models

Extract estimates from NONMEM ext file

Description

This function retrieves NONMEM estimates for use in the mrgsolve model when ⁠$NMEXT⁠ is invoked. See nmext().

Usage

read_nmext(
  run = NA_real_,
  project = getwd(),
  file = paste0(run, ".ext"),
  path = NULL,
  read_fun = c("data.table", "read.table"),
  index = "last"
)

Arguments

run	a run number or run identifier
project	the NONMEM project directory
file	the ext file name
path	full path and file name for ext file
read_fun	function to read the ext file; data.table::fread() will be used if available; otherwise utils::read.table() is used.
index	selects the table number whose results will be returned; use value "last" to select the last table in the .ext file; or pass an integer specifying the table number; in case there is exactly one table in the .ext file, pass the value "single" to bypass parsing the file to look for sub tables (this might be useful when BAYES analysis was performed as the only estimation method and there are 10000s of posterior samples in the file)

Value

A list with param, omega, and sigma in a format ready to be used to update a model object.

Examples

project <- system.file("nonmem", package = "mrgsolve")

est <- read_nmext(1005, project = project)

est$param

est$omega

est$sigma

est <- read_nmext(2005, project = project, index = 3)

metrumresearchgroup/mrgsolve documentation built on Feb. 13, 2024, 10:27 p.m.