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 = "Residue",
  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

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"))
#> Temporary DLL for differentials generated and loaded
fit_2 <- mkinfit(m_2, FOCUS_2006_D, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
pred_2 <- one_box(fit_2, ini = 1)
pred_2_saw <- sawtooth(pred_2, 2, 7)
plot(pred_2_saw, max_twa = 21, max_twa_var = "m1")