inst/doc/OtherPackagesGallery.R

In xtable: Export Tables to LaTeX or HTML

## ----set, include=FALSE-------------------------------------------------------
library(knitr)
opts_chunk$set(fig.path = 'Figures/other', debug = TRUE, echo = TRUE)
opts_chunk$set(out.width = '0.9\\textwidth')

## ----package, results='asis'------------------------------
library(xtable)
options(xtable.floating = FALSE)
options(xtable.timestamp = "")
options(width = 60)
set.seed(1234)

## ----dataspdep--------------------------------------------
library(spdep)
library(spatialreg, warn.conflicts = FALSE)
data("oldcol", package = "spdep")

data.in.sample <- COL.OLD[1:44,]
data.out.of.sample <- COL.OLD[45:49,]

listw.in.sample <- nb2listw(subset(COL.nb, !(1:49 %in% 45:49)))
listw.all.sample <- nb2listw(COL.nb)

COL.lag.eig <- lagsarlm(CRIME ~ INC + HOVAL, data = data.in.sample,
                        listw.in.sample)
class(COL.lag.eig)
COL.errW.GM <- GMerrorsar(CRIME ~ INC + HOVAL, data = data.in.sample,
                          listw.in.sample, returnHcov = TRUE)
class(COL.errW.GM)
COL.lag.stsls <- stsls(CRIME ~ INC + HOVAL, data = data.in.sample,
                       listw.in.sample)
class(COL.lag.stsls)

p1 <- predict(COL.lag.eig, newdata = data.out.of.sample,
              listw = listw.all.sample)
class(p1)
p2 <- predict(COL.lag.eig, newdata = data.out.of.sample,
              pred.type = "trend", type = "trend")
#type option for retrocompatibility with spdep 0.5-92
class(p2)

imp.exact <- impacts(COL.lag.eig, listw = listw.in.sample)
class(imp.exact)
imp.sim <- impacts(COL.lag.eig, listw = listw.in.sample, R = 200)
class(imp.sim)

## ----xtablesarlm, results = 'asis'------------------------
xtable(COL.lag.eig)

## ----xtablesarlmsumm, results = 'asis'--------------------
xtable(summary(COL.lag.eig, correlation = TRUE))

## ----xtablesarlmbooktabs, results = 'asis'----------------
print(xtable(COL.lag.eig), booktabs = TRUE)

## ----xtablegmsar, results = 'asis'------------------------
xtable(COL.errW.GM)

## ----xtablestsls, results = 'asis'------------------------
xtable(COL.lag.stsls)

## ----xtablesarlmpred, results = 'asis'--------------------
xtable(p1)

## ----xtablesarlmpred2, results = 'asis'-------------------
xtable(p2)

## ----xtablelagimpactexact, results = 'asis'---------------
xtable(imp.exact)

## ----xtablelagimpactmcmc, results = 'asis'----------------
xtable(imp.sim)

## ----minimalexample, results = 'hide'---------------------
library(spdep)
example(NY_data)
spautolmOBJECT <- spautolm(Z ~ PEXPOSURE + PCTAGE65P,data = nydata,
                           listw = listw_NY, family = "SAR",
                           method = "eigen", verbose = TRUE)
summary(spautolmOBJECT, Nagelkerke = TRUE)

## ----spautolmclass----------------------------------------
class(spautolmOBJECT)

## ----xtablespautolm, results = 'asis'---------------------
xtable(spautolmOBJECT,
       display = c("s",rep("f", 3), "e"), digits = 4)

## ----datasplm---------------------------------------------
library(splm)
data("Produc", package = "plm")
data("usaww",  package = "splm")
fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp
## the two standard specifications (SEM and SAR) one with FE
## and the other with RE:
## fixed effects panel with spatial errors
fespaterr <- spml(fm, data = Produc, listw = mat2listw(usaww),
                   model="within", spatial.error="b", Hess = FALSE)
summary(fespaterr)
## random effects panel with spatial lag
respatlag <- spml(fm, data = Produc, listw = mat2listw(usaww),
                   model="random", spatial.error="none", lag=TRUE)
summary(respatlag)
## calculate impact measures
impac1 <- splm:::impacts.splm(respatlag, 
                              listw = mat2listw(usaww, style = "W"), 
                              time = 17)
summary(impac1, zstats = TRUE, short = TRUE)

GM <- spgm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, data = Produc,
           listw = usaww, moments = "fullweights", spatial.error = TRUE)
class(GM)

## ----xtablesplm, results = 'asis'-------------------------
xtable(respatlag)

## ----xtablegm, results = 'asis'---------------------------
xtable(GM)

## ----xtablesplmimpacts, results = 'asis'------------------
xtable(impac1)

## ----datasphet--------------------------------------------
library(sphet)
data("columbus", package = "spdep")
listw <- nb2listw(col.gal.nb)
data("coldis", package = "sphet")
res.stsls <- stslshac(CRIME ~ HOVAL + INC, data = columbus, listw = listw,
                      distance = coldis, type = 'Triangular')
class(res.stsls)

## ----xtablesphet, results = 'asis'------------------------
xtable(res.stsls)

## ----zoo, results = 'asis'--------------------------------
library(zoo)
xDate <- as.Date("2003-02-01") + c(1, 3, 7, 9, 14) - 1
as.ts(xDate)
x <- zoo(rnorm(5), xDate)
xtable(x)

## ----zoots, results = 'asis'------------------------------
tempTs <- ts(cumsum(1 + round(rnorm(100), 0)),
              start = c(1954, 7), frequency = 12)
tempTable <- xtable(tempTs, digits = 0)
tempTable
tempZoo <- as.zoo(tempTs)
xtable(tempZoo, digits = 0)

## ----survival, results = 'asis'---------------------------
library(survival)
test1 <- list(time=c(4,3,1,1,2,2,3),
              status=c(1,1,1,0,1,1,0),
              x=c(0,2,1,1,1,0,0),
              sex=c(0,0,0,0,1,1,1))
coxFit <- coxph(Surv(time, status) ~ x + strata(sex), test1)
xtable(coxFit)