pp.plot2ci <-
function (x, segwd=2, segcol="steelblue", ...) 
{
	# placed in the public domain by Burns Statistics 2011

	# columns 1 and 4 are the estimates
	# columns 2 and 3 are the confidence limits for column 1
	# columns 5 and 6 are the confidence limits for column 4

	plot(x[,4], x[,1], xlim=range(x[,4:6]), ylim=range(x[,1:3]), ...)
	segments(x[,5], x[,2], y1=x[,3], lwd=segwd, col=segcol)
	segments(x[,6], x[,2], y1=x[,3], lwd=segwd, col=segcol)
	segments(x[,5], x[,2], x1=x[,6], lwd=segwd, col=segcol)
	segments(x[,5], x[,3], x1=x[,6], lwd=segwd, col=segcol)
}