https://github.com/cran/RandomFields
Raw File
Tip revision: a88578bf75299850a5f3ff3efba8d520c28d6a81 authored by Martin Schlather on 26 April 2005, 00:00:00 UTC
version 1.2.10
Tip revision: a88578b
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)

back to top