Function to apply local Bivariate Join Count statistics
Usage
st_local_bijoincount(
sfj,
varcol,
wt = NULL,
permutations = 999,
permutation_method = "complete",
significance_cutoff = 0.05,
cpu_threads = 6,
seed = 123456789
)
Arguments
- sfj
An sf (simple feature) object.
- varcol
The variable selected to calculate spatial lag, which is a character.
- wt
(optional) The spatial weights object,which can use
st_weights()
to construct,default is constructed byst_weights(sfj,'contiguity')
.- permutations
(optional) The number of permutations for the LISA computation.
- permutation_method
(optional) The permutation method used for the LISA computation. Options are 'complete', 'lookup'. Default is 'complete'.
- significance_cutoff
(optional) A cutoff value for significance p-values to filter not-significant clusters.
- cpu_threads
(optional) The number of cpu threads used for parallel LISA computation.
- seed
(optional) The seed for random number generator.
Author
Wenbo Lv lyu.geosocial@gmail.com
Examples
library(magrittr)
guerry_path = system.file("extdata", "Guerry.shp", package = "rgeoda")
guerry = sf::read_sf(guerry_path) %>% dplyr::mutate(InvCrm = 1 - TopCrm)
st_local_bijoincount(guerry,c("TopCrm", "InvCrm"))
#> [1] Not significant Significant Significant Not significant
#> [5] Not significant Not significant Significant Not significant
#> [9] Not significant Not significant Not significant Not significant
#> [13] Not significant Not significant Not significant Not significant
#> [17] Not significant Not significant Significant Significant
#> [21] Not significant Not significant Not significant Not significant
#> [25] Not significant Not significant Not significant Not significant
#> [29] Not significant Not significant Not significant Not significant
#> [33] Not significant Not significant Not significant Not significant
#> [37] Significant Not significant Not significant Not significant
#> [41] Not significant Not significant Not significant Not significant
#> [45] Not significant Not significant Not significant Not significant
#> [49] Not significant Not significant Significant Not significant
#> [53] Significant Not significant Not significant Not significant
#> [57] Not significant Not significant Not significant Not significant
#> [61] Not significant Not significant Not significant Not significant
#> [65] Not significant Not significant Not significant Not significant
#> [69] Significant Not significant Not significant Not significant
#> [73] Not significant Not significant Not significant Significant
#> [77] Not significant Not significant Not significant Not significant
#> [81] Not significant Not significant Not significant Not significant
#> [85] Not significant
#> Levels: Not significant Significant Undefined Isolated