Getting values from fitted variogram
Paolo, variogramLine can give you the semivariogram values from your variogram model for any distance: > vm = fit.variogram(vgm1, vgm(1,"Sph",300,1)) > variogramLine(vm, dist_vector = seq(100,1000,100)) Cheers, Jon
On 09-Jul-12 3:44, Ebrahim Jahanshiri wrote:
hi, if it helps, I think if you put vgm1 itself you will get a table like this :
vgm1
np dist gamma dir.hor dir.ver id
1 57 79.29244 0.1234479 0 0 var1
2 299 163.97367 0.2162185 0 0 var1
3 419 267.36483 0.3027859 0 0 var1
4 457 372.73542 0.4121448 0 0 var1
5 547 478.47670 0.4634128 0 0 var1
6 533 585.34058 0.5646933 0 0 var1
7 574 693.14526 0.5689683 0 0 var1
8 564 796.18365 0.6186769 0 0 var1
9 589 903.14650 0.6471479 0 0 var1
10 543 1011.29177 0.6915705 0 0 var1
11 500 1117.86235 0.7033984 0 0 var1
12 477 1221.32810 0.6038770 0 0 var1
13 452 1329.16407 0.6517158 0 0 var1
14 457 1437.25620 0.5665318 0 0 var1
15 415 1543.20248 0.5748227 0 0 var1
that contains the value of gamma (the average for each distance) and other
information. But these are exprimental semivariances. The models fit will
produce continuous gamma for the prediction purpose that I think are not as
important .
Date: Sat, 7 Jul 2012 07:14:37 -0430
From: Paolo Ramoni Perazzi <rpaolo1967 at gmail.com>
To: r-sig-geo at r-project.org
Subject: [R-sig-Geo] Getting values from fitted variogram
Message-ID:
<CANvsoOsB1kMPpL5r_5_aqTrgJHEKqdR==nznzuq1g5VgSRVqYA at mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1
I'm really sorry if this is a very simple question...
Using gstat and given the code:
data(meuse)
vgm1 <- variogram(log(zinc)~1, ~x+y, meuse)
fit.variogram(vgm1, vgm(1,"Sph",300,1))
I wonder if there is a simple way to get the value of "h"
corresponding to a given "gamma(h)"
Thanks a lot!
Paolo R.
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Jon Olav Sk?ien Joint Research Centre - European Commission Institute for Environment and Sustainability (IES) Land Resource Management Unit Via Fermi 2749, TP 440, I-21027 Ispra (VA), ITALY jon.skoien at jrc.ec.europa.eu Tel: +39 0332 789206 Disclaimer: Views expressed in this email are those of the individual and do not necessarily represent official views of the European Commission.