Thank you Johan for the clarification. Given this limitation with gstat, I was wondering if there is any package in R or so that supports kriging with anisotropy modelling using great circle distance. Or is gstat going to be developed to overcome this limitation? Thank you in advance! Kind regards, Jin
-----Original Message----- From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org] On Behalf Of Johan Van de Wauw Sent: Saturday, 20 August 2011 12:49 AM To: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Kriging wiht anisotropy modeling using great circle distance in gstat [SEC=UNCLASSIFIED] On Fri, Aug 19, 2011 at 3:37 AM, <jin...@ga.gov.au> wrote: > Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, : > gstat: value not allowed for: for long/lat data, anisotropy cannot be defined >> > > Is there anything wrong in my code or it is a bug? Any suggestions are > appreciated! > As you can read from the error message, nothing is wrong with your code or with gstat, but anisotropy is not supported when you use great circle distances. Sincerely Yours Johan _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo