search

r-spatial / spatialreg

spatialreg: spatial models estimation and testing
https://r-spatial.github.io/spatialreg/
41 stars 12 forks source link

na.action and precomputed eigenvalue bug #11

Closed rsbivand closed 4 years ago

rsbivand commented 4 years ago

See https://github.com/r-spatial/spatialreg/issues/10:

data(oldcol, package="spdep")
listw <- spdep::nb2listw(COL.nb, style="W")
ev <- eigenw(listw)
anyNA(COL.OLD$INC)
COL.OLD1 <- COL.OLD
is.na(COL.OLD1$INC) <- 10L
anyNA(COL.OLD1$INC)
res <- lagsarlm(CRIME ~ INC + HOVAL, data=COL.OLD1, listw, control=list(pre_eig=ev))
# Error in eigen_pre_setup(env, pre_eig = pre_eig, which = which) : 
#  length(pre_eig) == get("n", envir = env) is not TRUE