RFoptions(seed=0)
model <- RMgauss() + RMnugget(var=0.01)
z <- RFsimulate(model=model, 0:10, 0:10, grid=TRUE, n=4)
plot(z)
# gives the same model as
model <- RMplus(RMgauss(), RMnugget(var=0.01))
z <- RFsimulate(model=model, 0:10, 0:10, grid=TRUE, n=4)
plot(z)
RFoptions(seed=NA)
Run the code above in your browser using DataLab