lotriSep: Separate a lotri matrix into above and below lotri matrices

In nlmixrdevelopment/lotri: A Simple Way to Specify Symmetric, Block Diagonal Matrices

Description

This is used for creating nesting simulations in 'RxODE' and may not be useful for external function calls.

Usage

lotriSep(x, above, below, aboveStart = 1L, belowStart = 1L)

Arguments

x	lotri matrix
above	Named integer vector listing variability above the id level. Each element lists the number of population differences in the whole data-set (as integer)
below	Named integer vector listing variability below the id level. Each element lists the number of items below the individual level. For example with 3 occasions per indivdiual you could use 'c(occ=3L)'
aboveStart	Add the attribute of where THETA[#] will be added
belowStart	Add the attribute of where ETA[#] will be added

Value

List of two lotri matrices

Author(s)

Matthew Fidler

Examples

omega <- lotri(lotri(eta.Cl ~ 0.1,
                        eta.Ka ~ 0.1) | id(nu=100),
                  lotri(eye.Cl ~ 0.05,
                        eye.Ka ~ 0.05) | eye(nu=50),
                  lotri(iov.Cl ~ 0.01,
                        iov.Ka ~ 0.01) | occ(nu=200),
                  lotri(inv.Cl ~ 0.02,
                        inv.Ka ~ 0.02) | inv(nu=10))

lotriSep(omega, above=c(inv=10L), below=c(eye=2L, occ=4L))

nlmixrdevelopment/lotri documentation built on April 22, 2021, 2:21 a.m.