Hi all,
when I write
- Code: Select all
model.create(melem.name(2))
the theoretical range is equal to 1/3 times the practical range :
- Code: Select all
Covariance Part
---------------
- Exponential
Range = 1.000
Theo. Range = 0.334
Sill = 1.000
Total Sill = 1.000
But when the model is created from model.auto, it is the converse (and it seems false) :
- Code: Select all
vario=vario.grid(simtub(dbout=db.create(flag.grid=T,nx=c(100,100)),nbtuba=1000,model=model.create(melem.name(2),ndim=2,range=30)))
model.auto(vario,struct=melem.name(2))
Result :
- Code: Select all
Covariance Part
---------------
- Exponential
Sill = 0.749
Anisotropy = 1.000 1.009
Ranges = 22.238 22.036
Theo. Ranges= 66.619 66.014
Total Sill = 0.749
Furtheremore, the definition of the range parameter seems to have changed for the K-Bessel function.
Could you create a function allowing to compute the range from the scale and vice-versa for all the melem types.
Thanks a lot for your great package.