Applied patches supplied by Martin Gubri for spdep

[xtable.git] / pkg / vignettes / OtherPackagesGallery.Rnw

diff --git a/pkg/vignettes/OtherPackagesGallery.Rnw b/pkg/vignettes/OtherPackagesGallery.Rnw
index ffd5f9523acf3716b88bc8fefd331fe3917d650c..1676da4fe317913ebcb270faac27a9df2b833aec 100644 (file)
--- a/pkg/vignettes/OtherPackagesGallery.Rnw
+++ b/pkg/vignettes/OtherPackagesGallery.Rnw
@@ -1,4 +1,4 @@
-%\VignetteIndexEntry{xtable List of Tables Gallery}\r
+%\VignetteIndexEntry{xtable Other Packages Gallery}\r
 %\VignetteDepends{xtable, spdep, splm, sphet}\r
 %\VignetteKeywords{LaTeX, HTML, table}\r
 %\VignettePackage{xtable}\r
@@ -62,28 +62,34 @@ First load the package and create some objects.
 <<dataspdep>>=\r
 library(spdep)\r
 data("oldcol", package = "spdep")\r
-COL.lag.eig <- lagsarlm(CRIME ~ INC + HOVAL, data = COL.OLD[],\r
-                        nb2listw(COL.nb))\r
+\r
+data.in.sample <- COL.OLD[1:44,]\r
+data.out.of.sample <- COL.OLD[45:49,]\r
+\r
+listw.in.sample <- nb2listw(subset(COL.nb, !(1:49 %in% 45:49)))\r
+listw.all.sample <- nb2listw(COL.nb)\r
+\r
+COL.lag.eig <- lagsarlm(CRIME ~ INC + HOVAL, data = data.in.sample,\r
+                        listw.in.sample)\r
 class(COL.lag.eig)\r
-COL.errW.GM <- GMerrorsar(CRIME ~ INC + HOVAL, data = COL.OLD,\r
-                          nb2listw(COL.nb, style = "W"),\r
-                          returnHcov = TRUE)\r
+COL.errW.GM <- GMerrorsar(CRIME ~ INC + HOVAL, data = data.in.sample,\r
+                          listw.in.sample, returnHcov = TRUE)\r
 class(COL.errW.GM)\r
-COL.lag.stsls <- stsls(CRIME ~ INC + HOVAL, data = COL.OLD,\r
-                       nb2listw(COL.nb))\r
+COL.lag.stsls <- stsls(CRIME ~ INC + HOVAL, data = data.in.sample,\r
+                       listw.in.sample)\r
 class(COL.lag.stsls)\r
 \r
-p1 <- predict(COL.lag.eig, newdata = COL.OLD[45:49,],\r
-              listw = nb2listw(COL.nb))\r
+p1 <- predict(COL.lag.eig, newdata = data.out.of.sample,\r
+              listw = listw.all.sample)\r
 class(p1)\r
-p2 <- predict(COL.lag.eig, newdata = COL.OLD[45:49,],\r
+p2 <- predict(COL.lag.eig, newdata = data.out.of.sample,\r
               pred.type = "trend", type = "trend")\r
 #type option for retrocompatibility with spdep 0.5-92\r
 class(p2)\r
 \r
-imp.exact <- impacts(COL.lag.eig, listw = nb2listw(COL.nb))\r
+imp.exact <- impacts(COL.lag.eig, listw = listw.in.sample)\r
 class(imp.exact)\r
-imp.sim <- impacts(COL.lag.eig, listw = nb2listw(COL.nb), R = 200)\r
+imp.sim <- impacts(COL.lag.eig, listw = listw.in.sample, R = 200)\r
 class(imp.sim)\r
 @ %def\r
 \r