@@ -17,7 +17,7 @@ knitr::opts_chunk$set(
1717 fig.height = 6
1818)
1919knitr::opts_knit$set(global.par = TRUE)
20- pdfplots <- FALSE # default: FALSE; set this to TRUE only if you like pdf figures
20+ pdfplots <- TRUE # default: FALSE; set this to TRUE only if you like pdf figures
2121```
2222
2323``` {r, include = FALSE}
@@ -146,7 +146,7 @@ We plot the marginal posteriors (in blue) together with those under the improper
146146``` {r}
147147if (pdfplots) {
148148 pdf("6-4_2.pdf", width = 8, height = 3)
149- par(mar = c(2.5, 1.5, .1 , .1), mgp = c(1.6, .6, 0))
149+ par(mar = c(2.5, 1.5, 1.5 , .1), mgp = c(1.6, .6, 0))
150150}
151151par(mfrow = c(1, 3))
152152for (i in seq_len(nrow(beta.hat))) {
@@ -156,13 +156,13 @@ for (i in seq_len(nrow(beta.hat))) {
156156 curve(dt((x- res_conj1$bN[i])/post.sd.conj1[i], df=2* res_conj1$cN),
157157 from= res_conj1$bN[i]- 4*post.sd.conj1[i],
158158 to= res_conj1$bN[i]+ 4*post.sd.conj1[i],
159- add=TRUE, col=2,lty=2)
159+ add=TRUE, col=2,lty=2,lwd=2 )
160160 curve(dt((x- res_conj2$bN[i])/post.sd.conj2[i], df=2* res_conj2$cN),
161161 from= res_conj2$bN[i]- 4*post.sd.conj2[i],
162162 to= res_conj2$bN[i]+ 4*post.sd.conj2[i],
163- add=TRUE, col=3,lty=3)
163+ add=TRUE, col=3,lty=3,lwd=2 )
164164 legend("topright", c("improper", "B0=10", "B0=1"),
165- col = 1:3, lty = 1:3)
165+ col = 1:3, lty = 1:3,lwd=c(1,2,2) )
166166}
167167```
168168There is little difference to the improper prior for $B_0=10\textbf{I}$, however
0 commit comments