使用ff包进行克里格模拟

Kriging simulation using ff package

我正在尝试了解使用 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)

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