plot.gstatVariogram {gstat} | R Documentation |
Creates a variogram plot
plot.gstatVariogram(x, model = NULL, ylim, xlim, xlab = "distance", ylab = "semivariance", multipanel = TRUE, plot.numbers = FALSE, scales, ids = x$id, group.id = TRUE, skip, layout, ...) plot.variogramMap(x, np = FALSE, skip, ...)
x |
object of class "gstatVariogram", obtained from the function variogram, possibly containing directional or cross variograms |
model |
in case of a single variogram: a variogram model, as obtained from vgm or fit.variogram, to be drawn as a line in the variogram plot; in case of a set of variograms and cross variograms: a list with variogram models |
ylim |
numeric vector of length 2, limits of the y-axis |
xlim |
numeric vector of length 2, limits of the x-axis |
xlab |
x-axis label |
ylab |
y-axis label |
multipanel |
logical; if TRUE, directional variograms are plotted in different panels, if FALSE, directional variograms are plotted in the same graph, using color, colored lines and symbols to distinguish them |
plot.numbers |
logical; if TRUE, plot number of point pairs next to each plotted semivariance symbol |
scales |
optional argument that will be passed to xyplot in
case of the plotting of variograms and cross variograms |
ids |
ids of the data variables and variable pairs |
group.id |
logical; control for directional multivariate variograms: if TRUE, panels divide direction and colors indicate variables (ids), if FALSE panels divide variables/variable pairs and colors indicate direction |
skip |
logical; can be used to arrange panels, see xyplot |
layout |
integer vector; can be used to set panel layout: c(ncol,nrow) |
np |
logical; if TRUE, plot number of point pairs, if FALSE plot semivariances |
... |
any arguments that will be passed to the panel plotting functions
(such as auto.key in examples below) |
returns (or plots) the variogram plot
currently, plotting models and/or point pair numbers is not supported when a variogram is both directional and multivariable; also, three-dimensional directional variograms will probably not be displayed correctly.
Edzer J. Pebesma
variogram, fit.variogram, vgm variogram.line,
data(meuse) vgm1 <- variogram(log(zinc)~1, ~x+y, meuse) plot(vgm1) model.1 <- fit.variogram(vgm1,vgm(1,"Sph",300,1)) plot(vgm1, model=model.1) plot(vgm1, plot.numbers = TRUE, pch = "+") vgm2 <- variogram(log(zinc)~1, ~x+y, meuse, alpha=c(0,45,90,135)) plot(vgm2) # the following demonstrates plotting of directional models: model.2 <- vgm(.59,"Sph",926,.06,anis=c(0,0.3)) plot(vgm2, model=model.2) g = gstat(id="zinc < 200", form=I(zinc<200)~1,loc=~x+y,data=meuse) g = gstat(g, id="zinc < 400", form=I(zinc<400)~1,loc=~x+y,data=meuse) g = gstat(g, id="zinc < 800", form=I(zinc<800)~1,loc=~x+y,data=meuse) # calculate multivariable, directional variogram: v = variogram(g, alpha=c(0,45,90,135)) plot(v, group.id = FALSE, auto.key = TRUE) # id and id pairs panels plot(v, group.id = TRUE, auto.key = TRUE) # direction panels #if (require(sp)) { # plot(variogram(g, cutoff=1000, width=100, map=TRUE), # main = "(cross) semivariance maps") # plot(variogram(g, cutoff=1000, width=100, map=TRUE), np=TRUE, # main = "number of point pairs") #}