Computes distance-dependent estimates of Rao's quadratic entropy from a multivariate spatial point pattern in a simple (rectangular or circular) or complex sampling window. Computes optionally local confidence limits of the functions under the null hypothesis of either a random labelling or a species equivalence (see Details).
krfun(p, upto, by, nsim=0, dis = NULL, H0 = c("rl", "se"), alpha = 0.01)
A list of class "fads"
with essentially the following components:
a vector of regularly spaced out distances (seq(by,upto,by)
).
a data frame containing values of the function \(gr(r)\).
a data frame containing values of the function \(Kr(r)\).
Each component except r
is a data frame with the following variables:
a vector of estimated values for the observed point pattern.
a vector of theoretical values expected under the selected null hypothesis.
(optional) if nsim>0
a vector of the upper local confidence limits of a random distribution of the selected null hypothesis at a significant level \(\alpha\).
(optional) if nsim>0
a vector of the lower local confidence limits of a random distribution of the selected null hypothesis at a significant level \(\alpha\).
(optional) if nsim>0
a vector of local p-values of departure from the selected null hypothesis.
a "spp"
object defining a spatial point pattern in a given sampling window (see spp
).
maximum radius of the sample circles (see Details).
interval length between successive sample circles radii (see Details).
number of Monte Carlo simulations to estimate local confidence limits of the null hypothesis of a random allocation of species labels (see Details).
By default nsim = 0
, so that no confidence limits are computed.
(optional) a "dist"
object defining Euclidean distances between species. By default \(dis = NULL\) so that species are considered equidistant.
one of c("rl","se") to select either the null hypothesis of random labelling (H0 = "rl") or species equivalence (H0 = "se") (see Details). By default, the null hypothesis is random labelling.
if nsim>0
, significant level of the confidence limits. By default \(\alpha = 0.01\).
Function krfun
computes distance-dependent functions of Rao (1982) quadratic entropy (see divc
in package ade4
).
For a multivariate point pattern consisting of \(S\) species with intensity \(\lambda\)p, such functions can be estimated from the bivariate \(Kpq\)-functions between each pair of different species \(p\) and \(q\).
Function krfun
is thus a simple wrapper function of k12fun
and kfun
, standardized by Rao diversity coefficient (Pelissier & Goreaud 2014):
\(Kr(r) = sum(\lambda p * \lambda q * Kpq(r)*dpq) / (\lambda * \lambda * K(r) * HD)\).
\(gr(r) = sum(\lambda p * \lambda q * gpq(r)*dpq) / (\lambda * \lambda * g(r) * HD)\).
where \(dpq\) is the distance between species \(p\) and \(q\) defined by matrix dis
, typically a taxonomic, phylogenetic or functional distance, and \(HD=sum(Np*Nq*dpq/(N(N - 1)))\) is the unbiased version of Rao diversity coefficient (see Shimatani 2001). When dis = NULL
, species are considered each other equidistant and krfun
returns the same results than ksfun
.
The program introduces an edge effect correction term according to the method proposed by Ripley (1977)
and extended to circular and complex sampling windows by Goreaud & Pelissier (1999).
Theoretical values under the null hypothesis of either random labelling or species equivalence as well as local Monte Carlo confidence limits and p-values of departure from the null hypothesis (Besag & Diggle 1977) are estimated at each distance \(r\).
The random labelling hypothesis (H0 = "rl") is tested by reallocating species labels at random among all points of the pattern, keeping the point locations unchanged, so that expectations of \(gr(r)\) and \(Kr(r)\) are 1 for all \(r\).
The species equivalence hypothesis (H0 = "se") is tested by randomizing the between-species distances, keeping the point locations and distribution of species labels unchanged. The theoretical expectations of \(gr(r)\) and \(Kr(r)\) are thus \(gs(r)\) and \(Ks(r)\), respectively (see ksfun
).
Rao, C.R. 1982. Diversity and dissimilarity coefficient: a unified approach. Theoretical Population Biology, 21:24-43.
Shimatani, K. 2001. On the measurement of species diversity incorporating species differences. Oikos, 93, 135-147.
Goreaud F. & Pelissier R. 1999. On explicit formulas of edge effect correction for Ripley's K-function. Journal of Vegetation Science, 10:433-438.
Ripley B.D. 1977. Modelling spatial patterns. Journal of the Royal Statistical Society B, 39:172-192.
Pelissier, R. & Goreaud, F. 2014. ads package for R: A fast unbiased implementation of the k-function family for studying spatial point patterns in irregular-shaped sampling windows. Journal of Statistical Software, in press.
plot.fads
,
spp
,
ksfun
,
kdfun
,
divc
.