Read information from Nonmem ext files
NMreadExt(
file,
return,
as.fun,
modelname,
col.model,
auto.ext,
tableno = "max",
file.ext
)
Path to the ext file
The .ext file contains both final parameter
estimates and iterations of the estimates. If
return="pars"
(default) the final estimates are
returned in addition to what other parameter-level information
is found, like FIX, sd etc. as columns. If
return="iterations"
, the iterations are returned
(including objective function value). If return="obj"
objective function value at final estimate is returned. If
return="all"
, all er returned, though in separate
data.frames compiled in a list.
The default is to return data as a data.frame. Pass a function (say `tibble::as_tibble`) in as.fun to convert to something else. If `data.table`s are wanted, use `as.fun="data.table"`. The default can be configured using `NMdataConf()`.
See `?NMscanData`
See `?NMscanData`
If `TRUE` (default) the extension will automatically be modified using `NMdataConf()$file.ext`. This means `file` can be the path to an input or output control stream, and `NMreadExt` will still read the `.ext` file.
In case the ext file contains multiple tables, this argument controls which one to choose. The options are
"max" (default) Pick the table with the highest table number. This typically means the results from the last `$ESTIMATION` step are used.
"min" Pick results from the first table available.
"all" Keep all results. The tables can be distinguished by the `tableno` column.
an integer greater than 0, in which case the table with this table number will be picked.
Deprecated. Please use file
instead.
If return="all"
, a list with a final parameter
table and a table of the iterations. If return="pars"
,
only the parameter table, and if return="iterations"
only the iterations table. If you need both, it may be more
efficient to only read the file once and use
return="all"
. Often, only one of the two are needed,
and it more convenient to just extract one.
The parameter table returned if return="pars"
or return="all"
will contain columns based on the Nonmem 7.5 manual. It defines codes for different parameter-level values. They are:
-1e+09: se -1000000002: eigCor -1000000003: cond -1000000004: stdDevCor -1000000005: seStdDevCor -1000000006: FIX -1000000007: termStat -1000000008: partLik
The parameter name is in the parameter
column. The
"parameter type", like "THETA", "OMEGA", "SIGMA" are available in
the par.type
column. Counters are available in i
and
j
columns. j
will be NA
for
par.type=="THETA"
The objective function value is included as a parameter.
Notice that in case multiple tables are available in the `ext` file, the column names are taken from the first table. E.g., in case of SAEM/IMP estimation, the objective function values will be in the `SAEMOBJ` column, even for the IMP step. This may change in the future.