Обычный Кригинг в 2D и 3D: проблемы с дублированием?
Я хочу криге в 2-х мерной и 3-х мерной области (отдельно). У меня есть много моментов, которые были воспроизведены ниже.
НОТА:
Атрибут, который я хотел бы обозначить, - это Y. Они намеренно отрицательны, поскольку являются значениями экспонент.
Домен является трехмерным, поэтому с точки зрения xy некоторые точки кажутся совмещенными.
# define kriging domain
x <- sample(seq(0,300,by=5),10,replace=FALSE)
y <- sample(seq(0,300,by=5),10,replace=FALSE)
z <- sample(seq(0,30,by=0.5),100,replace=TRUE)
Y <- sample(seq(-3.0,-0.01,by=0.001),100*10*10,replace=TRUE)
(
mydata = data.frame(x=rep(x,each=10),y=rep(y,times=10),z=z,Y=Y)
grid = data.frame(x=rep(x,times=length(x)),y=rep(y,each=length(y)))
# 2_D Kriging
ok_2D <- krige(formula=Y~1, locations=~x+y, data=mydata, newdata=grid, model=fit.exp)
# 3_D Kriging
ok_3D <- krige(formula=Y~1, locations=~x+y+z, data=mydata, newdata=grid, model=fit.exp)
Ошибка:
> # 2_D Kriging
> ok_2D <- krige(formula=Y~1, locations=~x+y, data=mydata, newdata=grid, model=fit.exp)
[using ordinary kriging]
"chfactor.c", line 131: singular matrix in function LDLfactor()
Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, :
LDLfactor
>
> # 3_D Kriging
> ok_3D <- krige(formula=Y~1, locations=~x+y+z, data=mydata, newdata=grid, model=fit.exp)
Error in .subset(x, j) : only 0's may be mixed with negative subscripts
>
Мой вопрос таков: почему существует проблема с единственной матрицей? Я понимаю, что расположенные точки не должны быть проблемой при кригинге, так как веса будут присваивать нулевой вес одной из расположенных рядом точек.
Или я пропускаю еще одну проблему? Благодарю.