Plot model fits (observed and fitted) and the residuals for a row or column of an mmkin object

When x is a row selected from an mmkin object ([.mmkin), the same model fitted for at least one dataset is shown. When it is a column, the fit of at least one model to the same dataset is shown.

If the current plot device is a tikz device, then latex is being used for the formatting of the chi2 error level.

# S3 method for mmkin
plot(x, main = "auto", legends = 1,
  resplot = c("time", "errmod"), errmin_var = "All data", errmin_digits = 3,
  cex = 0.7, rel.height.middle = 0.9, ...)

Arguments

x

An object of class mmkin, with either one row or one column.
main

The main title placed on the outer margin of the plot.
legends

An index for the fits for which legends should be shown.
resplot

Should the residuals plotted against time, using mkinresplot, or as squared residuals against predicted values, with the error model, using mkinerrplot.
errmin_var

The variable for which the FOCUS chi2 error value should be shown.
errmin_digits

The number of significant digits for rounding the FOCUS chi2 error percentage.
cex

Passed to the plot functions and mtext.
rel.height.middle

The relative height of the middle plot, if more than two rows of plots are shown.
…

Further arguments passed to plot.mkinfit and mkinresplot.

Value

The function is called for its side effect.

Examples

  
  # Only use one core not to offend CRAN checks
  fits <- mmkin(c("FOMC", "HS"),
                list("FOCUS B" = FOCUS_2006_B, "FOCUS C" = FOCUS_2006_C), # named list for titles
                cores = 1, quiet = TRUE, error_model = "tc")
#> Warning: Optimisation did not converge:
#> iteration limit reached without convergence (10)
  plot(fits[, "FOCUS C"])
  plot(fits["FOMC", ])

  # We can also plot a single fit, if we like the way plot.mmkin works, but then the plot
  # height should be smaller than the plot width (this is not possible for the html pages
  # generated by pkgdown, as far as I know).
  plot(fits["FOMC", "FOCUS C"]) # same as plot(fits[1, 2])

  # Show the error models
  plot(fits["FOMC", ], resplot = "errmod")