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

## library(RandomFields, lib="~/TMP")


if (file.exists("source.R")) source("source.R")

x <- 1:1000
x <-  c(0, 1, 0.001)
model <- "fractalB"
param <- c(0,1,0,1, 1.70)
RFparameters(Storing=FALSE, Print=5)
#plot(seq(x[1], x[2], x[3]),
#     z<-GaussRF(x, model=model, param=param, gridtriple=TRUE))
e <- EmpiricalVariogram(x, data= GaussRF(x, model=model, param=param, n=1000,
                             gridtriple=TRUE),
                        bin=(0.5 + (0 : 20)) / 40, gridtriple=TRUE)
plot(e$c, e$e - param[2] * (e$c / param[4])^param[5])



plot(e$c, e$e / ( param[2] * (e$c / param[4])^param[5]))
back to top