https://github.com/cran/gstat
Tip revision: 1eda32e527c5c1386fbc8c54ec27917287ebdb8b authored by Edzer Pebesma on 27 February 2010, 00:00:00 UTC
version 0.9-70
version 0.9-70
Tip revision: 1eda32e
show.vgms.R
# $Id: show.vgms.q,v 1.6 2008-12-15 14:27:29 edzer Exp $
"show.vgms" <-
function(min = 1e-12 * max, max = 3, n = 50, sill = 1, range = 1,
models = as.character(vgm()$short[c(1:17)]), nugget = 0, kappa.range = 0.5,
plot = TRUE, ...)
{
zero.range.models = c("Nug", "Int", "Lin", "Err")
# print(models)
i = 0
if (length(kappa.range) > 1) { # loop over kappa values for Matern model:
if (missing(models))
models = "Mat"
data = matrix(NA, n * length(kappa.range), 2)
v.level = rep("", n * length(kappa.range))
for (kappa in kappa.range) {
v = vgm(sill, models, range, nugget = nugget, kappa = kappa)
x = variogramLine(v, 0, 1, 0)
data[(i*n+1), ] = as.matrix(x)
x = variogramLine(v, max, n - 1, min)
data[(i*n+2):((i+1)*n), ] = as.matrix(x)
m.name = paste("vgm(", sill, ",\"", models, "\",", range, sep = "")
if (nugget > 0)
m.name = paste(m.name, ",nugget=", nugget, sep = "")
m.name = paste(m.name, ",kappa=", kappa, ")", sep = "")
v.level[(i*n+1):((i+1)*n)] = rep(m.name, n)
i = i + 1
}
} else {
data = matrix(NA, n * length(models), 2)
v.level = rep("", n * length(models))
for (m in models) {
this.range = ifelse(!is.na(pmatch(m,zero.range.models)), 0, range)
v = vgm(sill, m, this.range, nugget = nugget, kappa = kappa.range)
x = variogramLine(v, 0, 1, 0)
data[(i*n+1), ] = as.matrix(x)
x = variogramLine(v, max, n - 1, min)
data[(i*n+2):((i+1)*n), ] = as.matrix(x)
m.name = paste("vgm(", sill, ",\"", m, "\",", this.range, sep = "")
if (nugget > 0)
m.name = paste(m.name, ",nugget=", nugget, sep = "")
m.name = paste(m.name, ")", sep = "")
v.level[(i*n+1):((i+1)*n)] = rep(m.name, n)
i = i + 1
}
}
dframe = data.frame(semivariance = data[,2], distance = data[,1],
model = factor(v.level, levels = unique(v.level)))
vgm.panel = function(x,y) {
n = length(x)
lpoints(x[1],y[1])
llines(x[2:n],y[2:n])
}
if (!plot)
dframe
else
xyplot(semivariance ~ distance | model, dframe,
panel = vgm.panel, as.table = TRUE, ...)
}