move Durbin factor functions to spdep

Author: rsbivand

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: spatialreg
 Version: 1.3-7
-Date: 2025-02-28
+Date: 2025-04-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")),
@@ -27,7 +27,7 @@ Authors@R: c(person("Roger", "Bivand", role = c("cre", "aut"), email = "Roger.Bi
              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.1), Matrix, sf
-Imports: spdep (>= 1.3-1), coda, methods, MASS, boot, splines, LearnBayes,
+Imports: spdep (>= 1.3-11), coda, methods, MASS, boot, splines, LearnBayes,
         nlme, multcomp
 Suggests: parallel, RSpectra, tmap, foreign, spam, knitr, lmtest, expm,
         sandwich, rmarkdown, igraph (>= 2.0.0), tinytest

NAMESPACE

@@ -16,7 +16,8 @@ importFrom("stats", "lm.fit", "as.formula", "coef", "lm", "model.extract",
 importFrom("spdep", "is.symmetric.glist", "lag.listw", "card",
  "nb2listw", "n.comp.nb", "listw2U", "listw2mat", "Szero", "listw2sn",
  "mat2listw", "is.symmetric.nb", "chkIDs", "get.spChkOption",
- "subset.listw", "spweights.constants")
+ "subset.listw", "spweights.constants", "have_factor_preds_mf",
+ "warn_factor_preds")
 
 import(Matrix#, except=c("expm")
 )

R/cat_durbin.R

@@ -2,30 +2,30 @@ library(spatialreg)
 data(oldcol, package="spdep")
 lw <- spdep::nb2listw(COL.nb)
 COL.OLD$fEW <- factor(COL.OLD$EW)
-COL.OLD$fDISCBD <- cut(COL.OLD$DISCBD, c(0, 1.5, 3, 4.5, 6))
+COL.OLD$fDISCBD <- ordered(cut(COL.OLD$DISCBD, c(0, 1.5, 3, 4.5, 6)))
 f <- formula(CRIME ~ INC + HOVAL + fDISCBD*fEW)
-expect_warning(COL.SLX <- lmSLX(f, data=COL.OLD, lw, Durbin=TRUE))
-expect_warning(COL.SLX <- lmSLX(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
-expect_warning(COL.SLX <- lmSLX(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fEW))
-expect_silent(COL.SLX <- lmSLX(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
-expect_warning(COL.err <- errorsarlm(f, data=COL.OLD, lw, Durbin=TRUE))
-expect_warning(COL.err <- errorsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
-expect_warning(COL.err <- errorsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
-expect_silent(COL.err <- errorsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
-expect_warning(COL.lag <- lagsarlm(f, data=COL.OLD, lw, Durbin=TRUE))
-expect_warning(COL.lag <- lagsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
-expect_warning(COL.lag <- lagsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
-expect_silent(COL.lag <- lagsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
-expect_warning(COL.sac <- sacsarlm(f, data=COL.OLD, lw, Durbin=TRUE))
-expect_warning(COL.sac <- sacsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
-expect_warning(COL.sac <- sacsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
-expect_silent(COL.sac <- sacsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
-expect_warning(COL.lag <- spBreg_lag(f, data=COL.OLD, lw, Durbin=TRUE))
-expect_warning(COL.lag <- spBreg_lag(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
-expect_warning(COL.lag <- spBreg_lag(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
-expect_silent(COL.lag <- spBreg_lag(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
-expect_warning(COL.err <- spBreg_err(f, data=COL.OLD, lw, Durbin=TRUE))
-expect_warning(COL.err <- spBreg_err(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
-expect_warning(COL.err <- spBreg_err(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
-expect_silent(COL.err <- spBreg_err(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
+expect_warning(COL.SLX0 <- lmSLX(f, data=COL.OLD, lw, Durbin=TRUE))
+expect_warning(COL.SLX1 <- lmSLX(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
+expect_warning(COL.SLX2 <- lmSLX(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fEW))
+expect_silent(COL.SLX3 <- lmSLX(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
+expect_warning(COL.err0 <- errorsarlm(f, data=COL.OLD, lw, Durbin=TRUE))
+expect_warning(COL.err1 <- errorsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
+expect_warning(COL.err2 <- errorsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
+expect_silent(COL.err3 <- errorsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
+expect_warning(COL.lag0 <- lagsarlm(f, data=COL.OLD, lw, Durbin=TRUE))
+expect_warning(COL.lag1 <- lagsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
+expect_warning(COL.lag2 <- lagsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
+expect_silent(COL.lag3 <- lagsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
+expect_warning(COL.sac0 <- sacsarlm(f, data=COL.OLD, lw, Durbin=TRUE))
+expect_warning(COL.sac1 <- sacsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
+expect_warning(COL.sac2 <- sacsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
+expect_silent(COL.sac3 <- sacsarlm(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
+expect_warning(COL.lag0 <- spBreg_lag(f, data=COL.OLD, lw, Durbin=TRUE))
+expect_warning(COL.lag1 <- spBreg_lag(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
+expect_warning(COL.lag2 <- spBreg_lag(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
+expect_silent(COL.lag3 <- spBreg_lag(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))
+expect_warning(COL.err0 <- spBreg_err(f, data=COL.OLD, lw, Durbin=TRUE))
+expect_warning(COL.err1 <- spBreg_err(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD*fEW))
+expect_warning(COL.err2 <- spBreg_err(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL + fDISCBD))
+expect_silent(COL.err3 <- spBreg_err(f, data=COL.OLD, lw, Durbin=~ INC + HOVAL))