spMC (version 0.2.2)

getlen: Estimation of Stratum Lengths for Embedded Markov Chain

Description

The function estimates the stratum lengths for a $d$-D spatial embedded Markov chain for a specified direction $\phi$.

Usage

getlen(data, coords, loc.id, direction, zero.allowed = FALSE)

Value

A list containing the following components:lengtha numerical vector with the stratum lengths along the given direction.categoriesa vector with the stratum categories.maxcensa vector with the maxima estimated censored lengths for each stratum.directionsa $d$-D numerical vector which represents the chosen direction.zerosa logical values which denotes the possible presence of zero lengths.

Rdversion

1.1

Details

Stratum lengths are the lengths occupied by the same $k$-th category along lines in the direction $\phi$.

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.

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

See Also

mlen, which.lines

Examples

data(ACM)
direction <- c(0,0,1)

# Compute the appertaining directional line for each location
loc.id <- which.lines(ACM[, 1:3], direction)

# Estimate stratum lengths
gl <- getlen(ACM$MAT5, ACM[, 1:3], loc.id, direction)