[R-sig-Geo] Ordinary co-kriging with gstat - computing time/syntax error?

Dear R-sig-Geo members, I want to perform co-kriging on the residuals on two regression models with gstat. I am following some tutorials (Rossiter, 2012; Malone et al., 2017), but I may have a mistake in the script. The computing time seems too long (I know it takes a lot of calculation time, but

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

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