[R-sig-Geo] Kriging wiht anisotropy modeling using great circle distance in gstat [SEC=UNCLASSIFIED]

[Jin.Li at ga.gov.au]

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-bounces at r-project.org [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Johan Van de Wauw
Sent: Saturday, 20 August 2011 12:49 AM
To: r-sig-geo at 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.Li at 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 at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo