Plot time series of decline data

Plot time series of decline data

# S3 method for one_box
plot(x, xlim = range(time(x)), ylim = c(0, max(x)),
  xlab = "Time", ylab = "Fraction of initial", max_twa = NULL,
  max_twa_var = dimnames(x)[[2]][1], ...)

Arguments

x

The object of type one_box to be plotted

xlim

Limits for the x axis

ylim

Limits for the y axis

xlab

Label for the x axis

ylab

Label for the y axis

max_twa

If a numeric value is given, the maximum time weighted average concentration(s) is/are shown in the graph.

max_twa_var

Variable for which the maximum time weighted average should be shown if max_twa is not NULL.

...

Further arguments passed to methods

See also

sawtooth

Examples

dfop_pred <- one_box("DFOP", parms = c(k1 = 0.2, k2 = 0.02, g = 0.7))
plot(dfop_pred)
plot(sawtooth(dfop_pred, 3, 7), max_twa = 21)

# Use a fitted mkinfit model
m_2 <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))
#> Successfully compiled differential equation model from auto-generated C code.
fit_2 <- mkinfit(m_2, FOCUS_2006_D, quiet = TRUE)
pred_2 <- one_box(fit_2)
pred_2_saw <- sawtooth(pred_2, 2, 7)
plot(pred_2_saw, max_twa = 21, max_twa_var = "m1")