vario.grid {RGeostats}R Documentation

Computes experimental variograms on a grid

Description

Computes various types of experimental variograms on a regular grid.

Usage

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

Arguments

db

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.

dirincr

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.

nlag

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

calcul

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

  • order4 : Order-4 Variogram

  • trans1 : Cross-to-Simple variogram G12/G1

  • trans2 : Cross-to-Simple variogram G12/G2

  • binormal : Cross-to-Simple variogram G12/sqrt(G1*G2)

opt.code

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

tolcode

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

dates

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

  • ...

means

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

vars

Array corresponding to the variance-covariance matrix

flag.indicator

When TRUE, the variable which corresponds to the first (and unique) Z-locator is converted beforehand into its indicators. The principle is described next.

The target variable is first converted into the lower closest integer value. Then the list of these different integer values is established and sorted in ascending order. Finally the initial variable is converted into indicators of these integer values.

vario.add

When defined, the current variogram computation is concatenated to the contents of the current variogram.

Value

An experimental vario-class structure.

See Also

db-class, vario-class,


[Package RGeostats version 14.0.10 Index]