example(NY_data)
## Not run:
# esar1f <- spautolm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, family="SAR", method="eigen")
# summary(esar1f)
# res <- MCMCsamp(esar1f, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# esar1fw <- spautolm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, weights=POP8, family="SAR", method="eigen")
# summary(esar1fw)
# res <- MCMCsamp(esar1fw, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# ecar1f <- spautolm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, family="CAR", method="eigen")
# summary(ecar1f)
# res <- MCMCsamp(ecar1f, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# esar1fw <- spautolm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, weights=POP8, family="SAR", method="eigen")
# summary(esar1fw)
# res <- MCMCsamp(esar1fw, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# ecar1fw <- spautolm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, weights=POP8, family="CAR", method="eigen")
# summary(ecar1fw)
# res <- MCMCsamp(ecar1fw, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# ## End(Not run)
esar0 <- errorsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
listw=listw_NY)
summary(esar0)
res <- MCMCsamp(esar0, mcmc=5000, burnin=500, listw=listw_NY)
summary(res)
## Not run:
# esar0w <- errorsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, weights=POP8)
# summary(esar0)
# res <- MCMCsamp(esar0w, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# esar1 <- errorsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, etype="emixed")
# summary(esar1)
# res <- MCMCsamp(esar1, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# lsar0 <- lagsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY)
# summary(lsar0)
# res <- MCMCsamp(lsar0, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# lsar1 <- lagsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, type="mixed")
# summary(lsar1)
# res <- MCMCsamp(lsar1, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# ssar0 <- sacsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY)
# summary(ssar0)
# res <- MCMCsamp(ssar0, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# ssar1 <- sacsarlm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data=nydata,
# listw=listw_NY, type="sacmixed")
# summary(ssar1)
# res <- MCMCsamp(ssar1, mcmc=5000, burnin=500, listw=listw_NY)
# summary(res)
# ## End(Not run)
Run the code above in your browser using DataLab