[R-sig-Geo] Generating uncertainty maps with Monte-Carlo simulation and large rasters

2016-10-11 Thread Mercedes Román
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

Re: [R-sig-Geo] Kriging with uncertain data

2016-10-11 Thread Santiago Beguería Portugués
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