https://github.com/cran/gstat
Tip revision: 1b2742a3d71abc7fbde73ca141b23e1375f185bf authored by Edzer Pebesma on 30 October 2011, 00:00:00 UTC
version 1.0-9
version 1.0-9
Tip revision: 1b2742a
hscat.R
hscat = function(formula, data, breaks, pch = 3, cex = .6, mirror = FALSE,
variogram.alpha = 0, ...) {
stopifnot(!missing(breaks))
x = variogram(formula, data, cloud = TRUE, cutoff = max(breaks),
alpha = variogram.alpha, ...)
x = as.data.frame(x)
x$class = cut(x$dist, breaks = breaks)
y = model.frame(formula, data)[[1]]
x$xx = y[x$left]
x$yy = y[x$right]
if (mirror)
x = data.frame(
xx = c(x$yy, y[x$left]),
yy = c(x$xx, y[x$right]),
class = c(x$class, x$class))
lab = as.character(formula)[2]
panel = function(x,y,subscripts, ...) {
xr = c(min(x),max(x))
llines(xr, xr)
lpoints(x,y,...)
ltext(min(x), max(y), paste("r =", signif(cor(x,y),3)), adj=c(0,0.5))
}
xyplot(xx~yy|class, x, panel = panel,
main = "lagged scatterplots", xlab = lab, ylab = lab, ...)
}