# ---- require(geoR) set.seed(11) a = grf(n=95, grid="reg", xlims=c(0,100), ylims=c(0,100), cov.model="matern", mean=50, cov.pars=c(100, 10), nugget=30, nsim=1, messages=TRUE) image(a) persp(a) points.geodata(a) # simulacao com grid irregular set.seed(12) b <- grf(n=95, grid="irreg", xlims=c(0,100), ylims=c(0,100), cov.model="matern", mean=50, cov.pars=c(100, 10), nugget=30, nsim=1, messages=TRUE) points.geodata(b) v <- variog(a, max.dist=70) X11() plot(v) ef = eyefit(v) likfit<-likfit(a, ini=c(87, 21)) likfit v <- variog(b, max.dist=70) X11() plot(v) ef = eyefit(v) likfit<-likfit(b, ini=c(86, 22)) likfit