[R-sig-Geo] Getting values from fitted variogram

Jon Olav Skoien jon.skoien at jrc.ec.europa.eu
Mon Jul 9 16:30:13 CEST 2012


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.



More information about the R-sig-Geo mailing list