## [SOLVED] units problem for variogram automatic fit

Here you can report the problems encountered when using a function relative to the variography:
- experimental variography (vario.calc, vario.grid, vmap.calc, ...)
- model fitting (model.auto)

### [SOLVED] units problem for variogram automatic fit

When making an automatic fit to a variogram with units say in e-07, the automatic fit provides non relevant outputs.
when I change the units of Z and thus those of the variogram to get "reasonnable" units, says between 0 and 1, then the automatic fit is correct.
bez

Posts: 18
Joined: Tue Oct 23, 2012 10:31 am

### Re: units problem for variogram automatic fit

Hi Nicolas,
Thank you for your feedback. When you talk about units, is it spatial units or variable units?

Nicolas
Nicolas Desassis

Posts: 196
Joined: Thu Sep 20, 2012 4:22 pm

### Re: units problem for variogram automatic fit

Hi Nicolas,
units of the variables
A+
Note : I had sent you an email with data and script to illustrate the problem.
bez

Posts: 18
Joined: Tue Oct 23, 2012 10:31 am

### Re: units problem for variogram automatic fit

When changing the scale of the variable of interest (say dividing it by 1000), the variance (and hence the units along the variogram axis) are divided by 1,000,000.
Then model.auto procedure must follow the same scale change.
If we consider 'diff' the gap between the experimental variogram and the model; during the iterations, if 'diff' becomes smaller than 'tolstop' the iterations stop. It is considered that the convergence has been reached.
Unfortunately this parameter 'tolstop' is not scaled as a function of the sill or the variance of the data. Then the user must adapt the value of 'tolstop' to the order of magnitude of the variance of the target variable.

In your case, due to very small variability of the data, the value of 'diff' was considered as smaller than 'tolstop' after a very limited number of iterations (say 2). The program stops while the fit is not correct.

The workaround is simply to reduce the value of 'tolcode'. Put it to 1.e-10 for example.
Didier Renard

Posts: 266
Joined: Thu Sep 20, 2012 4:22 pm

### Re: units problem for variogram automatic fit

Hi Didier
Fair enough. I understand the problem.
I did not pay attention to the fact that tolstop was in absolute distance and not relative to the variance.
Wouldn't it be more relevant to set the tolstop limit in terms of % of the variance ?
It would then be more robust to various situations.
Nicolas
bez

Posts: 18
Joined: Tue Oct 23, 2012 10:31 am

Return to Variography

### Who is online

Users browsing this forum: No registered users and 0 guests