Simulate: Geostatistical simulation

In italo-goncalves/geomod3D: Implicit 3D Geological Modeling and Geostatistics

Description

Generates a number of realizations of the modeled random field.

Usage

Simulate(object, ...)

## S4 method for signature 'GP'
Simulate(object, target, to = "value", Nsim,
  discount.noise = F)

## S4 method for signature 'SPGP'
Simulate(object, target, Nsim, to = "value",
  discount.noise = F, smooth = T, verbose = T, randnum = NULL)

Arguments

object	A GP object.
target	The spatial3DDataFrame object to receive the prediction.
to	The name of the column in which to write the prediction. Will be overwritten if it exists or created if not.
Nsim	The desired number of simulations.
discount.noise	Whether to sample from the GP latent function instead of the noisy values.
smooth	Whether to correct the output from the sparse GP.
verbose	Display status while running?

Details

Standard GP objects use the Cholesky method, which limits the number of data to be used and the number of test points in which to simulate. SPGP objects use a sequential method based on rank-one updates to the matrices involved, which only limits the number of pseudo-inputs that can be used.

Due to the nature of the sparse approximation of the covariance matrices, the output may appear noisy even if discount.noise = T. Set smooth = T to apply a correction to the generated samples.

Value

A 3D spatial object of the same class as target containing the simulations.

italo-goncalves/geomod3D documentation built on May 24, 2019, 2:49 p.m.