gstat包中的chol.default(A)krigeST中的R错误

gstat包中的chol.default(A)krigeST中的R错误,r,kriging,gstat,covariogram,R,Kriging,Gstat,Covariogram,我正在使用一个每小时一次的气温数据集,记录在一个相对较小的区域内约200个站点。我选择了一个时空变异函数(例如总和度量)来拟合我的数据,现在我正试图对我的相同台站进行预测,以填补NA(缺失值)缺口。当对每日聚合数据使用krigeST()函数时,一切似乎都很顺利,但当我以原始的每小时分辨率使用它时,总是会出现以下错误: Error in chol.default(A) the leading minor of order 68 is not positive definite 我在谷歌上搜索了一

我正在使用一个每小时一次的气温数据集,记录在一个相对较小的区域内约200个站点。我选择了一个时空变异函数(例如总和度量)来拟合我的数据,现在我正试图对我的相同台站进行预测,以填补NA(缺失值)缺口。当对每日聚合数据使用krigeST()函数时,一切似乎都很顺利,但当我以原始的每小时分辨率使用它时,总是会出现以下错误:

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

我在谷歌上搜索了一下,发现它与一个不完全正定的矩阵有关。然而,我不知道为什么会发生这种情况,我想知道你们中是否有人知道修复这种情况的方法(一种避免这种情况的解决方法)。

有几种可能导致奇异协方差矩阵。两个共同点:

  • 重复观测(相同的位置和时间戳)
  • 一种变异函数模型,不能充分区分观测值,导致观测值几乎完全相关