global {RGeostats}R Documentation

Geostatistical estimation variance of the arithmetic/kriging mean

Description

Geostatistical estimation variance of the arithmetic/kriging mean

Usage

global(dbin, model = model.input(), uc=c("1"), polygon = NA, mean = NA,
dbout=NA, calcul = "krige", ivar = 1, verbose = 1, ndisc=50, seed = 321441,
flag.regular=TRUE, flag.polin=TRUE, flag.wgt=FALSE, combine="set")

Arguments

dbin

The db-class structure containing the data file. Note that this file MUST be organized as a regular grid when using the transitive method.

model

The model-class structure containing the Model information

uc

The drift description. Use command uc.names for details.

polygon

The polygon-class structure which contains the polygon within which the mean value must be estimated.

mean

Array containing the mean of each variable, used in the case of Simple (Co-)Kriging

dbout

The db-class structure containing the discretization Grid File. If not provided, a grid is defaulted (the number of cells is given by the parameters 'ndisc').

calcul

Option for the calculation of the global estimation. It can be:

  • 'arith': for an estimation using the arithmetic mean

  • 'krige': for an estimation using kriging

  • 'trans': for an estimation using transitive method

  • 'mean' : for an estimation of the theoretical mean

ivar

Rank of the target variable (starting from 1). Note that the transitive global estimation is not available in the multivariate case.

verbose

1 for a verbose output; 0 otherwise

ndisc

Number of discretization steps used for:

  • the calculation of the mean covariogram over a cell (stratified case)

  • the calculation of the integral range (regular case)

seed

Seed used for the generation of random numbers. When 0, the seed is not initialized.

flag.regular

When TRUE, use the regular transitive method; otherwise use the stratified transitive method.

flag.polin

If TRUE, only the input data which belong to the Polygon are used for the estimation. Otherwise, all the data are used (even if they do not lie within the polygon). See also the comments on the argument 'combine'.

flag.wgt

When TRUE, a vector which contains the kriging weights is provided as a returned argument. This option is only available when calcul="krige".

combine

When a selection of active data is performed within a Polygon ('flag.polin' set to TRUE) and if a selection already exists, the 'combine' keyword indicates how the polygon selection must be combined with the previous one. For more information on the available keywords, check db.selcombine.

Value

A list containing the following information:

* zest : Global estimation

* sse : Standard deviation on the Global estimated Mean

* cv : CV

* wgt : a vector containing the kriging weights (of the whole data set, regardless of any selection) when calcul="krige", NA otherwise.


[Package RGeostats version 14.0.10 Index]