data(meuse)
coordinates(meuse) <- ~x+y
m <- vgm(.59, "Sph", 874, .04)
# leave-one-out cross validation:
out <- krige.cv(log(zinc)~1, meuse, m, nmax = 40)
criterio.cv(out)
# multiquadratic function
data(preci)
attach(preci)
# optimal sigma
tab <- rbf.tcv(sigma=1.488733, z=prec, coordinates=preci[,2:3],
n.neigh=9, func="M")
criterio.cv(tab)Run the code above in your browser using DataLab