# NOT RUN {
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_prp <- guerry_df['Crm_prp'][,1]
lit <- guerry_df['Litercy'][,1]
quantiles <- list(list(4,1,crm_prp), list(4,1, lit))
lisa <- local_multiquantilelisa(queen_w, quantiles)
clsts <- lisa_clusters(lisa)
clsts
# }
Run the code above in your browser using DataLab