"boa.plot.trace" <-
function(lnames, pname, annotate = boa.par("legend"))
#-------------------------------------------------------------------------------
#-------------------------------------------------------------------------------
{
work <- boa.chain("work")
ipname <- list()
ybar <- list()
xlim <- NULL
ylim <- NULL
key.names <- NULL
lnames <- intersect(names(work), lnames)
for(i in lnames) {
ipname[[i]] <- intersect(boa.pnames(work[[i]]), pname)
ylim <- range(ylim, work[[i]][, ipname[[i]]])
xlim <- range(xlim, boa.iter(work[[i]]))
key.names <- c(key.names, substring(i, first = 1, last = 16))
}
drawn <- is.vector(key.names)
if(drawn) {
val <- boa.par("par")
cex <- ifelse(is.null(val$cex), 1, val$cex)
lwd <- ifelse(is.null(val$lwd), 1, val$lwd)
plot(xlim, ylim, xlab = "Iteration", ylab = pname,
xlim = xlim, ylim = ylim, type = "n")
k <- 0
for(i in lnames) {
for(j in ipname[[i]]) {
k <- k + 1
lines(boa.iter(work[[i]]), work[[i]][, j], lty = k, lwd = lwd)
}
}
if(annotate)
legend(x = xlim[2], y = ylim[2], xjust = 1, yjust = 1,
legend = key.names, lty = 1:k, bty = "n",
cex = cex, lwd = lwd)
}
return(drawn)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.