simbipgs {RGeostats}R Documentation

Bi-Plurigaussian Conditional Simulations

Description

Bi-Plurigaussian conditional simulations

Usage

simbipgs(dbin = NA, dbout, dbprop=NA,
rule1 = rule.input(), rule2 = rule.input(),
model11 = NA, model12 = NA, model21=NA, model22=NA,
neigh=neigh.input(), props = NA, domain=0,
flag.gaus=FALSE, flag.prop=FALSE,
flag.check=FALSE, flag.show=FALSE, seed = 232131, 
nbsimu = 1, nbtuba = 100, nboot = 10, niter = 100, percent=0.1, 
toleps=1, radix = "SimBiPGS", modify.target = db.locmod())

Arguments

dbin

The db-class structure containing the data file. If not provided, non-conditional simulations are performed

dbout

The db-class structure which will contain the simulation outputs

dbprop

The db-class where the proportion variables must be located in the non-stationary case. See details in rule.check.

rule1

The rule-class structure which defines the Lithotype rule for the first PGS

rule2

The rule-class structure which defines the Lithotype rule for the second PGS

model11

The model-class structure containing the Model information for the first underlying gaussian random function of the first PGS. This model must be provided if the first gaussian random function is involved in the first Lithotype rule.

model12

The model-class structure containing the Model information for the first underlying gaussian random function of the first PGS. This model must be provided if the second gaussian random function is involved in the first Lithotype rule

model21

The model-class structure containing the Model information for the first underlying gaussian random function of the second PGS. This model must be provided if the first gaussian random function is involved in the second Lithotype rule.

model22

The model-class structure containing the Model information for the first underlying gaussian random function of the second PGS. This model must be provided if the second gaussian random function is involved in the second Lithotype rule

neigh

The neigh-class structure containing the Neighborhood information. Only the Unique or Bench neighborhoods are authorized.

props

Array giving the (constant) proportions of the different facies involved in both Lithotype rules. The number of elements of this array is the product of the number of facies involved in both Lithotype rules. See details in rule.check.

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.

flag.gaus

When FALSE, the simulation outcome is coded into facies. When TRUE, the (conditional) simulation is kept in the gaussian scale

flag.prop

When TRUE, the procedure returns the proportion of facies, rather than the simulations

flag.check

When TRUE (and if dbin is defined), the facies of conditioning data is compared to the facies of the closest grid node.

flag.show

When TRUE (and if dbin is defined), the grid node which coincides with a data is represented with the data facies (with a negative sign).

seed

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

nbsimu

Number of simulations

nbtuba

Number of turning bands

nboot

Number of iterations performed for bootstrap

niter

Maximum number of iterations for calculating the conditional expectation

percent

Amount of nugget effect added to the model if this initial model only contains Gaussian structures. If set to zero, the initial model is kept unchanged. This amount is defined as a percentage of the global variance.

toleps

Relative difference between the previous and the new value at a sample location below which the sample is considered as immobile

radix

Radix of the name given to the simulation outcomes stored in the output 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 data Db where the simulations outcomes have been added.


[Package RGeostats version 11.1.1 Index]