## source("xx.R")
q()

library(RandomFields)
source("/home/schlather/article/R/NEW.RF/RandomFields/R/ShowModels.R")

data(soil, package="geoR")
x <- soil$Linha
y <- soil$Coluna
z <- soil$pHAgua

#x <- 1:5
#y <- 1:8
#z <- runif(length(x) * length(y))
#xy <- as.matrix(expand.grid(x,y))
#

## xset b off
## gauss

bin <- c(-1, seq(0, 20, len=40))
(emp0 <- EmpiricalVariogram(x, y, data=z, grid=FALSE, bin=bin))
(emp1 <- EmpiricalVariogram(unique(y), unique(x), data=z, grid=TRUE, bin=bin))

by <- 0.15
RFparameters(Print=1, CE.force=TRUE, CE.trials=1, CE.userfft=TRUE)
ShowModels(seq(min(y), max(y), by=by), seq(min(x), max(x), by=by), emp=emp1,
           me="ci", fixed.rs=TRUE)