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

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)
#> 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")