実験的なバリオグラムを球状モデルとして選択したときに、どの方程式のセットを使用して実験的バリオグラムを適合させるのか疑問に思っていました。たとえば、球状モデルの次のパラメーターを使用しました。
model psill range
1 Nug 5384.346 0.0000
2 Sph 9042.061 944.9017
ここで、シル分散、c = 5384.34 + 9042.06 = 14426.41; および範囲、a = 944.9017。球体モデルは次の式で説明されます (書名: Geostatistics for Environmental Scientists)。
variogram = c*{(3h/2a)-0.5(h/a)^3} when, h<=a
variogram = c when, h>a
ここで、h = 距離。したがって、上記の式を使用して既知の距離 'h' を指定して値 'a' と 'c' をプラグインすると、座標 (h, バリオグラム) がモデルに収まるはずです。残念ながら、h の値が異なる場合はそうはなりません。「gstat」パッケージで異なるモデル (球面、指数など) の場合に使用される方程式のセットを知ることは可能ですか?