Thanks! Sorry, I had that problem because of the error message: "Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, : gstat: value not allowed for: for long/lat data, anisotropy cannot be defined"
Now I guess this message is displayed when using decimal degree as coordinates. Just later I understood that indeed metric coordinates are accepted. In fact it worked really smoothly after having converted to metres. Somewhere I should have skipped to read that metric coordinates should be used. Thanks again, Dario -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/gstat-kriging-anisotropy-tp7582312p7582316.html Sent from the R-sig-geo mailing list archive at Nabble.com. _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
