vario.grid {RGeostats}R Documentation

Computes experimental variograms on a grid


Computes various types of experimental variograms on a regular grid.


vario.grid(db, dirincr = NA, nlag = NA, calcul = "vg", 
           opt.code=0, tolcode=0, dates=NA, means=NA)



The db-class (organized as a grid) containing the data information used to calculate the experimental variogram. The variograms are calculated for the set of "z*" variables present in the DB. When present (defined in the Db), the field "w" gives the weights attached to each sample. It is compulsory only for transitive covariogram calculations.


A matrix giving the grid indices for the different directions in which the variograms must be calculated. When absent or in case of error, one direction is generated for each dimension of the space. When provided as a matrix, its first dimension must be equal to the space dimension, its second dimension gives the number of directions to be calculated. When provided as a list or a single element, as many directions are generated as there are elements in the list: the element in the list is considered as the angles of a rotation in the XoY plane.


A vector containing the number of lags per direction. If a constant is provided, it will serve for all the directions.


Character string giving the type of structure to compute:

  • vg : Variogram

  • cov : Covariance (Centered)

  • covnc : Non-centered Covariance

  • covg : Transitive Covariogram

  • mado : Madogram (Variogram of order 1)

  • rodo : Rodogram (Variogram of order 1/2)

  • poisson : Variogram of a Poisson-weighted variable

  • general1 : Generalized variogram of order 1

  • general2 : Generalized variogram of order 2

  • general3 : Generalized variogram of order 3


Array containing the option concerning the use of the code for each direction. This option defines the strategy for the code usage when considering a pair of samples:

  • 0 : a pair of samples is selected whatever their codes (if active)

  • 1 : a pair is retained if the code is active in the data base and the codes of the two samples are closer than tolcode

  • 2 : a pair is retained only if the codes of the two samples are different


Array giving the maximum distance between the codes of the two samples for each direction.


When defined and if the 'time' locator is defined, this value corresponds to the bounds of the Date Intervals. It is organized as follows:

  • The lower bound of the first Date Interval

  • The upper bound of the first Date Interval

  • The lower bound of the second Date Interval

  • The upper bound of the second Date Interval

  • ...


Array corresponding to the estimated means of the variable using for the Poisson variogram calculation


An experimental vario-class structure.

See Also

db-class, vario-class,

[Package RGeostats version 11.2.2 Index]