Dear all,
 
I have a very simple question about how I can include HPD confidence
lines in a Marginal posterior distribution scatterplot
I have the following code to draw the scatterplot(s) from an output
table;
 
par(mfrow=c(2,2))
 
                plot(temp1[,5],  temp1[,6],  xlab="log10(r)",
ylab="log10(tf)", pch=46)
                plot(temp2[,5],  temp2[,6],  xlab="log10(r)",
ylab="log10(tf)", pch=46)
                plot(temp3[,5],  temp3[,6],  xlab="log10(r)",
ylab="log10(tf)", pch=46)
                plot(temp4[,5],  temp4[,6],  xlab="log10(r)",
ylab="log10(tf)", pch=46)
 
                dev.print(device=pdf, height=12, width=10,
file="Marginal.pdf")
 
I would very much like to include say 0.1, 0.5 and 0.99 contours onto
the graphs, but I can't figure it out using any examples of how to do
it.  Can someone offer any suggestions?
 
Many thanks
Peter
 
 
	[[alternative HTML version deleted]]