Extract THETA Estimates from a NONMEM Listing File

Parses the THETA (fixed effects) parameter estimates and, if available, their standard errors and relative standard errors (RSE) from a NONMEM .lst file.

Usage

fetch_thetas(lst, digits = NA)

Arguments

lst

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

digits

Integer or NA. Number of decimal places to round RSE. Default NA (no rounding).

Value

A tibble with columns:

parameter

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

estimate

Numeric. Point estimate.

se

Numeric. Standard error. NA if no covariance step.

rse

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

Examples

path <- system.file("testdata", "full_cov.lst", package = "lstparsR")
lst  <- read_lst_file(path)
fetch_thetas(lst)
#> # A tibble: 12 × 4
#>    parameter     estimate       se     rse
#>    <chr>            <dbl>    <dbl>   <dbl>
#>  1 TH_1          34.1     3.37e+ 0 9.88e 0
#>  2 TH_2      387000       5.41e+ 7 1.40e 4
#>  3 TH_3           8.05    5.2 e- 1 6.46e 0
#>  4 TH_4           1.78    1.49e- 1 8.37e 0
#>  5 TH_5           0.498   2.78e- 2 5.58e 0
#>  6 TH_6           0.362   4.81e- 1 1.33e 2
#>  7 TH_7          -0.171   1.35e- 1 7.89e 1
#>  8 TH_8          -0.0314  2.69e- 1 8.57e 2
#>  9 TH_9          -0.163   2.78e- 1 1.71e 2
#> 10 TH10          -0.23    1.88e- 1 8.17e 1
#> 11 TH11           0.00539 1.98e- 2 3.67e 2
#> 12 TH12           0.2     1.90e+74 9.50e76