sessionInfo()
data(meuse)
coordinates(meuse) = ~x+y
meuse$value = log(meuse$zinc)
# set up intamap object:
krigingObject = createIntamapObject(
observations = meuse,
formulaString = as.formula('value~1'),class = "automap")
# do estimation steps:
krigingObject = estimateParameters(krigingObject)
krigingObject = methodParameters(krigingObject)
# Create a new object
krigingObject2 = createIntamapObject(observations = meuse,
formulaString = as.formula('value~1'),
params = list(methodParameters = krigingObject$methodParameters))
krigingObject$variogramModel
krigingObject2$variogramModel
Run the code above in your browser using DataLab