Create a kernel weights by specifying a bandwidth and a kernel method
kernel_weights(
sf_obj,
bandwidth,
kernel_method,
use_kernel_diagonals = FALSE,
power = 1,
is_inverse = FALSE,
is_arc = FALSE,
is_mile = TRUE
)An sf (simple feature) object
A positive numeric value of bandwidth
a string value, which has to be one of 'triangular', 'uniform', 'epanechnikov', 'quartic', 'gaussian'
(optional) FALSE (default) or TRUE, apply kernel on the diagonal of weights matrix
(optional) The power (or exponent) of a number says how many times to use the number in a multiplication.
(optional) FALSE (default) or TRUE, apply inverse on distance value
(optional) FALSE (default) or TRUE, compute arc distance between two observations
(optional) TRUE (default) or FALSE, convert distance unit from mile to km.
An instance of Weight-class
library(sf)
guerry_path <- system.file("extdata", "Guerry.shp", package = "rgeoda")
guerry <- st_read(guerry_path)
#> Reading layer `Guerry' from data source
#> `/Users/runner/work/_temp/Library/rgeoda/extdata/Guerry.shp'
#> using driver `ESRI Shapefile'
#> Simple feature collection with 85 features and 29 fields
#> Geometry type: MULTIPOLYGON
#> Dimension: XY
#> Bounding box: xmin: 47680 ymin: 1703258 xmax: 1031401 ymax: 2677441
#> Projected CRS: NTF (Paris) / Lambert zone II
bandwidth <- min_distthreshold(guerry)
kernel_w <- kernel_weights(guerry, bandwidth, kernel_method = "uniform")
summary(kernel_w)
#> name value
#> 1 number of observations: 85
#> 2 is symmetric: FALSE
#> 3 sparsity: 0.0434602076124567
#> 4 # min neighbors: 1
#> 5 # max neighbors: 7
#> 6 # mean neighbors: 3.69411764705882
#> 7 # median neighbors: 4
#> 8 has isolates: FALSE