require(geoR) data(s100) ls() class(s100) str(s100) plot(s100) plot(s100, low=T) data(parana) plot(parana, low=T) points(s100) points(s100, pt.div="rank") points(s100, pt.div="quar") points(s100, pt.div="quar", cex.min=1.3, cex.max=1.3) points(s100, pt.div="quint", cex.min=1.3, cex.max=1.3) args(plot.geodata) args(points.geodata) x11() plot(parana, low=T, trend="1st") require(MASS) boxcox(s100$data+5~1, lam=seq(0,4,leng=101)) summary(s100) v <- variog(s100, max.dist=1) plot(v) ef <- eyefit(v) ml <- likfit(s100, ini=ef) ml gr <- expand.grid(seq(0,1, len=100), seq(0,1, len=100)) points(s100) points(gr, pch=19, cex=0.25, col=2) kc <- krige.conv(s100, loc=gr, krige=krige.control(obj.model=ml)) image(kc, col=terrain.colors(21)) points(s100, add=T) image(kc, col=gray(seq(1, 0, len=21)))