## Not run:
# example(columbus)
# lmbase <- lm(CRIME ~ INC + HOVAL, data=columbus)
# lagcol <- SpatialFiltering(CRIME ~ 1, ~ INC + HOVAL, data=columbus,
# nb=col.gal.nb, style="W", alpha=0.1, verbose=TRUE)
# lagcol
# lmlag <- lm(CRIME ~ INC + HOVAL + fitted(lagcol), data=columbus)
# anova(lmlag)
# anova(lmbase, lmlag)
# set.seed(123)
# lagcol1 <- ME(CRIME ~ INC + HOVAL, data=columbus, family="gaussian",
# listw=nb2listw(col.gal.nb), alpha=0.1, verbose=TRUE)
# lagcol1
# lmlag1 <- lm(CRIME ~ INC + HOVAL + fitted(lagcol1), data=columbus)
# anova(lmlag1)
# anova(lmbase, lmlag1)
# set.seed(123)
# lagcol2 <- ME(CRIME ~ INC + HOVAL, data=columbus, family="gaussian",
# listw=nb2listw(col.gal.nb), alpha=0.1, stdev=TRUE, verbose=TRUE)
# lagcol2
# lmlag2 <- lm(CRIME ~ INC + HOVAL + fitted(lagcol2), data=columbus)
# anova(lmlag2)
# anova(lmbase, lmlag2)
# example(nc.sids)
# glmbase <- glm(SID74 ~ 1, data=nc.sids, offset=log(BIR74),
# family="poisson")
# set.seed(123)
# MEpois1 <- ME(SID74 ~ 1, data=nc.sids, offset=log(BIR74),
# family="poisson", listw=nb2listw(ncCR85_nb, style="B"), alpha=0.2, verbose=TRUE)
# MEpois1
# glmME <- glm(SID74 ~ 1 + fitted(MEpois1), data=nc.sids, offset=log(BIR74),
# family="poisson")
# anova(glmME, test="Chisq")
# anova(glmbase, glmME, test="Chisq")
# data(hopkins)
# hopkins_part <- hopkins[21:36,36:21]
# hopkins_part[which(hopkins_part > 0, arr.ind=TRUE)] <- 1
# hopkins.rook.nb <- cell2nb(16, 16, type="rook")
# glmbase <- glm(c(hopkins_part) ~ 1, family="binomial")
# set.seed(123)
# MEbinom1 <- ME(c(hopkins_part) ~ 1, family="binomial",
# listw=nb2listw(hopkins.rook.nb, style="B"), alpha=0.2, verbose=TRUE)
# glmME <- glm(c(hopkins_part) ~ 1 + fitted(MEbinom1), family="binomial")
# anova(glmME, test="Chisq")
# anova(glmbase, glmME, test="Chisq")
# ## End(Not run)
Run the code above in your browser using DataLab