1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
## ----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)
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
respatlag <- spml(fm, data = Produc, listw = mat2listw(usaww),
model="random", spatial.error="none", lag=TRUE)
class(respatlag)
GM <- spgm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, data = Produc,
listw = usaww, moments = "fullweights", spatial.error = TRUE)
class(GM)
imp.spml <- impacts(respatlag, listw = mat2listw(usaww, style = "W"), time = 17)
class(imp.spml)
## ----xtablesplm, results = 'asis'-------------------------
xtable(respatlag)
## ----xtablesplm1, results = 'asis'------------------------
xtable(GM)
## ----xtablesplmimpacts, results = 'asis'------------------
xtable(imp.spml)
## ----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)
|
