library(sp)
library(gstat)
data(meuse.all)
coordinates(meuse.all) <- ~ x + y
f <- function(x) -x ## log-prior for exponential distribution for the smoothing parameter
## Draw 100 samples from the posterior of eta given the data y.
OUTPUT <- DSSP(
formula = log(zinc) ~ 1, data = meuse.all, N = 100,
pars = c(0.001, 0.001), log_prior = f
)
plot(OUTPUT, contour_plots = FALSE)
Run the code above in your browser using DataLab