sim_path: Conditional Simulation Based on Path Algorithms

Description

The function simulates a random field through the Fixed Path algorithm or Random Path technique.

Usage

sim_path(x, data, coords, grid, radius, fixed = FALSE)

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.

radius

a numerical value that specifies a proper radius to search the nearest observed points within a $d$-D sphere.

fixed

a logical value; if TRUE, the fixed path algorithm is performed. The random path algorithm is performed by default.

Value

Details

These methods compute an 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.$ Path algorithms are based on Pickard random fields, so that the states of such chain at any unsampled location depends on the state of its nearest known neighbours in axial directions.

References

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.