require(geoR) ## variogramas data(ca20) summary(ca20) par(mfrow=c(2,3), mar=c(3,3,0,0), mgp=c(1.8, .7, 0)) plot(variog(ca20, max.dist=600), ylim=c(0, 150)) plot(variog(ca20, max.dist=600, trend=~area), ylim=c(0, 150)) plot(variog(ca20, max.dist=600, trend=~area+altitude), ylim=c(0, 150)) plot(variog(ca20, uvec=seq(0,600, l=11)), ylim=c(0, 150)) plot(variog(ca20, uvec=seq(0,600, l=21), trend=~area), ylim=c(0, 150)) plot(variog(ca20, uvec=seq(0,600, l=8), trend=~area+altitude), ylim=c(0, 150)) par(mfrow=c(1,1), mar=c(3,3,0,0), mgp=c(1.8, .7, 0)) plot(variog(ca20, trend=~area)) plot(variog(ca20, max.dist=600, trend=~area)) plot(variog(ca20, max.dist=800, trend=~area)) plot(variog(ca20, uvec=seq(0,600, l=6), trend=~area)) plot(variog(ca20, uvec=seq(0,600, l=15), trend=~area)) plot(variog(ca20, uvec=seq(0,600, l=21), trend=~area)) ca20.v <- variog(ca20, uvec=seq(0,600, l=15), trend=~area) ca20.v1 <- variog(ca20, uvec=seq(0,600, l=11), trend=~area) plot(ca20.v) plot(ca20.v, pts.range.cex=c(1, 3)) ca20.v ca20.ef <- eyefit(v) ca20.ef ca20.ef[[1]] ca20.vf <- variofit(ca20.v) ca20.vf1 <- variofit(ca20.v1) lines(ca20.vf, col=2) ca20.vf ca20.vf1 #args(variofit) ca20.vf2 <- variofit(ca20.v, cov.model="matern", kappa=2) ca20.vf2 ## passando valores iniciais ca20.vf <- variofit(ca20.v, ini=c(100, 40), nug=20) ca20.vf ca20.vf <- variofit(ca20.v, ini=ca20.ef[[1]]) ca20.vf