使用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个最接近的观察（或模拟）值用作条件数据。