plot(Rpost$supp, Rpost$pmf, col="black", type="l", xlim=c(0.5, 1.5), xlab="", ylab="")
abline(h=1/(20/0.01+1), col="red")
abline(v=res1$Rhat, col="blue")
abline(v=Rpost$supp[loc], col="purple")
plot(Rpost$supp, Rpost$pmf, col="black", type="l", xlim=c(0.5, 1.5), xlab="", ylab="")
abline(h=1/(20/0.01+1), col="red")
abline(v=res1$Rhat, col="blue")
abline(v=Rpost$supp[loc], col="purple")