diff --git a/DESCRIPTION b/DESCRIPTION
index ccc3d9a..8122b04 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: sphet
-Version: 2.0.6
+Version: 2.0.7
 Encoding: UTF-8
-Date: 2024-02-07
+Date: 2024-06-08
 Title: Estimation of Spatial Autoregressive Models with and without Heteroskedastic Innovations
 Authors@R: c(person("Gianfranco", "Piras", role = c("aut", "cre"), 
             email = "gpiras@mac.com", 
diff --git a/R/impacts.R b/R/impacts.R
index f0b5fe7..8d996c3 100644
--- a/R/impacts.R
+++ b/R/impacts.R
@@ -344,7 +344,7 @@ if((lambda > interval[2] ) | (lambda < interval[1])) warning("Value of the spati
 #' @examples
 #' require("sf", quietly=TRUE)
 #' library(coda)
-#' 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 <- spdep::read.gal(system.file("weights/columbus.gal", package="spData")[1])
 #' listw <- spdep::nb2listw(col.gal.nb)
 #' ev <- spatialreg::eigenw(listw)
@@ -689,7 +689,7 @@ impacts.stsls_sphet <- function(obj, ..., tr=NULL, R=NULL, listw=NULL,
 #' summary(impacts(lm.D3))
 #' 
 #' require("sf", quietly=TRUE)
-#' 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 <- spdep::read.gal(system.file("weights/columbus.gal", package="spData")[1])
 #' listw <- spdep::nb2listw(col.gal.nb)
 #' knear <- spdep::knearneigh(cbind(columbus$X, columbus$Y), 5)
@@ -933,7 +933,7 @@ impacts.ols_sphet <- function(obj, ..., tr=NULL, R=NULL, listw=NULL,
 #' @examples
 #' library(sphet)
 #' require("sf", quietly=TRUE)
-#' 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 <- spdep::read.gal(system.file("weights/columbus.gal", package="spData")[1])
 #' listw <- spdep::nb2listw(col.gal.nb)
 #' error1 <- spreg(CRIME ~ INC + HOVAL, columbus, listw, Durbin=TRUE,
diff --git a/man/impacts.error_sphet.Rd b/man/impacts.error_sphet.Rd
index 3470437..b3dcf87 100644
--- a/man/impacts.error_sphet.Rd
+++ b/man/impacts.error_sphet.Rd
@@ -44,7 +44,7 @@ Generate impacts for objects of class error_sphet created in sphet
 \examples{
 library(sphet)
 require("sf", quietly=TRUE)
-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 <- spdep::read.gal(system.file("weights/columbus.gal", package="spData")[1])
 listw <- spdep::nb2listw(col.gal.nb)
 error1 <- spreg(CRIME ~ INC + HOVAL, columbus, listw, Durbin=TRUE,
diff --git a/man/impacts.ols_sphet.Rd b/man/impacts.ols_sphet.Rd
index 2ecc2da..1413448 100644
--- a/man/impacts.ols_sphet.Rd
+++ b/man/impacts.ols_sphet.Rd
@@ -66,7 +66,7 @@ impacts(lm.D3)
 summary(impacts(lm.D3))
 
 require("sf", quietly=TRUE)
-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 <- spdep::read.gal(system.file("weights/columbus.gal", package="spData")[1])
 listw <- spdep::nb2listw(col.gal.nb)
 knear <- spdep::knearneigh(cbind(columbus$X, columbus$Y), 5)
diff --git a/man/impacts.stsls_sphet.Rd b/man/impacts.stsls_sphet.Rd
index 05c884c..5ff96e7 100644
--- a/man/impacts.stsls_sphet.Rd
+++ b/man/impacts.stsls_sphet.Rd
@@ -50,7 +50,7 @@ Generate impacts for objects of class lag_gmm created in sphet
 \examples{
 require("sf", quietly=TRUE)
 library(coda)
-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 <- spdep::read.gal(system.file("weights/columbus.gal", package="spData")[1])
 listw <- spdep::nb2listw(col.gal.nb)
 ev <- spatialreg::eigenw(listw)