Hi,
Thank you for your suggestions. I gave the gstat with weights argument a try,
but I’m not sure about the results / my implementation.
I have worked out an example with the meuse dataset so I can share it with you:
library(gstat)
library(sp)
data(meuse)
coordinates(meuse) <- ~x+y
data(meus
Dear list members,
I am trying to generate uncertainty maps of available water capacity
predicted with a pedotransfer function via Monte-Carlo simulation. I have
the raster layers with the mean and SD of the input variables, so I can
sample from their distributions (I ignore the correlations amo