R/plot_mcmctraces.R

In ddepn: Dynamic Deterministic Effects Propagation Networks: Infer signalling networks for timecourse RPPA data.

# Plot MCMC traces from inhibMCMC chains
# requires package coda
#
# Author: benderc
###############################################################################


plot_mcmctraces <- function(ret, thin=1) {
	samp <- ret$samplings
	ltr <- ret$ltraces
	lst <- list()
	for(i in 1:ncol(ltr)) {
		mcmcobj <- mcmc(ltr[,i],thin=thin)
		lst[[i]] <- mcmcobj
	}
	mcmclst <- mcmc.list(lst)
	plot(mcmclst, main=paste("Posterior traces,",ncol(ltr), "chains"),auto.layout=FALSE,density=FALSE)
}