Description Usage Arguments Details Value References See Also
Monte Carlo test of spatial segregation in a multivariate point process by simulating data from random re-labelling of the categorical marks.
1 | mcseg.test(pts, marks, h, stpts = NULL, ntest = 100, proc = TRUE)
|
pts |
matrix containing the |
marks |
numeric/character vector of the marked type labels of the point pattern. |
h |
numeric vector of the bandwidths at which to calculate the cross-validated likelihood function. |
stpts |
matrix containing the |
ntest |
integer with default 100, number of simulations for the Monte Carlo test. |
proc |
logical with default |
The null hypothesis is that the estimated risk surface is spatially
constant, i.e., the type-specific probabilities are
p_k(x)=p_k, for all k, see phat
. Each Monte Carlo
simulation is done by relabeling the data categorical marks at random
whilst preserving the observed number of cases of each type.
The segregation test can also be done pointwise, usually at a fine grid of points, to mark the areas where the estimated type-specific probabilities are significantly greater or smaller than the spatial average.
A list with components
numeric, p-value of the Monte Carlo test.
matrix, p-values of the test at each point in
stpts
(if stpts
is not NULL
), with each column corresponds
to one type
copy of the arguments pts, marks, h, stpts, ntest, proc
.
Kelsall, J. E. and Diggle, P. J. (1998) Spatial variation in risk: a nonparametric binary regression approach, Applied Statistics, 47, 559–573.
Diggle, P. J. and Zheng, P. and Durr, P. A. (2005) Nonparametric estimation of spatial segregation in a multivariate point process: bovine tuberculosis in Cornwall, UK. J. R. Stat. Soc. C, 54, 3, 645–658.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.