"boa.plot.acf" <-
function(lname, pname, annotate = boa.par("legend"))
#-------------------------------------------------------------------------------
#-------------------------------------------------------------------------------
{
drawn <- FALSE
parm <- boa.getparms(boa.chain("work")[[lname]], pname)
if(is.matrix(parm)) {
drawn <- TRUE
val <- boa.par("par")
cex <- ifelse(is.null(val$cex), 1, val$cex)
lwd <- ifelse(is.null(val$lwd), 1, val$lwd)
result <- acf(parm, plot = FALSE)
plot(result$lag, result$acf, xlab = "Lag", ylab = "Autocorrelation",
main = pname, ylim = c(-1, 1), type = "h", lwd = lwd)
abline(0, 0)
usr <- par("usr")
if (annotate)
legend(x = usr[2], y = 1, xjust = 1, yjust = 1, cex = cex,
legend = substring(lname, first = 1, last = 16), bty = "n")
}
return(drawn)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.