Nothing
## ---- fig.show='hold'---------------------------------------------------------
library(bayesAB)
A_binom <- rbinom(250, 1, .25)
B_binom <- rbinom(250, 1, .2)
## ---- fig.width=7.3-----------------------------------------------------------
plotBeta(100, 200) # looks a bit off
plotBeta(65, 200) # perfect
## ---- fig.show='hold'---------------------------------------------------------
AB1 <- bayesTest(A_binom, B_binom, priors = c('alpha' = 65, 'beta' = 200), n_samples = 1e5, distribution = 'bernoulli')
## ---- fig.width=7.3-----------------------------------------------------------
print(AB1)
summary(AB1)
plot(AB1)
## ---- fig.width=7.3-----------------------------------------------------------
A_pois <- rpois(250, 6.5)
B_pois <- rpois(250, 5.5)
plotGamma(30, 5) # 5-6 seem likely enough
## ---- fig.width=7.3-----------------------------------------------------------
AB2 <- bayesTest(A_pois, B_pois, priors = c('shape' = 30, 'rate' = 5), n_samples = 1e5, distribution = 'poisson')
print(AB2)
summary(AB2)
plot(AB2)
## -----------------------------------------------------------------------------
AB3 <- combine(AB1, AB2, f = `*`, params = c('Probability', 'Lambda'), newName = 'Expectation')
# also equivalent with %>% if you like piping
library(magrittr)
AB3 <- AB1 %>%
combine(AB2, f = `*`, params = c('Probability', 'Lambda'), newName = 'Expectation')
## ---- fig.width=7.3-----------------------------------------------------------
print(AB3)
summary(AB3)
plot(AB3)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.