来自gstat包的chol.default（A）krigeST中的R错误

Question

我正在处理一个每小时的气温数据集，记录在相对较小区域的~200个站点。我选择了一个时空变异函数（例如sum-metric）来拟合我的数据，现在我试图在同一个站点进行预测，以填补NA（缺失值）的空白。当在每日聚合数据上使用krigeST（）函数时，一切似乎都很顺利，但是当我在原始的每小时分辨率下使用它时，我总会得到以下错误：

Error in chol.default(A)
the leading minor of order 68 is not positive definite

我用谷歌搜索它，发现它与一个不完全正定的矩阵有关。但是，我不确定为什么会发生这种情况，并且想知道是否有人知道解决这个问题的方法（避免它的解决方法）。

Answer 1

有几种可能性导致奇异的协方差矩阵。两个常见的：

1. 重复观察（相同的位置和时间戳），
2. 变异函数模型，不能充分区分观察结果，导致近乎完全相关的观察结果。