potential.simulate {RGeostats}R Documentation

Perform simulations using the Potential Model

Description

Perform simulations using the Potential Model

Usage

potential.simulate(db, dbgrd, dbout, dbtgt = NA, model = model.input(),
                   neigh = neigh.input(), uc = c("1"),
                   nugget.grd = 0, nugget.tgt = 0, dist.tempere = NA,
                   flag.trans = TRUE, verbose = FALSE,
                   seed = 1332, nbsimu = 1, nbtuba = 100,
                   radix = "Potential", modify.target = db.locmod())

Arguments

db

The db-class which contains the iso-potential information. This file must be defined in 1, 2 or 3-D. It must contain a 'layer' attribute.

dbgrd

The db-class which contains the gradient information. This file must contain as many 'gradient' variables as the space dimension. This file is compulsory.

dbout

The db-class which will contain the simulated outputs. This object must have the same space dimension as the other Db files: then the potential will be simulated for each cell. But the space dimension of this output file can be one degree less: in that case, the program will find the coordinate of the different iso-potential values per cell. This latter case will be refered to as Interpolation mode.

dbtgt

The db-class which contains the tangent information. This file must contain as many 'tangent' variables as the space dimension. This file is optional.

model

The model-class which contains the structure describing the spatial characteristics of the Potential variable. The model is composed of one or several basic structures, which belong to a limited set of authorized structures:

  • Gaussian

  • Cubic

  • Spline (2nd-order)

An additional nugget effect is also authorized.

neigh

The neigh-class which contains the neighborhood description. Currently only the Unique Neighborhood is coded.

uc

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

nugget.grd

The value of the nugget effect on the Gradient component

nugget.tgt

The value of the nugget effect on the Tangent component

dist.tempere

Distance above which the simulations are tempered close enough to the kriged value.

flag.trans

When TRUE, the estimation is returned as the layer rank; otherwise the estimation is provided as the potential value.

verbose

Verbose option

seed

Seed for the random number generation

nbsimu

Number of simulations

nbtuba

Number of Turning Bands

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.

Details

When performing the simulations on 'dbout' target points using the Potential Method, the Potential values at Intercept (recentered to the value of this Potential at the first data point of the first Intercept) are dumped out. They can be accessed using the keypair mechanism as follows:

get.keypair("Potential.Simulation.xxx")

where 'xxx' stands for the rank of the simulation (starting from 1).

Value

The target Db where the simulated variables have been added.


[Package RGeostats version 14.0.10 Index]