VgmSemivariance: \(\gamma(0) = 0\), \(\gamma(h) = \sum psill\,g(h)\)
plus nuggets, with the unit-sill shapes of gstat vgm (Nug, Exp, Sph, Gau,
Exc, Mat, Ste, Cir, Lin, Bes, Pen, Per, Wav, Hol, Log, Pow) and Cub
(cubic), Cau (Cauchy), JBes (J-Bessel), Dmp (damped oscillation), Cos
(Chiles and Delfiner 2012). VgmCovariance and VgmCorrelogram
give \(C(0) - \gamma(h)\) and \(1 - \gamma(h)/C(0)\) for bounded models.
Identical to the Python arm morie.fn.vgmods.
Usage
VgmSemivariance(h, model)
VgmCovariance(h, model)
VgmCorrelogram(h, model)
Arguments
- h
Distance(s).
- model
Component list or list of nested components.
References
Chiles, J.-P. and Delfiner, P. (2012). Geostatistics: Modeling
Spatial Uncertainty, 2nd edn. Wiley.
Examples
VgmSemivariance(0.9, list(model = "Pen", psill = 1, range = 2))
#> [1] 0.7367636
VgmCovariance(0.9, list(model = "Cub", psill = 2, range = 2))
#> [1] 0.6361229
VgmCorrelogram(1, list(model = "Exp", psill = 3, range = 1))
#> [1] 0.3678794
