使用ff包进行Kriging仿真

Question

我正在尝试理解我可以使用ff包来克服错误的方法

"Error: cannot allocate vector of size 1.1 Mb"

使用克里金/高斯模拟时。我不知道如何更改输入数据。有什么想法可以帮我这么做吗？

我正在使用gstat包执行模拟，如下所示：

library(sp)
data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m, nsim=1000)

Answer 1

对于高斯模拟，您需要将邻域大小设置为某个值，否则内存消耗将无限增长。你可以这样做，例如通过设置nmax=50，将50个最近观察到的（或模拟的）值用作调节数据。