simsub {RGeostats}R Documentation

Substitution Simulation

Description

Simulation according to the Substitution Model

Usage

simsub(dbgrid, nfacies = 2, intensity = 0.1, nstates = NA, flag.direct = TRUE, 
       flag.coding = TRUE, flag.orient = FALSE, factor = 0, vector = NA, 
       trans = NA, seed = 232131, colfac=NA, colax=NA, colay=NA, colaz=NA,
       verbose = FALSE, domain = 0, radix = "Substitution", 
       modify.target = db.locmod())

Arguments

dbgrid

The db-class structure containing the target file. It must be organized as a grid.

nfacies

Number of facies.

intensity

Intensity of the Poisson point process used to implement the Poisson lines.

nstates

Number of states. If not specified, this number will be calculated from the number of Pisson planes generated.

flag.direct

When TRUE, the direction information simulation is performed internally. Otherwise, this outcome is read as the Z locator of the Db.

flag.coding

When TRUE, the coding process is performed internally. Otherwise, the result of the directing simulation is output.

flag.orient

When TRUE, the desorientation vector is taken into account

factor

This coefficient which must belong to [0,1] quantifies the desorientation strength: 0 corresponds to an isotropic outcome (no disorientation) whereas 1 gives a stratified output. If this factor should vary in space, use the colfac argument instead.

vector

Vector giving the orientation orthogonal to the layering of the disoriented outcome. Its dimension must match the space dimension. Otherwise the diagonal vector is defaulted. If this vector should vary in space, use the colax,colay,colaz arguments instead.

trans

The transition probability matrix. Its dimension must be nfacies * nfacies. Each one of its elements must lie bwteen 0 and 1. By row, the transition probabilities must add up to 1.

seed

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

colfac

Rank of the attribute containing a non-stationary desorientation factor. Note that this value must lie between 0 and 1. Set this argument to 0 (default value) if this non-stationary desorientation factor must not be used. Otherwise it will be considered only if the argument flag_orient is set to TRUE. The argument 'colfac' overwrites the argument 'factor' (if specified).

colax,colay,colaz

Rank of the attribute containing a non-stationary desorientation vector. Set these arguments to 0 (default value) if this non-stationary desorientation vector must not be used. Otherwise they will be considered only if the argument flag_orient is set to TRUE. The arguments 'colax,colay,colaz' overwrite the argument 'vector' (if defined).

verbose

When TRUE, the substitution simulation process issues control messages.

domain

When (strictly) positive, it corresponds to the reference domain value. If a Db attribute is defined as a Domain variable in a Db file, a sample will be selected only if the corresponding Domain variable is equal to the reference domain value. This Domain selection can serve as a secondary selection for the sample in the Input Db or as a mask in the Output Db.

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 target Db where the simulation outcome has been added.


[Package RGeostats version 11.0.6 Index]