[R-sig-Geo] How to draw two variogram curves in one plot
edzer.pebesma at uni-muenster.de
Thu Sep 4 10:18:33 CEST 2014
On 09/04/2014 08:06 AM, Edzer Pebesma wrote:
> On 09/04/2014 04:57 AM, lwwhdz wrote:
>> Dear all
>> I am doing some work on regression kriging in R. Now I want to plot the varigrom curves of the dependent variable (A) and the regression residuals (B) before kriging. With this plot, I can easily to compare the autocorrelation (e.g. sill, range & nugget) of A and B. I tried the following two commond, but it failed.
>> plot(variogram(A~1,data), vgm_a, pch = 19, col = "black", lwd=2, add=TRUE)
>> plot(variogram(B~1,data), vgm_b, pch = 19, col = "black", lwd=2)
>> Is there anyone have any ideas? Thanks very much!
> Yes, compose incrementally using base plot:
> a = variogram(A~1,data)
> b = variogram(B~1,data)
> plot(gamma ~ dist, a, pch = 19)
> points(gamma ~ dist, b, pch = 20) # use another pch!
> lines(variogramLine(vgm_a, max(a$gamma)))
> lines(variogramLine(vgm_b, max(a$gamma), col = 2))
this should of course have been max(a$dist) instead of max(a$gamma), and
even better the maximum of xlim, if you specify that.
> in the plot() command you probably want to specify xlim and ylim in
> order to have (0,0) as the plot origin.
>> Wang Li
>> Phd Candidate in University of Chinese Academy of Sciences (UCAS) and Institute of Remote Sensing and Digital Earth, CAS,The State Key Laboratory of Remote Sensing Science
>> CAS Olympic S&T Park, No. 20 at Datun Road, Chaoyang, Beijing, P.R.C. (P. O. Box 9718)
>> lwwhdz at sina.com
>> lwwhdz at gmail.com
>> [[alternative HTML version deleted]]
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251
83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 490 bytes
Desc: OpenPGP digital signature
More information about the R-sig-Geo