Raw File
hana.R
## check, hana

if (EXTENDED.TESTING <- file.exists("source.R")) source("source.R")

x <-  c(0, 1, if (EXTENDED.TESTING) 0.001 else 0.05)
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