spatstat (version 1.59-0)

linearKdot: Multitype K Function (Dot-type) for Linear Point Pattern

Description

For a multitype point pattern on a linear network, estimate the multitype $K$ function which counts the expected number of points (of any type) within a given distance of a point of type $i$.

Usage

linearKdot(X, i, r=NULL, …, correction="Ang")

Arguments

X

The observed point pattern, from which an estimate of the dot type $K$ function $K_{i\bullet }(r)$ will be computed. An object of class "lpp" which must be a multitype point pattern (a marked point pattern whose marks are a factor).

i

Number or character string identifying the type (mark value) of the points in X from which distances are measured. Defaults to the first level of marks(X).

r

numeric vector. The values of the argument $r$ at which the $K$-function $K_{i\bullet }(r)$ should be evaluated. There is a sensible default. First-time users are strongly advised not to specify this argument. See below for important conditions on $r$.

correction

Geometry correction. Either "none" or "Ang". See Details.

…

Ignored.

Value

An object of class "fv" (see fv.object).

Warnings

The argument i is interpreted as a level of the factor marks(X). Beware of the usual trap with factors: numerical values are not interpreted in the same way as character values.

Details

This is a counterpart of the function Kdot for a point pattern on a linear network (object of class "lpp").

The argument i will be interpreted as levels of the factor marks(X). If i is missing, it defaults to the first level of the marks factor.

The argument r is the vector of values for the distance $r$ at which $K_{i\bullet }(r)$ should be evaluated. The values of $r$ must be increasing nonnegative numbers and the maximum $r$ value must not exceed the radius of the largest disc contained in the window.

References

Baddeley, A, Jammalamadaka, A. and Nair, G. (to appear) Multitype point process analysis of spines on the dendrite network of a neuron. Applied Statistics (Journal of the Royal Statistical Society, Series C), In press.

See Also

Kdot, linearKcross, linearK.

Examples

# NOT RUN {
   data(chicago)
   K <- linearKdot(chicago, "assault")
# }