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]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.