R/plotConvergence.R

In birte: Bayesian Inference of Regulatory Influence on Expression (biRte)

plotConvergence = function(res, title="") {
	if(! "log_lik_trace" %in% names(res))
		stop("Need entry log_lik_trace!")
	if(! "nburnin" %in% names(res))
		stop("Need entry nburnin!")

	take = which(res$log_lik_trace != Inf)
	if(length(take) > 0){
		plot(take, res$log_lik_trace[take], type="l", main=paste("log-likelihood along MCMC sampling", title), xlab="step", ylab="log likelihood")
		abline(v=res$nburnin, lty=3)		
	}
}