R/plot.DS_GF_micro.R

In BayesGOF: Bayesian Modeling via Frequentist Goodness-of-Fit

plot.DS_GF_micro <-
function(x, ...){
	if(dim(x$post.fit)[2]==2){
	plot(x$post.fit$theta.vals, x$post.fit$parm.pos,  
		type = "l", lwd = 2, col = "blue", lty = "dashed",
		xlab = expression(theta), ylab = "", font.main = 1,
		cex.lab=1.45, cex.axis=1.5, cex.main=1.75, cex.sub=1.5, ...)
	title(ylab = expression(paste(hat(pi)(theta[i]~"|"~y[i]))), line = 2.3, cex.lab=1.45)
	} else {
	plot(x$post.fit$theta.vals, x$post.fit$ds.pos,  
		type = "l", lwd = 2, col = "red3",
		xlab = expression(theta), ylab = "", font.main = 1,
		cex.lab=1.45, cex.axis=1.5, cex.main=1.75, cex.sub=1.5, ...)
	title(ylab = expression(paste(hat(pi)(theta[i]~"|"~y[i]))), line = 2.3, cex.lab=1.45)
	lines(x$post.fit$theta.vals, x$post.fit$parm.pos,
		col = "blue", lty = "dashed", lwd = 2)
	}
	}