I am looking for someone with experience in how to handle the anisotropy-parameters in the
likfit() function which is part of the GeoR package in R.
I am using
likfit() to generate the
obj.m-parameter for the function
krige.conv(). The data I use consists of scattered points (5 to 30) on a ~50*50 grid.
likfit(), I want the parameters
psiR to be estimated from the available data. This works fine by stating:
fix.psiA = FALSE, fix.psiR = FALSE
However, the estimated values show quite a huge range for
psiR, given that the datasets I used are not fundamentaly different from each other. (It is a set of soil moisture measurements with values from 0 to ~45.)
psiR ranges from 1 to about 8000... can that be right? Most of the values are in a range from 1 to 10, but I can not tell why some datasets produce these very large values for
I was unable to find further information on what this parameter exactly does. I do understand that this parameter regulates the dependency of an estimated value according to it's location relative to the ambient measured values. But I do not know in which way this is accomplished.
I am sorry for the possibly bad english, not a native speaker. I am also sorry for not posting my code, but it is quite long and I felt that it is not necessary in behalf of my question.
Thank you very much for your interest. I will post additional details, if required.