[SOLVED] Extracting theoretical range value

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] Extracting theoretical range value

Postby Thibault Cariou » Fri Apr 12, 2019 8:55 am

Hi everyone,

I was wondering if there was a way to extract the theoretical range that appears in the output of model.auto().
When I look at the R object only the practical range is given.

I understand that I could try to calculate the theoretical range but I'm using a Matérn variogram model which makes it difficult.

Thanks!

Thibault.
Last edited by Thibault Cariou on Fri Apr 19, 2019 1:34 pm, edited 1 time in total.
Thibault Cariou
 
Posts: 3
Joined: Thu Apr 11, 2019 2:54 pm

Re: Extracting theoretical range value

Postby Didier Renard » Tue Apr 16, 2019 8:30 am

Hello

As you noticed, the contents of the R object which characterizes the Model (of RGeostats) only contains the practical range.
This was a initial decision of the programming team. As geostatisticians, we also know that we need to access
to the theoretical range ... sometimes.
As the information is not stored in the internal structure, the access to theoretical range has been added in the following way:

1) During the creation phase: we can define the structure either by practical or theoretical range:
Code: Select all
qq1 = model.create(vartype="K-Bessel",param=1.2,sill=4,scale=2)


gives the same result as:

Code: Select all
qq2 = model.create(vartype="K-Bessel",param=1.2,sill=4,range=7.9843)


To double check my statement, just print qq1 and qq2.

2) When the structure exists, created using either one or the other definition sentences, we can use the assessors to get either the practical range:

Code: Select all
qq1[1]$range


or
Code: Select all
qq1[1]$scale


Hope this will help.
Didier Renard
 
Posts: 337
Joined: Thu Sep 20, 2012 4:22 pm

Re: Extracting theoretical range value

Postby Thibault Cariou » Fri Apr 19, 2019 1:33 pm

Hi,

Thank you for your answer.
It has given me the information that I wanted!

Best regards,

Thibault
Thibault Cariou
 
Posts: 3
Joined: Thu Apr 11, 2019 2:54 pm


Return to Variography

Who is online

Users browsing this forum: No registered users and 3 guests

cron