library(sp) library(gstat) #library(RColorBrewer) #library(TeachingDemos) library(maps) #library(mapproj) library(geoR) library(geoRglm) ######################################## help(coordinates) help(spplot) help(variogram) data(meuse) meuse coordinates(meuse) <- ~x+y spplot(meuse,c("cadmium")) mypalette_green<-brewer.pal(7,"Greens") spplot(meuse,c("cadmium"),col.regions=mypalette_green) meuse_var<-variogram(cadmium~1, data = meuse) plot(meuse_var)