## Not run:
# ### load data
# data(meuse,meuse.blocks)
#
# ### plot blocks
# plot(meuse.blocks)
#
# ### compute the approximated block variance of each block in
# ### meuse.blocks without the definition of neighbours blocks (default)
# preCK_1 <- preCKrige(newdata = meuse.blocks,
# model = covmodel("exponential", 0.05, 0.15, scale = 192.5),
# pwidth = 75, pheight = 75)
#
# ### plot block approximation of block 59
# plot(preCK_1, 59)
#
#
# ### define neighbours
# if(require(spdep))
# {
# neighbours <- poly2nb(meuse.blocks)
# class(neighbours)
# ### neighbours should be an object of the class "list"
# class(neighbours) <- "list"
# ### compute the approximated block variance-covariance matrices of each block in
# ### meuse.blocks without the defined block neighbours
# preCK_2 <- preCKrige(newdata = meuse.blocks, neighbours = neighbours,
# model = covmodel("exponential", 0.05, 0.15, scale = 192.5),
# pwidth = 75, pheight = 75)
#
# ### plot block approximation of block 59 and its
# ### block neighbours
# plot(preCK_1, 59)
# }
# if(!require(spdep))
# {
# cat("Please, install the package spdep to excute this example.\n")
# }
#
# ## End(Not run)
Run the code above in your browser using DataLab