From c6079a807e2b400fe0c772603392aeacd887da2f Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 8 May 2019 20:57:48 +0200 Subject: Add functionality to plot the error model by plotting squared residuals against predicted values, and showing the variance function used in the fitted error model. Rebuild docs --- docs/reference/plot.mmkin.html | 19 ++++++++++++++----- 1 file changed, 14 insertions(+), 5 deletions(-) (limited to 'docs/reference/plot.mmkin.html') diff --git a/docs/reference/plot.mmkin.html b/docs/reference/plot.mmkin.html index b3f2de66..3037ca24 100644 --- a/docs/reference/plot.mmkin.html +++ b/docs/reference/plot.mmkin.html @@ -141,8 +141,9 @@ If the current plot device is a tikz device, 
# S3 method for mmkin
-plot(x, main = "auto", legends = 1, errmin_var = "All data", errmin_digits = 3,
-              cex = 0.7, rel.height.middle = 0.9, ...)
+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

@@ -159,6 +160,12 @@ If the current plot device is a tikz device, + + + + @@ -187,15 +194,17 @@ If the current plot device is a tikz device,

Examples

- 
  # Only use one core not to offend CRAN checks
+    
  
  # 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)
+                cores = 1, quiet = TRUE, error_model = "tc")
   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])

+  plot(fits["FOMC", "FOCUS C"]) # same as plot(fits[1, 2])

+  # Show the error models
+  plot(fits["FOMC", ], resplot = "errmod")
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.