vardiag (version 0.2-2)

PlotDiag.varobj: Diagnostic Plots for Variograms

Description

Produces diagnostic plots for a variogram object.

Usage

PlotDiag.varobj(v, region = NULL, xyi = 0, zmv = 0)

Value

No value is returned. The function is called for its side effects, i.e. creating plots.

Arguments

v

A variogram object

region

a n by 2 matrix of vertices of a polygon, which represents the region of interest.

xyi

the index of the pair of observations, which shall be marked in the plots.

zmv

the value of zmv= determines which values will be plotted as third dimension (via circles with radius corresponding to the value) in the map view. The following values are possible:

0

no third dimension,

1

values of the original observations,

2

Cook's distances,

3

reduction in lack of fit.

Author

Ernst Glatzer

Details

Four plots are created on the current graphics device.

These plots are a map view, a square root cloud with the fitted model, a leave-one-out plot and a plot of decorrelated residuals. If the argument region is given, this region is plotted as a polygon in the map view. If the argument xyi is given, the corresponding pair of observations is marked as a blue and a red point in the map view and connected by a magenta line. Furthermore all the points in the square root cloud, which contain the first (second) observation in the pair are marked blue (red). The one point that contains both observations of the selected pair is marked magenta. The points in the leave-one-out plot are marked analogously.

References

E. Glatzer, W.G. Mueller, 2004. Residual diagnostics for variogram fitting. Computers and Geosciences 30, pp 859-866.

See Also

interact.varobj

Examples

data(tulln)
PlotDiag.varobj(vs50,region=tu1)
PlotDiag.varobj(vs50,region=tu1,xyi=144)