Cálculos de quantis distribuição beta

beta.bayes.R

###
# Brincando de Estatística Bayesiana
#
require(tidyverse)
n1 <- 9.5
n2 <- 0.5

n<-68
n.inocente<- 68
n.culpado <- n-n.inocente


# priori
qbeta(c(0.025,0.975), n1,n2)
#posteriori
qbeta(c(0.025,0.975), n1+n.culpado,n2+n.inocente)

x<-rbeta(1000, n1+n.culpado,n2+n.inocente)
y<-rbeta(1000, n1,n2)

sim.beta <-
  data.frame(
    data=
      c(y,x),
    kind=
      factor(
        c(
          rep('Priori',length(y)),
          rep('Posteriori',length(x))
        )
        )
  )

ggplot(data=sim.beta, aes(data, fill=kind)) +
  geom_density()