gstat包中的chol.default(A)krigeST中的R错误
我正在使用一个每小时一次的气温数据集,记录在一个相对较小的区域内约200个站点。我选择了一个时空变异函数(例如总和度量)来拟合我的数据,现在我正试图对我的相同台站进行预测,以填补NA(缺失值)缺口。当对每日聚合数据使用krigeST()函数时,一切似乎都很顺利,但当我以原始的每小时分辨率使用它时,总是会出现以下错误: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 我在谷歌上搜索了一
Error in chol.default(A)
the leading minor of order 68 is not positive definite
我在谷歌上搜索了一下,发现它与一个不完全正定的矩阵有关。然而,我不知道为什么会发生这种情况,我想知道你们中是否有人知道修复这种情况的方法(一种避免这种情况的解决方法)。有几种可能导致奇异协方差矩阵。两个共同点: