rm(list = ls())
setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
load("beta-samples.Rda")
library(ggplot2)
# library(boot)  #this was part of Fletch's program
# but I don't see it being used anywhere
# "boot" is a bootstrapping library
ggplot( beta_frame ) +
  geom_density2d( aes(x=prior_beta1, y=prior_beta2),
  color="#000066", linetype=2 ) +
  geom_density2d( aes(x=posterior_beta1, y=posterior_beta2),
  color="#0000cc" ) +
  theme_bw() +
  xlab( expression( beta[0] ) ) +
  ylab( expression( beta[1] ) ) +
  ggtitle( expression(
  paste( "Approximate prior and posterior for ",
  beta[0], ", ", beta[1] ) ) )