From 8e368868df629c1bf754340a5b7e614357c968be Mon Sep 17 00:00:00 2001 From: Roger Bivand Date: Mon, 10 Jun 2024 13:13:08 +0200 Subject: [PATCH] fixes #50 --- DESCRIPTION | 4 ++-- vignettes/SpatialFiltering.Rmd | 22 ++++++++++++---------- 2 files changed, 14 insertions(+), 12 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index 8e874c5..6207cf9 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,6 +1,6 @@ Package: spatialreg Version: 1.3-4 -Date: 2024-06-08 +Date: 2024-06-10 Title: Spatial Regression Analysis Encoding: UTF-8 Authors@R: c(person("Roger", "Bivand", role = c("cre", "aut"), email = "Roger.Bivand@nhh.no", comment=c(ORCID="0000-0003-2392-6140")), @@ -25,7 +25,7 @@ Authors@R: c(person("Roger", "Bivand", role = c("cre", "aut"), email = "Roger.Bi person(given = "Mauricio", family = "Sarrias", role = c("ctb"), email = "mauricio.sarrias@ucn.cl"), person(given = "JuanTomas", family = "Sayago", role = c("ctb"), email = "juantomas.sayago@gmail.com"), person("Michael", "Tiefelsdorf", role = "ctb")) -Depends: R (>= 3.3.0), spData (>= 2.3.2), Matrix, sf +Depends: R (>= 3.3.0), spData (>= 2.3.1), Matrix, sf Imports: spdep (>= 1.3-1), coda, methods, MASS, boot, splines, LearnBayes, nlme, multcomp Suggests: parallel, RSpectra, tmap, foreign, spam, knitr, lmtest, expm, diff --git a/vignettes/SpatialFiltering.Rmd b/vignettes/SpatialFiltering.Rmd index e755306..11d9534 100644 --- a/vignettes/SpatialFiltering.Rmd +++ b/vignettes/SpatialFiltering.Rmd @@ -54,7 +54,13 @@ by Pedro Peres-Neto and is discussed in @dray+legendre+peres-neto:06 and ```{r echo=TRUE} library(spdep) require("sf", quietly=TRUE) -NY8 <- st_read(system.file("shapes/NY8_utm18.gpkg", package="spData")) +if (packageVersion("spData") >= "2.3.2") { + NY8 <- sf::st_read(system.file("shapes/NY8_utm18.gpkg", package="spData")) +} else { + NY8 <- sf::st_read(system.file("shapes/NY8_bna_utm18.gpkg", package="spData")) + sf::st_crs(NY8) <- "EPSG:32618" + NY8$Cases <- NY8$TRACTCAS +} NY_nb <- read.gal(system.file("weights/NY_nb.gal", package="spData"), override.id=TRUE) ``` @@ -89,22 +95,18 @@ set.seed(111) ``` ```{r, echo=TRUE} -nyME <- ME(Cases~PEXPOSURE+PCTAGE65P+PCTOWNHOME, data=NY8, offset=log(POP8), family="poisson", listw=NYlistwW, alpha=0.44) +nyME <- ME(Cases~PEXPOSURE+PCTAGE65P+PCTOWNHOME, data=NY8, offset=log(POP8), family="poisson", listw=NYlistwW, alpha=0.46) ``` ```{r, echo=TRUE} nyME -NY8$eigen_24 <- fitted(nyME)[,1] -NY8$eigen_223 <- fitted(nyME)[,2] +NY8$eigen_1 <- fitted(nyME)[,1] +NY8$eigen_2 <- fitted(nyME)[,2] ``` -```{r, results='asis',echo=FALSE, eval=FALSE} -.iwidth <- 6 -.iheight <- 4 -.ipointsize <- 10 -library(RColorBrewer) +```{r, echo=TRUE} #gry <- brewer.pal(9, "Greys")[-1] -spplot(NY8, c("eigen_24", "eigen_223"), col.regions=grey.colors(6, 0.95, 0.55, 2.2), cuts=5) +plot(NY8[,c("eigen_1", "eigen_2")]) ``` ```{r, echo=TRUE}