Introducción al pensamiento Bayesiano

Hay una probabilidad del 90% de que la proporción de interés se encuentre entre 0.256 y 0.513.Las medidas de resumen anteriores son exactas ya que son obtenidas usando las funciones en Rpara la densidad posterior beta. Un método alternativo para el cálculo de estas medidas deresumen se basa en la simulación. Es posible simular una gran cantidad de valores a partir de ladensidad posterior beta y luego calcular las medidas de resumen sobre estos valores. Usando elcomando rbeta se simulan 1000 valores para p:> ps = rbeta(1000, alpha+x, beta+n-x)y graficando la distribución posterior a partir de un histograma obtenido sobre los valoressimulados:> hist(ps, xlab="p", main="")