I am trying to compute the conventional benefit with RGeostats.
For that, I carry out an anamorphosis of my data and build a program that uses the development of the conventional benefit on the Hermit polynomials as done in the non-linear lesson of Jacques Rivoirard (cours C-139, pp.49-50, 1991). In this lesson, he uses the Rodrigues formulation of the Hermit polynomials that introduces a normalization term.
Is it this formulation used in Rgeostats and if not, is the value of the field psi (coefficients of the anamorphosis) normalized by the factorial of the index i in the sum or not ?
It seems that the selectivity function enables to compute the conventional benefit for a defined dataset but is it possible to estimate the conventional benefit by disjunctive kriging or do I have to implement the method described in Rivoirard,1991 ?
I am going further than my previous post: if I effectively want to estimate the conventional benefit by disjunctive kriging, I will need the coefficients of the anamorphosis.
In my previous version of RGeostats, I have access to these coefficients via the slot psi of the object returned by anam.fit but it does no longer exist in the last version of the library. Has it been replaced by the slot hermite ? If not, is there a way to recover the coeffcients ?