@@ -44,10 +44,10 @@ Now we can plot the data and its empirical autocorrelation function.
4444
4545``` {r, echo = -c(1:2)}
4646if (pdfplots) {
47- pdf("7-2_1.pdf", width = 8 , height = 5)
48- par(mar = c(2.5 , 1.5, 1.5, .1), mgp = c(1.5, .5, 0), lwd = 2 )
47+ pdf("7-2_1.pdf", width = 12 , height = 5)
48+ par(mar = c(2.6 , 1.5, 1.5, .1), mgp = c(1.5, .5, 0), lwd = 1.5 )
4949}
50- par(mfrow = c(1,2))
50+ par(mfrow = c(1, 2))
5151ts.plot(logret, main = "U.S. GDP log returns")
5252acf(logret, lag = 8, main = "")
5353title("Empirical autocorrelation function")
@@ -184,7 +184,7 @@ corresponding to the highest lag in each of the models.
184184if (pdfplots) {
185185 pdf("7-2_2.pdf", width = 8, height = 5)
186186}
187- par(mfrow = c(2, 2), mar = c(2.5, 1.5, 1.5, .1), mgp = c(1.5, .5, 0))
187+ par(mfrow = c(2, 2), mar = c(2.5, 1.5, 1.5, .1), mgp = c(1.5, .5, 0), lwd = 1.5 )
188188for (p in 1:4) {
189189 hist(res[[p]]$betas[, p + 1], freq = FALSE, main = bquote(AR(.(p))),
190190 xlab = bquote(phi[.(p)]), ylab = "", breaks = seq(-.75, .75, .02))
@@ -222,7 +222,7 @@ We now visualize the posterior of the intercept under all those priors.
222222
223223``` {r, echo = -c(1:2)}
224224if (pdfplots) {
225- pdf("7-2_3.pdf", width = 8 , height = 5)
225+ pdf("7-2_3.pdf", width = 12 , height = 5, lwd = 1. 5)
226226}
227227par(mfrow = c(1, 2), mar = c(2.5, 1.5, 1.5, .1), mgp = c(1.5, .5, 0))
228228dens_improper <- density(res_improper$betas[, 1], bw = "SJ", adj = 2)
@@ -254,9 +254,9 @@ variance.
254254
255255``` {r, echo = -c(1:2)}
256256if (pdfplots) {
257- pdf("7-2_4.pdf", width = 8 , height = 5)
257+ pdf("7-2_4.pdf", width = 9 , height = 5)
258258}
259- par(mfrow = c(4, 2), mar = c(2.5, 1.5, 1.5, .1), mgp = c(1.5, .5, 0))
259+ par(mfrow = c(4, 2), mar = c(2.5, 1.5, 1.5, .1), mgp = c(1.5, .5, 0), lwd = 1.5 )
260260for (p in 1:3) {
261261 dens_improper <- density(res_improper$betas[, p + 1], bw = "SJ", adj = 2)
262262 dens_semi_1 <- density(res_semi_1$betas[, p + 1], bw = "SJ", adj = 2)
@@ -317,7 +317,7 @@ stationarity for $p = 1, 2, 3$.
317317if (pdfplots) {
318318 pdf("7-2_5.pdf", width = 9, height = 3)
319319}
320- par(mfrow = c(1, 3), mar = c(2.5, 2.5, 1.2, .1), mgp = c(1.5, .5, 0))
320+ par(mfrow = c(1, 3), mar = c(2.5, 2.5, 1.2, .1), mgp = c(1.5, .5, 0), lwd = 1.5 )
321321hist(res[[1]]$betas[,2], breaks = 20, freq = FALSE, main = "AR(1)",
322322 xlab = bquote(phi), ylab = "")
323323
@@ -349,8 +349,8 @@ First, we plot the data and its empirical autocorrelation function.
349349
350350``` {r, echo = -c(1:2)}
351351if (pdfplots) {
352- pdf("7-2_6.pdf", width = 8 , height = 5)
353- par(mar = c(2.5 , 1.5, 1.5, .1), mgp = c(1.5, .5, 0))
352+ pdf("7-2_6.pdf", width = 12 , height = 5)
353+ par(mar = c(2.6 , 1.5, 1.5, .1), mgp = c(1.5, .5, 0), lwd = 1.5 )
354354}
355355par(mfrow = c(1,2))
356356ts.plot(inflation, main = "EU Inflation")
@@ -376,7 +376,7 @@ for (p in 1:3) {
376376if (pdfplots) {
377377 pdf("7-2_7.pdf", width = 9, height = 3)
378378}
379- par(mfrow = c(1, 3), mar = c(2.5, 2.5, 1.2, .1), mgp = c(1.5, .5, 0))
379+ par(mfrow = c(1, 3), mar = c(2.5, 2.5, 1.2, .1), mgp = c(1.5, .5, 0), lwd = 1.5 )
380380hist(res[[1]]$betas[,2], breaks = 20, freq = FALSE, main = "AR(1)",
381381 xlab = bquote(phi), ylab = "")
382382abline(v = 1, col = 2)
0 commit comments