Returns the Getis Ord G*i in standardized form.
gstar_i(shp, adj, wts, spatial_mat, epsg = 3857)sf data frame. Optional if adj or spatial_mat provided.
zero indexed adjacency list. Optional if shp or spatial_mat provided.
Required. Numeric vector with weights to use for Moran's I.
matrix of spatial weights. Optional if shp or adj provided.
numeric EPSG code to planarize to. Default is 3857.
vector of G*i scores
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