if (requireNamespace("spdep", quietly = TRUE)) {
library(spdep)
data(used.cars)
moran.test(used.cars$price.1960, nb2listw(usa48.nb))
moran.plot(used.cars$price.1960, nb2listw(usa48.nb),
labels=rownames(used.cars))
uc.lm <- lm(price.1960 ~ tax.charges, data=used.cars)
summary(uc.lm)
lm.morantest(uc.lm, nb2listw(usa48.nb))
lm.morantest.sad(uc.lm, nb2listw(usa48.nb))
lm.LMtests(uc.lm, nb2listw(usa48.nb))
#if (requireNamespace("spatialreg", quietly = TRUE)) {
# library(spatialreg)
# uc.err <- errorsarlm(price.1960 ~ tax.charges, data=used.cars,
# nb2listw(usa48.nb), tol.solve=1.0e-13,
# control=list(tol.opt=.Machine$double.eps^0.3))
# summary(uc.err)
# uc.lag <- lagsarlm(price.1960 ~ tax.charges, data=used.cars,
# nb2listw(usa48.nb), tol.solve=1.0e-13,
# control=list(tol.opt=.Machine$double.eps^0.3))
# summary(uc.lag)
# uc.lag1 <- lagsarlm(price.1960 ~ 1, data=used.cars,
# nb2listw(usa48.nb), tol.solve=1.0e-13,
# control=list(tol.opt=.Machine$double.eps^0.3))
# summary(uc.lag1)
# uc.err1 <- errorsarlm(price.1960 ~ 1, data=used.cars,
# nb2listw(usa48.nb), tol.solve=1.0e-13,
# control=list(tol.opt=.Machine$double.eps^0.3))
# summary(uc.err1)
#}
}
Run the code above in your browser using DataLab