sim: Random Field Simulation

Description

The function simulates a random field. The simulation methods available are based on Indicator Kriging techniques (IK and CK), Fixed and Random Path (PATH) and Multinomial Categorical Simulation (MCS).

Usage

sim(x, data, coords, grid, method = "ik", ...)

Arguments

an object of the class multi_tpfit, typically with the output of the function multi_tpfit.

data

a categorical data vector of length $n$.

coords

an $n x d$ matrix where each row denotes the $d$-D coordinates of data locations.

grid

an $m x d$ matrix where each row denotes the $d$-D coordinates in the simulation grid.

method

a character object specifying the method to simulate the random field. Possible choises are "ik" (by default) for the indicator Kriging, "ck" for the indicator coKriging, "path" for the fixed and random path and "mcs" for the multinomial categorical simulation method.

...

other arguments to pass to the functions sim_ik, sim_ck, sim_path or sim_mcs.

Value

Details

The methods implemented compute the approximation of posterior probabilities $$\Pr\left(Z(\mathbf{s}_0) = z_k \left\vert \bigcap_{i = 1}^n Z(\mathbf{s}_i) = z(\mathbf{s}_i)\right.\right).$$ $where i = 1, ..., n.$

Once the probabilities are calculated for all the points in the simulation grid, the predictions (based on most probable category) and simulations are returned.

References

Allard, D., D'Or, D., Froidevaux, R. (2011) An efficient maximum entropy approach for categorical variable prediction. European Journal of Soil Science, 62(3), 381-393.

Carle, S. F., Fogg, G. E. (1996) Transition Probability-Based Indicator Geostatistics. Mathematical Geosciences, 28(4), 453-476.

Carle, S. F. (1999) T-PROGS: Transition Probability Geostatistical Software. University of California, Davis.

Li, W. (2007) A Fixed-Path Markov Chain Algorithm for Conditional Simulation of Discrete Spatial Variables. Mathematical Geology, 39(2), 159-176.

Li, W. (2007) Markov Chain Random Fields for Estimation of Categorical Variables. Mathematical Geology, 39(June), 321-335.

Pickard, D. K. (1980) Unilateral Markov Fields. Advances in Applied Probability, 12(3), 655-671.

Sartore, L. (2010) Geostatistical models for 3-D data. M.Phil. thesis, Ca' Foscari University of Venice.

Weise, T. (2009) Global Optimization Algorithms - Theory and Application. http://www.it-weise.de/.