plot.nlmixr2AugPred: Plot a nlmixr2 augPred object

In nlmixr2plot: Nonlinear Mixed Effects Models in Population PK/PD, Plot Functions

Plot a nlmixr2 augPred object

Description

Plot a nlmixr2 augPred object

Usage

## S3 method for class 'nlmixr2AugPred'
plot(x, y, ...)

Arguments

x	augPred object
y	ignored, used to mach plot generic
...	Other arguments (ignored)

Value

Nothing called for side effects

Examples

library(nlmixr2est)
## The basic model consiss of an ini block that has initial estimates
one.compartment <- function() {
  ini({
    tka <- 0.45 # Log Ka
    tcl <- 1 # Log Cl
    tv <- 3.45    # Log V
    eta.ka ~ 0.6
    eta.cl ~ 0.3
    eta.v ~ 0.1
    add.sd <- 0.7
  })
  # and a model block with the error sppecification and model specification
  model({
    ka <- exp(tka + eta.ka)
    cl <- exp(tcl + eta.cl)
    v <- exp(tv + eta.v)
    d/dt(depot) = -ka * depot
    d/dt(center) = ka * depot - cl / v * center
    cp = center / v
    cp ~ add(add.sd)
  })
}

## The fit is performed by the function nlmixr/nlmix2 specifying the model, data and estimate
fit <- nlmixr2est::nlmixr2(one.compartment, theo_sd,  est="saem", saemControl(print=0))

# augPred shows more points for the fit:

a <- nlmixr2est::augPred(fit)

# you can plot it with plot(augPred object)
plot(a)

nlmixr2plot documentation built on Oct. 20, 2022, 9:06 a.m.