spMC (version 0.2.0)

imgMultiTransiogram: Images with Pseudo-empirical Multidimensional Transiograms

Description

The function plots $2$-D sections of a pseudo-empirical multidimensional transiogram computed without any ellipsoidal interpolation.

Usage

imgMultiTransiogram(data, coords, mpoints, which.dire, max.dist,
                    tolerance = pi/8, main, mar, ask = TRUE, ...,
                    nlevels = 10, contour = FALSE, rotation = NULL,
                    mle = FALSE)

Value

An image is produced on the current graphics device. No values are returned.

Rdversion

1.1

Details

A multidimensional transiogram is a diagram which shows the transition probabilities for a single pair of categories. The probability is computed for any lag vector $h$ through $$\mbox{expm} (\Vert h \Vert R_h),$$ where entries of $R_h$ are not ellipsoidally interpolated, but they are estimated for the direction specified by the vector $h$.

In particular cases, some entries of the estimated matrix $R_h$ might be not finite, so that the exponential matrix is computable and the resulting transition probabilities are set to be NaN. If mle = TRUE, this problem may be partially solved.

The exponential matrix is evaluated by the scaling and squaring algorithm.

References

Carle, S. F., Fogg, G. E. (1997) Modelling Spatial Variability with One and Multidimensional Continuous-Lag Markov Chains. Mathematical Geology, 29(7), 891-918.

Higham, N. J. (2008) Functions of Matrices: Theory and Computation. Society for Industrial and Applied Mathematics.

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

See Also

image.multi.tpfit, image, plot.transiogram

Examples

data(ACM)

# Plot 2-D pseudo-empirical sections of
# a multidimensional transiogram
imgMultiTransiogram(ACM$MAT3, ACM[, 1:3], 10, 
    max.dist=c(20, 10, 5), which.dire=c(1, 3),
    col = rev(heat.colors(500)), breaks = 0:500 / 500,
    mar = .7, mle = TRUE)