Plot time series of decline data

Source: R/twa.R

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