R/set.R In gstat: Spatial and Spatio-Temporal Geostatistical Modelling, Prediction and Simulation

```# \$Id: set.q,v 1.7 2008-12-16 14:59:22 edzer Exp \$
gstat.set <- function(set) {
if(!is.list(set))
stop("set should be a list")
if (length(set) == 0)
return(NULL)
ret = NULL
n = names(set)
for (i in (1:length(set))) {
val = set[[i]]
if (n[i] == "method")
str = paste("method: ", val, ";", sep="")
else {
if (is.character(val))
val = paste("'", val, "'", sep = "")
str = paste("set ", n[i], " = ", val, ";", sep="")
}
ret = c(ret, str)
}
ret
}

if(!is.list(set))
stop("set should be a list")
if (length(set) == 0)
return(NULL)
n = names(set)
for (i in (1:length(set))) {
if (n[i] == "method")
ret = .Call(gstat_set_method, set[[i]])
else
ret = .Call(gstat_set_set, n[i], set[[i]])
}
invisible(ret)
}

if (is.character(obj\$merge) && length(obj\$merge) == 2)
obj\$merge = list(c(obj\$merge[1], 1, obj\$merge[2], 1))

if (!is.list(obj\$merge))
stop("merge argument should be list or character vector of lenght 2")

ret = NULL
for (i in 1:length(obj\$merge)) {
m = obj\$merge[[i]]
if (is.character(m) && length(m) == 4) {
id = match(m[c(1,3)], names(obj\$data)) - 1 # name ->> id
if (any(is.na(id)))
stop(paste("could not match all ids:", m[c(1,3)]))
col = as.integer(m[c(2,4)]) - 1
if (any(is.na(col)) || any(col < 0))
stop("merge: parameters should be positive integers")
ret = .Call(gstat_set_merge, id[1], col[1], id[2], col[2])
} else stop(
"list elements of merge should be lenght 4 character vectors")
}
ret
}
```

Try the gstat package in your browser

Any scripts or data that you put into this service are public.

gstat documentation built on April 6, 2023, 5:21 p.m.