rm(list = ls())
setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
load("post-c-samp.Rda")
# When running without case 18
# comment previous "load" and uncomment next line
# load("post-c-samp-18.Rda")
iterates=length(beta.c[1,])
a <- c( 1, .577 )
b <- c( .577, 1 )
# define logit function
logit <- function(p)
{
lg = log(p/(1-p))
return(lg)}

prior_p1 <- rbeta(iterates, a[1], b[1])
prior_p2 <- rbeta(iterates, a[2], b[2])
prior_beta_c_1 <- 0.2717391 * logit( prior_p1 ) +
                     0.7282609 * logit( prior_p2 )
prior_beta_c_2 <- -0.05 * logit( prior_p1 ) +
                          0.05 * logit( prior_p2 )

beta_c_frame <- data.frame( cbind( prior_beta_c_1,
        prior_beta_c_2, beta.c[1,], beta.c[2,] ) )

# When running without case 18
# comment out next 2 commands
names( beta_c_frame ) <- c( "prior_beta_c_1",
"prior_beta_c_2","posterior_beta_c_1", "posterior_beta_c_2" )
save(beta_c_frame,file="beta-c-samp.Rda")
# When running without case 18 uncomment next 2 commands (4 lines)
#names( beta_c_frame ) <- c( "prior_beta_c_18_1",
#"prior_beta_c_18_2","posterior_beta_c_18_1",
# "posterior_beta_c_18_2" )
#save(beta_c_frame,file="beta-c-samp-18.Rda")