Nothing
#' Plotting all parameters
#'
#' Plots all parameter values, including tree topologies (see makeplot.topology) over the length of the MCMC chain
#'
#' @param chains A set of rwty.chain objects
#' @param burnin The number of trees to omit as burnin.
#' @param facet Boolean denoting whether to make a facet plot.
#' @param free_y TRUE/FALSE to turn free y scales on the facetted plots on or off (default FALSE). Only works if facet = TRUE.
#' @param strip Number indicating which column to strip off (i.e., strip=1 removes first column, which is necessary for most MCMC outputs in which the first column is just the generation).
#' You can skip multiple columns by passing a vector of columns to skip, e.g., strip=c(1,4,6).
#'
#' @return param.plot Returns a list of ggplot objects.
#'
#' @keywords parameter, plot, convergence, mcmc, phylogenetics
#'
#' @export makeplot.all.params
#' @examples
#' \dontrun{
#' data(fungus)
#' makeplot.all.params(fungus, burnin=20)
#' }
makeplot.all.params <- function(chains, burnin = 0, facet=TRUE, free_y=FALSE, strip = 1){
chains = check.chains(chains)
chain = chains[[1]]
if(is.null(chain$ptable)) stop("No parameters associated with your chains")
params <- names(chain$ptable)[-strip]
param.plots <- lapply(params, FUN = function(x) makeplot.param(parameter = x, burnin = burnin, chains = chains, facet = facet))
t.plot = makeplot.topology(chains, burnin = burnin, facet = facet)
param.plots[[length(param.plots)+1]] <- t.plot
plot.names = c(paste(params, ".trace", sep=""), "topology.trace.plot")
names(param.plots) <- plot.names
return(param.plots)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.