[R-sig-Geo] distances reported by gstat when cokriging

Tom Gottfried tom.gottfried at tum.de
Sat Mar 17 09:01:37 CET 2012


Am 16.03.2012 16:13, schrieb Edzer Pebesma:
>>
>>> When you say, "Is there a way to get the distances for the second
>>> covariate?", do you mean, as R object?
>>
>> Well, that would be luxurious. Currently I copy-paste those data into
>> Calc for further analysis. Of course I can compute the distances somehow
>> outside from gstat, but is there a way to get gstat doing this?
>
> Of course there is (modifying the c code) but my advice would be to take
> a shorter path -- sp::spDists? rgeos::gDistance?

The original goal of my dilettantish copy-pasting was to compare kriging 
weights vs. distance. I think I can solve that more elegantly with your 
code from `demo(weights)' and e.g. spDists.

Thanks!
Tom

-- 
Technische Universität München
Department für Pflanzenwissenschaften
Lehrstuhl für Grünlandlehre
Alte Akademie 12
85350 Freising / Germany
Phone: ++49 (0)8161 715324
Fax:   ++49 (0)8161 713243
email: tom.gottfried at wzw.tum.de
http://www.wzw.tum.de/gruenland



More information about the R-sig-Geo mailing list