block-omega-sigma: Create or remove $OMEGA/$SIGMA BLOCKs

In NMproject: Script Based 'NONMEM' Model Development

Create or remove $OMEGA/$SIGMA BLOCKs

Description

Manipulate $OMEGA (and $SIGMA) BLOCKs to introduce or remove correlations.

Usage

block(iomega, eta_numbers = NA, diag_init = 0.01)

unblock(iomega, eta_numbers)

Arguments

iomega	A tibble output from init_omega() or init_sigma().
eta_numbers	Numeric vector. ETA numbers to put into a block or unblock for block() and unblock(), respectively. Must be contiguous.
diag_init	Numeric. Default value for off diagonal elements.

Value

An nm object with modified ctl_contents field.

See Also

init_theta(), init_omega(), init_sigma()

Examples

# create example object m1 from package demo files
exdir <- system.file("extdata", "examples", "theopp", package = "NMproject")
m1 <- new_nm(run_id = "m1", 
             based_on = file.path(exdir, "Models", "ADVAN2.mod"),
             data_path = file.path(exdir, "SourceData", "THEOPP.csv"))

io <- m1 %>% init_omega()
io <- io %>% block(c(2, 3))
m1 <- m1 %>% init_omega(io)
m1 %>% dollar("OMEGA") ## to display $OMEGA


# create example object m1 from package demo files
exdir <- system.file("extdata", "examples", "theopp", package = "NMproject")
m1 <- new_nm(run_id = "m1", 
             based_on = file.path(exdir, "Models", "ADVAN2.mod"),
             data_path = file.path(exdir, "SourceData", "THEOPP.csv"))

## first create a block
io <- m1 %>% init_omega()
io <- io %>% block(c(2, 3))
m1 <- m1 %>% init_omega(io)
m1 %>% dollar("OMEGA") ## to display $OMEGA

## now unblock
io <- io %>% unblock(c(2, 3))
m1 <- m1 %>% init_omega(io)
m1 %>% dollar("OMEGA") ## to display $OMEGA

NMproject documentation built on Sept. 30, 2022, 1:06 a.m.