Extract OMEGA (ETA) Estimates from a NONMEM Listing File

Parses the OMEGA covariance matrix diagonal (variance of random effects) and, if available, standard errors, RSE, and ETA shrinkage from a NONMEM .lst file.

Usage

fetch_etas(lst, digits = NA, shk_digits = NA)

Arguments

lst

An object of class "lst" from read_lst_file().

digits

Integer or NA. Rounding for rse. Default NA.

shk_digits

Integer or NA. Rounding for shrinkage. Default NA.

Value

A tibble with columns:

parameter

Character. Parameter label (e.g. "ETA1").

estimate

Numeric. Diagonal variance estimate.

se

Numeric. Standard error. NA if no covariance step.

rse

Numeric. Relative standard error (%). NA if no covariance step.

shrinkage

Numeric. ETA shrinkage (%). NA if not reported.

Examples

path <- system.file("testdata", "full_cov.lst", package = "lstparsR")
lst  <- read_lst_file(path)
fetch_etas(lst)
#> # A tibble: 7 × 5
#>   parameter estimate        se      rse shrinkage
#>   <chr>        <dbl>     <dbl>    <dbl>     <dbl>
#> 1 ETA1       0.2      3.89e- 2  1.94e 1      5.08
#> 2 ETA2       0.255    4.63e- 2  1.82e 1      1.62
#> 3 ETA3       0.1      2   e+73  2   e76    100   
#> 4 ETA4       0.1      1.99e+73  1.99e76    100   
#> 5 ETA5       0.00001  6.85e+68  6.85e75     97.6 
#> 6 ETA6       0.00001 NA        NA           99.0 
#> 7 ETA7       0.00001 NA        NA           97.0