Estimate the intensity threshold for the reconstruction of spatial structures

Usage

findIntensityThreshold(ppp, mark_select = NULL, bndw = NULL, dim, steps = 250)

Arguments

ppp: point pattern object of class ppp
mark_select: character; name of mark that is to be selected for the reconstruction
bndw: numeric; bandwith of the sigma parameter in the density estimation, if no value is given the bandwith is estimated using cross validation with the bw.diggle function.
dim: numeric; x dimension of the final reconstruction.
steps: numeric; value used to filter the density estimates, where only densities greater than the maximum value divided by threshold are considered. Default is 250.

Value

numeric; estimated intensity threshold

Examples

spe <- imcdatasets::Damond_2019_Pancreas("spe", full_dataset = FALSE)
#> see ?imcdatasets and browseVignettes('imcdatasets') for documentation
#> loading from cache
ppp <- SPE2ppp(spe, marks = "cell_category", image_col = "image_name", image_id = "E04")
findIntensityThreshold(ppp, mark_select = "islet", dim = 250)
#> [1] 0.003948522