Я пытаюсь провести пространственно-временной анализ некоторых данных, используя krigST из пакета gstat в R. Модель, которую я использую:
sm <- vgmST(stModel="sumMetric",
space=vgm(80,"Exp",30,0),
time=vgm(80,"Exp",30,0),
joint=vgm(80,"Exp",30,0),
stAni=0.00003,method="L-BFGS-B")
После подбора STvariogram, результат:
sill.s : 7.453354e+01 range.s : 6.546394e+01 nugget.s : 3.292897e+01
sill.t : 0.000000e+00 range.t : 3.062291e+01 nugget.t : 2.179542e+00
sill.st : 7.453226e+01 range.st : 6.546647e+01 nugget.st : 0.000000e+00 1
anis : .490116e-08
Однако, когда я выполняю krigST, используя вышеуказанную модель, я получаю следующую ошибку:
Error in chol.default(A) :
the leading minor of order 186 is not positive definite
Может кто-нибудь помочь мне решить проблему, пожалуйста?
Заранее спасибо.