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")
-#> Warning: Optimisation did not converge:
-#> iteration limit reached without convergence (10)
-  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")
-
-  # }
-
-

- -