如果删除数据帧列，autoKrige会给出不同的结果

Question

我在这里可能做错了什么，但是我发现如果我通过删除不相关的列来简化我的数据框架，自动化库中的autoKrige函数会给出不同的结果。我在自动程序库中用meuse数据再现了这个问题。

library(automap)
data(meuse)
colnames(meuse)
 [1] "x"       "y"       "cadmium" "copper"  "lead"    "zinc"    "elev"   
 [8] "dist"    "om"      "ffreq"   "soil"    "lime"    "landuse" "dist.m" 

coordinates(meuse) =~ x+y
data(meuse.grid)
gridded(meuse.grid) =~ x+y
kriging_result_01 = autoKrige(zinc~1, meuse)
plot(kriging_result_01)

meuse <- NULL
data(meuse)
meuse <- meuse[, c(1,2,6)]

coordinates(meuse) =~ x+y
data(meuse.grid)
gridded(meuse.grid) =~ x+y
kriging_result_02 = autoKrige(zinc~1, meuse)
plot(kriging_result_02)

identical(kriging_result_01, kriging_result_02)
[1] FALSE

情节也各不相同。

这是预期的行为吗？

谢谢，比尔

Answer 1

问题不是移除的柱子，而是网格上应用的拟合模型。使用您的示例，您可以看到：

library(automap)
set.seed(42)

data(meuse)
coordinates(meuse) =~ x+y
data(meuse.grid)
gridded(meuse.grid) =~ x+y
kriging_result_01 = autoKrige(zinc~1, meuse)

meuse <- NULL
data(meuse)
meuse <- meuse[, c(1,2,6)]
coordinates(meuse) =~ x+y

kriging_result_02 = autoKrige(zinc~1, meuse)


kriging_result_01$krige_output@grid
                            x1           x2
cellcentre.offset 178635.12844 329744.86186
cellsize              32.93423     32.93423
cells.dim             84.00000    118.00000

和

kriging_result_02$krige_output@grid
                            x1           x2
cellcentre.offset 178614.42379 329741.35016
cellsize              32.93423     32.93423
cells.dim             85.00000    118.00000

如您所见，网格略有不同，分别为4999和4996个网格点。

如果您使用

kriging_result_01 = autoKrige(zinc~1, meuse, meuse.grid)
kriging_result_02 = autoKrige(zinc~1, meuse, meuse.grid)

kriging_result_01$krige_output@grid
                       x      y
cellcentre.offset 178460 329620
cellsize              40     40
cells.dim             78    104

kriging_result_02$krige_output@grid
                       x      y
cellcentre.offset 178460 329620
cellsize              40     40
cells.dim             78    104

网格是相似的。