Spatial

Macam-macam model Variogram

library(gstat)
vgm()
##    short                                      long
## 1    Nug                              Nug (nugget)
## 2    Exp                         Exp (exponential)
## 3    Sph                           Sph (spherical)
## 4    Gau                            Gau (gaussian)
## 5    Exc        Exclass (Exponential class/stable)
## 6    Mat                              Mat (Matern)
## 7    Ste Mat (Matern, M. Stein's parameterization)
## 8    Cir                            Cir (circular)
## 9    Lin                              Lin (linear)
## 10   Bes                              Bes (bessel)
## 11   Pen                      Pen (pentaspherical)
## 12   Per                            Per (periodic)
## 13   Wav                                Wav (wave)
## 14   Hol                                Hol (hole)
## 15   Log                         Log (logarithmic)
## 16   Pow                               Pow (power)
## 17   Spl                              Spl (spline)
## 18   Leg                            Leg (Legendre)
## 19   Err                   Err (Measurement error)
## 20   Int                           Int (Intercept)
show.vgms()

show.vgms(models = c("Exp", "Mat", "Gau"), nugget = 0.1)

# show a set of Matern models with different smoothness:
show.vgms(kappa.range = c(.1, .2, .5, 1, 2, 5, 10), max = 10)
## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

# show a set of Exponential class models with different shape parameter:
show.vgms(kappa.range = c(.05, .1, .2, .5, 1, 1.5, 1.8, 1.9, 2), models = "Exc", max = 10)
## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

# show a set of models with different shape parameter of M. Stein's representation of the Matern:
show.vgms(kappa.range = c(.01, .02, .05, .1, .2, .5, 1, 2, 5, 1000), models = "Ste", max = 2)
## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length

## Warning in v.level[(i * n + 1):((i + 1) * n)] <- rep(m.name, n): number of items
## to replace is not a multiple of replacement length