估计空间模型的各向异性参数的问题

机器算法验证 r 空间的
2022-04-02 17:47:23

我正在寻找在如何处理函数中的各向异性参数方面有经验的人,该likfit()函数是 R 中GeoR包的一部分。我likfit()用来生成obj.m函数的 -parameter krige.conv()我使用的数据由 ~50*50 网格上的散点(5 到 30)组成。

likfit()中,我希望参数psiApsiR从可用数据中进行估计。通过说明以下内容可以正常工作:

fix.psiA = FALSE, fix.psiR = FALSE

psiR然而,考虑到我使用的数据集彼此之间并没有根本不同,估计值显示出相当大的范围。(它是一组土壤水分测量值,值从 0 到 ~45。)

psiR范围从 1 到大约 8000... 可以吗?大多数值都在 1 到 10 的范围内,但我不知道为什么有些数据集会为psiR.

我无法找到有关此参数确切作用的更多信息。我确实了解此参数根据估计值相对于环境测量值的位置来调节估计值的依赖性。但我不知道这是通过什么方式完成的。

我很抱歉没有发布我的代码,但它很长,我觉得没有必要代表我的问题。

非常感谢您的关注。如果需要,我将发布更多详细信息。

1个回答

简而言之,用这些稀疏的数据识别各向异性是没有希望的

有问题的两个参数psiApsiR描述了各向异性(“几何各向异性”的角度和比率分别为:有关详细信息,请咨询GSLIBJournel & Huijbregts,因为Diggle & Ribeiro JrgeoR中的文档在各向异性方面确实不足)。对于相对较少的数据点,很有可能——事实上,对于土壤数据(众所周知,它是可变的)——在某些方向上几乎没有检测到空间相关性,而在其他方向上似乎存在一些相关性。这可能导致近乎无限的比率。此外,如果只有一个方向的趋势并且没有被消除,那么这种趋势将产生强烈的各向异性。

您的问题是,对于任何类型的参数估计,五个点都太少30仍然太少而无法可靠地识别各向异性。文献中的经验法则建议您至少需要30100点只是为了开始估计参数和计算预测(即克里金法)。(所有的经验法则都有例外,但听起来这些数据不足以证明是合格的。)如果你不假设一个各向同性模型,你需要在至少四个基本方向上探索方向变异函数,每个这样的变异函数都是基于大约510每个点,这又太小了。要识别各向异性,需要关于100点。

解决方法是施加各向同性变异函数(或从独立于数据的考虑确定各向异性)并希望获得最好的结果。预计预测误差会很大。