traceplot | R Documentation |
Produce trace-plot for fit if applicable
traceplot(x, ...) ## S3 method for class 'nlmixr2FitCore' traceplot(x, ...)
x |
fit object |
... |
other parameters |
Fit traceplot or nothing.
Rik Schoemaker, Wenping Wang & Matthew L. Fidler
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 <- nlmixr2(one.compartment, theo_sd, est="saem", saemControl(print=0)) # This shows the traceplot of the fit (useful for saem) traceplot(fit)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.