Exploratory Geostatistical Plots

This function produces a \(2 \times 2\) display with the following plots: the first indicates the spatial locations assign different colors to data in different quartiles, the next two shows data against the X and Y coordinates and the last is an histogram of the data values or optionally, a 3-D plot with spatial locations and associated data values.

spatial, dplot
# S3 method for geodata
plot(x, coords=x$coords, data = x$data,
             borders, trend="cte", lambda = 1, = 1,
             weights.divide = "units.m", lowess = FALSE, scatter3d = FALSE,
             density = TRUE, rug = TRUE, qt.col, …)

a list containing elements coords and data described next. Typically an object of the class "geodata" - a geoR data-set. If not provided the arguments coords and data must be provided instead.


an \(n \times 2\) matrix containing in each row Euclidean coordinates of the n data locations. By default it takes the element coords of the argument geodata.


a vector with data values. By default it takes the element data of the argument geodata.


If an \(n \times 2\) matrix or data-frame with the borders of the area is provided, the borders are included in the first plot. By default it searches for a element named "borders" in the geodata object.


specifies the mean part of the model. The options are: "cte" (constant mean - default option), "1st" (a first order polynomial on the coordinates), "2nd" (a second order polynomial on the coordinates), or a formula of the type ~X where X is a matrix with the covariates (external trend). If provided the trend is "removed" using the function lm and the residuals are plotted.


value of the Box-Cox transformation parameter. Two particular cases are \(\lambda = 1\) which corresponds to no transformation and \(\lambda = 0\) corresponding to the log-transformation.

indicates the column number for the data to be plotted. Only valid if more than one data-set is available i.e., if the argument data is a matrix.


if a vector of weights with the same length as the data is provided each data is divided by the corresponding element in this vector. Defaults divides the data by the element units.m in the data object, if present, otherwise no action is taken and original data is used. The usage of units.m is common for data objects to be analysed using the package geoRglm.


logical. Indicates whether the function lowess should be used in the plots of the data against the coordinates.


logical. If TRUE the last plot is produced by scatterplot3d showing a 3d plot with data locations and corresponding values.


logical. If TRUE (default) a line with density estimation is added to the histogram.


logical. If TRUE a rug plot is added to the histogram.


colors for the quartiles in the first plot. If missing defaults to blue, green, yellow and red.

further arguments to be passed to the function hist or scatterplot3d.


A plot is produced on the graphics device. No values are returned.


Further information on the package geoR can be found at:

See Also

points.geodata, scatterplot3d, lowess, density, rug.

  • plot.geodata
plot(s100, scatter3d=TRUE)
plot(s100, qt.col=1)

plot(ca20)                       # original data
plot(ca20, trend=~altitude+area) # residuals from an external trend
plot(ca20, trend='1st')          # residuals from a polynomial trend

plot(sic.100, bor=sic.borders)           # original data
plot(sic.100, bor=sic.borders, lambda=0) # logarithm of the data
# }
Documentation reproduced from package geoR, version 1.8-1, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.