Revision b2f6358d2679f8a415bd4f01a706fc38b04bc0c5 authored by Martin Schlather on 26 April 2005, 00:00:00 UTC, committed by Gabor Csardi on 26 April 2005, 00:00:00 UTC
1 parent dec8043
xx.R
## 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)
Computing file changes ...
