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.

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

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.
ylab: Label for the y axis.
standardized: Should the residuals be standardized? This option is passed to mkinresplot, it only takes effect if resplot = "time".
show_errmin: Should the chi2 error level be shown on top of the plots to the left?
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.
ymax: Maximum y axis value for plot.mkinfit.
...: Further arguments passed to plot.mkinfit and mkinresplot.

Value

The function is called for its side effect.

Details

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

Author

Johannes Ranke

Examples


  # \dontrun{
  # 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")
  plot(fits[, "FOCUS C"])

  plot(fits["FOMC", ])

  plot(fits["FOMC", ], show_errmin = FALSE)


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

  # }