{RGeostats}R Documentation

Automatic Model Fitting


Automatic Model Fitting

Usage, struct =,4,5,2,3)), constraints=NA,
      auth.aniso = TRUE, auth.rotation = TRUE, 
      auth.locksame = FALSE, auth.lock2d = FALSE, flag.goulard=TRUE,
      flag.noreduce = FALSE, flag.norm.sill=FALSE,flag.keep.intstr=FALSE,
      param = NA, lower = NA, upper = NA, equal = NA, properties=NA,
      draw = TRUE, wmode=2, maxiter = 1000, verbose=0, tolstop = 1.e-5, 
      epsdelta = 1.e-5, tolsigma = 5, delta=1, ...)



The vario-class containing the set of experimental variograms


List of basic structures to be used. It can be given as a vector of character strings or as a vector of indexes. In the latter case, the basic structures will be the result of


Vector giving the constraints on the cumulative sills of the model. This vector must be dimensionned to the number of variables in the variogram. All the defined terms must be positive (or undefined). Otherwise, no constraints is applied on the Sill estimation. This is incompatible with the 'flag.norm.sill' flag.


When TRUE, the anisotropy will be checked


When TRUE, the fit will look for rotation search in anisotropy. However note that an initial rotation is always defined (based on the calculation directions of the experimental variogram). Therefore when this switch is OFF, the initial rotation is maintained.


When TRUE and if an anisotropy is allowed (auth.aniso), all the basic structures should share the same rotation.


When TRUE, the anisotropy rotation will only consider a rotation around the Z-axis. The Dip and Plunge angles are kept equal to zero.


When TRUE, the Goulard algorithm is used for the determination of the Sills (quicker). Otherwise Sills are infered using the standard Foxleg algorithm.


When TRUE, the useless basic structures must not be discarded.


When TRUE, the automatic fit must fulfill the constraint that the sum of the sills must always be equal to 1. This flag is incompatible with the definition of 'constraints'.


When TRUE, the fitting procedure must always maintain at least one intrinsic basic structure. If no intrinsic basic structure is provided in the initial template, the automatic fitting procedure fails.


List of initial values for the parameters to be fitted. For details see model.param.define.


List of lower bounds for the parameters to be fitted For details see model.param.define.


List of upper bounds for the parameters to be fitted For details see model.param.define.


List of constant values assigned to the parameters to be fitted For details see model.param.define.


This argument defines the transformation to be applied to the model. For more details, see properties.define.


When TRUE, both the experimental variograms and the Model are represented graphically. Otherwise, no graphic is produced.


Type of the weighting function used in the fitting procedure. This function is defined in the case of several directional experimental variograms, calculated in a multivariate case:

  • 0: The weight is constant

  • 1: The weight is proportional to the number of pairs

  • 2: The weight is proportional to the number of pairs and inverse proportional to the distance

  • 3: The weight is inverse proportional to the number of lags for each direction


Maximum number of iterations


Verbose option

  • 1 Indications are given regarding the convergence of the algorithm used for automatic model fitting.

  • 0 No message is printed except the warning when convergence is not reached

  • -1 Silent mode


Tolerance for the distance between the experiment and the model.


Tolerance for the direction increment used in the gradient search


Percentage of the variance below which a basic structure will be discarded


Delta value


Arguments passed to vario.plot and model.plot.


The model-class with fitted parameters

When needed, you can access to each matrix of sills by using the keypair mechanism using the keyword "Fitted_Sill_'i'" where 'i' stands for the rank of the each basic structure (starting from 1).

You can also access to the eigen decomposition of the matrix of Sills of rank 'i' where 'i' stands for the rank of each basic structure (starting from 1), using:


Function originally funded by Geovariances and ANR CRISCO2.


Nicolas Desassis and Didier Renard (Geosciences/Mines-ParisTech)


Desassis, N., and Renard, D. (2013) Automatic Variogram Modeling by Iterative Least Squares: Univariate and Multivariate Cases. Mathematical Geosciences 45 (4) pp 453-470.


	title="Experimental Variogram in 2-D space to be Fitted")
model <-,maxiter=100,
        title="Experimental variogram in 2-D space")

[Package RGeostats version 11.0.6 Index]