https://github.com/cran/fields
Tip revision: 8eab500c3dad2103092ff68706417414fe53e16b authored by Doug Nychka on 22 September 2009, 20:23:49 UTC
version 6.01
version 6.01
Tip revision: 8eab500
sim.krig.image.R
# fields, Tools for spatial data
# Copyright 2004-2007, Institute for Mathematics Applied Geosciences
# University Corporation for Atmospheric Research
# Licensed under the GPL -- www.gpl.org/licenses/gpl.html
"sim.krig.image" <- function(out, nreps = 10) {
temp1 <- list()
temp2 <- as.list(1:nreps)
class(temp1) <- "sim.krig.image"
temp1$call <- out$call
temp1$grid <- out$grid
hold2 <- look <- out$surface$z
temp <- krig.image.parameters(out)
rho <- temp$rho
shat.MLE <- temp$shat.MLE
for (k in 1:nreps) {
look <- sqrt(rho) * sim.rf(out$cov.obj)
Y <- c(look[out$index]) + rnorm(out$N) * shat.MLE
hold2 <- krig.image(out$x, Y, lambda = out$lambda, cov.function = out$cov.function,
cov.obj = out$cov.obj, kmax = 100, start = NULL)$surface$z
temp2[[k]] <- look - hold2 + out$surface$z
NULL
}
temp1$out <- temp2
return(temp1)
}