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]]