# NOT RUN {
# Summarize Data
summary(ContrivedData)
# Empirical semivariance for variable y
raw.var <- semivariance(x=ContrivedData$y,coords = cbind(ContrivedData$s.1,
ContrivedData$s.2))
as.matrix(raw.var)
# Estimation using metropolis.krige()
#' # Set seed
set.seed(1241060320)
M <- 100
#M<-10000
contrived.run <- metropolis.krige(y ~ x.1 + x.2, coords = c("s.1","s.2"),
data = ContrivedData, n.iter = M, range.tol = 0.05)
# Parametric semivariance
parametric.var <- semivariance(contrived.run)
as.matrix(parametric.var)
as.data.frame(parametric.var)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab