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.