melem.init {RGeostats}R Documentation

Initialize a Basic Structure

Description

Initialize a new melem-class

Usage

melem.init(ndim=2, nvar=1, vartype=melem.name(3), range=NA, scale=NA, sill=NA,
           param=1, order=0, flag.aniso=FALSE, aniso.rotmat=NA,
           aniso.angles=NA, flag.check.sill=FALSE)

Arguments

vartype

Type of the Basic Structure to be created. By default, this type corresponds to the Nugget Effect. The list of types available is:

  • Nugget effect

  • Exponential

  • Spherical

  • Cubic

  • Gaussian

  • Cardinal Sine

  • J-Bessel

  • K-Bessel

  • Gamma

  • Cauchy

  • Stable

  • Linear

  • Cosinus (1D)

  • Triangle (1D)

  • Cosine Exponential

  • Exp2dfact

  • Expfact

  • Reg1D (1D)

  • Pentamodel

  • Power

  • Order-1 GC

Some are restricted to 1-D case. The two last types should be discarded when building a covariance Model

ndim

Space dimension.

nvar

Number of variables.

range

This parameter provides the range that should be attributed to the new Basic Structure (if required). Its contents depends on the isotropy / anisotropy characteristics defined in argument 'flag.aniso'.

  • In the isotropic case, this argument corresponds to the isotropic range of the basic structure.

  • In the anisotropic case, this argument must be a vector (dimension equal to 'ndim'): the i-th element is the extension of the anisotropic ellipsoid in the i-th direction (after rotation). Otherwise an error is issued.

scale

Theoretical range(s) of the new basic structure. If not defined (but needed for this type of basic structures), range can be defined instead.

Depending on the isotropic / anisotropic case, the argument 'scale' must be either a single value or a vector (see 'range' for details).

sill

This parameter provides the matrix of sills that should be attributed to the new Basic Structure (if required). The dimension of the matrix is 'nvar'.

param

This parameter provides the third parameter that should be attributed to the new Basic Structure (if required).

order

Parameter giving the maximum order of the IRF which corresponds to the model defined. If mode=-1, only covariances can be entered; if mode=0, all covariances and variograms are allowed.

flag.aniso

This parameter is used to specify that the new Basic Structure if anisotropic.

aniso.rotmat

This parameter provides the rotation matrix which defines the main axes of the rotation ellipsoid. It may not be used together with argument 'aniso.angles'.

aniso.angles

A vector giving the set of rotations. See util.ang2mat for details. It may not be used together with argument 'aniso.rotmat'.

flag.check.sill

When TRUE, the matrix of sill for each covariance component is checked against definite positiveness. This constraint is too severe for the whole model to be valid: therefore it is only optional.

Value

The new melem-class Basic Structure.


[Package RGeostats version 14.0.10 Index]