plot.xvalid {geoR}R Documentation

Plot Cross-Validation Results

Description

This function produces ten plots with the results produced by the cross-validation function xvalid.

Usage

plot(x, coords, borders = NULL, ask = TRUE,
            error = TRUE, std.error = TRUE, data.predicted = TRUE,
            pp = TRUE, map = TRUE, histogram = TRUE,
            error.predicted = TRUE, error.data = TRUE, ...)

Arguments

x an object of the class "xvalid", typically an output from the function xvalid.
coords an n x 2 object containing coordinates of the (cross-)validation locations.
borders optional. Takes a two column matrix or data-frame with coordinates of the borders. If provided the borders are included in the errors maps.
ask logical. Defines whether or not the user is prompted before each plot is produced.
error logical. Defines whether the plots for the errors (error = data - predicted) will be produced.
std.error logical. Defines whether the plots for the standardised errors will be produced.
data.predicted logical defining whether a plot of data versus predicted should be displayed. Defaults to TRUE.
pp logical defining whether a pp plot should be displayed. Defaults to TRUE.
map logical defining whether a map of the errors should be displayed. Defaults to TRUE.
histogram logical defining whether a histogram of the errors should be displayed. Defaults to TRUE.
error.predicted logical defining whether a plot of errors versus predicted should be displayed. Defaults to TRUE.
error.data logical defining whether a plot of errors versus data should be displayed. Defaults to TRUE.
... other arguments to be passed to the function plot.

Value

No value returned. Plots are produced on the current graphics device.

See Also

xvalid for the cross-validation computations.

Examples

if(is.R()) data(s100)
wls <- variofit(variog(s100, max.dist = 1), ini = c(.5, .5), fix.n = TRUE)
xvl <- xvalid(s100, model = wls)
#
op <- par(no.readonly = TRUE)
par(mfcol = c(3,2))
par(mar = c(3,3,0,1))
par(mgp = c(2,1,0))
plot(xvl, error = FALSE, ask = FALSE)
plot(xvl, std.err = FALSE, ask = FALSE)
par(op)


[Package Contents]