mlayer.krig {RGeostats}R Documentation

Multi-Layers estimation

Description

Multi-Layers estimation

Usage

mlayer.krig(dbin, dbout=NA, model = model.input(), neigh = neigh.input(), 
  	    flag.z = TRUE, flag.vel = FALSE, flag.cumul=FALSE,flag.ext = FALSE, 
	    flag.std=FALSE, flag.bayes=FALSE, irf.rank = 0, match.time = FALSE, 
            prior.mean=NA, prior.vars=NA,
	    colrefd = NA, colreft = NA, colrefb = NA, verbose=FALSE,
	    radix = "MLayers", modify.target = db.locmod())

Arguments

dbin

The db-class structure containing the 2-D data file. The file must be 2-D and contain one Z-variable corresponding to the sample depth and a locator "Layer".

dbout

The db-class structure containing the 2-D ouput file where the additional variables must be defined and where the results will be stored.

model

The model-class structure containing the Model information. The number of variables defines the number of layers

neigh

The neigh-class structure containing the Neighborhood information. The neighborhood is not used yet.

flag.z

When TRUE, the results are converted back into depth surfaces. When FALSE, the results are provided either in thickness (if flag.vel=FALSE) or in interval velocity (if flag.vel=TRUE)

flag.vel

Should be TRUE if the estimation must be performed using velocities (the Model should correspond to interval velocities) or FALSE when working with depth (the Model should correspond to interval thickness). If flag.vel is TRUE, the layer time maps should be defined in the grid output file (locator "time" or "f" if match_time is TRUE).

flag.cumul

Should be TRUE if the calculations must be performed in Depth rather than in Thickness. This is essential when considering the standard deviation of estimation errors.

flag.ext

Should be TRUE if External Drift should be used. Then the grid output file should contain the corresponding maps (locator "f").

flag.std

When TRUE, the standard deviation of the estimation error must be calculated.

flag.bayes

When TRUE, the drift coefficients are tuned using Bayesian hypothesis. The user must then define a vector of prior means and a vertor of prior variances for each drift coefficient

irf.rank

Rank of the Intrisic Random Function. Should be -1 (for strict stationarity), 0 (for intrinsic) or 1 (for a first order drift).

match.time

When TRUE, the external drift maps and the time maps coincide in the output grid file. They are both defined using the locator "f".

prior.mean

Vector giving the prior guess for the mean of the drift coefficients. The dimension of this vector must be equal to the number of drift functions multiplied by the number of layers.

prior.vars

Vector giving the prior guess for the variance of the drift coefficients (the cross-covariance of the priors is assumed to be zero). The dimension of this vector must be equal to the number of drift functions multiplied by the number of layers.

colrefd

Rank of the attribute containing the Reference Depth map. If 0, the reference depth is set to 0

colreft

Rank of the attribute containing the Reference Time map. If 0, the reference time is set to 0.

colrefb

Rank of the attribute containing the Bottom Depth map (if defined).

verbose

Verbose option which recapitulates the main options of the multi-layer estimation.

radix

Radix of the name given to the variables storing the results in the target Db.

modify.target

Decides whether or not the newly created variables will have their locator defined or not. For more information, see db.locmod.

Value

The output Db grid file where the estimations have been added:


[Package RGeostats version 11.1.2 Index]