getlen: Estimation of Stratum Lengths for Embedded Markov Chain

In spMC: Continuous-Lag Spatial Markov Chains

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)

Arguments

data	a categorical data vector of length n.
coords	an n \times d matrix where each row denotes the d-D coordinates of data locations.
loc.id	a vector of n values which indicates the directional line of each location. It is usually the output of the function which_lines.
direction	a d-D numerical vector (or versor) which represents the chosen direction.
zero.allowed	a logical value which allows to return zero stratum lengths. It is FALSE by default.

Details

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

Value

A list containing the following components:

length	a numerical vector with the stratum lengths along the given direction.
categories	a vector with the stratum categories.
maxcens	a vector with the maxima estimated censored lengths for each stratum.
directions	a d-D numerical vector which represents the chosen direction.
zeros	a logical values which denotes the possible presence of zero lengths.

Author(s)

Luca Sartore [e-mail]

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)

spMC documentation built on May 3, 2023, 9:13 a.m.