Simple computation of spatial statistic functions of distance to characterize the spatial structures of mapped objects, following Marcon, Traissac, Puech, and Lang (2015)
dbmss is an R package for simple computation of spatial statistic functions of distance to characterize the spatial structures of mapped objects, including classical ones (Ripley's K and others) and more recent ones used by spatial economists (Duranton and Overman's Kd, Marcon and Puech's M). Relies on spatstat for some core calculation.
A quick introduction is in vignette("docs", "dbmss")
.
A full documentation is available in section section of the same vignette (Reference Guide). It is a continuous update of the paper published in the Journal of Statistical Software (Marcon et al., 2015).
Marcon, E., Traissac, S., Puech, F. and Lang, G. (2015). Tools to Characterize Point Patterns: dbmss for R. Journal of Statistical Software. 67(3): 1-15.
init.c
no longer necessary.wmppp()
accepts characters in PointType
. A PointName
column is searched to set the row names of the marks.wmppp()
failed if some point laid outside the window.Mhat
and mhat
with argument Individual
.wmppp()
. They are preserved in the fv
object returned by individual Mhat
and mhat
to identify points more easily.[.wmppp
function added to spatstat genericsexportPattern("^[[:alpha:]]+")
.adaptIntegrate()
.Mhat
, mhat
and Kdhat
and their envelopes.Kd
now includes reflection: the estimation of Kd
is more accurate than before, but values may vary from previous versions close to the minimum distance.Kdhat
(bug introduced in v.2.2-4). Fixed.Mhat
and mhat
C++ code is now parallelized thanks to RcppParallel.Kdhat
(around 5% faster).r
where ignored when using the approximate estimation of Kd
, so Kd
was overestimated when r
was smaller than usual (much less than the default values). Corrected.Kdhat
. Fixed.m
function added.Controls=NULL
for Dhat
: controls are by default all points except for cases.Kd
and m
are both computed by default up to one third of the diameter of the window, other options are added.Original
and Adjust
arguments allow to change the default bandwith used by Kd
and m
.FillEnveloppe
corrected: FillEnvelope
is the new name.Khat
and Kenvelope
is the same as reference type instead of "".Kd
estimation now accepts "" as reference type to use all points.Kd
on big data sets with little RAM.FillEnvelope
returned wrong quantiles for local confidence intervals. Fixed.NA
values are found in the simulations.Kd
is computed by default up to the median (instead of mean) distance between points.Kdhat
with Reference and Neighbor points of different types was calculated with only half of point pairs and sometimes crashed R. Fixed.Kdhat
and M
are computed through C++ loops instead of spatstat pairdist()
. Much faster, and uses far less memory.Kdhat
and M
propose a default value for argument r
.Kdhat
probability density estimation is cut at the lowest distance between points instead of 0.Kdhat
with Reference and Neighbor points of the same type and Weighted=TRUE
returned an error. Fixed.X.r
became Xhat.wmppp
.fv
, envelopes of class envelope
.Kd
is that of Duranton and Overman (2005). An optional parameter is added to used improved bandwith selection.Kd
was run twice. Fixed.Kinhom
function added.