Skip to content

Instantly share code, notes, and snippets.

@toyeiei
Created January 30, 2019 03:28
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save toyeiei/a4528a373dbf7e25ed4048d0e4ec4000 to your computer and use it in GitHub Desktop.
Save toyeiei/a4528a373dbf7e25ed4048d0e4ec4000 to your computer and use it in GitHub Desktop.
# create a prior probability vector
prior <- 0:10 * 0.1
# writing what we know
beta <- 0.20
alpha <- 0.05
power <- 1 - beta
# bayes theorem
posterior <- power * prior / (power * prior + alpha * (1-prior))
round(posterior, 2)
# plot posterior
plot(prior, posterior, pch = 16, type = "b")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment