local_quantilelisa: Quantile LISA statistics
In lixun910/rgeoda: R Library for Spatial Data Analysis

Description Usage Arguments Value Examples

The function to apply quantile LISA statistics

local_quantilelisa(
  w,
  k,
  q,
  data,
  permutations = 999,
  significance_cutoff = 0.05,
  cpu_threads = 6,
  seed = 123456789
)

`w`	An instance of Weight object
`k`	A value indicates the number of quantiles. Value range e.g. [1, 10]
`q`	A value indicates which quantile or interval used in local join count statistics. Value stars from 1.
`data`	A numeric array of selected variable
`permutations`	The number of permutations for the LISA computation
`significance_cutoff`	A cutoff value for significance p-values to filter not-significant clusters
`cpu_threads`	The number of cpu threads used for parallel LISA computation
`seed`	The seed for random number generator

lisa_obj An instance of LISA (LocalSpatialAutocorrelation) object

guerry_path <- system.file("extdata", "Guerry.shp", package = "rgeoda")
guerry <- geoda_open(guerry_path)
queen_w <- queen_weights(guerry)
guerry_df <- as.data.frame(guerry) # use as data.frame
crm_prs <- guerry_df['Crm_prs'][,1]
lisa <- local_quantilelisa(queen_w, k=4, q=1, crm_prs)
clsts <- lisa_clusters(lisa)
clsts