Baseado numa questão de um post no G+ do meu amigo Thiago Santos:

Considere um supermercado de classe média, como um Pão de Açúcar (ou um colégio particular, ou um clube, ou um restaurante). Assuma como prior, que o pais é justo e que há iguais chances de você ver brancos e negros entre os clientes do supermercado e que a população afrodescendente é de ao menos 50% na nossa população. Agora, tome uma observação: suponha que você observou 200 clientes no supermercado e dentre eles apenas 10 eram afrodescendentes. Como isso afeta seu posterior, isto é, como você revê sua crença na igualdade considerando tais dados?



In [1]:

    
%matplotlib inline
from IPython.core.pylabtools import figsize
from IPython.html.widgets import interact
import numpy as np
from matplotlib import pyplot as plt
figsize(11, 9)

import scipy.stats as stats

dist = stats.beta
x = np.linspace(0, 1, 100)

def plot_posterior(**kwargs):
    populacao = kwargs['populacao']
    num_ad = kwargs['afrodescendentes']
    
    fig, ax = plt.subplots(figsize=(4, 3),
                           subplot_kw={'axisbg':'#EEEEEE',
                                       'axisbelow':True})
    ax.grid(color='w', linewidth=2, linestyle='solid')
    
    y = dist.pdf(x, 1 + num_ad, 1 + populacao - num_ad)
    ax.plot(x, y, label="Universo de %d pessoas,\n %d afrodescendentes vistos" % (populacao, num_ad))
    ax.fill_between(x, 0, y, color="#348ABD", alpha=0.4)
    ax.vlines(0.5, 0, 4, color="k", linestyles="--", lw=1)

    leg = plt.legend()
    leg.get_frame().set_alpha(0.4)
    ax.autoscale(tight=True)

i = interact(plot_posterior, populacao=(1,400), afrodescendentes=(1,200))



In [ ]: