Kmmhat: Estimation of the Kmm function
In dbmss: Distance-Based Measures of Spatial Structures

Kmmhat

R Documentation

Estimation of the Kmm function

Description

Estimates of the Kmm function

Usage

Kmmhat(X, r = NULL, ReferenceType = "", CheckArguments = TRUE)

Arguments

`X`	A weighted, marked, planar point pattern (`wmppp.object`).
`r`	A vector of distances. If `NULL`, a sensible default value is chosen (512 intervals, from 0 to half the diameter of the window) following spatstat.
`ReferenceType`	One of the point types. Others are ignored. Default is all point types.
`CheckArguments`	Logical; if `TRUE`, the function arguments are verified. Should be set to `FALSE` to save time in simulations for example, when the arguments have been checked elsewhere.

Details

The Kmm function is used to test the independence of marks.

Value

An object of class fv, see fv.object, which can be plotted directly using plot.fv.

Note

The function is computed using markcorrint in spatstat.

References

Penttinen, A., Stoyan, D. and Henttonen, H. M. (1992). Marked Point Processes in Forest Statistics. Forest Science 38(4): 806-824.

Penttinen, A. (2006). Statistics for Marked Point Patterns. in The Yearbook of the Finnish Statistical Society. The Finnish Statistical Society, Helsinki: 70-91.

Examples

data(paracou16)
# Keep only 50% of points to run this example
X <- as.wmppp(rthin(paracou16, 0.5))
autoplot(X, 
  labelSize = expression("Basal area (" ~cm^2~ ")"), 
  labelColor = "Species")

# Calculate Kmm
r <- seq(0, 30, 2)
(Paracou <- Kmmhat(X, r))

# Plot
autoplot(Paracou, ./(pi*r^2) ~ r)

dbmss documentation built on Sept. 11, 2024, 9:19 p.m.