dbmss (version 2.0.3)

KdEnvelope: Estimation of the confidence envelope of the Kd function under its null hypothesis

Description

Simulates point patterns according to the null hypothesis and returns the envelope of Kd according to the confidence level.

Usage

KdEnvelope(X, r, NumberOfSimulations = 100, Alpha = 0.05, ReferenceType, 
           NeighborType = ReferenceType, Weighted = FALSE, Original = TRUE, 
           SimulationType = "RandomLocation", Global = FALSE)

Arguments

X

A point pattern (wmppp.object).

r

A vector of distances. No default value is available.

NumberOfSimulations

The number of simulations to run.

Alpha

The risk level.

ReferenceType

One of the point types. By default, the same as reference type.

NeighborType

One of the point types.

Weighted

Logical; if TRUE, estimates the Kemp function.

Original

Logical; if TRUE (by default), the original bandwidth selection by Duranton and Overman (2005) following Silverman (2006: eq 3.31) is used. If FALSE, it is calculated following Sheather and Jones (1991), i.e. the state o

SimulationType

A string describing the null hypothesis to simulate. The null hypothesis may be "RandomLocation": points are redistributed on the actual locations; "RandomLabeling": randomizes point types, keeping locations and weights unchang

Global

Logical; if TRUE, a global envelope sensu Duranton and Overman (2005) is calculated.

Value

• An envelope object (envelope). There are methods for print and plot for this class. The fv contains the observed value of the function, its average simulated value and the confidence enveloppe.

Details

This envelope is local by default, that is to say it is computed separately at each distance. See Loosmore and Ford (2006) for a discussion. The global enveloppe is calculated by iteration: the simulations reaching one of the upper or lower values at any distance are eliminated at each step. The process is repeated until Alpha / Number of simulations simulations are dropped. The remaining upper and lower bounds at all distances constitute the global envelope. Interpolation is used if the exact ratio cannot be reached.

References

Duranton, G. and Overman, H. G. (2005). Testing for Localisation Using Micro-Geographic Data. Review of Economic Studies 72(4): 1077-1106. Kenkel, N. C. (1988). Pattern of Self-Thinning in Jack Pine: Testing the Random Mortality Hypothesis. Ecology 69(4): 1017-1024. Loosmore, N. B. and Ford, E. D. (2006). Statistical inference using the G or K point pattern spatial statistics. Ecology 87(8): 1925-1931. Marcon, E. and F. Puech (2012). A typology of distance-based measures of spatial concentration. HAL SHS. 00679993. Silverman, B. W. (1986). Density estimation for statistics and data analysis. Chapman and Hall, London.

See Also

Kdhat

Examples

data(paracou16)
plot(paracou16[paracou16$marks$PointType=="Q. Rosea"])

# Calculate confidence envelope
r <- 0:100
plot(KdEnvelope(paracou16, r, ReferenceType="Q. Rosea", Global=TRUE))