A plotting utile for quick visualization of the quality of an ABC
analysis from an object of class "abc" generated with methods
"loclinear" or "neuralnet" (see abc for
details). Four plots are currently available: a density plot of the
prior distribution, a density plot of the posterior distribution, a
scatter plot of the Euclidean distances as a function of the parameter
values, and a Normal Q-Q plot of the residuals from the
regression.
# S3 method for abc
plot(x, param, subsample = 1000, true = NULL, file = NULL,
postscript = FALSE, onefile = TRUE, ask =
!is.null(deviceIsInteractive()), ...)
an object of class "abc" generated with methods
"loclinear" or "neuralnet" (see abc for
details).
a vector or matrix of parameter values from the simulations that
were used in the original call to abc.
the number of rows (simulations) to be plotted. Rows are randomly
selected from param.
a vector of true parameter values, if known. Vertical lines are drawn at these values.
a character string giving the name of the file. See
postscript for details on accepted file names. If
NULL (the default) plots are printed to the null device
(e.g. X11). If not NULL plots are printed on a
pdf device. See also postscript.
logical; if FALSE (default) plots are printed on a
pdf device, if TRUE on a postscript device.
logical, if TRUE (the default) allow multiple figures in one
file. If FALSE, generate a file name containing the page
number for each page. See postscript for further details.
logical; if TRUE (the default), the user is asked before each plot, see
par(ask=.).
other parameters passed to plot.
In order to use this function, one of the regression correction
methods had to be used in the original call to abc,
i.e. "loclinear" or "neuralnet" (see abc
for details). Four plots are printed for each parameter. (i) A density
plot of the prior distribution. (ii) A density plot of the posterior
distribution using the regression correction (red thick lines) and,
for reference, using the simple rejection method (black fine
lines). The prior distribution (in the posterior distributions' range)
is also displayed (dashed lines). (iii) A scatter plot of the log
Euclidean distances as a function of the true parameter values. Points
corresponding to the accepted simulations are displayed in red. (iv) A
Normal Q-Q plot of the residuals from the regression, thus from
lsfit when method was "loclinear", and from
nnet when method was "neuralnet" in the original
abc.
For plots (i) and (iii) not the whole data but a subsample is used,
the size of which can be is given by subsample. This is to
avoid plots that may take too much time to print.
If a parameter transformation was applied in the original call to
abc, the same transformations are applied to the
parameters for plotting (on plots (i)-(iii)).
abc, hist.abc, summary.abc