如何在顺序高斯模拟中控制采样?
例如,在以下代码中,我如何保证蒙特卡罗样本不是双重的?
library(sp)
library(gstat)
data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(zinc~1, meuse, meuse.grid, model = m, nsim=10, debug.level=-1)
我知道示例函数中的默认设置是“replace = FALSE”而不是替换生成的样本。 SGS是一样的吗?