two_proportions.R
In TeachBayes: Teaching Bayesian Inference

## ----global_options, include=FALSE---------------------------------------
knitr::opts_chunk$set(fig.width=6, 
                      fig.height=4, fig.path='Figs/',
                      echo=TRUE, warning=FALSE, message=FALSE)

## ------------------------------------------------------------------------
library(TeachBayes)

## ------------------------------------------------------------------------
prior <- testing_prior(.05, .95, 10,
                       uniform=TRUE)

## ------------------------------------------------------------------------
draw_two_p(prior)

## ------------------------------------------------------------------------
(diff_dist <- two_p_summarize(prior))
prob_plot(diff_dist)

## ------------------------------------------------------------------------
y1n1 <- c(10, 20)
y2n2 <- c(8, 24)

## ------------------------------------------------------------------------
post <- two_p_update(prior, y1n1, y2n2)

## ------------------------------------------------------------------------
draw_two_p(post)

## ------------------------------------------------------------------------
prob_plot(two_p_summarize(post))

## ------------------------------------------------------------------------
prior <- testing_prior(.05, .95, 10, pequal=0.5)

## ------------------------------------------------------------------------
draw_two_p(prior)

## ------------------------------------------------------------------------
prob_plot(two_p_summarize(prior))

## ------------------------------------------------------------------------
y1n1 <- c(10, 20)
y2n2 <- c(8, 24)

## ------------------------------------------------------------------------
post <- two_p_update(prior, y1n1, y2n2)

## ------------------------------------------------------------------------
draw_two_p(post)

## ------------------------------------------------------------------------
prob_plot(two_p_summarize(post))