r-spatial · rsbivand · Jun 10, 2024 · Jun 8, 2024 · Jun 10, 2024
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: spdep
 Version: 1.3-5
-Date: 2024-06-02
+Date: 2024-06-10
 Title: Spatial Dependence: Weighting Schemes, Statistics
 Encoding: UTF-8
 Authors@R: c(person("Roger", "Bivand", role = c("cre", "aut"), 
@@ -38,7 +38,7 @@ Authors@R: c(person("Roger", "Bivand", role = c("cre", "aut"),
 	person("René", "Westerholt", role="ctb"),
 	person("Levi", "Wolf", role = "ctb"),
 	person("Danlin", "Yu", role = "ctb"))
-Depends: R (>= 3.3.0), methods, spData (>= 0.2.6.0), sf 
+Depends: R (>= 3.3.0), methods, spData (>= 2.3.1), sf 
 Imports: stats, deldir, boot (>= 1.3-1), graphics, utils, grDevices, units, s2, e1071, sp (>= 1.0)
 Suggests: spatialreg (>= 1.2-1), Matrix, parallel, dbscan, RColorBrewer, lattice, xtable, foreign, igraph, RSpectra, knitr, classInt, tmap, spam, ggplot2, rmarkdown, tinytest, rgeoda
 URL: https://github.com/r-spatial/spdep/, https://r-spatial.github.io/spdep/

diff --git a/NEWS.md b/NEWS.md
@@ -1,4 +1,6 @@
-# Version 1.3-4 (development)
+# Version 1.3-5 (development)
+
+* #157 migrate ESRI Shapefile to GPKG files; convert bhicv.shp to GPKG 
 
 * #155 Throw error if `hotspot` despatched on object without a `"quadr"` attribute
 

diff --git a/inst/etc/shapes/bhicv.dbf b/inst/etc/shapes/bhicv.dbf
diff --git a/inst/etc/shapes/bhicv.gpkg b/inst/etc/shapes/bhicv.gpkg
diff --git a/inst/etc/shapes/bhicv.shp b/inst/etc/shapes/bhicv.shp
diff --git a/inst/etc/shapes/bhicv.shx b/inst/etc/shapes/bhicv.shx
diff --git a/inst/etc/shapes/columbus.dbf b/inst/etc/shapes/columbus.dbf
diff --git a/inst/etc/shapes/columbus.shp b/inst/etc/shapes/columbus.shp
diff --git a/inst/etc/shapes/columbus.shx b/inst/etc/shapes/columbus.shx
diff --git a/inst/etc/shapes/eire.dbf b/inst/etc/shapes/eire.dbf
diff --git a/inst/etc/shapes/eire.shp b/inst/etc/shapes/eire.shp
diff --git a/inst/etc/shapes/eire.shx b/inst/etc/shapes/eire.shx
diff --git a/man/EBImoran.mc.Rd b/man/EBImoran.mc.Rd
@@ -60,7 +60,7 @@ following components:
 \seealso{\code{\link{moran}}, \code{\link{moran.mc}}, \code{\link{EBest}}}
 
 \examples{
-nc.sids <- st_read(system.file("shapes/sids.shp", package="spData")[1], quiet=TRUE)
+nc.sids <- st_read(system.file("shapes/sids.gpkg", package="spData")[1], quiet=TRUE)
 rn <- as.character(nc.sids$FIPS)
 ncCC89_nb <- read.gal(system.file("weights/ncCC89.gal", package="spData")[1],
  region.id=rn)

diff --git a/man/EBest.Rd b/man/EBest.Rd
@@ -32,7 +32,7 @@ and a \code{parameters} attribute list with components:
 \seealso{\code{\link{EBlocal}}, \code{\link{probmap}}, \code{\link{EBImoran.mc}}}
 
 \examples{
-auckland <- st_read(system.file("shapes/auckland.shp", package="spData")[1], quiet=TRUE)
+auckland <- st_read(system.file("shapes/auckland.gpkg", package="spData")[1], quiet=TRUE)
 res <- EBest(auckland$M77_85, 9*auckland$Und5_81)
 attr(res, "parameters")
 auckland$estmm000 <- res$estmm*1000

diff --git a/man/EBlocal.Rd b/man/EBlocal.Rd
@@ -37,7 +37,7 @@ and a \code{parameters} attribute list with components (if both are zero, the es
 \seealso{\code{\link{EBest}}, \code{\link{probmap}}}
 
 \examples{
-auckland <- st_read(system.file("shapes/auckland.shp", package="spData")[1], quiet=TRUE)
+auckland <- st_read(system.file("shapes/auckland.gpkg", package="spData")[1], quiet=TRUE)
 auckland.nb <- poly2nb(auckland)
 res <- EBlocal(auckland$M77_85,  9*auckland$Und5_81, auckland.nb)
 auckland$est000 <- res$est*1000

diff --git a/man/SD.RStests.Rd b/man/SD.RStests.Rd
@@ -42,7 +42,7 @@ The results in the example below agree with those in Table 3, p. 22 in Koley and
 \code{\link{lm}}, \code{\link{lm.RStests}}
 }
 \examples{
-columbus <- sf::st_read(system.file("shapes/columbus.shp", package="spData")[1])
+columbus <- sf::st_read(system.file("shapes/columbus.gpkg", package="spData")[1])
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 col.listw <- nb2listw(col.gal.nb, style="W")
 lm_obj <- lm(CRIME ~ INC + HOVAL, data=columbus)

diff --git a/man/autocov_dist.Rd b/man/autocov_dist.Rd
@@ -35,7 +35,7 @@ for details.}
 
 \seealso{\code{\link{nb2listw}}}
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 #xy <- cbind(columbus$X, columbus$Y)
 xy <- st_coordinates(st_centroid(st_geometry(columbus),
  of_largest_polygon=TRUE))

diff --git a/man/bhicv.Rd b/man/bhicv.Rd
@@ -24,7 +24,7 @@
 %%\format{}
 %%\details{}
 \examples{
-bh <- st_read(system.file("etc/shapes/bhicv.shp",
+bh <- st_read(system.file("etc/shapes/bhicv.gpkg",
       package="spdep")[1])
 }
 \keyword{data}% at least one, from doc/KEYWORDS
diff --git a/man/choynowski.Rd b/man/choynowski.Rd
@@ -31,7 +31,7 @@ Interactive Spatial Data Analysis, Harlow: Longman, pp. 300--303.}
 
 \seealso{\code{\link{probmap}}}
 \examples{
-auckland <- st_read(system.file("shapes/auckland.shp", package="spData")[1], quiet=TRUE)
+auckland <- st_read(system.file("shapes/auckland.gpkg", package="spData")[1], quiet=TRUE)
 auckland.nb <- poly2nb(auckland)
 res <- choynowski(auckland$M77_85, 9*auckland$Und5_81)
 resl <- choynowski(auckland$M77_85, 9*auckland$Und5_81, legacy=TRUE)

diff --git a/man/columbus.Rd b/man/columbus.Rd
@@ -13,7 +13,7 @@
 The data set is now part of the spData package
 }
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 }
 

diff --git a/man/compon.Rd b/man/compon.Rd
@@ -22,7 +22,7 @@ A list of:
 \author{Nicholas Lewin-Koh \email{[email protected]}}
 \seealso{\code{\link{plot.nb}}}
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", 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)))
 plot(col.gal.nb, coords, col="grey")

diff --git a/man/diffnb.Rd b/man/diffnb.Rd
@@ -21,7 +21,7 @@ A neighbours list with class \code{nb}
 \author{Roger Bivand \email{[email protected]}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 coords <- st_centroid(st_geometry(columbus), of_largest_polygon=TRUE)
 rn <- row.names(columbus)
 knn1 <- knearneigh(coords, 1)

diff --git a/man/dnearneigh.Rd b/man/dnearneigh.Rd
@@ -32,7 +32,7 @@ for neighbours satisfying the distance criteria. See \code{\link{card}} for deta
 \seealso{\code{\link{knearneigh}}, \code{\link{card}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 coords <- st_centroid(st_geometry(columbus), of_largest_polygon=TRUE)
 rn <- row.names(columbus)
 k1 <- knn2nb(knearneigh(coords))

diff --git a/man/edit.nb.Rd b/man/edit.nb.Rd
@@ -31,7 +31,7 @@ The function returns an object of class \code{nb} with the edited list of intege
 
 \examples{
 \dontrun{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 class(columbus)
 if (FALSE) nnb1 <- edit.nb(col.gal.nb, polys=as(columbus, "Spatial"))
 }

diff --git a/man/globalG.test.Rd b/man/globalG.test.Rd
@@ -42,7 +42,7 @@ p. 276; Bivand RS, Wong DWS 2018 Comparing implementations of global and local i
 \seealso{\code{\link{localG}}}
 
 \examples{
-nc.sids <- st_read(system.file("shapes/sids.shp", package="spData")[1], quiet=TRUE)
+nc.sids <- st_read(system.file("shapes/sids.gpkg", package="spData")[1], quiet=TRUE)
 sidsrate79 <- (1000*nc.sids$SID79)/nc.sids$BIR79
 dists <- c(10, 20, 30, 33, 40, 50, 60, 70, 80, 90, 100)
 ndists <- length(dists)

diff --git a/man/graphneigh.Rd b/man/graphneigh.Rd
@@ -103,7 +103,7 @@ The helper functions return an \code{nb} object with a list of integer
 \code{\link{knn2nb}}, \code{\link{card}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 sf_obj <- st_centroid(st_geometry(columbus), of_largest_polygon)
 sp_obj <- as(sf_obj, "Spatial")
 coords <- st_coordinates(sf_obj)

diff --git a/man/include.self.Rd b/man/include.self.Rd
@@ -21,7 +21,7 @@ The function returns an object of class \code{nb} with a list of integer vectors
 \seealso{\code{\link{summary.nb}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 coords <- st_coordinates(st_centroid(columbus))
 summary(col.gal.nb, coords)

diff --git a/man/joincount.multi.Rd b/man/joincount.multi.Rd
@@ -34,7 +34,7 @@ A matrix with class \code{jcmulti} with row and column names for observed and ex
 \seealso{\code{\link{joincount.test}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 HICRIME <- cut(columbus$CRIME, breaks=c(0,35,80), labels=c("low","high"))
 (nb <- poly2nb(columbus))
 lw <- nb2listw(nb, style="B")

diff --git a/man/joincount.test.Rd b/man/joincount.test.Rd
@@ -62,7 +62,7 @@ print(is.symmetric.nb(COL.k4.nb))
 joincount.test(HICRIME, nb2listw(COL.k4.nb, style="B"))
 cat("Note non-symmetric weights matrix - use listw2U()\n")
 joincount.test(HICRIME, listw2U(nb2listw(COL.k4.nb, style="B")))
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 HICRIME <- cut(columbus$CRIME, breaks=c(0,35,80), labels=c("low","high"))
 (nb <- poly2nb(columbus))
 lw <- nb2listw(nb, style="B")

diff --git a/man/knearneigh.Rd b/man/knearneigh.Rd
@@ -32,7 +32,7 @@ A list of class \code{knn}
 \code{\link{knn2nb}}, \code{\link[dbscan]{kNN}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 coords <- st_centroid(st_geometry(columbus), of_largest_polygon=TRUE)
 col.knn <- knearneigh(coords, k=4)
 plot(st_geometry(columbus), border="grey")

diff --git a/man/knn2nb.Rd b/man/knn2nb.Rd
@@ -24,7 +24,7 @@ The function returns an object of class \code{nb} with a list of integer vectors
 \seealso{\code{\link{knearneigh}}, \code{\link{card}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 coords <- st_coordinates(st_centroid(columbus))
 col.knn <- knearneigh(coords, k=4)
 plot(st_geometry(columbus), border="grey")

diff --git a/man/listw2sn.Rd b/man/listw2sn.Rd
@@ -34,7 +34,7 @@ sn2listw(sn, style = NULL, zero.policy = NULL, from_mat2listw=FALSE)
 \seealso{\code{\link{nb2listw}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 col.listw <- nb2listw(col.gal.nb)
 col.listw$neighbours[[1]]

diff --git a/man/lm.morantest.exact.Rd b/man/lm.morantest.exact.Rd
@@ -52,7 +52,7 @@ must be one of greater (default), less or two.sided.}
 
 \seealso{\code{\link{lm.morantest.sad}}}
 \examples{
-eire <- st_read(system.file("shapes/eire.shp", package="spData")[1])
+eire <- st_read(system.file("shapes/eire.gpkg", package="spData")[1])
 row.names(eire) <- as.character(eire$names)
 st_crs(eire) <- "+proj=utm +zone=30 +ellps=airy +units=km"
 eire.nb <- poly2nb(eire)

diff --git a/man/lm.morantest.sad.Rd b/man/lm.morantest.sad.Rd
@@ -68,7 +68,7 @@ Geographical Analysis, 34, pp. 187--206; Bivand RS, Wong DWS 2018 Comparing impl
 \seealso{\code{\link{lm.morantest}}}
 
 \examples{
-eire <- st_read(system.file("shapes/eire.shp", package="spData")[1])
+eire <- st_read(system.file("shapes/eire.gpkg", package="spData")[1])
 row.names(eire) <- as.character(eire$names)
 st_crs(eire) <- "+proj=utm +zone=30 +ellps=airy +units=km"
 eire.nb <- poly2nb(eire)

diff --git a/man/localGS.Rd b/man/localGS.Rd
@@ -38,7 +38,7 @@ Sonea, A. and Westerholt, R. (2021): Geographic and temporal access to basic ban
 
 \examples{
 \donttest{
-boston.tr <- sf::st_read(system.file("shapes/boston_tracts.shp", package="spData")[1])
+boston.tr <- sf::st_read(system.file("shapes/boston_tracts.gpkg", package="spData")[1])
 boston.tr_utm <- st_transform(boston.tr, 32619) #26786
 
 boston_listw1 <- nb2listwdist(dnearneigh(st_centroid(boston.tr_utm), 1, 2000),

diff --git a/man/localmoran.exact.Rd b/man/localmoran.exact.Rd
@@ -72,7 +72,7 @@ When the alternative is being tested, a list of left and right M products in att
 
 \seealso{\code{\link{lm.morantest.exact}}, \code{\link{localmoran.sad}}}
 \examples{
-eire <- st_read(system.file("shapes/eire.shp", package="spData")[1])
+eire <- st_read(system.file("shapes/eire.gpkg", package="spData")[1])
 row.names(eire) <- as.character(eire$names)
 st_crs(eire) <- "+proj=utm +zone=30 +ellps=airy +units=km"
 eire.nb <- poly2nb(eire)

diff --git a/man/localmoran.sad.Rd b/man/localmoran.sad.Rd
@@ -105,7 +105,7 @@ Geographical Analysis, 34, pp. 187--206.}
  \code{\link{lm.morantest.sad}}, \code{\link{errorsarlm}}}
 
 \examples{
-eire <- st_read(system.file("shapes/eire.shp", package="spData")[1])
+eire <- st_read(system.file("shapes/eire.gpkg", package="spData")[1])
 row.names(eire) <- as.character(eire$names)
 st_crs(eire) <- "+proj=utm +zone=30 +ellps=airy +units=km"
 eire.nb <- poly2nb(eire)

diff --git a/man/localmoran_bv.Rd b/man/localmoran_bv.Rd
@@ -37,7 +37,7 @@ I_i^B= cx_i\Sigma_j{w_{ij}y_j}
 }
 \examples{
 # load columbus data
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData"))
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData"))
 nb <- poly2nb(columbus)
 listw <- nb2listw(nb)
 set.seed(1)

diff --git a/man/mat2listw.Rd b/man/mat2listw.Rd
@@ -33,7 +33,7 @@ A \code{listw} object with the following members:
 
 \seealso{\code{\link{nb2listw}}, \code{\link{nb2mat}}}
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col005 <- dnearneigh(st_coordinates(st_centroid(st_geometry(columbus),
  of_largest_polygon=TRUE)), 0, 0.5, as.character(columbus$NEIGNO))
 summary(col005)

diff --git a/man/moran.test.Rd b/man/moran.test.Rd
@@ -83,7 +83,7 @@ moran.test(crime, nb2listw(COL.nb, style="W"), zero.policy=TRUE,
 moran.test(crime, nb2listw(COL.nb, style="W"), zero.policy=TRUE,
  na.action=na.exclude)
 moran.test(crime, nb2listw(COL.nb, style="W"), na.action=na.pass)
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col_geoms <- st_geometry(columbus)
 col_geoms[1] <- st_buffer(col_geoms[1], dist=-0.05)
 st_geometry(columbus) <- col_geoms

diff --git a/man/mstree.Rd b/man/mstree.Rd
@@ -45,7 +45,7 @@ mstree(nbw, ini = NULL)
 %%\seealso{ ~~objects to See Also as \code{\link{help}}, ~~~ }
 \examples{
 ### loading data
-bh <- st_read(system.file("etc/shapes/bhicv.shp",
+bh <- st_read(system.file("etc/shapes/bhicv.gpkg",
       package="spdep")[1], quiet=TRUE)
 st_crs(bh) <- "OGC:CRS84"
 ### data padronized

diff --git a/man/nb2lines.Rd b/man/nb2lines.Rd
@@ -40,7 +40,7 @@ df2sn(df, i="i", i_ID="i_ID", j="j", wt="wt")
 \seealso{\code{\link{sn2listw}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 res <- listw2lines(nb2listw(col.gal.nb), st_geometry(columbus))
 summary(res)

diff --git a/man/nb2listw.Rd b/man/nb2listw.Rd
@@ -42,7 +42,7 @@ and attributes:
 \seealso{\code{\link{summary.nb}}, \code{\link{read.gal}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 coords <- st_coordinates(st_centroid(columbus))
 cards <- card(col.gal.nb)

diff --git a/man/nb2mat.Rd b/man/nb2mat.Rd
@@ -31,7 +31,7 @@ An n by n matrix, where n=length(neighbours)
 \seealso{\code{\link{nb2listw}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col005 <- dnearneigh(st_coordinates(st_centroid(st_geometry(columbus),
  of_largest_polygon=TRUE)), 0, 0.5, as.character(columbus$NEIGNO))
 summary(col005)

diff --git a/man/nbdists.Rd b/man/nbdists.Rd
@@ -24,7 +24,7 @@ A list with class \code{nbdist}
 \seealso{\code{\link{summary.nb}}, \code{\link{nb2listw}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 coords <- st_coordinates(st_centroid(columbus))
 dlist <- nbdists(col.gal.nb, coords)

diff --git a/man/nblag.Rd b/man/nblag.Rd
@@ -23,7 +23,7 @@ returns a list of lagged neighbours lists each with class \code{nb}
 \seealso{\code{\link{summary.nb}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 coords <- st_coordinates(st_centroid(columbus))
 summary(col.gal.nb, coords)

diff --git a/man/nboperations.Rd b/man/nboperations.Rd
@@ -39,7 +39,7 @@ complement.nb(nb.obj)
 \code{\link{setdiff.nb}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 coords <- st_coordinates(st_centroid(columbus))
 col.tri.nb <- tri2nb(coords)

diff --git a/man/plot.nb.Rd b/man/plot.nb.Rd
@@ -30,7 +30,7 @@ neighbours lists}
 \seealso{\code{\link{summary.nb}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", package="spData")[1], quiet=TRUE)
 col.gal.nb <- read.gal(system.file("weights/columbus.gal", package="spData")[1])
 plot(col.gal.nb, st_geometry(columbus))
 title(main="GAL order 1 links with first nearest neighbours in red", cex.main=0.6)

diff --git a/man/poly2nb.Rd b/man/poly2nb.Rd
@@ -32,7 +32,7 @@ From 1.1-9 with \pkg{sf} 1.0-0, \pkg{s2} is used in bounding box indexing intern
 \seealso{\code{\link{summary.nb}}, \code{\link{card}}}
 
 \examples{
-columbus <- st_read(system.file("shapes/columbus.shp", package="spData")[1], quiet=TRUE)
+columbus <- st_read(system.file("shapes/columbus.gpkg", 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)))
 xx <- poly2nb(as(columbus, "Spatial"))
@@ -62,7 +62,7 @@ fg[maxconts] <- "red"
 fg[xx[[maxconts]]] <- "green"
 plot(st_geometry(columbus), col=fg)
 title(main="Region with largest number of contiguities", cex.main=0.6)
-nc.sids <- st_read(system.file("shapes/sids.shp", package="spData")[1], quiet=TRUE)
+nc.sids <- st_read(system.file("shapes/sids.gpkg", package="spData")[1], quiet=TRUE)
 system.time(xxnb <- poly2nb(nc.sids))
 system.time(xxnb <- poly2nb(as(nc.sids, "Spatial")))
 plot(st_geometry(nc.sids))

diff --git a/man/probmap.Rd b/man/probmap.Rd
@@ -29,7 +29,7 @@ probmap(n, x, row.names=NULL, alternative="less")
 \seealso{\code{\link{EBest}}, \code{\link{EBlocal}}, \code{\link{ppois}}}
 
 \examples{
-auckland <- st_read(system.file("shapes/auckland.shp", package="spData")[1], quiet=TRUE)
+auckland <- st_read(system.file("shapes/auckland.gpkg", package="spData")[1], quiet=TRUE)
 res <- probmap(auckland$M77_85, 9*auckland$Und5_81)
 rt <- sum(auckland$M77_85)/sum(9*auckland$Und5_81)
 ppois_pmap <- numeric(length(auckland$Und5_81))

diff --git a/man/skater.Rd b/man/skater.Rd
@@ -72,7 +72,7 @@ skater(edges, data, ncuts, crit, vec.crit, method = c("euclidean",
 \seealso{See Also as \code{\link{mstree}}}
 \examples{
 ### loading data
-bh <- st_read(system.file("etc/shapes/bhicv.shp",
+bh <- st_read(system.file("etc/shapes/bhicv.gpkg",
       package="spdep")[1], quiet=TRUE)
 st_crs(bh) <- "OGC:CRS84"
 ### data standardized 

diff --git a/man/sp.correlogram.Rd b/man/sp.correlogram.Rd
@@ -52,7 +52,7 @@ of the Institute of British Geographers}, 66, 95--118.}
 \seealso{\code{\link{nblag}}, \code{\link{moran}}, \code{\link[stats]{p.adjust}}}
 
 \examples{
-nc.sids <- st_read(system.file("shapes/sids.shp", package="spData")[1], quiet=TRUE)
+nc.sids <- st_read(system.file("shapes/sids.gpkg", package="spData")[1], quiet=TRUE)
 rn <- as.character(nc.sids$FIPS)
 ncCC89_nb <- read.gal(system.file("weights/ncCC89.gal", package="spData")[1],
  region.id=rn)