Returns the Getis Ord G*i in standardized form.

gstar_i(shp, adj, wts, spatial_mat, epsg = 3857)

Arguments

shp

sf data frame. Optional if adj or spatial_mat provided.

adj

zero indexed adjacency list. Optional if shp or spatial_mat provided.

wts

Required. Numeric vector with weights to use for Moran's I.

spatial_mat

matrix of spatial weights. Optional if shp or adj provided.

epsg

numeric EPSG code to planarize to. Default is 3857.

Value

vector of G*i scores

Examples

library(dplyr)
data('checkerboard')
checkerboard <- checkerboard %>% mutate(m = as.numeric((id + i) %% 2 == 0))
gstar_i(shp = checkerboard, wts = checkerboard$m)
#> Warning: Planarizing skipped. `x` missing CRS.
#>  [1]  1.425573 -1.760216  1.760216 -1.760216  1.760216 -1.760216  1.760216
#>  [8] -1.425573 -1.760216  2.049390 -2.049390  2.049390 -2.049390  2.049390
#> [15] -2.049390  1.760216  1.760216 -2.049390  2.049390 -2.049390  2.049390
#> [22] -2.049390  2.049390 -1.760216 -1.760216  2.049390 -2.049390  2.049390
#> [29] -2.049390  2.049390 -2.049390  1.760216  1.760216 -2.049390  2.049390
#> [36] -2.049390  2.049390 -2.049390  2.049390 -1.760216 -1.760216  2.049390
#> [43] -2.049390  2.049390 -2.049390  2.049390 -2.049390  1.760216  1.760216
#> [50] -2.049390  2.049390 -2.049390  2.049390 -2.049390  2.049390 -1.760216
#> [57] -1.425573  1.760216 -1.760216  1.760216 -1.760216  1.760216 -1.760216
#> [64]  1.425573