spdep
spdep copied to clipboard
Example in `n.comp.nb()` not reproducible
I am trying to find out the best way to identify the disconnected subgraphs from a nb
object. In trying to explore how to convert an nb
object to an igraph
object I ran into the following issue.
library(sf)
#> Warning: package 'sf' was built under R version 4.3.1
#> Linking to GEOS 3.11.0, GDAL 3.5.3, PROJ 9.1.0; sf_use_s2() is TRUE
library(spdep)
#> Loading required package: spData
#> To access larger datasets in this package, install the spDataLarge
#> package with: `install.packages('spDataLarge',
#> repos='https://nowosad.github.io/drat/', type='source')`
columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
coords <- st_coordinates(st_centroid(st_geometry(columbus)))
col2 <- droplinks(col.gal.nb, 21)
res <- n.comp.nb(col2)
B <- as(nb2listw(col2, style="B", zero.policy=TRUE), "CsparseMatrix")
#> Error in as(nb2listw(col2, style = "B", zero.policy = TRUE), "CsparseMatrix"): no method or default for coercing "listw" to "CsparseMatrix"
Created on 2024-03-31 with reprex v2.0.2
Coercion method is in spatialreg
.