plot.wa {analogue}R Documentation

Plot diagnostics for a weighted averaging model

Description

Two plots (selectable by which) are currently available: a plot of estimated against observed values, a plot of residuals against estimated values.

Usage

## S3 method for class 'wa':
plot(x,
     which = 1:2,
     caption = c("Inferred vs Observed", "Residuals vs Fitted"),
     max.bias = TRUE,
     n.bias = 10,
     sub.caption = NULL,
     main = "",
     ask = prod(par("mfcol")) < length(which) &&
                                  dev.interactive(),
     ...,
     panel = if (add.smooth) panel.smooth else points,
     add.smooth = getOption("add.smooth"))

Arguments

x an object of class "wa".
which which aspects of the "wa" object to plot if a subset of the plots is required, specify a subset of the numbers 1:2.
caption captions to appear above the plots.
max.bias logical, should max bias lines be added to residuals.
n.bias numeric, number of sections to calculate maximum bias for.
sub.caption common title-above figures if there are multiple; used as ‘sub’ (s.‘title’) otherwise. If NULL, as by default, a possibly shortened version of deparse(x$call) is used.
main title to each plot-in addition to the above caption.
ask logical; if TRUE, the user is asked before each plot, see par(ask=.).
... graphical arguments passed to other graphics functions.
panel panel function. The useful alternative to points, panel.smooth, can be chosen by add.smooth = TRUE.
add.smooth logical indicating if a smoother should be added to fitted & residuals plots; see also panel above.

Details

This plotting function is modelled closely on plot.lm and many of the conventions and defaults for that function are replicated here.

sub.caption - by default the function call - is shown as a subtitle (under the x-axis title) on each plot when plots are on separate pages, or as a subtitle in the outer margin (if any) when there are multiple plots per page.

Value

One or more plots, drawn on the current device.

Author(s)

Gavin L. Simpson. Code borrows heavily from plot.lm.

See Also

mat

Examples

## continue the RLGH example from ?wa
example(wa)

## diagnostics for the WA model
par(mfrow = c(1,2))
plot(mod)
par(mfrow = c(1,1))


[Package analogue version 0.6-6 Index]