## 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]))