R/plot.DS_FB_obj.R

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

plot.DS_FB_obj <-
function(x, plot.type = "posterior",...){

if(dim(x$post.fit)[2] == 4){	
	if(plot.type == "posterior"){
	par(mar=c(5,5,4,2)+0.3) #changes left margin to make large labels fit
	plot(x$post.fit$theta.vals, x$post.fit$ds.post, 
		type = "l", yaxs = "i", 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.post, lty = "dashed", col = "blue", lwd = 2)
	lines(x$post.fit$theta.vals, x$post.fit$finite.post, col = "darkgreen", lwd = 2)
	} else {
		par(mar=c(5,5,4,2)+0.3) #changes left margin to make large labels fit
		plot(x$prior.fit$theta.vals, x$prior.fit$ds.prior, 
			type = "l", yaxs = "i", 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$prior.fit$theta.vals, x$prior.fit$parm.prior, col = "blue", lty = "dashed", lwd = 2)
		lines(x$prior.fit$theta.vals, x$prior.fit$finite.prior, col = "darkgreen", lwd = 2)
	}
} else {
	if(plot.type == "posterior"){
	par(mar=c(5,5,4,2)+0.3) #changes left margin to make large labels fit
	plot(x$post.fit$theta.vals, x$post.fit$parm.post, 
		type = "l", yaxs = "i", lwd = 2, col = "blue",  
		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$finite.post, col = "darkgreen", lwd = 2)
	} else {
		par(mar=c(5,5,4,2)+0.3) #changes left margin to make large labels fit
		plot(x$prior.fit$theta.vals, x$prior.fit$parm.prior, 
			type = "l", yaxs = "i", lwd = 2, col = "blue",
			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$prior.fit$theta.vals, x$prior.fit$finite.prior, col = "darkgreen", lwd = 2)
	}
}
}